Trabecular bone score (TBS) the new parameter of 2D texture analysis for the evaluation of 3D bone micro architecture status

X-ray is a technology that is used for numerous applications in the medical field. The process of X-ray projection gives a 2-dimension (2D) grey-level texture from a 3- dimension (3D) object. Until now no clear demonstration or correlation has positioned the 2D texture analysis as a valid indirect evaluation of the 3D microarchitecture. TBS is a new texture parameter based on the measure of the experimental variogram. TBS evaluates the variation between 2D image grey-levels. The aim of this study was to evaluate existing correlations between 3D bone microarchitecture parameters - evaluated from μCT reconstructions - and the TBS value, calculated on 2D projected images. 30 dried human cadaveric vertebrae were acquired on a micro-scanner (eXplorer Locus, GE) at isotropic resolution of 93 μm. 3D vertebral body models were used. The following 3D microarchitecture parameters were used: Bone volume fraction (BV/TV), Trabecular thickness (TbTh), trabecular space (TbSp), trabecular number (TbN) and connectivity density (ConnD). 3D/2D projections has been done by taking into account the Beer-Lambert Law at X-ray energy of 50, 100, 150 KeV. TBS was assessed on 2D projected images. Correlations between TBS and the 3D microarchitecture parameters were evaluated using a linear regression analysis. Paired T-test is used to assess the X-ray energy effects on TBS. Multiple linear regressions (backward) were used to evaluate relationships between TBS and 3D microarchitecture parameters using a bootstrap process. BV/TV of the sample ranged from 18.5 to 37.6% with an average value at 28.8%. Correlations' analysis showedthat TBSwere strongly correlatedwith ConnD(0.856≤r≤0.862; p<0.001),with TbN (0.805≤r≤0.810; p<0.001) and negatively with TbSp (−0.714≤r≤−0.726; p<0.001), regardless X-ray energy. Results show that lower TBS values are related to "degraded" microarchitecture, with low ConnD, low TbN and a high TbSp. The opposite is also true. X-ray energy has no effect onTBS neither on the correlations betweenTBS and the 3Dmicroarchitecture parameters. In this study, we demonstrated that TBS was significantly correlated with 3D microarchitecture parameters ConnD and TbN, and negatively with TbSp, no matter what X-ray energy has been used. This article is part of a Special Issue entitled ECTS 2011. Disclosure of interest: None declared.

3D Automated Lung Nodule Segmentation in HRCT

A fully-automated 3D image analysis method is proposed to segment lung nodules in HRCT. A specific gray-level mathematical morphology operator, the SMDC-connection cost, acting in the 3D space of the thorax volume is defined in order to discriminate lung nodules from other dense (vascular) structures. Applied to clinical data concerning patients with pulmonary carcinoma, the proposed method detects isolated, juxtavascular and peripheral nodules with sizes ranging from 2 to 20 mm diameter. The segmentation accuracy was objectively evaluated on real and simulated nodules. The method showed a sensitivity and a specificity ranging from 85% to 97% and from 90% to 98%, respectively.

JULIDE: a software tool for 3D reconstruction and statistical analysis of autoradiographic mouse brain sections.

In this article we introduce JULIDE, a software toolkit developed to perform the 3D reconstruction, intensity normalization, volume standardization by 3D image registration and voxel-wise statistical analysis of autoradiographs of mouse brain sections. This software tool has been developed in the open-source ITK software framework and is freely available under a GPL license. The article presents the complete image processing chain from raw data acquisition to 3D statistical group analysis. Results of the group comparison in the context of a study on spatial learning are shown as an illustration of the data that can be obtained with this tool.

Functional MRI Mapping of the Human Auditory Cortex

