Application of MFBD Algorithms to Image Reconstruction Under Anisoplanatic Conditions

All optical systems that operate in or through the atmosphere suffer from turbulence induced image blur. Both military and civilian surveillance, gun-sighting, and target identification systems are interested in terrestrial imaging over very long horizontal paths, but atmospheric turbulence can blur the resulting images beyond usefulness. My dissertation explores the performance of a multi-frame-blind-deconvolution technique applied under anisoplanatic conditions for both Gaussian and Poisson noise model assumptions. The technique is evaluated for use in reconstructing images of scenes corrupted by turbulence in long horizontal-path imaging scenarios and compared to other speckle imaging techniques. Performance is evaluated via the reconstruction of a common object from three sets of simulated turbulence degraded imagery representing low, moderate and severe turbulence conditions. Each set consisted of 1000 simulated, turbulence degraded images. The MSE performance of the estimator is evaluated as a function of the number of images, and the number of Zernike polynomial terms used to characterize the point spread function. I will compare the mean-square-error (MSE) performance of speckle imaging methods and a maximum-likelihood, multi-frame blind deconvolution (MFBD) method applied to long-path horizontal imaging scenarios. Both methods are used to reconstruct a scene from simulated imagery featuring anisoplanatic turbulence induced aberrations. This comparison is performed over three sets of 1000 simulated images each for low, moderate and severe turbulence-induced image degradation. The comparison shows that speckle-imaging techniques reduce the MSE 46 percent, 42 percent and 47 percent on average for low, moderate, and severe cases, respectively using 15 input frames under daytime conditions and moderate frame rates. Similarly, the MFBD method provides, 40 percent, 29 percent, and 36 percent improvements in MSE on average under the same conditions. The comparison is repeated under low light conditions (less than 100 photons per pixel) where improvements of 39 percent, 29 percent and 27 percent are available using speckle imaging methods and 25 input frames and 38 percent, 34 percent and 33 percent respectively for the MFBD method and 150 input frames. The MFBD estimator is applied to three sets of field data and the results presented. Finally, a combined Bispectrum-MFBD Hybrid estimator is proposed and investigated. This technique consistently provides a lower MSE and smaller variance in the estimate under all three simulated turbulence conditions.

Valve reconstruction or replacement for long-term biopsy-induced tricuspid regurgitation following heart transplantation

Tricuspid regurgitation following heart transplantation can become a severe problem in a subset of patients, where medical therapy fails. Operative findings are described and results of subsequent results with surgical intervention including repair and replacement are analysed. Although follow-up is short, tricuspid replacement seems superior to reconstruction following heart transplantation. Best results are obtained, if replacement is performed, before right ventricular function deteriorates.

A system for 3-D reconstruction of a patient-specific surface model from calibrated X-ray images

This paper presents a system for 3-D reconstruction of a patient-specific surface model from calibrated X-ray images. Our system requires two X-ray images of a patient with one acquired from the anterior-posterior direction and the other from the axial direction. A custom-designed cage is utilized in our system to calibrate both images. Starting from bone contours that are interactively identified from the X-ray images, our system constructs a patient-specific surface model of the proximal femur based on a statistical model based 2D/3D reconstruction algorithm. In this paper, we present the design and validation of the system with 25 bones. An average reconstruction error of 0.95 mm was observed.

Quantitative assessment of oral orbicular muscle deformation after cleft lip reconstruction: an ultrasound elastography study

Reconstruction of a cleft lip leads inevitably to scar tissue formation. Scar tissue within the restored oral orbicular muscle might be assessed by quantification of the local contractility of this muscle. Furthermore, information about the contraction capability of the oral orbicular muscle is crucial for planning the revision surgery of an individual patient. We used ultrasound elastography to determine the local deformation (strain) of the upper lip and to differentiate contracting muscle from passive scar tissue. Raw ultrasound data (radio-frequency format; rf-) were acquired, while the lips were brought from normal state into a pout condition and back in normal state, in three patients and three normal individuals. During this movement, the oral orbicular muscle contracts and, consequently, thickens in contrast to scar tissue that will not contract, or even expand. An iterative coarse-to-fine strain estimation method was used to calculate the local tissue strain. Analysis of the raw ultrasound data allows estimation of tissue strain with a high precision. The minimum strain that can be assessed reproducibly is 0.1%. In normal individuals, strain of the orbicular oral muscle was in the order of 20%. Also, a uniform strain distribution in the oral orbicular muscle was found. However, in patients deviating values were found in the region of the reconstruction and the muscle tissue surrounding that. In two patients with a successful reconstruction, strain was reduced by 6% in the reconstructed region with respect to the normal parts of the muscle (from 22% to 16% and from 25% to 19%). In a patient with severe aesthetical and functional disability, strain decreased from 30% in the normal region to 5% in the reconstructed region. With ultrasound elastography, the strain of the oral orbicular muscle can be quantified. In healthy subjects, the strain profiles and maximum strain values in all parts of the muscle were similar. The maximum strain of the muscle during pout was 20% +/- 1%. In surgically repaired cleft lips, decreased deformation was observed.

