3D noise power spectrum applied on clinical mdct scanners: effects of reconstruction algorithms and reconstruction filters

The noise power spectrum (NPS) is the reference metric for understanding the noise content in computed tomography (CT) images. To evaluate the noise properties of clinical multidetector (MDCT) scanners, local 2D and 3D NPSs were computed for different acquisition reconstruction parameters.A 64- and a 128-MDCT scanners were employed. Measurements were performed on a water phantom in axial and helical acquisition modes. CT dose index was identical for both installations. Influence of parameters such as the pitch, the reconstruction filter (soft, standard and bone) and the reconstruction algorithm (filtered-back projection (FBP), adaptive statistical iterative reconstruction (ASIR)) were investigated. Images were also reconstructed in the coronal plane using a reformat process. Then 2D and 3D NPS methods were computed.In axial acquisition mode, the 2D axial NPS showed an important magnitude variation as a function of the z-direction when measured at the phantom center. In helical mode, a directional dependency with lobular shape was observed while the magnitude of the NPS was kept constant. Important effects of the reconstruction filter, pitch and reconstruction algorithm were observed on 3D NPS results for both MDCTs. With ASIR, a reduction of the NPS magnitude and a shift of the NPS peak to the low frequency range were visible. 2D coronal NPS obtained from the reformat images was impacted by the interpolation when compared to 2D coronal NPS obtained from 3D measurements.The noise properties of volume measured in last generation MDCTs was studied using local 3D NPS metric. However, impact of the non-stationarity noise effect may need further investigations.

Shoulder pathology: estimating dominant upper-limb activity using 3D kinematic sensors.

Introduction. Quantification of daily upper-limb activity is a key determinant in evaluation of shoulder surgery. For a number of shoulder diseases, problem in performing daily activities have been expressed in terms of upper-limb usage and non-usage. Many instruments measure upper-limb movement but do not focus on the differentiations between the use of left or right shoulder. Several methods have been used to measure it using only accelerometers, pressure sensors or video-based analysis. However, there is no standard or widely used objective measure for upper-limb movement. We report here on an objective method to measure the movement of upper-limb and we examined the use of 3D accelerometers and 3D gyroscopes for that purpose. Methods. We studied 8 subjects with unilateral pathological shoulder (8 rotator cuff disease: 53 years old ± 8) and compared them to 18 control subjects (10 right handed, 8 left handed: 32 years old ± 8, younger than the patient group to be almost sure they don_t have any unrecognized shoulder pathology). The Simple Shoulder Test (SST) and Disabilities of the Arm and Shoulder Score (DASH) questionnaires were completed by each subject. Two modules with 3 miniature capacitive gyroscopes and 3 miniature accelerometers were fixed by a patch on the dorsal side of the distal humerus, and one module with 3 gyroscopes and 3 accelerometers were fixed on the thorax. The subject wore the system during one day (8 hours), at home or wherever he/she went. We used a technique based on the 3D acceleration and the 3D angular velocities from the modules attached on the humerus. Results. As expected, we observed that for the stand and sit postures the right side is more used than the left side for a healthy right-handed person(idem on the left side for a healthy left-handed person). Subjects used their dominant upper-limb 18% more than the non-dominant upper-limb. The measurements on patients in daily life have shown that the patient has used more his non affected and non dominant side during daily activity if the dominant side = affected shoulder. If the dominant side affected shoulder, the difference can be showed only during walking period. Discussion-Conclusion. The technique developed and used allowed the quantification of the difference between dominant and non dominant side, affected and unaffected upper-limb activity. These results were encouraging for future evaluation of patients with shoulder injuries, before and after surgery. The feasibility and patient acceptability of the method using body fixed sensors for ambulatory evaluation of upper limbs kinematics was shown.

Estimating traveltimes and groundwater flow patterns using 3D time-lapse crosshole ERT imaging of electrical resistivity fluctuations induced by infiltrating river water

The infiltration of river water into aquifers is of high relevance to drinking-water production and is a key driver of biogeochemical processes in the hyporheic and riparian zone, but the distribution and quantification of the infiltrating water are difficult to determine using conventional hydrological methods (e.g., borehole logging and tracer tests). By time-lapse inverting crosshole ERT (electrical resistivity tomography) monitoring data, we imaged groundwater flow patterns driven by river water infiltrating a perialpine gravel aquifer in northeastern Switzerland. This was possible because the electrical resistivity of the infiltrating water changed during rainfall-runoff events. Our time-lapse resistivity models indicated rather complex flow patterns as a result of spatially heterogeneous bank filtration and aquifer heterogeneity. The upper part of the aquifer was most affected by the river infiltrate, and the highest groundwater velocities and possible preferential flow occurred at shallow to intermediate depths. Time series of the reconstructed resistivity models matched groundwater electrical resistivity data recorded on borehole loggers in the upper and middle parts of the aquifer, whereas the resistivity models displayed smaller variations and delayed responses with respect to the logging data. in the lower part. This study demonstrated that crosshole ERT monitoring of natural electrical resistivity variations of river infiltrate could be used to image and quantify 3D bank filtration and aquifer dynamics at a high spatial resolution.