Mapping the human auditory cortex with standard functional imaging techniques is difficult because of its small size and angular position along the Sylvian fissure. As a result, the exact number and location of auditory cortex areas in the human remains unknown. In a first experiment, we measured the two largest tonotopic areas of primary auditory cortex (PAC, Al and R) using high-resolution functional MRI at 7 Tesla relative to the underlying anatomy of Heschl's gyrus (HG). The data reveals a clear anatomical- functional relationship that indicates the location of PAC across the range of common morphological variants of HG (single gyri, partial duplication and complete duplication). Human PAC tonotopic areas are oriented along an oblique posterior-to-anterior axis with mirror-symmetric frequency gradients perpendicular to HG, as in the macaque. In a second experiment, we tested whether these primary frequency-tuned units were modulated by selective attention to preferred vs. non-preferred sound frequencies in the dynamic manner needed to account for human listening abilities in noisy environments, such as cocktail parties or busy streets. We used a dual-stream selective attention experiment where subjects attended to one of two competing tonal streams presented simultaneously to different ears. Attention to low-frequency tones (250 Hz) enhanced neural responses within low-frequency-tuned voxels relative to high (4000 Hz), and vice versa when at-tention switched from high to low. Human PAC is able to tune into attended frequency channels and can switch frequencies on demand, like a radio. In a third experiment, we investigated repetition suppression effects to environmental sounds within primary and non-primary early-stage auditory areas, identified with the tonotopic mapping design. Repeated presentations of sounds from the same sources, as compared to different sources, gave repetition suppression effects within posterior and medial non-primary areas of the right hemisphere, reflecting their potential involvement in semantic representations. These three studies were conducted at 7 Tesla with high-resolution imaging. However, 7 Tesla scanners are, for the moment, not yet used for clinical diagnosis and mostly reside in institutions external to hospitals. Thus, hospital-based clinical functional and structural studies are mainly performed using lower field systems (1.5 or 3 Tesla). In a fourth experiment, we acquired tonotopic maps at 3 and 7 Tesla and evaluated the consistency of a tonotopic mapping paradigm between scanners. Mirror-symmetric gradients within PAC were highly similar at 7 and 3 Tesla across renderings at different spatial resolutions. We concluded that the tonotopic mapping paradigm is robust and suitable for definition of primary tonotopic areas, also at 3 Tesla. Finally, in a fifth study, we considered whether focal brain lesions alter tonotopic representations in the intact ipsi- and contralesional primary auditory cortex in three patients with hemispheric or cerebellar lesions, without and with auditory complaints. We found evidence for tonotopic reorganisation at the level of the primary auditory cortex in cases of brain lesions independently of auditory complaints. Overall, these results reflect a certain degree of plasticity within primary auditory cortex in different populations of subjects, assessed at different field strengths. - La cartographie du cortex auditif chez l'humain est difficile à réaliser avec des techniques d'imagerie fonctionnelle standard, étant donné sa petite taille et position angulaire le long de la fissure sylvienne. En conséquence, le nombre et l'emplacement exacts des différentes aires du cortex auditif restent inconnus chez l'homme. Lors d'une première expérience, nous avons mesuré, avec de l'imagerie par résonance magnétique à haute intensité (IRMf à 7 Tesla) chez des sujets humains sains, deux larges aires au sein du cortex auditif primaire (PAC; Al et R) avec une représentation spécifique des fréquences pures préférées - ou tonotopie. Nos résultats ont démontré une relation anatomico- fonctionnelle qui définit clairement la position du PAC à travers toutes les variantes du gyrus d'Heschl's (HG). Les aires tonotopiques du PAC humain sont orientées le long d'un axe postéro-antérieur oblique avec des gradients de fréquences spécifiques perpendiculaires à HG, d'une manière similaire à celles mesurées chez le singe. Dans une deuxième expérience, nous avons testé si ces aires primaires pouvaient être modulées, de façon dynamique, par une attention sélective pour des fréquences préférées par rapport à celles non-préférées. Cette modulation est primordiale lors d'interactions sociales chez l'humain en présence de bruits distracteurs tels que d'autres discussions ou un environnement sonore nuisible (comme par exemple, dans la circulation routière). Dans cette étude, nous avons utilisé une expérience d'attention sélective où le sujet devait être attentif à une des deux voies sonores présentées simultanément à chaque oreille. Lorsque le sujet portait était attentif aux sons de basses fréquences (250 Hz), la réponse neuronale relative à ces fréquences augmentait par rapport à celle des hautes fréquences (4000 Hz), et vice versa lorsque l'attention passait des hautes aux basses fréquences. De ce fait, nous pouvons dire que PAC est capable de focaliser sur la fréquence attendue et de changer de canal selon la demande, comme une radio. Lors d'une troisième expérience, nous avons étudié les effets de suppression due à la répétition de sons environnementaux dans les aires auditives primaires et non-primaires, d'abord identifiées via le protocole de la première étude. La présentation répétée de sons provenant de la même source sonore, par rapport à de sons de différentes sources sonores, a induit un effet de suppression dans les aires postérieures et médiales auditives non-primaires de l'hémisphère droite, reflétant une implication de ces aires dans la représentation de la catégorie sémantique. Ces trois études ont été réalisées avec de l'imagerie à haute résolution à 7 Tesla. Cependant, les scanners 7 Tesla ne sont pour le moment utilisés que pour de la recherche fondamentale, principalement dans des institutions externes, parfois proches du patient mais pas directement à son chevet. L'imagerie fonctionnelle et structurelle clinique se fait actuellement principalement avec des infrastructures cliniques à 1.5 ou 3 Tesla. Dans le cadre dune quatrième expérience, nous avons avons évalués la cohérence du paradigme de cartographie tonotopique à travers différents scanners (3 et 7 Tesla) chez les mêmes sujets. Nos résultats démontrent des gradients de fréquences définissant PAC très similaires à 3 et 7 Tesla. De ce fait, notre paradigme de définition des aires primaires auditives est robuste et applicable cliniquement. Finalement, nous avons évalués l'impact de lésions focales sur les représentations tonotopiques des aires auditives primaires des hémisphères intactes contralésionales et ipsilésionales chez trois patients avec des lésions hémisphériques ou cérébélleuses avec ou sans plaintes auditives. Nous avons trouvé l'évidence d'une certaine réorganisation des représentations topographiques au niveau de PAC dans le cas de lésions cérébrales indépendamment des plaintes auditives. En conclusion, nos résultats démontrent une certaine plasticité du cortex auditif primaire avec différentes populations de sujets et différents champs magnétiques.

