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.

Hierarchical Markov Random Fields Applied to Model Soft Tissue Deformations on Graphics Hardware

Many methodologies dealing with prediction or simulation of soft tissue deformations on medical image data require preprocessing of the data in order to produce a different shape representation that complies with standard methodologies, such as mass–spring networks, finite element method s (FEM). On the other hand, methodologies working directly on the image space normally do not take into account mechanical behavior of tissues and tend to lack physics foundations driving soft tissue deformations. This chapter presents a method to simulate soft tissue deformations based on coupled concepts from image analysis and mechanics theory. The proposed methodology is based on a robust stochastic approach that takes into account material properties retrieved directly from the image, concepts from continuum mechanics and FEM. The optimization framework is solved within a hierarchical Markov random field (HMRF) which is implemented on the graphics processor unit (GPU See Graphics processing unit ).

3D Ultrasound imaging--a useful non-invasive tool to detect AV fistulas in transplanted kidneys

BACKGROUND: A precise, non-invasive, non-toxic, repeatable, convenient and inexpensive follow-up of renal transplants, especially following biopsies, is in the interest of nephrologists. Formerly, the rate of biopsies leading to AV fistulas had been underestimated. Imaging procedures suited to a detailed judgement of these vascular malformations are to be assessed. METHODS: Three-dimensional (3D) reconstruction techniques of ultrasound flow-directed and non-flow-directed energy mode pictures were compared with a standard procedure, gadolinium-enhanced nuclear magnetic resonance imaging angiography (MRA) using the phase contrast technique. RESULTS: Using B-mode and conventional duplex information, AV fistulas were localized in the upper pole of the kidney transplant of the index patient. The 3D reconstruction provided information about the exact localization and orientation of the fistula in relation to other vascular structures, and the flow along the fistula. The MRA provided localization and orientation information, but less functional information. Flow-directed and non-flow-directed energy mode pictures could be reconstructed to provide 3D information about vascular malformations in transplanted kidneys. CONCLUSION: In transplanted kidneys, 3D-ultrasound angiography may be equally as effective as MRA in localizing and identifying AV malformations. Advantages of the ultrasound method are that it is cheaper, non-toxic, non-invasive, more widely availability and that it even provides more functional information. Future prospective studies will be necessary to evaluate the two techniques further.

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.

Algorithms For Automatic And Robust Registration Of 3D Head Scans

wo methods for registering laser-scans of human heads and transforming them to a new semantically consistent topology defined by a user-provided template mesh are described. Both algorithms are stated within the Iterative Closest Point framework. The first method is based on finding landmark correspondences by iteratively registering the vicinity of a landmark with a re-weighted error function. Thin-plate spline interpolation is then used to deform the template mesh and finally the scan is resampled in the topology of the deformed template. The second algorithm employs a morphable shape model, which can be computed from a database of laser-scans using the first algorithm. It directly optimizes pose and shape of the morphable model. The use of the algorithm with PCA mixture models, where the shape is split up into regions each described by an individual subspace, is addressed. Mixture models require either blending or regularization strategies, both of which are described in detail. For both algorithms, strategies for filling in missing geometry for incomplete laser-scans are described. While an interpolation-based approach can be used to fill in small or smooth regions, the model-driven algorithm is capable of fitting a plausible complete head mesh to arbitrarily small geometry, which is known as "shape completion". The importance of regularization in the case of extreme shape completion is shown.

A 950 yr temperature reconstruction from Duckhole Lake, southern Tasmania, Australia

A lack of quantitative high resolution paleoclimate data from the Southern Hemisphere limits the ability to examine current trends within the context of long-term natural climate variability. This study presents a temperature reconstruction for southern Tasmania based on analyses of a sediment core from Duckhole Lake (43.365°S, 146.875°E). The relationship between non-destructive whole core scanning reflectance spectroscopy measurements in the visible spectrum (380–730 nm) and the instrumental temperature record (ad 1911–2000) was used to develop a calibration-in-time reflectance spectroscopy-based temperature model. Results showed that a trough in reflectance from 650 to 700 nm, which represents chlorophyll and its derivatives, was significantly correlated to annual mean temperature. A calibration model was developed (R = 0.56, p auto < 0.05, root mean squared error of prediction (RMSEP) = 0.21°C, five-year filtered data, calibration period 1911–2000) and applied down-core to reconstruct annual mean temperatures in southern Tasmania over the last c. 950 years. This indicated that temperatures were initially cool c. ad 1050, but steadily increased until the late ad 1100s. After a brief cool period in the ad 1200s, temperatures again increased. Temperatures steadily decreased during the ad 1600s and remained relatively stable until the start of the 20th century when they rapidly decreased, before increasing from ad 1960s onwards. Comparisons with high resolution temperature records from western Tasmania, New Zealand and South America revealed some similarities, but also highlighted differences in temperature variability across the mid-latitudes of the Southern Hemisphere. These are likely due to a combination of factors including the spatial variability in climate between and within regions, and differences between records that document seasonal (i.e. warm season/late summer) versus annual temperature variability. This highlights the need for further records from the mid-latitudes of the Southern Hemisphere in order to constrain past natural spatial and seasonal/annual temperature variability in the region, and to accurately identify and attribute changes to natural variability and/or anthropogenic activities.

