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.

USING AUTO- AND CROSS-CORRELATIONS FROM SEISMIC NOISE TO MONITOR VELOCITY CHANGES AT VILLARRICA VOLCANO, CHILE

We used the Green's functions from auto-correlations and cross-correlations of seismic ambient noise to monitor temporal velocity changes in the subsurface at Villarrica volcano in the Southern Andes of Chile. Campaigns were conducted from March to October 2010 and February to April 2011 with 8 broadband and 6 short-period stations, respectively. We prepared the data by removing the instrument response, normalizing with a root-mean-square method, whitening the spectra, and filtering from 1 to 10 Hz. This frequency band was chosen based on the relatively high background noise level in that range. Hour-long auto- and cross-correlations were computed and the Green's functions stacked by day and total time. To track the temporal velocity changes we stretched a 24 hour moving window of correlation functions from 90% to 110% of the original and cross correlated them with the total stack. All of the stations' auto-correlations detected what is interpreted as an increase in velocity in 2010, with an average increase of 0.13%. Cross-correlations from station V01, near the summit, to the other stations show comparable changes that are also interpreted as increases in velocity. We attribute this change to the closing of cracks in the subsurface due either to seasonal snow loading or regional tectonics. In addition to the common increase in velocity across the stations, there are excursions in velocity on the same order lasting several days. Amplitude decreases as the station's distance from the vent increases suggesting these excursions may be attributed to changes within the volcanic edifice. In at least two occurrences the amplitudes at stations V06 and V07, the stations farthest from the vent, are smaller. Similar short temporal excursions were seen in the auto-correlations from 2011, however, there was little to no increase in the overall velocity.

Individualization of 5-fluorocytosine therapy

Using a convenient and fast HPLC procedure we determined serum concentrations of the fungistatic agent 5-fluorocytosine (5-FC) in 375 samples from 60 patients treated with this drug. The mean trough concentration (n = 127) was 64.3 mg/l (range: 11.8-208.0 mg/l), the mean peak concentration (n = 122) was 99.9 mg/l (range: 25.6-263.8 mg/l), the mean nonpeak/nontrough concentration (n = 126) was 80.1 mg/l (range: 10.5-268.0 mg/l). Totally 134 (35.7%) samples were outside the therapeutic range (25-100 mg/l), 108 (28.8%) being too high, 26 (6.9%) being too low. Forty-four (73%) patients showed 5-FC serum concentrations outside the therapeutic range at least once during the treatment course. In a prospective study we performed 65 dosage predictions on 30 patients by use of a 3-point method previously developed for aminoglycoside dosage adaptation. The mean absolute prediction error of the dosage adaptation was +0.7 mg/l (range: -26.0 to +28.0 mg/l). The root mean square prediction error was 10.7 mg/l. The mean predicted concentration (65.3 mg/l) agreed very well with the mean measured concentration (64.6 mg/l). The frequency distribution of 5-FC serum concentrations indicates that 5-FC monitoring is important. The applied pharmacokinetic method allows individual adaptations of 5-FC dosage with a clinically acceptable prediction error.

Landmark-based augmented reality system for paranasal and transnasal endoscopic surgeries

BACKGROUND: In this paper we present a landmark-based augmented reality (AR) endoscope system for endoscopic paranasal and transnasal surgeries along with fast and automatic calibration and registration procedures for the endoscope. METHODS: Preoperatively the surgeon selects natural landmarks or can define new landmarks in CT volume. These landmarks are overlaid, after proper registration of preoperative CT to the patient, on the endoscopic video stream. The specified name of the landmark, along with selected colour and its distance from the endoscope tip, is also augmented. The endoscope optics are calibrated and registered by fast and automatic methods. Accuracy of the system is evaluated in a metallic grid and cadaver set-up. RESULTS: Root mean square (RMS) error of the system is 0.8 mm in a controlled laboratory set-up (metallic grid) and was 2.25 mm during cadaver studies. CONCLUSIONS: A novel landmark-based AR endoscope system is implemented and its accuracy is evaluated. Augmented landmarks will help the surgeon to orientate and navigate the surgical field. Studies prove the capability of the system for the proposed application. Further clinical studies are planned in near future.

Nighttime vagal cardiac control and plasma fibrinogen levels in a population of working men and women