Coronary MR imaging using free-breathing 3D steady-state free precession with radial k-space sampling.

PURPOSE: To investigate the potential of free-breathing 3D steady-state free precession (SSFP) imaging with radial k-space sampling for coronary MR-angiography (MRA), coronary projection MR-angiography and coronary vessel wall imaging. MATERIALS AND METHODS: A navigator-gated free-breathing T2-prepared 3D SSFP sequence (TR = 6.1 ms, TE = 3.0 ms, flip angle = 120 degrees, field-of-view = 360 mm(2)) with radial k-space sampling (384 radials) was implemented for coronary MRA. For projection coronary MRA, this sequence was combined with a 2D selective aortic spin tagging pulse. Coronary vessel wall imaging was performed using a high-resolution inversion-recovery black-blood 3D radial SSFP sequence (384 radials, TR = 5.3 ms, TE = 2.7 ms, flip angle = 55 degrees, reconstructed resolution 0.35 x 0.35 x 1.2 mm(3)) and a local re-inversion pulse. Six healthy volunteers (two for each sequence) were investigated. Motion artifact level was assessed by two radiologists. Results: In coronary MRA, the coronary lumen was displayed with a high signal and high contrast to the surrounding lumen. Projection coronary MRA demonstrated selective visualization of the coronary lumen while surrounding tissue was almost completely suppressed. In coronary vessel wall imaging, the vessel wall was displayed with a high signal when compared to the blood pool and the surrounding tissue. No visible motion artifacts were seen. Conclusion: 3D radial SSFP imaging enables coronary MRA, coronary projection MRA and coronary vessel wall imaging with a low motion artifact level.

Assessment of myocardial edema and area-at-risk in acute myocardial infarction by CMR: Evaluation of a novel T2-mapping method