The Advantages of Parametric Modeling for the Reconstruction of Historic Buildings

Consecrated in 1297 as the monastery church of the four years earlier founded St. Catherine’s monastery, the Gothic Church of St. Catherine was largely destroyed in a devastating bombing raid on January 2nd 1945. To counteract the process of disintegration, the departments of geo-information and lower monument protection authority of the City of Nuremburg decided to getting done a three dimensional building model of the Church of St. Catherine’s. A heterogeneous set of data was used for preparation of a parametric architectural model. In effect the modeling of historic buildings can profit from the so called BIM method (Building Information Modeling), as the necessary structuring of the basic data renders it into very sustainable information. The resulting model is perfectly suited to deliver a vivid impression of the interior and exterior of this former mendicant orders’ church to present observers.

Nonviral Gene Delivery of Growth and Differentiation Factor 5 to Human Mesenchymal Stem Cells Injected into a 3D Bovine Intervertebral Disc Organ Culture System

Intervertebral disc (IVD) cell therapy with unconditioned 2D expanded mesenchymal stem cells (MSC) is a promising concept yet challenging to realize. Differentiation of MSCs by nonviral gene delivery of growth and differentiation factor 5 (GDF5) by electroporation mediated gene transfer could be an excellent source for cell transplantation. Human MSCs were harvested from bone marrow aspirate and GDF5 gene transfer was achieved by in vitro electroporation. Transfected cells were cultured as monolayers and as 3D cultures in 1.2% alginate bead culture. MSC expressed GDF5 efficiently for up to 21 days. The combination of GDF5 gene transfer and 3D culture in alginate showed an upregulation of aggrecan and SOX9, two markers for chondrogenesis, and KRT19 as a marker for discogenesis compared to untransfected cells. The cells encapsulated in alginate produced more proteoglycans expressed in GAG/DNA ratio. Furthermore, GDF5 transfected MCS injected into an IVD papain degeneration organ culture model showed a partial recovery of the GAG/DNA ratio after 7 days. In this study we demonstrate the potential of GDF5 transfected MSC as a promising approach for clinical translation for disc regeneration.

Beyond the tree of texts: Building an empirical model of scribal variation through graph analysis of texts and stemmata

Stemmatology, or the reconstruction of the transmission history of texts, is a field that stands particularly to gain from digital methods. Many scholars already take stemmatic approaches that rely heavily on computational analysis of the collated text (e.g. Robinson and O’Hara 1996; Salemans 2000; Heikkilä 2005; Windram et al. 2008 among many others). Although there is great value in computationally assisted stemmatology, providing as it does a reproducible result and allowing access to the relevant methodological process in related fields such as evolutionary biology, computational stemmatics is not without its critics. The current state-of-the-art effectively forces scholars to choose between a preconceived judgment of the significance of textual differences (the Lachmannian or neo-Lachmannian approach, and the weighted phylogenetic approach) or to make no judgment at all (the unweighted phylogenetic approach). Some basis for judgment of the significance of variation is sorely needed for medieval text criticism in particular. By this, we mean that there is a need for a statistical empirical profile of the text-genealogical significance of the different sorts of variation in different sorts of medieval texts. The rules that apply to copies of Greek and Latin classics may not apply to copies of medieval Dutch story collections; the practices of copying authoritative texts such as the Bible will most likely have been different from the practices of copying the Lives of local saints and other commonly adapted texts. It is nevertheless imperative that we have a consistent, flexible, and analytically tractable model for capturing these phenomena of transmission. In this article, we present a computational model that captures most of the phenomena of text variation, and a method for analysis of one or more stemma hypotheses against the variation model. We apply this method to three ‘artificial traditions’ (i.e. texts copied under laboratory conditions by scholars to study the properties of text variation) and four genuine medieval traditions whose transmission history is known or deduced in varying degrees. Although our findings are necessarily limited by the small number of texts at our disposal, we demonstrate here some of the wide variety of calculations that can be made using our model. Certain of our results call sharply into question the utility of excluding ‘trivial’ variation such as orthographic and spelling changes from stemmatic analysis.