Patterns and transitions in substance use among young Swiss men: a latent transition analysis approach

This study investigates the potential stages of drug use. Data from the longitudinal Cohort Study on Substance Use Risk Factors were used (N = 5,116). Drug use (alcohol, tobacco, and 16 illicit drugs) over the previous 12 months was assessed at two time points. Patterns and trajectories of drug use were studied using latent transition analysis (LTA). This study's substantive contributions are twofold. First, the pattern of drug use displayed the well-known sequence of drug involvement (licit drugs to cannabis to other illicit drugs), but with an added distinction between two kinds of illicit drugs ("middle-stage" drugs: uppers, hallucinogens, inhaled drugs; and "final-stage" drugs: heroin, ketamine, GHB/GBL, research chemicals, crystal meth, and spice). Second, subgroup membership was stable over time, as the most likely transition was remaining in the same latent class.

Comparison of 3D segmented gradient-echo and steady-state free precession coronary MRI sequences in patients with coronary artery disease.

OBJECTIVE: Our objective was to compare two state-of-the-art coronary MRI (CMRI) sequences with regard to image quality and diagnostic accuracy for the detection of coronary artery disease (CAD). SUBJECTS AND METHODS: Twenty patients with known CAD were examined with a navigator-gated and corrected free-breathing 3D segmented gradient-echo (turbo field-echo) CMRI sequence and a steady-state free precession sequence (balanced turbo field-echo). CMRI was performed in a transverse plane for the left coronary artery and a double-oblique plane for the right coronary artery system. Subjective image quality (1- to 4-point scale, with 1 indicating excellent quality) and objective image quality parameters were independently determined for both sequences. Sensitivity, specificity, and accuracy for the detection of significant (> or = 50% diameter) coronary artery stenoses were determined as defined in invasive catheter X-ray coronary angiography. RESULTS: Subjective image quality was superior for the balanced turbo field-echo approach (1.8 +/- 0.9 vs 2.3 +/- 1.0 for turbo field-echo; p < 0.001). Vessel sharpness, signal-to-noise ratio, and contrast-to-noise ratio were all superior for the balanced turbo field-echo approach (p < 0.01 for signal-to-noise ratio and contrast-to-noise ratio). Of the 103 segments, 18% of turbo field-echo segments and 9% of balanced turbo field-echo segments had to be excluded from disease evaluation because of insufficient image quality. Sensitivity, specificity, and accuracy for the detection of significant coronary artery stenoses in the evaluated segments were 92%, 67%, 85%, respectively, for turbo field-echo and 82%, 82%, 81%, respectively, for balanced turbo field-echo. CONCLUSION: Balanced turbo field-echo offers improved image quality with significantly fewer nondiagnostic segments when compared with turbo field-echo. For the detection of CAD, both sequences showed comparable accuracy for the visualized segments.

Pore-scale modeling of two-phase flow instabilities in porous media