BACKGROUND: Elevated plasma fibrinogen levels have prospectively been associated with an increased risk of coronary artery disease in different populations. Plasma fibrinogen is a measure of systemic inflammation crucially involved in atherosclerosis. The vagus nerve curtails inflammation via a cholinergic antiinflammatory pathway. We hypothesized that lower vagal control of the heart relates to higher plasma fibrinogen levels. METHODS: Study participants were 559 employees (age 17-63 years; 89% men) of an airplane manufacturing plant in southern Germany. All subjects underwent medical examination, blood sampling, and 24-hour ambulatory heart rate recording while kept on their work routine. The root mean square of successive differences in RR intervals during the night period (nighttime RMSSD) was computed as the heart rate variability index of vagal function. RESULTS: After controlling for demographic, lifestyle, and medical factors, nighttime RMSSD explained 1.7% (P = 0.001), 0.8% (P = 0.033), and 7.8% (P = 0.007), respectively, of the variance in fibrinogen levels in all subjects, men, and women. Nighttime RMSSD and fibrinogen levels were stronger correlated in women than in men. In all workers, men, and women, respectively, there was a mean +/- SEM increase of 0.41 +/- 0.13 mg/dL, 0.28 +/- 0.13 mg/dL, and 1.16 +/- 0.41 mg/dL fibrinogen for each millisecond decrease in nighttime RMSSD. CONCLUSIONS: Reduced vagal outflow to the heart correlated with elevated plasma fibrinogen levels independent of the established cardiovascular risk factors. This relationship seemed comparably stronger in women than men. Such an autonomic mechanism might contribute to the atherosclerotic process and its thrombotic complications.

Prediction of soil organic carbon for Ethiopian highlands using soil spectroscopy

Soil spectroscopy was applied for predicting soil organic carbon (SOC) in the highlands of Ethiopia. Soil samples were acquired from Ethiopia’s National Soil Testing Centre and direct field sampling. The reflectance of samples was measured using a FieldSpec 3 diffuse reflectance spectrometer. Outliers and sample relation were evaluated using principal component analysis (PCA) and models were developed through partial least square regression (PLSR). For nine watersheds sampled, 20% of the samples were set aside to test prediction and 80% were used to develop calibration models. Depending on the number of samples per watershed, cross validation or independent validation were used.The stability of models was evaluated using coefficient of determination (R2), root mean square error (RMSE), and the ratio performance deviation (RPD). The R2 (%), RMSE (%), and RPD, respectively, for validation were Anjeni (88, 0.44, 3.05), Bale (86, 0.52, 2.7), Basketo (89, 0.57, 3.0), Benishangul (91, 0.30, 3.4), Kersa (82, 0.44, 2.4), Kola tembien (75, 0.44, 1.9),Maybar (84. 0.57, 2.5),Megech (85, 0.15, 2.6), andWondoGenet (86, 0.52, 2.7) indicating that themodels were stable. Models performed better for areas with high SOC values than areas with lower SOC values. Overall, soil spectroscopy performance ranged from very good to good.

Biogeochemical variability during the past 3.6 million years recorded by FTIR spectroscopy in the sediment record of Lake El'gygytgyn, Far East Russian Arctic

Abstract. A number of studies have shown that Fourier transform infrared spectroscopy (FTIRS) can be applied to quantitatively assess lacustrine sediment constituents. In this study, we developed calibration models based on FTIRS for the quantitative determination of biogenic silica (BSi; n = 420; gradient: 0.9–56.5 %), total organic carbon (TOC; n = 309; gradient: 0–2.9 %), and total inorganic carbon (TIC; n = 152; gradient: 0–0.4 %) in a 318 m-long sediment record with a basal age of 3.6 million years from Lake El’gygytgyn, Far East Russian Arctic. The developed partial least squares (PLS) regression models yield high cross-validated (CV) R2 CV = 0.86–0.91 and low root mean square error of crossvalidation (RMSECV) (3.1–7.0% of the gradient for the different properties). By applying these models to 6771 samples from the entire sediment record, we obtained detailed insight into bioproductivity variations in Lake El’gygytgyn throughout the middle to late Pliocene and Quaternary. High accumulation rates of BSi indicate a productivity maximum during the middle Pliocene (3.6–3.3 Ma), followed by gradually decreasing rates during the late Pliocene and Quaternary. The average BSi accumulation during the middle Pliocene was �3 times higher than maximum accumulation rates during the past 1.5 million years. The indicated progressive deterioration of environmental and climatic conditions in the Siberian Arctic starting at ca. 3.3 Ma is consistent with the first occurrence of glacial periods and the finally complete establishment of glacial–interglacial cycles during the Quaternary.

