Kernel-based Mapping of Meteorological Fields in Complex Orography

Due to the advances in sensor networks and remote sensing technologies, the acquisition and storage rates of meteorological and climatological data increases every day and ask for novel and efficient processing algorithms. A fundamental problem of data analysis and modeling is the spatial prediction of meteorological variables in complex orography, which serves among others to extended climatological analyses, for the assimilation of data into numerical weather prediction models, for preparing inputs to hydrological models and for real time monitoring and short-term forecasting of weather.In this thesis, a new framework for spatial estimation is proposed by taking advantage of a class of algorithms emerging from the statistical learning theory. Nonparametric kernel-based methods for nonlinear data classification, regression and target detection, known as support vector machines (SVM), are adapted for mapping of meteorological variables in complex orography.With the advent of high resolution digital elevation models, the field of spatial prediction met new horizons. In fact, by exploiting image processing tools along with physical heuristics, an incredible number of terrain features which account for the topographic conditions at multiple spatial scales can be extracted. Such features are highly relevant for the mapping of meteorological variables because they control a considerable part of the spatial variability of meteorological fields in the complex Alpine orography. For instance, patterns of orographic rainfall, wind speed and cold air pools are known to be correlated with particular terrain forms, e.g. convex/concave surfaces and upwind sides of mountain slopes.Kernel-based methods are employed to learn the nonlinear statistical dependence which links the multidimensional space of geographical and topographic explanatory variables to the variable of interest, that is the wind speed as measured at the weather stations or the occurrence of orographic rainfall patterns as extracted from sequences of radar images. Compared to low dimensional models integrating only the geographical coordinates, the proposed framework opens a way to regionalize meteorological variables which are multidimensional in nature and rarely show spatial auto-correlation in the original space making the use of classical geostatistics tangled.The challenges which are explored during the thesis are manifolds. First, the complexity of models is optimized to impose appropriate smoothness properties and reduce the impact of noisy measurements. Secondly, a multiple kernel extension of SVM is considered to select the multiscale features which explain most of the spatial variability of wind speed. Then, SVM target detection methods are implemented to describe the orographic conditions which cause persistent and stationary rainfall patterns. Finally, the optimal splitting of the data is studied to estimate realistic performances and confidence intervals characterizing the uncertainty of predictions.The resulting maps of average wind speeds find applications within renewable resources assessment and opens a route to decrease the temporal scale of analysis to meet hydrological requirements. Furthermore, the maps depicting the susceptibility to orographic rainfall enhancement can be used to improve current radar-based quantitative precipitation estimation and forecasting systems and to generate stochastic ensembles of precipitation fields conditioned upon the orography.

Volume-targeted and whole-heart coronary magnetic resonance angiography using an intravascular contrast agent.

PURPOSE: To compare volume-targeted and whole-heart coronary magnetic resonance angiography (MRA) after the administration of an intravascular contrast agent. MATERIALS AND METHODS: Six healthy adult subjects underwent a navigator-gated and -corrected (NAV) free breathing volume-targeted cardiac-triggered inversion recovery (IR) 3D steady-state free precession (SSFP) coronary MRA sequence (t-CMRA) (spatial resolution = 1 x 1 x 3 mm(3)) and high spatial resolution IR 3D SSFP whole-heart coronary MRA (WH-CMRA) (spatial resolution = 1 x 1 x 2 mm(3)) after the administration of an intravascular contrast agent B-22956. Subjective and objective image quality parameters including maximal visible vessel length, vessel sharpness, and visibility of coronary side branches were evaluated for both t-CMRA and WH-CMRA. RESULTS: No significant differences (P = NS) in image quality were observed between contrast-enhanced t-CMRA and WH-CMRA. However, using an intravascular contrast agent, significantly longer vessel segments were measured on WH-CMRA vs. t-CMRA (right coronary artery [RCA] 13.5 +/- 0.7 cm vs. 12.5 +/- 0.2 cm; P < 0.05; and left circumflex coronary artery [LCX] 11.9 +/- 2.2 cm vs. 6.9 +/- 2.4 cm; P < 0.05). Significantly more side branches (13.3 +/- 1.2 vs. 8.7 +/- 1.2; P < 0.05) were visible for the left anterior descending coronary artery (LAD) on WH-CMRA vs. t-CMRA. Scanning time and navigator efficiency were similar for both techniques (t-CMRA: 6.05 min; 49% vs. WH-CMRA: 5.51 min; 54%, both P = NS). CONCLUSION: Both WH-CMRA and t-CMRA using SSFP are useful techniques for coronary MRA after the injection of an intravascular blood-pool agent. However, the vessel conspicuity for high spatial resolution WH-CMRA is not inferior to t-CMRA, while visible vessel length and the number of visible smaller-diameter vessels and side-branches are improved.