A papain-induced disc degeneration model for the assessment of thermo-reversible hydrogel-cells therapeutic approach

Nucleus pulposus (NP) regeneration by the application of injectable cell-embedded hydrogels is an appealing approach for tissue engineering. We investigated a thermo-reversible hydrogel (TR-HG), based on a modified polysaccharide with a thermo-reversible polyamide [poly(N-isopropylacrylamide), pNIPAM], which is made to behave as a liquid at room temperature and hardens at > 32 °C. In order to test the hydrogel, a papain-induced bovine caudal disc degeneration model (PDDM), creating a cavity in the NP, was employed. Human mesenchymal stem cells (hMSCs) or autologous bovine NP cells (bNPCs) were seeded in TR-HG; hMSCs were additionally preconditioned with rhGDF-5 for 7 days. Then, TR-HG was reversed to a fluid and the cell suspension injected into the PDDM and kept under static loading for 7 days. Experimental design was: (D1) fresh disc control + PBS injection; (D2) PDDM + PBS injection; (D3) PDDM + TR-HG (material control); (D4) PDDM + TR-HG + bNPCs; (D5) PDDM + TR-HG + hMSCs. Magnetic resonance imaging performed before and after loading, on days 9 and 16, allowed imaging of the hydrogel-filled PDDM and assessment of disc height and volume changes. In gel-injected discs the NP region showed a major drop in volume and disc height during culture under static load. The RT–PCR results of injected hMSCs showed significant upregulation of ACAN, COL2A1, VCAN and SOX9 during culture in the disc cavity, whereas the gene expression profile of NP cells remained unchanged. The cell viability of injected cells (NPCs or hMSCs) was maintained at over 86% in 3D culture and dropped to ~72% after organ culture. Our results underline the need for load-bearing hydrogels that are also cyto-compatible.

A prehistoric thoroughfare between the Ganges and the Himalayas

The greater Himalayan region, including the Tibetan plateau in the north and the Gangetic plain in the south, served as the principal prehistoric thoroughfare for the peopling of East and Southeast Asia. The descendants of ancient migrants through this region ultimately settled lands as far away as New Zealand, Madagascar and the Americas. Several of the keys to understanding the ethnogenesis of human diversity in Asia include the Father Tongue correlation, possible refugia during the Last Glacial Maximum and the hypothesis that language families may have arisen as the result of demographic bottlenecks in prehistory. Ethnolinguistically informed inferences based on Asian Y chromosomal phylogeography permit a reconstruction of episodes of ethnolinguistic prehistory which lie beyond the linguistic event horizon, i.e. beyond the time depth empirically accessible to historical linguistics. The origins of the language families which make up the hypothetical Uralo-Siberian and East Asian linguistic phyla are argued to have lain in the northeastern corner of the Indian subcontinent. Several other Asian language families are shown to be tied to the subcontinent. The Centripetal Migration model, which assumes that migrations in quest of a better life unfolded in both centrifugal and centripetal directions with respect to technologically more advanced centres of civilisation, is opposed to the Farming Language Dispersal theory, which assumes that all linguistic dispersals were driven by agricultural centrifugal migration.

A new methodological protocol for the use of dendrogeomorphological data in flood risk analysis

Dendrogeomorphology uses information sources recorded in the roots, trunks and branches of trees and bushes located in the fluvial system to complement (or sometimes even replace) systematic and palaeohydrological records of past floods. The application of dendrogeomorphic data sources and methods to palaeoflood analysis over nearly 40 years has allowed improvements to be made in frequency and magnitude estimations of past floods. Nevertheless, research carried out so far has shown that the dendrogeomorphic indicators traditionally used (mainly scar evidence), and their use to infer frequency and magnitude, have been restricted to a small, limited set of applications. New possibilities with enormous potential remain unexplored. New insights in future research of palaeoflood frequency and magnitude using dendrogeomorphic data sources should: (1) test the application of isotopic indicators (16O/18O ratio) to discover the meteorological origin of past floods; (2) use different dendrogeomorphic indicators to estimate peak flows with 2D (and 3D) hydraulic models and study how they relate to other palaeostage indicators; (3) investigate improved calibration of 2D hydraulic model parameters (roughness); and (4) apply statistics-based cost–benefit analysis to select optimal mitigation measures. This paper presents an overview of these innovative methodologies, with a focus on their capabilities and limitations in the reconstruction of recent floods and palaeofloods.

