Application of a new laser Doppler imaging system in planning and monitoring of surgical flaps

There is a demand for technologies able to assess the perfusion of surgical flaps quantitatively and reliably to avoid ischemic complications. The aim of this study is to test a new high-speed high-definition laser Doppler imaging (LDI) system (FluxEXPLORER, Microvascular Imaging, Lausanne, Switzerland) in terms of preoperative mapping of the vascular supply (perforator vessels) and postoperative flow monitoring. The FluxEXPLORER performs perfusion mapping of an area 9 x 9 cm with a resolution of 256 x 256 pixels within 6 s in high-definition imaging mode. The sensitivity and predictability to localize perforators is expressed by the coincidence of preoperatively assessed LDI high flow spots with intraoperatively verified perforators in nine patients. 18 free flaps are monitored before, during, and after total ischemia. 63% of all verified perforators correspond to a high flow spot, and 38% of all high flow spots correspond to a verified perforator (positive predictive value). All perfused flaps reveal a value of above 221 perfusion units (PUs), and all values obtained in the ischemic flaps are beneath 187 PU. In summary, we conclude that the present LDI system can serve as a reliable, fast, and easy-to-handle tool to detect ischemia in free flaps, whereas perforator vessels cannot be detected appropriately.

Assessment of melanocytic skin lesions with a high-definition laser Doppler imaging system

Early detection is a major goal in the management of malignant melanoma. Besides clinical assessment many noninvasive technologies such as dermoscopy, digital dermoscopy and in vivo laser scanner microscopy are used as additional methods. Herein we tested a system to assess lesional perfusion as a tool for early melanoma detection.

A Compensating Anastigmatic Submillimeter Array Imaging System for STEAMR

In this paper, we present a novel technique for the removal of astigmatism in submillimeter-wave optical systems through employment of a specific combination of so-called astigmatic off-axis reflectors. This technique treats an orthogonally astigmatic beam using skew Gaussian beam analysis, from which an anastigmatic imaging network is derived. The resultant beam is considered truly stigmatic, with all Gaussian beam parameters in the orthogonal directions being matched. This is thus considered an improvement over previous techniques wherein a beam corrected for astigmatism has only the orthogonal beam amplitude radii matched, with phase shift and phase radius of curvature not considered. This technique is computationally efficient, negating the requirement for computationally intensive numerical analysis of shaped reflector surfaces. The required optical surfaces are also relatively simple to implement compared to such numerically optimized shaped surfaces. This technique is implemented in this work as part of the complete optics train for the STEAMR antenna. The STEAMR instrument is envisaged as a mutli-beam limb sounding instrument operating at submillimeter wavelengths. The antenna optics arrangement for this instrument uses multiple off-axis reflectors to control the incident radiation and couple them to their corresponding receiver feeds. An anastigmatic imaging network is successfully implemented into an optical model of this antenna, and the resultant design ensures optimal imaging of the beams to the corresponding feed horns. This example also addresses the challenges of imaging in multi-beam antenna systems.

Transconvolution and the virtual positron emission tomograph--a new method for cross calibration in quantitative PET∕CT imaging

