Effect of irradiation distance on image contrast in epi-optoacoustic imaging of human volunteers

In combined clinical optoacoustic (OA) and ultrasound (US) imaging, epi-mode irradiation and detection integrated into one single probe offers flexible imaging of the human body. The imaging depth in epi-illumination is, however, strongly affected by clutter. As shown in previous phantom experiments, the location of irradiation plays an important role in clutter generation. We investigated the influence of the irradiation geometry on the local image contrast of clinical images, by varying the separation distance between the irradiated area and the acoustic imaging plane of a linear ultrasound transducer in an automated scanning setup. The results for different volunteers show that the image contrast can be enhanced on average by 25% and locally by more than a factor of two, when the irradiated area is slightly separated from the probe. Our findings have an important impact on the design of future optoacoustic probes for clinical application.

Image Space Adaptive Rendering

In this thesis, we develop an adaptive framework for Monte Carlo rendering, and more specifically for Monte Carlo Path Tracing (MCPT) and its derivatives. MCPT is attractive because it can handle a wide variety of light transport effects, such as depth of field, motion blur, indirect illumination, participating media, and others, in an elegant and unified framework. However, MCPT is a sampling-based approach, and is only guaranteed to converge in the limit, as the sampling rate grows to infinity. At finite sampling rates, MCPT renderings are often plagued by noise artifacts that can be visually distracting. The adaptive framework developed in this thesis leverages two core strategies to address noise artifacts in renderings: adaptive sampling and adaptive reconstruction. Adaptive sampling consists in increasing the sampling rate on a per pixel basis, to ensure that each pixel value is below a predefined error threshold. Adaptive reconstruction leverages the available samples on a per pixel basis, in an attempt to have an optimal trade-off between minimizing the residual noise artifacts and preserving the edges in the image. In our framework, we greedily minimize the relative Mean Squared Error (rMSE) of the rendering by iterating over sampling and reconstruction steps. Given an initial set of samples, the reconstruction step aims at producing the rendering with the lowest rMSE on a per pixel basis, and the next sampling step then further reduces the rMSE by distributing additional samples according to the magnitude of the residual rMSE of the reconstruction. This iterative approach tightly couples the adaptive sampling and adaptive reconstruction strategies, by ensuring that we only sample densely regions of the image where adaptive reconstruction cannot properly resolve the noise. In a first implementation of our framework, we demonstrate the usefulness of our greedy error minimization using a simple reconstruction scheme leveraging a filterbank of isotropic Gaussian filters. In a second implementation, we integrate a powerful edge aware filter that can adapt to the anisotropy of the image. Finally, in a third implementation, we leverage auxiliary feature buffers that encode scene information (such as surface normals, position, or texture), to improve the robustness of the reconstruction in the presence of strong noise.

Optimized aperture shapes for depth estimation

The finite depth of field of a real camera can be used to estimate the depth structure of a scene. The distance of an object from the plane in focus determines the defocus blur size. The shape of the blur depends on the shape of the aperture. The blur shape can be designed by masking the main lens aperture. In fact, aperture shapes different from the standard circular aperture give improved accuracy of depth estimation from defocus blur. We introduce an intuitive criterion to design aperture patterns for depth from defocus. The criterion is independent of a specific depth estimation algorithm. We formulate our design criterion by imposing constraints directly in the data domain and optimize the amount of depth information carried by blurred images. Our criterion is a quadratic function of the aperture transmission values. As such, it can be numerically evaluated to estimate optimized aperture patterns quickly. The proposed mask optimization procedure is applicable to different depth estimation scenarios. We use it for depth estimation from two images with different focus settings, for depth estimation from two images with different aperture shapes as well as for depth estimation from a single coded aperture image. In this work we show masks obtained with this new evaluation criterion and test their depth discrimination capability using a state-of-the-art depth estimation algorithm.

Enhanced-depth optical coherence tomography for imaging horizontal rectus muscles in Graves' orbitopathy.

PURPOSE Graves' orbitopathy (GO) is an extraocular eye disease with symptoms ranging from minor discomfort from dry eyes to strabismus and visual loss. One of the hallmarks of active GO is visible hyperemia at the insertion of the extraocular muscles. The aim of the present study was to evaluate the use of enhanced-depth imaging spectral domain anterior segment optical coherence tomography (EDI SD AS-OCT) for detecting pathological changes in horizontal recti muscles of patients with GO. METHODS Prospective cross sectional study of 27 eyes. Only women were included. EDI AS-OCT was used to measure the thickness of the tendons of the horizontal recti muscles in a predefined area in patients with GO and healthy controls. RESULTS EDI AS-OCT was able to image the tendons of the horizontal recti muscles in both healthy controls and patients suffering from GO. The mean thickness of the medial rectus muscle (MR) tendon was 256.4 μm [±17.13 μm standard deviation (SD)] in the GO group and, therefore, significantly thicker (p = 0.046) than in the healthy group which had a mean thickness of 214.7 μm (±5.516 μm SD). There was no significant difference in the mean thickness of the tendon of the lateral recti muscles (LRs) between these groups. CONCLUSION This is the first report showing that EDI AS-OCT is suitable to detect swelling at the insertion site of the MR muscle in GO. MR tendon thickness may be a useful parameter to monitor activity in these patients.