Antinociceptive effect of buprenorphine and evaluation of the nociceptive withdrawal reflex in foals.

OBJECTIVE To elicit and evaluate the NWR (nociceptive withdrawal reflex) in 2 and 11 day old foals, to investigate if buprenorphine causes antinociception and determine if the NWR response changes with increasing age. The effect of buprenorphine on behaviour was also evaluated. STUDY DESIGN Prospective, experimental cross-over trial. ANIMALS Nine Norwegian Fjord research foals. METHODS Buprenorphine, 10 μg kg(-1) was administered intramuscularly (IM) to the same foal at 2 days and at 11 days of age. The NWR and the effect of buprenorphine were evaluated by electromyograms recorded from the left deltoid muscle following electrical stimulation of the left lateral palmar nerve at the level of the pastern. Mentation, locomotor activity and respiratory rate were recorded before and after buprenorphine administration. RESULTS We were able to evoke the NWR and temporal summation in foals using this model. Buprenorphine decreased the root mean square amplitude following single electrical stimulation (p < 0.001) in both age groups, and increased the NWR threshold following single electrical stimulation in 2 day old foals (p = 0.0012). Repeated electrical stimulation at 2 Hz was more effective to elicit temporal summation compared to 5 Hz (p < 0.001). No effect of age upon the NWR threshold was found (p = 0.34). Sedation when left undisturbed (11 occasions), increased locomotor activity when handled (9 occasions) and tachypnea (13 occasions) were common side-effects of buprenorphine. CONCLUSION AND CLINICAL RELEVANCE These findings indicate that buprenorphine has antinociceptive effect in foals. Opioid side effects often recognized in adult horses also occur in foals.

Image registration for triggered and non-triggered DTI of the human kidney: Reduced variability of diffusion parameter estimation

BACKGROUND: To investigate if non-rigid image-registration reduces motion artifacts in triggered and non-triggered diffusion tensor imaging (DTI) of native kidneys. A secondary aim was to determine, if improvements through registration allow for omitting respiratory-triggering. METHODS: Twenty volunteers underwent coronal DTI of the kidneys with nine b-values (10-700 s/mm2 ) at 3 Tesla. Image-registration was performed using a multimodal nonrigid registration algorithm. Data processing yielded the apparent diffusion coefficient (ADC), the contribution of perfusion (FP ), and the fractional anisotropy (FA). For comparison of the data stability, the root mean square error (RMSE) of the fitting and the standard deviations within the regions of interest (SDROI ) were evaluated. RESULTS: RMSEs decreased significantly after registration for triggered and also for non-triggered scans (P < 0.05). SDROI for ADC, FA, and FP were significantly lower after registration in both medulla and cortex of triggered scans (P < 0.01). Similarly the SDROI of FA and FP decreased significantly in non-triggered scans after registration (P < 0.05). RMSEs were significantly lower in triggered than in non-triggered scans, both with and without registration (P < 0.05). CONCLUSION: Respiratory motion correction by registration of individual echo-planar images leads to clearly reduced signal variations in renal DTI for both triggered and particularly non-triggered scans. Secondarily, the results suggest that respiratory-triggering still seems advantageous.J. Magn. Reson. Imaging 2014. (c) 2014 Wiley Periodicals, Inc.

Repeated electrical stimulations as a tool to evoke temporal summation of nociceptive inputs in healthy, non-medicated experimental sheep.

The nociceptive withdrawal reflex (NWR) model is used in animal pain research to quantify nociception. The aim of this study was to evaluate the NWR evoked by repeated stimulations in healthy, non-medicated standing sheep. Repeated electrical stimulations were applied at 5Hz for 2s to the digital nerves of the right thoracic and the pelvic limbs of 25 standing sheep. The stimulation intensities applied were fractions (0.5, 0.6, 0.7, 0.8, 0.9 and 1) of the individual previously determined nociceptive threshold (It) after single stimulation. Surface-electromyographic activity (EMG) was recorded from the deltoid, the femoral biceps or the peroneus tertius muscles. The repeated stimulation threshold (RS It) was reached if at least one stimulus in the train was followed by a reflex with a minimal root-mean-square-amplitude (RMSA) of 20μV. The behavioural reaction following each series of stimulations was scored on a scale from 0 (no reaction) to 5 (vigorous whole-body reaction). For the deltoid muscle, RS It was 2.3mA (1.6-3mA) with a reaction score of 2 (1-2) and at a fraction of 0.6 (0.5-0.8)×It. For the biceps femoris muscle, RS It was 2.9mA (2.6-4mA) with a reaction score of 1 (1-2) at a fraction of and 0.55 (0.4-0.7)×It while for the peroneus tertius muscle RS It was 3mA (2.8-3.5mA) with a reaction score of 1 (1-2) and at a fraction of 0.8 (0.8-0.95)×It. Both, RMSA and reaction scores increased significantly with increasing stimulation intensities in all muscles (p<0.001). The repeated application of electrical stimuli led to temporal summation of nociceptive inputs and therefore a reduction of the stimulus intensity evoking a withdrawal reaction in healthy, standing sheep. Data achieved in this study can now serve as reference for further clinical or experimental applications of the model in this species.