PURPOSE Positron emission tomography (PET)∕computed tomography (CT) measurements on small lesions are impaired by the partial volume effect, which is intrinsically tied to the point spread function of the actual imaging system, including the reconstruction algorithms. The variability resulting from different point spread functions hinders the assessment of quantitative measurements in clinical routine and especially degrades comparability within multicenter trials. To improve quantitative comparability there is a need for methods to match different PET∕CT systems through elimination of this systemic variability. Consequently, a new method was developed and tested that transforms the image of an object as produced by one tomograph to another image of the same object as it would have been seen by a different tomograph. The proposed new method, termed Transconvolution, compensates for differing imaging properties of different tomographs and particularly aims at quantitative comparability of PET∕CT in the context of multicenter trials. METHODS To solve the problem of image normalization, the theory of Transconvolution was mathematically established together with new methods to handle point spread functions of different PET∕CT systems. Knowing the point spread functions of two different imaging systems allows determining a Transconvolution function to convert one image into the other. This function is calculated by convolving one point spread function with the inverse of the other point spread function which, when adhering to certain boundary conditions such as the use of linear acquisition and image reconstruction methods, is a numerically accessible operation. For reliable measurement of such point spread functions characterizing different PET∕CT systems, a dedicated solid-state phantom incorporating (68)Ge∕(68)Ga filled spheres was developed. To iteratively determine and represent such point spread functions, exponential density functions in combination with a Gaussian distribution were introduced. Furthermore, simulation of a virtual PET system provided a standard imaging system with clearly defined properties to which the real PET systems were to be matched. A Hann window served as the modulation transfer function for the virtual PET. The Hann's apodization properties suppressed high spatial frequencies above a certain critical frequency, thereby fulfilling the above-mentioned boundary conditions. The determined point spread functions were subsequently used by the novel Transconvolution algorithm to match different PET∕CT systems onto the virtual PET system. Finally, the theoretically elaborated Transconvolution method was validated transforming phantom images acquired on two different PET systems to nearly identical data sets, as they would be imaged by the virtual PET system. RESULTS The proposed Transconvolution method matched different PET∕CT-systems for an improved and reproducible determination of a normalized activity concentration. The highest difference in measured activity concentration between the two different PET systems of 18.2% was found in spheres of 2 ml volume. Transconvolution reduced this difference down to 1.6%. In addition to reestablishing comparability the new method with its parameterization of point spread functions allowed a full characterization of imaging properties of the examined tomographs. CONCLUSIONS By matching different tomographs to a virtual standardized imaging system, Transconvolution opens a new comprehensive method for cross calibration in quantitative PET imaging. The use of a virtual PET system restores comparability between data sets from different PET systems by exerting a common, reproducible, and defined partial volume effect.

Ability of dGEMRIC and T2 mapping to evaluate cartilage repair after microfracture: a goat study

OBJECTIVE: To investigate the ability of delayed gadolinium-enhanced magnetic resonance (MR) imaging of cartilage (dGEMRIC) and T2 mapping to evaluate the quality of repair tissue after microfracture. DESIGN: Twelve knees from 12 goats were studied. An osteochondral defect (diameter, 6mm; depth, 3mm) with microfracture was created in the weight-bearing aspect of both the medial and lateral femoral condyles. Goats were euthanized at 24 weeks (n=6) and 48 weeks (n=6) postsurgery. Pre-contrast R1 (R1pre) and post-contrast R1 (R1post) measurements for dGEMRIC and a pre-contrast T2 measurement for T2 mapping were performed with a 3T MR imaging system. MR imaging findings were compared with histological and biochemical assessments. RESULTS: In native cartilage, significant correlations were observed between the R1post and the glycosaminoglycan (GAG) concentration, as well as DeltaR1 (difference between the R1pre and R1post) and the GAG concentration (P<0.05). In repair tissue, a significant correlation was observed between DeltaR1 and the GAG concentration (P<0.05), but not between the R1post and the GAG concentration. In both repair tissue and native cartilage, no correlation was observed between T2 and the water concentration or between T2 and the hydroxyproline (HP) concentration. A zonal variation of T2 and a clear dependence of T2 on the angles relative to B0 were observed in native cartilage, but not in repair tissue. CONCLUSION: dGEMRIC with DeltaR1 measurement might be useful for the evaluation of the GAG concentration in repair tissue after microfracture. T2 mapping might be useful for the differentiation of repair tissue after microfracture from native cartilage; however, its potential to assess the specific biochemical markers in native cartilage as well as repair tissue may be limited.

Modifying peripheral IV catheters with side holes and side slits results in favorable changes in fluid dynamic properties during the injection of iodinated contrast material