Les problèmes d'écoulements multiphasiques en média poreux sont d'un grand intérêt pour de nombreuses applications scientifiques et techniques ; comme la séquestration de C02, l'extraction de pétrole et la dépollution des aquifères. La complexité intrinsèque des systèmes multiphasiques et l'hétérogénéité des formations géologiques sur des échelles multiples représentent un challenge majeur pour comprendre et modéliser les déplacements immiscibles dans les milieux poreux. Les descriptions à l'échelle supérieure basées sur la généralisation de l'équation de Darcy sont largement utilisées, mais ces méthodes sont sujettes à limitations pour les écoulements présentant de l'hystérèse. Les avancées récentes en terme de performances computationnelles et le développement de méthodes précises pour caractériser l'espace interstitiel ainsi que la distribution des phases ont favorisé l'utilisation de modèles qui permettent une résolution fine à l'échelle du pore. Ces modèles offrent un aperçu des caractéristiques de l'écoulement qui ne peuvent pas être facilement observées en laboratoire et peuvent être utilisé pour expliquer la différence entre les processus physiques et les modèles à l'échelle macroscopique existants. L'objet premier de la thèse se porte sur la simulation numérique directe : les équations de Navier-Stokes sont résolues dans l'espace interstitiel et la méthode du volume de fluide (VOF) est employée pour suivre l'évolution de l'interface. Dans VOF, la distribution des phases est décrite par une fonction fluide pour l'ensemble du domaine et des conditions aux bords particulières permettent la prise en compte des propriétés de mouillage du milieu poreux. Dans la première partie de la thèse, nous simulons le drainage dans une cellule Hele-Shaw 2D avec des obstacles cylindriques. Nous montrons que l'approche proposée est applicable même pour des ratios de densité et de viscosité très importants et permet de modéliser la transition entre déplacement stable et digitation visqueuse. Nous intéressons ensuite à l'interprétation de la pression capillaire à l'échelle macroscopique. Nous montrons que les techniques basées sur la moyenne spatiale de la pression présentent plusieurs limitations et sont imprécises en présence d'effets visqueux et de piégeage. Au contraire, une définition basée sur l'énergie permet de séparer les contributions capillaires des effets visqueux. La seconde partie de la thèse est consacrée à l'investigation des effets d'inertie associés aux reconfigurations irréversibles du ménisque causé par l'interface des instabilités. Comme prototype pour ces phénomènes, nous étudions d'abord la dynamique d'un ménisque dans un pore angulaire. Nous montrons que, dans un réseau de pores cubiques, les sauts et reconfigurations sont si fréquents que les effets d'inertie mènent à différentes configurations des fluides. A cause de la non-linéarité du problème, la distribution des fluides influence le travail des forces de pression, qui, à son tour, provoque une chute de pression dans la loi de Darcy. Cela suggère que ces phénomènes devraient être pris en compte lorsque que l'on décrit l'écoulement multiphasique en média poreux à l'échelle macroscopique. La dernière partie de la thèse s'attache à démontrer la validité de notre approche par une comparaison avec des expériences en laboratoire : un drainage instable dans un milieu poreux quasi 2D (une cellule Hele-Shaw avec des obstacles cylindriques). Plusieurs simulations sont tournées sous différentes conditions aux bords et en utilisant différents modèles (modèle intégré 2D et modèle 3D) afin de comparer certaines quantités macroscopiques avec les observations au laboratoire correspondantes. Malgré le challenge de modéliser des déplacements instables, où, par définition, de petites perturbations peuvent grandir sans fin, notre approche numérique apporte de résultats satisfaisants pour tous les cas étudiés. - Problems involving multiphase flow in porous media are of great interest in many scientific and engineering applications including Carbon Capture and Storage, oil recovery and groundwater remediation. The intrinsic complexity of multiphase systems and the multi scale heterogeneity of geological formations represent the major challenges to understand and model immiscible displacement in porous media. Upscaled descriptions based on generalization of Darcy's law are widely used, but they are subject to several limitations for flow that exhibit hysteric and history- dependent behaviors. Recent advances in high performance computing and the development of accurate methods to characterize pore space and phase distribution have fostered the use of models that allow sub-pore resolution. These models provide an insight on flow characteristics that cannot be easily achieved by laboratory experiments and can be used to explain the gap between physical processes and existing macro-scale models. We focus on direct numerical simulations: we solve the Navier-Stokes equations for mass and momentum conservation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. In the VOF the distribution of the phases is described by a fluid function (whole-domain formulation) and special boundary conditions account for the wetting properties of the porous medium. In the first part of this thesis we simulate drainage in a 2-D Hele-Shaw cell filled with cylindrical obstacles. We show that the proposed approach can handle very large density and viscosity ratios and it is able to model the transition from stable displacement to viscous fingering. We then focus on the interpretation of the macroscopic capillary pressure showing that pressure average techniques are subject to several limitations and they are not accurate in presence of viscous effects and trapping. On the contrary an energy-based definition allows separating viscous and capillary contributions. In the second part of the thesis we investigate inertia effects associated with abrupt and irreversible reconfigurations of the menisci caused by interface instabilities. As a prototype of these phenomena we first consider the dynamics of a meniscus in an angular pore. We show that in a network of cubic pores, jumps and reconfigurations are so frequent that inertia effects lead to different fluid configurations. Due to the non-linearity of the problem, the distribution of the fluids influences the work done by pressure forces, which is in turn related to the pressure drop in Darcy's law. This suggests that these phenomena should be taken into account when upscaling multiphase flow in porous media. The last part of the thesis is devoted to proving the accuracy of the numerical approach by validation with experiments of unstable primary drainage in a quasi-2D porous medium (i.e., Hele-Shaw cell filled with cylindrical obstacles). We perform simulations under different boundary conditions and using different models (2-D integrated and full 3-D) and we compare several macroscopic quantities with the corresponding experiment. Despite the intrinsic challenges of modeling unstable displacement, where by definition small perturbations can grow without bounds, the numerical method gives satisfactory results for all the cases studied.