Design, Calibration, and Evaluation of Depth-of-Interaction-Capable PET Detector Modules

High-resolution, small-bore PET systems suffer from a tradeoff between system sensitivity, and image quality degradation. In these systems long crystals allow mispositioning of the line of response due to parallax error and this mispositioning causes resolution blurring, but long crystals are necessary for high system sensitivity. One means to allow long crystals without introducing parallax errors is to determine the depth of interaction (DOI) of the gamma ray interaction within the detector module. While DOI has been investigated previously, newly available solid state photomultipliers (SSPMs) well-suited to PET applications and allow new modules for investigation. Depth of interaction in full modules is a relatively new field, and so even if high performance DOI capable modules were available, the appropriate means to characterize and calibrate the modules are not. This work presents an investigation of DOI capable arrays and techniques for characterizing and calibrating those modules. The methods introduced here accurately and reliably characterize and calibrate energy, timing, and event interaction positioning. Additionally presented is a characterization of the spatial resolution of DOI capable modules and a measurement of DOI effects for different angles between detector modules. These arrays have been built into a prototype PET system that delivers better than 2.0 mm resolution with a single-sided-stopping-power in excess of 95% for 511 keV g's. The noise properties of SSPMs scale with the active area of the detector face, and so the best signal-to-noise ratio is possible with parallel readout of each SSPM photodetector pixel rather than multiplexing signals together. This work additionally investigates several algorithms for improving timing performance using timing information from multiple SSPM pixels when light is distributed among several photodetectors.

Faecal pellet production of copepoda collected in water depth up to 30 meters in the eastern Mediterranean Sea in Oktober 2008 during SES_GR2

Mesozooplankton is collected by vertical tows within the Black sea water body mass layer in the NE Aegean, using a WP-2 200 µm net equipped with a large non-filtering cod-end (10 l). Macrozooplankton organisms are removed using a 2000 µm net. A few unsorted animals (approximately 100) are placed inside several glass beaker of 250 ml filled with GF/F or 0.2 µm Nucleopore filtered seawater and with a 100 µm net placed 1 cm above the beaker bottom. Beakers are then placed in an incubator at natural light and maintaining the in situ temperature. After 1 hour pellets are separated from animals and placed in separated flasks and preserved with formalin. Pellets are counted and measured using an inverted microscope. Animals are scanned and counted using an image analysis system. Carbon- Specific faecal pellet production is calculated from a) faecal pellet production, b) individual carbon: Animals are scanned and their body area is measured using an image analysis system. Body volume is then calculated as an ellipsoid using the major and minor axis of an ellipse of same area as the body. Individual carbon is calculated from a carbon- total body volume of organisms (relationship obtained for the Mediterranean Sea by Alcaraz et al. (2003) divided by the total number of individuals scanned and c) faecal pellet carbon: Faecal pellet length and width is measured using an inverted microscope. Faecal pellet volume is calculated from length and width assuming cylindrical shape. Conversion of faecal pellet volume to carbon is done using values obtained in the Mediterranean from: a) faecal pellet density 1,29 g cm**3 (or pg µm**3) from Komar et al. (1981); b) faecal pellet DW/WW=0,23 from Elder and Fowler (1977) and c) faecal pellet C%DW=25,5 Marty et al. (1994).

Faecal pellet production of copepoda collected in water depth up to 100 meters in the eastern Mediterranean Sea in March and April 2008 during SES_GR1