OBJECTIVE: The purpose of this study was to compare a standard peripheral end-hole angiocatheter with those modified with side holes or side slits using experimental optical techniques to qualitatively compare the contrast material exit jets and using numeric techniques to provide flow visualization and quantitative comparisons. MATERIALS AND METHODS: A Schlieren imaging system was used to visualize the angiocatheter exit jet fluid dynamics at two different flow rates. Catheters were modified by drilling through-and-through side holes or by cutting slits into the catheters. A commercial computational fluid dynamics package was used to calculate numeric results for various vessel diameters and catheter orientations. RESULTS: Experimental images showed that modifying standard peripheral IV angiocatheters with side holes or side slits qualitatively changed the overall flow field and caused the exiting jet to become less well defined. Numeric calculations showed that the addition of side holes or slits resulted in a 9-30% reduction of the velocity of contrast material exiting the end hole of the angiocatheter. With the catheter tip directed obliquely to the wall, the maximum wall shear stress was always highest for the unmodified catheter and was always lowest for the four-side-slit catheter. CONCLUSION: Modified angiocatheters may have the potential to reduce extravasation events in patients by reducing vessel wall shear stress.

Evaluation of a new six degrees of freedom couch for radiation therapy

Purpose: The aim of this work is to evaluate the geometric accuracy of a prerelease version of a new six degrees of freedom (6DoF) couch. Additionally, a quality assurance method for 6DoF couches is proposed. Methods: The main principle of the performance tests was to request a known shift for the 6DoF couch and to compare this requested shift with the actually applied shift by independently measuring the applied shift using different methods (graph paper, laser, inclinometer, and imaging system). The performance of each of the six axes was tested separately as well as in combination with the other axes. Functional cases as well as realistic clinical cases were analyzed. The tests were performed without a couch load and with a couch load of up to 200 kg and shifts in the range between −4 and +4 cm for the translational axes and between −3° and +3° for the rotational axes were applied. The quality assurance method of the new 6DoF couch was performed using a simple cube phantom and the imaging system. Results: The deviations (mean ± one standard deviation) accumulated over all performance tests between the requested shifts and the measurements of the applied shifts were −0.01 ± 0.02, 0.01 ± 0.02, and 0.01 ± 0.02 cm for the longitudinal, lateral, and vertical axes, respectively. The corresponding values for the three rotational axes couch rotation, pitch, and roll were 0.03° ± 0.06°, −0.04° ± 0.12°, and −0.01° ± 0.08°, respectively. There was no difference found between the tests with and without a couch load of up to 200 kg. Conclusions: The new 6DoF couch is able to apply requested shifts with high accuracy. It has the potential to be used for treatment techniques with the highest demands in patient setup accuracy such as those needed in stereotactic treatments. Shifts can be applied efficiently and automatically. Daily quality assurance of the 6DoF couch can be performed in an easy and efficient way. Long-term stability has to be evaluated in further tests.

67P/Churyumov-Gerasimenko: Activity between March and June 2014 as observed from Rosetta/OSIRIS

Aims. 67P/Churyumov-Gerasimenko is the target comet of the ESA’s Rosetta mission. After commissioning at the end of March 2014, the Optical, Spectroscopic, and Infrared Remote Imaging System (OSIRIS) onboard Rosetta, started imaging the comet and its dust environment to investigate how they change and evolve while approaching the Sun. Methods. We focused our work on Narrow Angle Camera (NAC) orange images and Wide Angle Camera (WAC) red and visible-610 images acquired between 2014 March 23 and June 24 when the nucleus of 67P was unresolved and moving from approximately 4.3 AU to 3.8 AU inbound. During this period the 67P – Rosetta distance decreased from 5 million to 120 thousand km. Results. Through aperture photometry, we investigated how the comet brightness varies with heliocentric distance. 67P was likely already weakly active at the end of March 2014, with excess flux above that expected for the nucleus. The comet’s brightness was mostly constant during the three months of approach observations, apart from one outburst that occurred around April 30 and a second increase in flux after June 20. Coma was resolved in the profiles from mid-April. Analysis of the coma morphology suggests that most of the activity comes from a source towards the celestial north pole of the comet, but the outburst that occurred on April 30 released material in a different direction.

The rotation state of 67P/Churyumov-Gerasimenko from approach observations with the OSIRIS cameras on Rosetta