Multi-Model Data Fusion to Improve an Early Warning System for Hypo-/Hyperglycemic Events

Correct predictions of future blood glucose levels in individuals with Type 1 Diabetes (T1D) can be used to provide early warning of upcoming hypo-/hyperglycemic events and thus to improve the patient's safety. To increase prediction accuracy and efficiency, various approaches have been proposed which combine multiple predictors to produce superior results compared to single predictors. Three methods for model fusion are presented and comparatively assessed. Data from 23 T1D subjects under sensor-augmented pump (SAP) therapy were used in two adaptive data-driven models (an autoregressive model with output correction - cARX, and a recurrent neural network - RNN). Data fusion techniques based on i) Dempster-Shafer Evidential Theory (DST), ii) Genetic Algorithms (GA), and iii) Genetic Programming (GP) were used to merge the complimentary performances of the prediction models. The fused output is used in a warning algorithm to issue alarms of upcoming hypo-/hyperglycemic events. The fusion schemes showed improved performance with lower root mean square errors, lower time lags, and higher correlation. In the warning algorithm, median daily false alarms (DFA) of 0.25%, and 100% correct alarms (CA) were obtained for both event types. The detection times (DT) before occurrence of events were 13.0 and 12.1 min respectively for hypo-/hyperglycemic events. Compared to the cARX and RNN models, and a linear fusion of the two, the proposed fusion schemes represents a significant improvement.

Trends of mean and extreme temperature indices since 1874 at low-elevation sites in the Southern Alps

We describe the recovery of three daily meteorological records for the southern Alps (Domodossola, Riva del Garda, and Rovereto), all starting in the second half of the nineteenth century. We use these new data, along with additional records, to study regional changes in the mean temperature and extreme indices of heat waves and cold spells frequency and duration over the period 1874–2015. The records are homogenized using subdaily cloud cover observations as a constraint for the statistical model, an approach that has never been applied before in the literature. A case study based on a record of parallel observations between a traditional meteorological window and a modern screen shows that the use of cloud cover can reduce the root-mean-square error of the homogenization by up to 30% in comparison to an unaided statistical correction. We find that mean temperature in the southern Alps has increased by 1.4°C per century over the analyzed period, with larger increases in daily minimum temperatures than maximum temperatures. The number of hot days in summer has more than tripled, and a similar increase is observed in duration of heat waves. Cold days in winter have dropped at a similar rate. These trends are mainly caused by climate change over the last few decades.

Diffusion tensor imaging of the human kidney: Does image registration permit scanning without respiratory triggering?

PURPOSE To investigate if image registration of diffusion tensor imaging (DTI) allows omitting respiratory triggering for both transplanted and native kidneys MATERIALS AND METHODS: Nine kidney transplant recipients and eight healthy volunteers underwent renal DTI on a 3T scanner with and without respiratory triggering. DTI images were registered using a multimodal nonrigid registration algorithm. Apparent diffusion coefficient (ADC), the contribution of perfusion (FP ), and the fractional anisotropy (FA) were determined. Relative root mean square errors (RMSE) of the fitting and the standard deviations of the derived parameters within the regions of interest (SDROI ) were evaluated as quality criteria. RESULTS Registration significantly reduced RMSE in all DTI-derived parameters of triggered and nontriggered measurements in cortex and medulla of both transplanted and native kidneys (P < 0.05 for all). In addition, SDROI values were lower with registration for all 16 parameters in transplanted kidneys (14 of 16 SDROI values were significantly reduced, P < 0.04) and for 15 of 16 parameters in native kidneys (9 of 16 SDROI values were significantly reduced, P < 0.05). Comparing triggered versus nontriggered DTI in transplanted kidneys revealed no significant difference for RMSE (P > 0.14) and for SDROI (P > 0.13) of all parameters. In contrast, in native kidneys relative RMSE from triggered scans were significantly lower than those from nontriggered scans (P < 0.02), while SDROI was slightly higher in triggered compared to nontriggered measurements in 15 out of 16 comparisons (significantly for two, P < 0.05). CONCLUSION Registration improves the quality of DTI in native and transplanted kidneys. Diffusion parameters in renal allografts can be measured without respiratory triggering. In native kidneys, respiratory triggering appears advantageous. J. Magn. Reson. Imaging 2016.