High-resolution summer temperature reconstruction from Lake Silvaplana based on in-situ reflectance spectroscopy

Annually laminated (varved) sediments of proglacial Lake Silvaplana (46 ̊27’N, 9 ̊48’E, 1791 m a.s.l., Engadine, eastern Swiss Alps) provide an excellent archive for quantitative high-resolution (seasonal – annual) reconstruction of high- and lowfrequency climate signals back to AD 1580. The chronology of the core is based on varve counting, Cs-137, Pb-210 and event stratigraphy. In this study we present a reconstruction based on in-situ reflectance spectroscopy. In situ reflectance spectroscopy is known as a cost- and time-effective non destructtive method for semi-quantitative analysis of pigments (e.g., chlorines and carotenoids) and of lithoclastic sediment fractions. Reflectance-dependent absorption (RDA) was measured with a Gretac Macbeth spectrolino at 2 mm resolution. The spectral coverage ranges from 380 nm to 730 nm at 10 nm band resolution. In proglacial Lake Silvaplana, 99% of the sediment is lithoclastic prior to AD 1950. Therefore, we concentrate on absorption features that are characteristic for lithoclastic sediment fractions. In Lake Silvaplana, two significant correlations that are stable in time were found between RDA typical for lithoclastics and meteorological data: (1) the time series R 570 /R 630 (ratio between RDA at 570 nm and 630 nm) of varves in Lake Silvaplana and May to October temperatures at nearby station of Sils correlate highly significantly (calibration period AD 1864 – 1951, r = 0.74, p < 0.01 for 5ptsmoothed series; RMSE is 0.28 ̊C, RE = 0.41 and CE = 0.38), and (2) the minimum reflectance within the 690nm band (min690) data correlate with May to October (calibration period AD 1864 – 1951, r = 0.68, p < 0.01 for 5pt-smoothed series; RMSE = 0.22 ̊C, RE = 0.5, CE = 0.31). Both proxy series (min690nm and R 570 /R 630 values) are internally highly consistent (r = 0.8, p < 0.001). In proglacial Lake Silvaplana the largest amount of sediment is transported by glacial meltwater. The melting season spans approximately from May to October, which gives us a good understanding of the geophysical processes explaining the correlations between lithoclastic proxies and the meteorological data. The reconstructions were extended back to AD 1580 and show a broad corresponddence with fully independent reconstructions from tree rings and documentary data.

Comparison of climate field reconstruction techniques: Application to Europe

This paper presents a comparison of principal component (PC) regression and regularized expectation maximization (RegEM) to reconstruct European summer and winter surface air temperature over the past millennium. Reconstruction is performed within a surrogate climate using the National Center for Atmospheric Research (NCAR) Climate System Model (CSM) 1.4 and the climate model ECHO-G 4, assuming different white and red noise scenarios to define the distortion of pseudoproxy series. We show how sensitivity tests lead to valuable “a priori” information that provides a basis for improving real world proxy reconstructions. Our results emphasize the need to carefully test and evaluate reconstruction techniques with respect to the temporal resolution and the spatial scale they are applied to. Furthermore, we demonstrate that uncertainties inherent to the predictand and predictor data have to be more rigorously taken into account. The comparison of the two statistical techniques, in the specific experimental setting presented here, indicates that more skilful results are achieved with RegEM as low frequency variability is better preserved. We further detect seasonal differences in reconstruction skill for the continental scale, as e.g. the target temperature average is more adequately reconstructed for summer than for winter. For the specific predictor network given in this paper, both techniques underestimate the target temperature variations to an increasing extent as more noise is added to the signal, albeit RegEM less than with PC regression. We conclude that climate field reconstruction techniques can be improved and need to be further optimized in future applications.