Aims. Approach observations with the Optical, Spectroscopic, and Infrared Remote Imaging System (OSIRIS) experiment onboard Rosetta are used to determine the rotation period, the direction of the spin axis, and the state of rotation of comet 67P’s nucleus. Methods. Photometric time series of 67P have been acquired by OSIRIS since the post wake-up commissioning of the payload in March 2014. Fourier analysis and convex shape inversion methods have been applied to the Rosetta data as well to the available ground-based observations. Results. Evidence is found that the rotation rate of 67P has significantly changed near the time of its 2009 perihelion passage, probably due to sublimation-induced torque. We find that the sidereal rotation periods P1 = 12.76129 ± 0.00005 h and P2 = 12.4043 ± 0.0007 h for the apparitions before and after the 2009 perihelion, respectively, provide the best fit to the observations. No signs of multiple periodicity are found in the light curves down to the noise level, which implies that the comet is presently in a simple rotation state around its axis of largest moment of inertia. We derive a prograde rotation model with spin vector J2000 ecliptic coordinates λ = 65° ± 15°, β = + 59° ± 15°, corresponding to equatorial coordinates RA = 22°, Dec = + 76°. However, we find that the mirror solution, also prograde, at λ = 275° ± 15°, β = + 50° ± 15° (or RA = 274°, Dec = + 27°), is also possible at the same confidence level, due to the intrinsic ambiguity of the photometric problem for observations performed close to the ecliptic plane.

Fractures on comet 67P/Churyumov-Gerasimenko observed by Rosetta/OSIRIS

The Optical, Spectroscopic, and Infrared Remote Imaging System (OSIRIS) experiment onboard the Rosetta spacecraft currently orbiting comet 67P/Churyumov-Gerasimenko has yielded unprecedented views of a comet's nucleus. We present here the first ever observations of meter-scale fractures on the surface of a comet. Some of these fractures form polygonal networks. We present an initial assessment of their morphology, topology, and regional distribution. Fractures are ubiquitous on the surface of the comet's nucleus. Furthermore, they occur in various settings and show different topologies suggesting numerous formation mechanisms, which include thermal insulation weathering, orbital-induced stresses, and possibly seasonal thermal contraction. However, we conclude that thermal insolation weathering is responsible for creating most of the observed fractures based on their morphology and setting in addition to thermal models that indicate diurnal temperature ranges exceeding 200K and thermal gradients of similar to 15K/min at perihelion are possible. Finally, we suggest that fractures could be a facilitator in surface evolution and long-term erosion.

Spectrophotometric properties of the nucleus of comet 67P/Churyumov-Gerasimenko from the OSIRIS instrument onboard the ROSETTA spacecraft

Context. The Rosetta mission of the European Space Agency has been orbiting the comet 67P/Churyumov-Gerasimenko (67P) since August 2014 and is now in its escort phase. A large complement of scientific experiments designed to complete the most detailed study of a comet ever attempted are onboard Rosetta. Aims. We present results for the photometric and spectrophotometric properties of the nucleus of 67P derived from the OSIRIS imaging system, which consists of a Wide Angle Camera (WAC) and a Narrow Angle Camera (NAC). The observations presented here were performed during July and the beginning of August 2014, during the approach phase, when OSIRIS was mapping the surface of the comet with several filters at different phase angles (1.3 degrees-54 degrees). The resolution reached up to 2.1 m/px. Methods. The OSIRIS images were processed with the OSIRIS standard pipeline, then converted into I/F. radiance factors and corrected for the illumination conditions at each pixel using the Lommel-Seeliger disk law. Color cubes of the surface were produced by stacking registered and illumination-corrected images. Furthermore, photometric analysis was performed both on disk-averaged photometry in several filters and on disk-resolved images acquired with the NAC orange filter, centered at 649 ran, using Hapke modeling. Results. The disk-averaged phase function of the nucleus of 67P shows a strong opposition surge with a G parameter value of -0.13 +/- 0.01 in the HG system formalism and an absolute magnitude H-v(1, 1, 0) = 15.74 +/- 0.02 mag. The integrated spectrophotometry in 20 filters covering the 250-1000 nm wavelength range shows a red spectral behavior, without clear absorption bands except for a potential absorption centered at similar to 290 rim that is possibly due to SO2 ice. The nucleus shows strong phase reddening, with disk-averaged spectral slopes increasing from 11%/( 100 nm) to 16%/(100 nm) in the 1.3 degrees-54 degrees phase angle range. The geometric albedo of the comet is 6.5 +/- 0.2% at 649 nm, with local variations of up to similar to 16% in the Hapi region. From the disk-resolved images we computed the spectral slope together with local spectrophotometry and identified three distinct groups of regions (blue, moderately red, and red). The Hapi region is the brightest, the bluest in term of spectral slope, and the most active surface on the comet. Local spectrophotometry shows an enhancement of the flux in the 700-750 nm that is associated with coma emissions.