Northern Hemisphere temperature reconstruction during the last millennium using multiple annual proxies

Previous studies have either exclusively used annual tree-ring data or have combined tree-ring series with other, lower temporal resolution proxy series. Both approaches can lead to significant uncertainties, as tree-rings may underestimate the amplitude of past temperature variations, and the validity of non-annual records cannot be clearly assessed. In this study, we assembled 45 published Northern Hemisphere (NH) temperature proxy records covering the past millennium, each of which satisfied 3 essential criteria: the series must be of annual resolution, span at least a thousand years, and represent an explicit temperature signal. Suitable climate archives included ice cores, varved lake sediments, tree-rings and speleothems. We reconstructed the average annual land temperature series for the NH over the last millennium by applying 3 different reconstruction techniques: (1) principal components (PC) plus second-order autoregressive model (AR2), (2) composite plus scale (CPS) and (3) regularized errors-in-variables approach (EIV). Our reconstruction is in excellent agreement with 6 climate model simulations (including the first 5 models derived from the fifth phase of the Coupled Model Intercomparison Project (CMIP5) and an earth system model of intermediate complexity (LOVECLIM), showing similar temperatures at multi-decadal timescales; however, all simulations appear to underestimate the temperature during the Medieval Warm Period (MWP). A comparison with other NH reconstructions shows that our results are consistent with earlier studies. These results indicate that well-validated annual proxy series should be used to minimize proxy-based artifacts, and that these proxy series contain sufficient information to reconstruct the low-frequency climate variability over the past millennium.

A model-data comparison of Holocene timberline changes in the Swiss Alps reveals past and future drivers of mountain forest dynamics

Mountain vegetation is strongly affected by temperature and is expected to shift upwards with climate change. Dynamic vegetation models are often used to assess the impact of climate on vegetation and model output can be compared with paleobotanical data as a reality check. Recent paleoecological studies have revealed regional variation in the upward shift of timberlines in the Northern and Central European Alps in response to rapid warming at the Younger Dryas/Preboreal transition ca. 11700years ago, probably caused by a climatic gradient across the Alps. This contrasts with previous studies that successfully simulated the early Holocene afforestation in the (warmer) Central Alps with a chironomid-inferred temperature reconstruction from the (colder) Northern Alps. We use LandClim, a dynamic landscape vegetation model to simulate mountain forests under different temperature, soil and precipitation scenarios around Iffigsee (2065m a.s.l.) a lake in the Northwestern Swiss Alps, and compare the model output with the paleobotanical records. The model clearly overestimates the upward shift of timberline in a climate scenario that applies chironomid-inferred July-temperature anomalies to all months. However, forest establishment at 9800 cal. BP at Iffigsee is successfully simulated with lower moisture availability and monthly temperatures corrected for stronger seasonality during the early Holocene. The model-data comparison reveals a contraction in the realized niche of Abies alba due to the prominent role of anthropogenic disturbance after ca. 5000 cal. BP, which has important implications for species distribution models (SDMs) that rely on equilibrium with climate and niche stability. Under future climate projections, LandClim indicates a rapid upward shift of mountain vegetation belts by ca. 500m and treeline positions of ca. 2500m a.s.l. by the end of this century. Resulting biodiversity losses in the alpine vegetation belt might be mitigated with low-impact pastoralism to preserve species-rich alpine meadows.

Temperature reconstruction from 10 to 120 kyr b2k from the NGRIP ice core

In order to reconstruct the temperature of the North Greenland Ice Core Project (NGRIP) site, new measurements of δ15N have been performed covering the time period from the beginning of the Holocene to Dansgaard–Oeschger (DO) event 8. Together with previously measured and mostly published δ15N data, we present for the first time a NGRIP temperature reconstruction for the whole last glacial period from 10 to 120 kyr b2k (thousand years before 2000 AD) including every DO event based on δ15N isotope measurements combined with a firn densification and heat diffusion model. The detected temperature rises at the onset of DO events range from 5 °C (DO 25) up to 16.5 °C (DO 11) with an uncertainty of ±3 °C. To bring measured and modelled data into agreement, we had to reduce the accumulation rate given by the NGRIP ss09sea06bm timescale in some periods by 30 to 35%, especially during the last glacial maximum. A comparison between reconstructed temperature and δ18Oice data confirms that the isotopic composition of the stadial was strongly influenced by seasonality. We evidence an anticorrelation between the variations of the δ18Oice sensitivity to temperature (referred to as α) and obliquity in agreement with a simple Rayleigh distillation model. Finally, we suggest that α might be influenced by the Northern Hemisphere ice sheet volume.