The SES_UNLUATA_GR1-Mesozooplankton faecal pellet production rates dataset is based on samples taken during March and April 2008 in the Northern Libyan Sea, Southern Aegean Sea and in the North-Eastern Aegean Sea. Mesozooplankton is collected by vertical tows within the 0-100 m layer or within the Black sea water body mass layer in the case of the NE Aegean, using a WP-2 200 µm net equipped with a large non-filtering cod-end (10 l). Macrozooplankton organisms are removed using a 2000 µm net. A few unsorted animals (approximately 100) are placed inside several glass beaker of 250 ml filled with GF/F or 0.2 µm Nucleopore filtered seawater and with a 100 µm net placed 1 cm above the beaker bottom. Beakers are then placed in an incubator at natural light and maintaining the in situ temperature. After 1 hour pellets are separated from animals and placed in separated flasks and preserved with formalin. Pellets and are counted and measured using an inverted microscope. Animals are scanned and counted using an image analysis system. Carbon- Specific faecal pellet production is calculated from a) faecal pellet production, b) individual carbon: Animals are scanned and their body area is measured using an image analysis system. Body volume is then calculated as an ellipsoid using the major and minor axis of an ellipse of same area as the body. Individual carbon is calculated from a carbon- total body volume of organisms (relationship obtained for the Mediterranean Sea by Alcaraz et al. (2003) divided by the total number of individuals scanned and c) faecal pellet carbon: Faecal pellet length and width is measured using an inverted microscope. Faecal pellet volume is calculated from length and width assuming cylindrical shape. Conversion of faecal pellet volume to carbon is done using values obtained in the Mediterranean from: a) faecal pellet density 1,29 g cm**3 (or pg µm**3) from Komar et al. (1981); b) faecal pellet DW/WW=0,23 from Elder and Fowler (1977) and c) faecal pellet C%DW=25,5 Marty et al. (1994).

Faecal pellet production of copepoda collected at different water depth in the eastern Mediterranean Sea during SES_GR2

The SES_GR2-Mesozooplankton faecal pellet production rates dataset is based on samples taken during August and September 2008 in the Northern Libyan Sea, Southern Aegean Sea and the North-Eastern Aegean Sea. Mesozooplankton is collected by vertical tows within the 0-100 m layer or within the Black sea water body mass layer in the case of the NE Aegean, using a WP-2 200 µm net equipped with a large non-filtering cod-end (10 l). Macrozooplankton organisms are removed using a 2000 µm net. A few unsorted animals (approximately 100) are placed inside several glass beaker of 250 ml filled with GF/F or 0.2 µm Nucleopore filtered seawater and with a 100 µm net placed 1 cm above the beaker bottom. Beakers are then placed in an incubator at natural light and maintaining the in situ temperature. After 1 hour pellets are separated from animals and placed in separated flasks and preserved with formalin. Pellets are counted and measured using an inverted microscope. Animals are scanned and counted using an image analysis system. Carbon- Specific faecal pellet production is calculated from a) faecal pellet production, b) individual carbon: Animals are scanned and their body area is measured using an image analysis system. Body volume is then calculated as an ellipsoid using the major and minor axis of an ellipse of same area as the body. Individual carbon is calculated from a carbon- total body volume of organisms (relationship obtained for the Mediterranean Sea by Alcaraz et al. (2003) divided by the total number of individuals scanned and c) faecal pellet carbon: Faecal pellet length and width is measured using an inverted microscope. Faecal pellet volume is calculated from length and width assuming cylindrical shape. Conversion of faecal pellet volume to carbon is done using values obtained in the Mediterranean from: a) faecal pellet density 1,29 g cm**3 (or pg µm**3) from Komar et al. (1981); b) faecal pellet DW/WW=0,23 from Elder and Fowler (1977) and c) faecal pellet C%DW=25,5 Marty et al. (1994).

Faecal pellet production of copepoda collected in water depth up to 20 meters in the eastern Mediterranean Sea in March and April 2008 during SES_GR1

The SES_GR1-Mesozooplankton faecal pellet production rates dataset is based on samples taken during April 2008 in the North-Eastern Aegean Sea. Mesozooplankton is collected by vertical tows within the Black sea water body mass layer in the NE Aegean, using a WP-2 200 µm net equipped with a large non-filtering cod-end (10 l). Macrozooplankton organisms are removed using a 2000 µm net. A few unsorted animals (approximately 100) are placed inside several glass beaker of 250 ml filled with GF/F or 0.2 µm Nucleopore filtered seawater and with a 100 µm net placed 1 cm above the beaker bottom. Beakers are then placed in an incubator at natural light and maintaining the in situ temperature. After 1 hour pellets are separated from animals and placed in separated flasks and preserved with formalin. Pellets are counted and measured using an inverted microscope. Animals are scanned and counted using an image analysis system. Carbon- Specific faecal pellet production is calculated from a) faecal pellet production, b) individual carbon: Animals are scanned and their body area is measured using an image analysis system. Body volume is then calculated as an ellipsoid using the major and minor axis of an ellipse of same area as the body. Individual carbon is calculated from a carbon- total body volume of organisms (relationship obtained for the Mediterranean Sea by Alcaraz et al. (2003) divided by the total number of individuals scanned and c) faecal pellet carbon: Faecal pellet length and width is measured using an inverted microscope. Faecal pellet volume is calculated from length and width assuming cylindrical shape. Conversion of faecal pellet volume to carbon is done using values obtained in the Mediterranean from: a) faecal pellet density 1,29 g cm**3 (or pg µm**3) from Komar et al. (1981); b) faecal pellet DW/WW=0,23 from Elder and Fowler (1977) and c) faecal pellet C%DW=25,5 Marty et al. (1994).