Baseline Results of an Acute Ischemic Stroke Registry Including Acute Multimodal Imaging (ASTRAL)

Background: The DEFUSE (n_74) and EPITHET (n_101) studies have in common that a baseline MRI was obtained prior to treatment (tPA in DEFUSE; tPA or placebo in EPITHET) in the 3-6 hour time-window. There were however important methodological differences between the studies. A standardized reanalysis of pooled data was undertaken to determine the effect of these differences on baseline characteristics and study outcomes. Methods: To standardize the studies 1) the DWI and PWI source images were reprocessed and segmented using automated image processing software (RAPID); 2) patients were categorized according to their baseline MRI profile as either Target Mismatch (PWITmax_6/DWI ratio_ 1.8 and an absolute mismatch _15mL), Malignant (DWI or PWITmax_10 lesion _ 100 mL), or No Mismatch. 3) favorable clinical response was defined as NIHSS score of 0-1 or a _8 points improvement on the NIHSSS at day 90. Results: Prior to standardization there was no difference in the proportion of Target Mismatch patients between EPITHET and DEFUSE (54% vs 49%, p_0.6), but the EPITHET study had more patients with the Malignant profile than DEFUSE (35% vs 9%, p_0.01) and fewer patients that had No Mismatch (11% vs 42%, p_0.01). These differences in baseline MRI profiles between EPITHET and DEFUSE were largely eliminated by standardized processing of PWI and DWI images with RAPID software (Target Mismatch 49% vs 48%; Malignant 15% vs 8%; No Mismatch 36% vs 25%; p_NS for all comparisons) Reperfusion was strongly associated with a favorable clinical response in mismatch patients (figure). This relationship was not affected by the standardization procedures (pooled odds ratio of 8.8 based on original data and 6.6 based on standardized data). Conclusion: Standardization of image analyses procedures in acute stroke is important as non-standardized techniques introduce significant variability in DWI and PWI imaging characteristics. Despite methodological differences, the DEFUSE and EPITHET studies show a consistent and robust association between reperfusion and favorable clinical response in Target Mismatch patients regardless of standardization. These data support an RCT of iv tPA in the 3-6 hour time-window for Target Mismatch patients identified using RAPID.

Occlusive arterial disease of abdominal aorta and lower extremities: comparison of helical CT angiography with transcatheter angiography.

The purpose of this study was to evaluate helical CT angiography in the assessment of occlusive arterial disease of abdominal aorta and the lower extremities. Sixteen patients underwent both transcatheter angiography and helical CT. Helical CT was inconclusive in 6.2% of segments whereas angiography was inconclusive in 5%. The overall sensitivity of helical CT was 91% and specificity 93%. Segmental analysis found a sensitivity of 43% in infrapopliteal arteries, and a specificity of 86%.

Sparse image reconstruction on the sphere: implications of a new sampling theorem.

We study the impact of sampling theorems on the fidelity of sparse image reconstruction on the sphere. We discuss how a reduction in the number of samples required to represent all information content of a band-limited signal acts to improve the fidelity of sparse image reconstruction, through both the dimensionality and sparsity of signals. To demonstrate this result, we consider a simple inpainting problem on the sphere and consider images sparse in the magnitude of their gradient. We develop a framework for total variation inpainting on the sphere, including fast methods to render the inpainting problem computationally feasible at high resolution. Recently a new sampling theorem on the sphere was developed, reducing the required number of samples by a factor of two for equiangular sampling schemes. Through numerical simulations, we verify the enhanced fidelity of sparse image reconstruction due to the more efficient sampling of the sphere provided by the new sampling theorem.

An efficient total variation algorithm for super-resolution in fetal brain MRI with adaptive regularization.

Although fetal anatomy can be adequately viewed in new multi-slice MR images, many critical limitations remain for quantitative data analysis. To this end, several research groups have recently developed advanced image processing methods, often denoted by super-resolution (SR) techniques, to reconstruct from a set of clinical low-resolution (LR) images, a high-resolution (HR) motion-free volume. It is usually modeled as an inverse problem where the regularization term plays a central role in the reconstruction quality. Literature has been quite attracted by Total Variation energies because of their ability in edge preserving but only standard explicit steepest gradient techniques have been applied for optimization. In a preliminary work, it has been shown that novel fast convex optimization techniques could be successfully applied to design an efficient Total Variation optimization algorithm for the super-resolution problem. In this work, two major contributions are presented. Firstly, we will briefly review the Bayesian and Variational dual formulations of current state-of-the-art methods dedicated to fetal MRI reconstruction. Secondly, we present an extensive quantitative evaluation of our SR algorithm previously introduced on both simulated fetal and real clinical data (with both normal and pathological subjects). Specifically, we study the robustness of regularization terms in front of residual registration errors and we also present a novel strategy for automatically select the weight of the regularization as regards the data fidelity term. Our results show that our TV implementation is highly robust in front of motion artifacts and that it offers the best trade-off between speed and accuracy for fetal MRI recovery as in comparison with state-of-the art methods.