Size-frequency distribution of boulders ≥7 m on comet 67P/Churyumov-Gerasimenko

Aims. We derive for the first time the size-frequency distribution of boulders on a comet, 67P/Churyumov-Gerasimenko (67P), computed from the images taken by the Rosetta/OSIRIS imaging system. We highlight the possible physical processes that lead to these boulder size distributions. Methods. We used images acquired by the OSIRIS Narrow Angle Camera, NAC, on 5 and 6 August 2014. The scale of these images (2.44−2.03 m/px) is such that boulders ≥7 m can be identified and manually extracted from the datasets with the software ArcGIS. We derived both global and localized size-frequency distributions. The three-pixel sampling detection, coupled with the favorable shadowing of the surface (observation phase angle ranging from 48° to 53°), enables unequivocally detecting boulders scattered all over the illuminated side of 67P. Results. We identify 3546 boulders larger than 7 m on the imaged surface (36.4 km2), with a global number density of nearly 100/km2 and a cumulative size-frequency distribution represented by a power-law with index of −3.6 +0.2/−0.3. The two lobes of 67P appear to have slightly different distributions, with an index of −3.5 +0.2/−0.3 for the main lobe (body) and −4.0 +0.3/−0.2 for the small lobe (head). The steeper distribution of the small lobe might be due to a more pervasive fracturing. The difference of the distribution for the connecting region (neck) is much more significant, with an index value of −2.2 +0.2/−0.2. We propose that the boulder field located in the neck area is the result of blocks falling from the contiguous Hathor cliff. The lower slope of the size-frequency distribution we see today in the neck area might be due to the concurrent processes acting on the smallest boulders, such as i) disintegration or fragmentation and vanishing through sublimation; ii) uplifting by gas drag and consequent redistribution; and iii) burial beneath a debris blanket. We also derived the cumulative size-frequency distribution per km2 of localized areas on 67P. By comparing the cumulative size-frequency distributions of similar geomorphological settings, we derived similar power-law index values. This suggests that despite the selected locations on different and often opposite sides of the comet, similar sublimation or activity processes, pit formation or collapses, as well as thermal stresses or fracturing events occurred on multiple areas of the comet, shaping its surface into the appearance we see today.

Dust measurements in the coma of comet 67P/Churyumov-Gerasimenko inbound to the Sun

Critical measurements for understanding accretion and the dust/gas ratio in the solar nebula, where planets were forming 4.5 billion years ago, are being obtained by the GIADA (Grain Impact Analyser and Dust Accumulator) experiment on the European Space Agency's Rosetta spacecraft orbiting comet 67P/Churyumov-Gerasimenko. Between 3.6 and 3.4 astronomical units inbound, GIADA and OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) detected 35 outflowing grains of mass 10(-10) to 10(-7) kilograms, and 48 grains of mass 10(-5) to 10(-2) kilograms, respectively. Combined with gas data from the MIRO (Microwave Instrument for the Rosetta Orbiter) and ROSINA (Rosetta Orbiter Spectrometer for Ion and Neutral Analysis) instruments, we find a dust/gas mass ratio of 4 +/- 2 averaged over the sunlit nucleus surface. A cloud of larger grains also encircles the nucleus in bound orbits from the previous perihelion. The largest orbiting clumps are meter-sized, confirming the dust/gas ratio of 3 inferred at perihelion from models of dust comae and trails.