Objectives. The goal of this study is to evaluate a T2-mapping sequence by: (i) measuring the reproducibility intra- and inter-observer variability in healthy volunteers in two separate scanning session with a T2 reference phantom; (2) measuring the mean T2 relaxation times by T2-mapping in infarcted myocardium in patients with subacute MI and compare it with patient's the gold standard X-ray coronary angiography and healthy volunteers results. Background. Myocardial edema is a consequence of an inflammation of the tissue, as seen in myocardial infarct (MI). It can be visualized by cardiovascular magnetic resonance (CMR) imaging using the T2 relaxation time. T2-mapping is a quantitative methodology that has the potential to address the limitation of the conventional T2-weighted (T2W) imaging. Methods. The T2-mapping protocol used for all MRI scans consisted in a radial gradient echo acquisition with a lung-liver navigator for free-breathing acquisition and affine image registration. Mid-basal short axis slices were acquired.T2-maps analyses: 2 observers semi- automatically segmented the left ventricle in 6 segments accordingly to the AHA standards. 8 healthy volunteers (age: 27 ± 4 years; 62.5% male) were scanned in 2 separate sessions. 17 patients (age : 61.9 ± 13.9 years; 82.4% male) with subacute STEMI (70.6%) and NSTEMI underwent a T2-mapping scanning session. Results. In healthy volunteers, the mean inter- and intra-observer variability over the entire short axis slice (segment 1 to 6) was 0.1 ms (95% confidence interval (CI): -0.4 to 0.5, p = 0.62) and 0.2 ms (95% CI: -2.8 to 3.2, p = 0.94, respectively. T2 relaxation time measurements with and without the correction of the phantom yielded an average difference of 3.0 ± 1.1 % and 3.1 ± 2.1 % (p = 0.828), respectively. In patients, the inter-observer variability in the entire short axis slice (S1-S6), was 0.3 ms (95% CI: -1.8 to 2.4, p = 0.85). Edema location as determined through the T2-mapping and the coronary artery occlusion as determined on X-ray coronary angiography correlated in 78.6%, but only in 60% in apical infarcts. All except one of the maximal T2 values in infarct patients were greater than the upper limit of the 95% confidence interval for normal myocardium. Conclusions. The T2-mapping methodology is accurate in detecting infarcted, i.e. edematous tissue in patients with subacute infarcts. This study further demonstrated that this T2-mapping technique is reproducible and robust enough to be used on a segmental basis for edema detection without the need of a phantom to yield a T2 correction factor. This new quantitative T2-mapping technique is promising and is likely to allow for serial follow-up studies in patients to improve our knowledge on infarct pathophysiology, on infarct healing, and for the assessment of novel treatment strategies for acute infarctions.

SwissTargetPrediction: a web server for target prediction of bioactive small molecules.

Bioactive small molecules, such as drugs or metabolites, bind to proteins or other macro-molecular targets to modulate their activity, which in turn results in the observed phenotypic effects. For this reason, mapping the targets of bioactive small molecules is a key step toward unraveling the molecular mechanisms underlying their bioactivity and predicting potential side effects or cross-reactivity. Recently, large datasets of protein-small molecule interactions have become available, providing a unique source of information for the development of knowledge-based approaches to computationally identify new targets for uncharacterized molecules or secondary targets for known molecules. Here, we introduce SwissTargetPrediction, a web server to accurately predict the targets of bioactive molecules based on a combination of 2D and 3D similarity measures with known ligands. Predictions can be carried out in five different organisms, and mapping predictions by homology within and between different species is enabled for close paralogs and orthologs. SwissTargetPrediction is accessible free of charge and without login requirement at http://www.swisstargetprediction.ch.

Functional calibration procedure for 3D knee joint angle description using inertial sensors.

Measurement of three-dimensional (3D) knee joint angle outside a laboratory is of benefit in clinical examination and therapeutic treatment comparison. Although several motion capture devices exist, there is a need for an ambulatory system that could be used in routine practice. Up-to-date, inertial measurement units (IMUs) have proven to be suitable for unconstrained measurement of knee joint differential orientation. Nevertheless, this differential orientation should be converted into three reliable and clinically interpretable angles. Thus, the aim of this study was to propose a new calibration procedure adapted for the joint coordinate system (JCS), which required only IMUs data. The repeatability of the calibration procedure, as well as the errors in the measurement of 3D knee angle during gait in comparison to a reference system were assessed on eight healthy subjects. The new procedure relying on active and passive movements reported a high repeatability of the mean values (offset<1 degrees) and angular patterns (SD<0.3 degrees and CMC>0.9). In comparison to the reference system, this functional procedure showed high precision (SD<2 degrees and CC>0.75) and moderate accuracy (between 4.0 degrees and 8.1 degrees) for the three knee angle. The combination of the inertial-based system with the functional calibration procedure proposed here resulted in a promising tool for the measurement of 3D knee joint angle. Moreover, this method could be adapted to measure other complex joint, such as ankle or elbow.

Detailed 3D seismic imaging of a fault zone beneath Lake Geneva, Switzerland