Comparison of Inoculation with the InoqulA and WASP Automated Systems with Manual Inoculation.

The quality of sample inoculation is critical for achieving an optimal yield of discrete colonies in both monomicrobial and polymicrobial samples to perform identification and antibiotic susceptibility testing. Consequently, we compared the performance between the InoqulA (BD Kiestra), the WASP (Copan), and manual inoculation methods. Defined mono- and polymicrobial samples of 4 bacterial species and cloudy urine specimens were inoculated on chromogenic agar by the InoqulA, the WASP, and manual methods. Images taken with ImagA (BD Kiestra) were analyzed with the VisionLab version 3.43 image analysis software to assess the quality of growth and to prevent subjective interpretation of the data. A 3- to 10-fold higher yield of discrete colonies was observed following automated inoculation with both the InoqulA and WASP systems than that with manual inoculation. The difference in performance between automated and manual inoculation was mainly observed at concentrations of >10(6) bacteria/ml. Inoculation with the InoqulA system allowed us to obtain significantly more discrete colonies than the WASP system at concentrations of >10(7) bacteria/ml. However, the level of difference observed was bacterial species dependent. Discrete colonies of bacteria present in 100- to 1,000-fold lower concentrations than the most concentrated populations in defined polymicrobial samples were not reproducibly recovered, even with the automated systems. The analysis of cloudy urine specimens showed that InoqulA inoculation provided a statistically significantly higher number of discrete colonies than that with WASP and manual inoculation. Consequently, the automated InoqulA inoculation greatly decreased the requirement for bacterial subculture and thus resulted in a significant reduction in the time to results, laboratory workload, and laboratory costs.

Accelerated Microstructure Imaging via Convex Optimisation for regions with multiple fibres (AMICOx)

This paper reviews and extends our previous work to enable fast axonal diameter mapping from diffusion MRI data in the presence of multiple fibre populations within a voxel. Most of the existing mi-crostructure imaging techniques use non-linear algorithms to fit their data models and consequently, they are computationally expensive and usually slow. Moreover, most of them assume a single axon orientation while numerous regions of the brain actually present more complex configurations, e.g. fiber crossing. We present a flexible framework, based on convex optimisation, that enables fast and accurate reconstructions of the microstructure organisation, not limited to areas where the white matter is coherently oriented. We show through numerical simulations the ability of our method to correctly estimate the microstructure features (mean axon diameter and intra-cellular volume fraction) in crossing regions.

Does abdominoplasty with liposuction of the love handles yield a shorter scar? An analysis with abdominal 3D laser scanning.

The aim of this study was to evaluate the combination of abdominoplasty with liposuction of both flanks with regards to length of scar, complications, and patient's satisfaction. A retrospective analysis of 35 patients who underwent esthetic abdominoplasty at our institution between 2002 and 2004 was performed. Thirteen patients underwent abdominoplasty with liposuction of both flanks, 22 patients underwent conventional abdominoplasty. Liposuction of the flanks did not increase the rate of complications of the abdominoplasty procedures. We found a tendency toward shorter scars in patients who underwent abdominoplasty combined with liposuction of the flanks. Implementation of 3-dimensional laser surface scanning to objectify the postoperative outcomes, documented a comparable degree of flatness of the achieved body contouring in both procedures. 3-dimensional laser surface scanning can be a valuable tool to objectify assessment of postoperative results.

Chronic exercise preserves lean muscle mass in masters athletes.

Aging is commonly associated with a loss of muscle mass and strength, resulting in falls, functional decline, and the subjective feeling of weakness. Exercise modulates the morbidities of muscle aging. Most studies, however, have examined muscle-loss changes in sedentary aging adults. This leaves the question of whether the changes that are commonly associated with muscle aging reflect the true physiology of muscle aging or whether they reflect disuse atrophy. This study evaluated whether high levels of chronic exercise prevents the loss of lean muscle mass and strength experienced in sedentary aging adults. A cross-section of 40 high-level recreational athletes ("masters athletes") who were aged 40 to 81 years and trained 4 to 5 times per week underwent tests of health/activity, body composition, quadriceps peak torque (PT), and magnetic resonance imaging of bilateral quadriceps. Mid-thigh muscle area, quadriceps area (QA), subcutaneous adipose tissue, and intramuscular adipose tissue were quantified in magnetic resonance imaging using medical image processing, analysis, and visualization software. One-way analysis of variance was used to examine age group differences. Relationships were evaluated using Spearman correlations. Mid-thigh muscle area (P = 0.31) and lean mass (P = 0.15) did not increase with age and were significantly related to retention of mid-thigh muscle area (P < 0.0001). This occurred despite an increase in total body fat percentage (P = 0.003) with age. Mid-thigh muscle area (P = 0.12), QA (P = 0.17), and quadriceps PT did not decline with age. Specific strength (strength per QA) did not decline significantly with age (P = 0.06). As muscle area increased, PT increased significantly (P = 0.008). There was not a significant relationship between intramuscular adipose tissue (P = 0.71) or lean mass (P = 0.4) and PT. This study contradicts the common observation that muscle mass and strength decline as a function of aging alone. Instead, these declines may signal the effect of chronic disuse rather than muscle aging. Evaluation of masters athletes removes disuse as a confounding variable in the study of lower-extremity function and loss of lean muscle mass. This maintenance of muscle mass and strength may decrease or eliminate the falls, functional decline, and loss of independence that are commonly seen in aging adults.