Geochemistry of volcanic rocks from the Hikurangi and Manihiki Plateaus

Here we present the first radiometric age data and a comprehensive geochemical data set (including major and trace element and Sr-Nd-Pb-Hf isotope ratios) for samples from the Hikurangi Plateau basement and seamounts on and adjacent to the plateau obtained during the R/V Sonne 168 cruise, in addition to age and geochemical data from DSDP Site 317 on the Manihiki Plateau. The 40Ar/39Ar age and geochemical data show that the Hikurangi basement lavas (118-96 Ma) have surprisingly similar major and trace element and isotopic characteristics to the Ontong Java Plateau lavas (ca. 120 and 90 Ma), primarily the Kwaimbaita-type composition, whereas the Manihiki DSDP Site 317 lavas (117 Ma) have similar compositions to the Singgalo lavas on the Ontong Java Plateau. Alkalic, incompatible-element-enriched seamount lavas (99-87 Ma and 67 Ma) on the Hikurangi Plateau and adjacent to it (Kiore Seamount), however, were derived from a distinct high time-integrated U/Pb (HIMU)-type mantle source. The seamount lavas are similar in composition to similar-aged alkalic volcanism on New Zealand, indicating a second wide-spread event from a distinct source beginning ca. 20 Ma after the plateau-forming event. Tholeiitic lavas from two Osbourn seamounts on the abyssal plain adjacent to the northeast Hikurangi Plateau margin have extremely depleted incompatible element compositions, but incompatible element characteristics similar to the Hikurangi and Ontong Java Plateau lavas and enriched isotopic compositions intermediate between normal mid-ocean-ridge basalt (N-MORB) and the plateau basement. These younger (~52 Ma) seamounts may have formed through remelting of mafic cumulate rocks associated with the plateau formation. The similarity in age and geochemistry of the Hikurangi, Ontong Java and Manihiki Plateaus suggest derivation from a common mantle source. We propose that the Greater Ontong Java Event, during which ?1% of the Earth's surface was covered with volcanism, resulted from a thermo-chemical superplume/dome that stalled at the transition zone, similar to but larger than the structure imaged presently beneath the South Pacific superswell. The later alkalic volcanism on the Hikurangi Plateau and the Zealandia micro-continent may have been part of a second large-scale volcanic event that may have also triggered the final breakup stage of Gondwana, which resulted in the separation of Zealandia fragments from West Antarctica.

Geochemical studies of the Cretaceous-Tertiary boundary in ODP Holes 113-689B and 113-690C

In a study of ODP Hole 689B no iridium (Ir) anomaly was found in Sections 1 through 6 of Core 25X or in Core 26X from the top down to section 2, 3-12 cm. The background Ir abundance averaged 11 parts per trillion (ppt) and a clay-enriched region had nearly the same average, 26 ± 12 ppt. If the Cretaceous-Tertiary (K-T) contact is in the region studied, then sedimentation was not continuous, and the K-T boundary was probably either not deposited or it was eroded away. In a study of Cores 15X and 16X of ODP Hole 690C, an iridium peak with a maximum abundance of 1566 ± 222 ppt was found in Section 4 of Core 15X at 39-40 cm with a half-width of 6.6 cm. Background abundances were ~15 ppt and distinctly higher Ir abundances were observed from 119 cm below to 72 cm above the main peak. The Ir distribution below the main peak is attributed to bioturbation by organisms with burrows extending at least 0.4 m. The Ir distribution above the main peak may be due to the same cause but other explanations may be significant. There are variable enrichments of clay in the mainly CaCO3 sediment of Core 15X, and the stratigraphically lowest part of the most abundant clay deposits is found (within 2 cm) in the same position as the main Ir peak. The clay deposit, which is estimated to be about 50% of the sediment, extends upward ~19 cm and then slowly decreases to a background level of 10% over 1 m. The degree of homogeneity of the clay-rich interval suggests it was not due to episodic volcanism but may have been due to a decrease of the CaCO3 deposition rate which was possibly triggered by the impact of a large asteroid or comet on the Earth.