Free-breathing whole-heart coronary MRA with 3D radial SSFP and self-navigated image reconstruction.

Respiratory motion is a major source of artifacts in cardiac magnetic resonance imaging (MRI). Free-breathing techniques with pencil-beam navigators efficiently suppress respiratory motion and minimize the need for patient cooperation. However, the correlation between the measured navigator position and the actual position of the heart may be adversely affected by hysteretic effects, navigator position, and temporal delays between the navigators and the image acquisition. In addition, irregular breathing patterns during navigator-gated scanning may result in low scan efficiency and prolonged scan time. The purpose of this study was to develop and implement a self-navigated, free-breathing, whole-heart 3D coronary MRI technique that would overcome these shortcomings and improve the ease-of-use of coronary MRI. A signal synchronous with respiration was extracted directly from the echoes acquired for imaging, and the motion information was used for retrospective, rigid-body, through-plane motion correction. The images obtained from the self-navigated reconstruction were compared with the results from conventional, prospective, pencil-beam navigator tracking. Image quality was improved in phantom studies using self-navigation, while equivalent results were obtained with both techniques in preliminary in vivo studies.

Effects of computing parameters and measurement locations on the estimation of 3D NPS in non-stationary MDCT images.

The goal of this study was to investigate the impact of computing parameters and the location of volumes of interest (VOI) on the calculation of 3D noise power spectrum (NPS) in order to determine an optimal set of computing parameters and propose a robust method for evaluating the noise properties of imaging systems. Noise stationarity in noise volumes acquired with a water phantom on a 128-MDCT and a 320-MDCT scanner were analyzed in the spatial domain in order to define locally stationary VOIs. The influence of the computing parameters in the 3D NPS measurement: the sampling distances bx,y,z and the VOI lengths Lx,y,z, the number of VOIs NVOI and the structured noise were investigated to minimize measurement errors. The effect of the VOI locations on the NPS was also investigated. Results showed that the noise (standard deviation) varies more in the r-direction (phantom radius) than z-direction plane. A 25 × 25 × 40 mm(3) VOI associated with DFOV = 200 mm (Lx,y,z = 64, bx,y = 0.391 mm with 512 × 512 matrix) and a first-order detrending method to reduce structured noise led to an accurate NPS estimation. NPS estimated from off centered small VOIs had a directional dependency contrary to NPS obtained from large VOIs located in the center of the volume or from small VOIs located on a concentric circle. This showed that the VOI size and location play a major role in the determination of NPS when images are not stationary. This study emphasizes the need for consistent measurement methods to assess and compare image quality in CT.

3D culture of postnatal retinal stem cells using self assembling polymers

PURPOSE We have previously shown that retinal stem cells (RSCs) can be isolated from the radial glia population of the newborn mouse retina (Angénieux et al., 2006). These RSCs have a great capacity to renew and to generate a large number of neurons including cells differentiated towards the photoreceptor lineage (Mehri-Soussi et al., 2006). However, recent published results from our lab revealed that such cells have a poor integration and survival rate after grafting. The uncontrolled environment of a retina seems to prevent good integration and survival after grafting in vivo. To bypass this problem, we are evaluating the possibility of generating in vitro a hemi-retinal tissue before transplantation. METHODS RSC were expanded and cells passaged <10 were seeded in a solution containing poly-ethylene-glycol (PEG) polymer based hydrogels crosslinked with peptides that are chosen to be substrates for matrix metalloproteinases. Various doses of cross linkers peptides allowing connections between PEG polymers were tested. Different growth factors were studied to stimulate cell proliferation and differentiation. RESULTS Cells survived only in the presence of EGF and FGF-2 and generated colonies with a sphere shape. No cells migrated within the gel. To improve the migration and the repartition of the cells in the gels, the integrin ligand RGDSP was added into the gel. In the presence of FGF-2 and EGF, newly formed cell clusters appeared by cell proliferation within several days, but again no outspreading of cells was observed. No difference was even seen when the stiffness of the hydrogels or the concentration of the integrin ligand RGDSP were changed. However, our preliminary results show that RSCs still form spheres when laminin is entrapped in the gel, but they started to spread out having a neuronal morphology after around 2 weeks. The neuronal population was assessed by the presence of the neuronal marker b-tubulin-III. This differentiation was achieved after successive steps of stimulations including FGF-2 and EGF, and then only FGF-2. Glial cells were also present. Further characterizations are under process. CONCLUSIONS RSC can be grown in 3D. Preliminary results show that neuronal cell phenotype acquisition can be instructed by exogenous stimulations and factors linked to the gel. Further developments are necessary to form a homogenous tissue containing retinal cells.