Three-dimensional cholangio-spiral CT demonstration of a post-traumatic bile leak in a child.

Bilioma is a rare complication of traumatic liver injury, and the precise site of bile leak is often difficult to demonstrate with a non-invasive technique. We report a case of post-traumatic bile leak in a 15-year-old girl in whom spiral CT after intravenous cholangiography allowed excellent preoperative demonstration of the extent of the liver rupture and an exact location of the bile leak. We think that spiral-CT cholangiography could be an accurate, non-invasive technique to investigate the biliary system in cases of paediatric liver trauma.

Reliable acetabular cup orientation with a new gravity-assisted guidance system.

Acetabular cup orientation is a key factor determining hip stability, and standard mechanical guides have shown little help in improving alignment. An in vitro study was carried out to compare the accuracy and precision of a new gravity-assisted guidance system with a standard mechanical guide. Three hundred ten cups were impacted by 5 surgeons, and the final cup orientation was measured. With the new guide, the average error in anteversion was 0.4 degrees , compared with 10.4 degrees with the standard guide and 0.3 degrees and -4.7 degrees , respectively, for abduction angles. The average time required for orienting the cups was similar for both guides. The accuracy and reproducibility obtained with the new guide were better (P < .0001). These good results would require a clinical validation.

Segmentation of brain structures in presence of a space-occupying lesion.

Brain deformations induced by space-occupying lesions may result in unpredictable position and shape of functionally important brain structures. The aim of this study is to propose a method for segmentation of brain structures by deformation of a segmented brain atlas in presence of a space-occupying lesion. Our approach is based on an a priori model of lesion growth (MLG) that assumes radial expansion from a seeding point and involves three steps: first, an affine registration bringing the atlas and the patient into global correspondence; then, the seeding of a synthetic tumor into the brain atlas providing a template for the lesion; finally, the deformation of the seeded atlas, combining a method derived from optical flow principles and a model of lesion growth. The method was applied on two meningiomas inducing a pure displacement of the underlying brain structures, and segmentation accuracy of ventricles and basal ganglia was assessed. Results show that the segmented structures were consistent with the patient's anatomy and that the deformation accuracy of surrounding brain structures was highly dependent on the accurate placement of the tumor seeding point. Further improvements of the method will optimize the segmentation accuracy. Visualization of brain structures provides useful information for therapeutic consideration of space-occupying lesions, including surgical, radiosurgical, and radiotherapeutic planning, in order to increase treatment efficiency and prevent neurological damage.

Double-oblique free-breathing high resolution three-dimensional coronary magnetic resonance angiography.

OBJECTIVES: The goal of the present study was to develop a strategy for three-dimensional (3D) volume acquisition along the major axes of the coronary arteries. BACKGROUND: For high-resolution 3D free-breathing coronary magnetic resonance angiography (MRA), coverage of the coronary artery tree may be limited due to excessive measurement times associated with large volume acquisitions. Planning the 3D volume along the major axis of the coronary vessels may help to overcome such limitations. METHODS: Fifteen healthy adult volunteers and seven patients with X-ray angiographically confirmed coronary artery disease underwent free-breathing navigator-gated and corrected 3D coronary MRA. For an accurate volume targeting of the high resolution scans, a three-point planscan software tool was applied. RESULTS: The average length of contiguously visualized left main and left anterior descending coronary artery was 81.8 +/- 13.9 mm in the healthy volunteers and 76.2 +/- 16.5 mm in the patients (p = NS). For the right coronary artery, a total length of 111.7 +/- 27.7 mm was found in the healthy volunteers and 79.3 +/- 4.6 mm in the patients (p = NS). Comparing coronary MRA and X-ray angiography, a good agreement of anatomy and pathology was found in the patients. CONCLUSIONS: Double-oblique submillimeter free-breathing coronary MRA allows depiction of extensive parts of the native coronary arteries. The results obtained in patients suggest that the method has the potential to be applied in broader prospective multicenter studies where coronary MRA is compared with X-ray angiography.