An efficient high-resolution (HR) three-dimensional (3D) seismic reflection system for small-scale targets in lacustrine settings was developed. In Lake Geneva, near the city of Lausanne, Switzerland, the offshore extension of a complex fault zone well mapped on land was chosen for testing our system. A preliminary two-dimensional seismic survey indicated structures that include a thin (<40 m) layer of subhorizontal Quaternary sediments that unconformably overlie south-east-dipping Tertiary Molasse beds and a major fault zone (Paudeze Fault Zone) that separates Plateau and Subalpine Molasse (SM) units. A 3D survey was conducted over this test site using a newly developed three-streamer system. It provided high-quality data with a penetration to depths of 300 m below the water bottom of non-aliased signal for dips up to 30degrees and with a maximum vertical resolution of 1.1 m. The data were subjected to a conventional 3D processing sequence that included post-stack time migration. Tests with 3D pre-stack depth migration showed that such techniques can be applied to HR seismic surveys. Delineation of several horizons and fault surfaces reveals the potential for small-scale geologic and tectonic interpretation in three dimensions. Five major seismic facies and their detailed 3D geometries can be distinguished. Three fault surfaces and the top of a molasse surface were mapped in 3D. Analysis of the geometry of these surfaces and their relative orientation suggests that pre-existing structures within the Plateau Molasse (PM) unit influenced later faulting between the Plateau and SM. In particular, a change in strike of the PM bed dip may indicate a fold formed by a regional stress regime, the orientation of which was different from the one responsible for the creation of the Paudeze Fault Zone. This structure might have later influenced the local stress regime and caused the curved shape of the Paudeze Fault in our surveyed area.

Debris flow susceptibility mapping at a regional scale along the National Road N7, Argentina

A first assessment of debris flow susceptibility at a large scale was performed along the National Road N7, Argentina. Numerous catchments are prone to debris flows and likely to endanger the road-users. A 1:50,000 susceptibility map was created. The use of a DEM (grid 30 m) associated to three complementary criteria (slope, contributing area, curvature) allowed the identification of potential source areas. The debris flow spreading was estimated using a process- and GISbased model (Flow-R) based on basic probabilistic and energy calculations. The best-fit values for the coefficient of friction and the mass-to-drag ratio of the PCM model were found to be ? = 0.02 and M/D = 180 and the resulting propagation on one of the calibration site was validated using the Coulomb friction model. The results are realistic and will be useful to determine which areas need to be prioritized for detailed studies.

Evolution of a debris flow channel monitored using a 3D terrestrial laser scanner

Like numerous torrents in mountainous regions, the Illgraben creek (canton of Wallis, SW Switzerland) produces almost every year several debris flows. The total area of the active catchment is only 4.7 km², but large events ranging from 50'000 to 400'000 m³ are common (Zimmermann 2000). Consequently, the pathway of the main channel often changes suddenly. One single event can for instance fill the whole river bed and dig new several-meters-deep channels somewhere else (Bardou et al. 2003). The quantification of both, the rhythm and the magnitude of these changes, is very important to assess the variability of the bed's cross section and long profile. These parameters are indispensable for numerical modelling, as they should be considered as initial conditions. To monitor the channel evolution an Optech ILRIS 3D terrestrial laser scanner (LIDAR) was used. LIDAR permits to make a complete high precision 3D model of the channel and its surroundings by scanning it from different view points. The 3D data are treated and interpreted with the software Polyworks from Innovmetric Software Inc. Sequential 3D models allow for the determination of the variation in the bed's cross section and long profile. These data will afterwards be used to quantify the erosion and the deposition in the torrent reaches. To complete the chronological evolution of the landforms, precise digital terrain models, obtained by high resolution photogrammetry based on old aerial photographs, will be used. A 500 m long section of the Illgraben channel was scanned on 18th of August 2005 and on 7th of April 2006. These two data sets permit identifying the changes of the channel that occurred during the winter season. An upcoming scanning campaign in September 2006 will allow for the determination of the changes during this summer. Preliminary results show huge variations in the pathway of the Illgraben channel, as well as important vertical and lateral erosion of the river bed. Here we present the results of a river bank on the left (north-western) flank of the channel (Figure 1). For the August 2005 model the scans from 3 viewpoints were superposed, whereas the April 2006 3D image was obtained by combining 5 separate scans. The bank was eroded. The bank got eroded essentially on its left part (up to 6.3 m), where it is hit by the river and the debris flows (Figures 2 and 3). A debris cone has also formed (Figure 3), which suggests that a part of the bank erosion is due to shallow landslides. They probably occur when the river erosion creates an undercut slope. These geometrical data allow for the monitoring of the alluvial dynamics (i.e. aggradation and degradation) on different time scales and the influence of debris flows occurrence on these changes. Finally, the resistance against erosion of the bed's cross section and long profile will be analysed to assess the variability of these two key parameters. This information may then be used in debris flow simulation.