Oxidation of methane at the CH4/H2O-(CO2) transition zone in the external part of the Central Alps, Switzerland: Evidence from stable isotope investigations

With the aim of understanding the mechanisms that control the metamorphic transition from the CH4- to the H2O-(CO2)-dominated fluid zone in the Helvetic domain of the Central Alps of Switzerland, fluid inclusions in quartz, illite ``crystallinity'' index, vitrinite reflectance, and the stable isotope compositions of vein and whole rock minerals and fluids trapped in quartz were investigated along four cross-sections. Increasing temperature during prograde metamorphism led to the formation of dry gas by hydrocarbon cracking in the CH4-zone. Fluid immiscibility in the H2O-CH4-(CO2)-NaCl system resulted in cogenetic, CH4- and H2O-dominated fluid inclusions. In the CH4-zone, fluids were trapped at temperatures <= 270 +/- 5 degrees C. The end of the CH4-zone is markedby a sudden increase of CO2 content in the gas phase of fluid inclusions. At temperatures > 270 +/- 5 degrees C, in the H2O-zone, the total amount of volatiles within the fluid decreased below 1 mol% with no immiscibility. This resulted m total homogenization temperatures of H2O-(CO2-CH4)-NaCl inclusions below 180 degrees C. Hydrogen isotope compositions of methane in fluid inclusion have delta D values of less than -100 parts per thousand in the CH4-zone, typical for an origin through cracking of higher hydrocarbons, but where the methane has not equilibrated with the pore water. delta D values of fluid inclusion water are around -40 parts per thousand., in isotopic equilibrium with phyllosilicates of the whole rocks. Within the CH4 to H2O(CO2) transition zone, delta D(H2O) values in fluid inclusions decrease to -130 parts per thousand interpreted to reflect the contribution of deuterium depleted water from methane oxidation. In the H2O-zone, delta D(H2O) values increase again towards an average of -30 parts per thousand which is again consistent with isotopic equilibrium with host-rock phyllosilicates. delta C-13 values of methane in fluid inclusions from the CH4-zone are around -27 parts per thousand in isotopic equilibrium with calcite in veins and whole rocks. The delta C-13(CH4) values decrease to less than -35 parts per thousand at the transition to the H2O-zone and are no longer in equilibrium with the carbonates in the whole rocks. delta C-13 values of CO, are variable but too low to be in equilibrium with the wall rock fluids, compatible with a contribution of CO2 from closed system oxidation of methane. Differences in isotopic composition between host-rock and Alpine fissure carbonate are generally small, suggesting that the amount of CO2 produced by oxidation of methane was small compared to the C-budget in the rocks and local pore fluids were buffered by the wall rocks during precipitation of calcite within the fissures. (c) 2006 Elsevier B.V. All rights reserved.

Sea-level changes and ammonite faunal turnover during the Lias/Dogger transition in the western Tethys

The aim of this paper was to investigate the possible connections between ammonite faunal turnover and the eustatic events recorded in Tethyan sequences during the middle Toarcian/early Bajocian time interval. For this we have analysed the biostratigraphic ranges, at the subzone level, of approximately 600 ammonite species belonging to 160 genera from several selected sections of the western Tethys (Mediterranean and Submediterranean provinces). The analysis of taxon ranges enabled us to plot curves for ammonite faunal turnovers, inter-subzonal distance, and diversity. Comparing the mentioned curves with Tethyan sequences [Hardenbol et al., 19981, we find that sea-level changes correlate well with origination and extinction events and faunal diversity. Most of the faunal turnovers correlate with stratigraphic events. Extinction events with their corresponding decrease in diversity correlate with regressive intervals and with major or minor sequence boundaries. Origination events and their corresponding increase in diversity were clearly connected with transgressions in Tethyan sequences. In several cases, the major sequence boundary and the subsequent transgressive phase correlate with major ammonite faunal turnover, whereas minor or medium sequence boundaries generally gave rise to minor or medium turnovers.