Comparison of central macular thickness (CMT) measurement between time domain and spectral domain (SD) optical coherence tomography (OCT) and reproducibility of SD OCT in active neovascular AMD

Purpose: To evaluate the reproducibility of Cirrus-SD OCT measurements and to compare central macular thickness (CMT) measurements between TD-Stratus and SD-Cirrus OCT in patients with active exudative AMD. Methods: Consecutive case series of patients with active exudative AMD seen in the Medical Retina Department. Patients underwent 1 scan with Stratus (macular thickness map protocol) and 5 scans with Cirrus (Macular Cube protocol) at the same visit by the same experienced examiner. To be included, patients best-corrected visual acuity (BCVA) had to be >20/200 while all scans had to be of sufficient quality, well-centered and at least one Cirrus scan with CMT >300 microns. The repeatability of the SD Cirrus was estimated by using all 5 CMT measurements and the mean of the Cirrus measurements was compared with the CMT obtained by TD Stratus. Results: Cirrus OCT demonstrated high intraobserver repeatability at the central foveal region (ICC 96%). The mean of the CMT measurements was 321microns for Stratus and 387 microns for Cirrus. The average difference was 65m (SD=30). The coefficient of concordance between Stratus and Cirrus CMT measurements was rho=0,749 with a high precision and a moderate accuracy. The equation of the line of regression between Stratus and meanCirrus is given by the following: M_stratus = 0,848 x m_cirrus - 4,496 (1).Conclusions: The Cirrus macular cube protocol allows reproducible CMT measurements in patients with active exudative AMD. In cases of upgrading from TD to SD use and vice versa, there is the possibility to predict the measurements by using the equation (1). These real life data and conclusions can help in improving our clinical management of patients with neovascular AMD.

Inversion of crosshole seismic data in heterogenous environments : Comparison of waveform and ray-based approaches

High-resolution tomographic imaging of the shallow subsurface is becoming increasingly important for a wide range of environmental, hydrological and engineering applications. Because of their superior resolution power, their sensitivity to pertinent petrophysical parameters, and their far reaching complementarities, both seismic and georadar crosshole imaging are of particular importance. To date, corresponding approaches have largely relied on asymptotic, ray-based approaches, which only account for a very small part of the observed wavefields, inherently suffer from a limited resolution, and in complex environments may prove to be inadequate. These problems can potentially be alleviated through waveform inversion. We have developed an acoustic waveform inversion approach for crosshole seismic data whose kernel is based on a finite-difference time-domain (FDTD) solution of the 2-D acoustic wave equations. This algorithm is tested on and applied to synthetic data from seismic velocity models of increasing complexity and realism and the results are compared to those obtained using state-of-the-art ray-based traveltime tomography. Regardless of the heterogeneity of the underlying models, the waveform inversion approach has the potential of reliably resolving both the geometry and the acoustic properties of features of the size of less than half a dominant wavelength. Our results do, however, also indicate that, within their inherent resolution limits, ray-based approaches provide an effective and efficient means to obtain satisfactory tomographic reconstructions of the seismic velocity structure in the presence of mild to moderate heterogeneity and in absence of strong scattering. Conversely, the excess effort of waveform inversion provides the greatest benefits for the most heterogeneous, and arguably most realistic, environments where multiple scattering effects tend to be prevalent and ray-based methods lose most of their effectiveness.

Source wavelet estimation for waveform inversion of crosshole georadar data

AbstractFor a wide range of environmental, hydrological, and engineering applications there is a fast growing need for high-resolution imaging. In this context, waveform tomographic imaging of crosshole georadar data is a powerful method able to provide images of pertinent electrical properties in near-surface environments with unprecedented spatial resolution. In contrast, conventional ray-based tomographic methods, which consider only a very limited part of the recorded signal (first-arrival traveltimes and maximum first-cycle amplitudes), suffer from inherent limitations in resolution and may prove to be inadequate in complex environments. For a typical crosshole georadar survey the potential improvement in resolution when using waveform-based approaches instead of ray-based approaches is in the range of one order-of- magnitude. Moreover, the spatial resolution of waveform-based inversions is comparable to that of common logging methods. While in exploration seismology waveform tomographic imaging has become well established over the past two decades, it is comparably still underdeveloped in the georadar domain despite corresponding needs. Recently, different groups have presented finite-difference time-domain waveform inversion schemes for crosshole georadar data, which are adaptations and extensions of Tarantola's seminal nonlinear generalized least-squares approach developed for the seismic case. First applications of these new crosshole georadar waveform inversion schemes on synthetic and field data have shown promising results. However, there is little known about the limits and performance of such schemes in complex environments. To this end, the general motivation of my thesis is the evaluation of the robustness and limitations of waveform inversion algorithms for crosshole georadar data in order to apply such schemes to a wide range of real world problems.One crucial issue to making applicable and effective any waveform scheme to real-world crosshole georadar problems is the accurate estimation of the source wavelet, which is unknown in reality. Waveform inversion schemes for crosshole georadar data require forward simulations of the wavefield in order to iteratively solve the inverse problem. Therefore, accurate knowledge of the source wavelet is critically important for successful application of such schemes. Relatively small differences in the estimated source wavelet shape can lead to large differences in the resulting tomograms. In the first part of my thesis, I explore the viability and robustness of a relatively simple iterative deconvolution technique that incorporates the estimation of the source wavelet into the waveform inversion procedure rather than adding additional model parameters into the inversion problem. Extensive tests indicate that this source wavelet estimation technique is simple yet effective, and is able to provide remarkably accurate and robust estimates of the source wavelet in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity as well as significant ambient noise in the recorded data. Furthermore, our tests also indicate that the approach is insensitive to the phase characteristics of the starting wavelet, which is not the case when directly incorporating the wavelet estimation into the inverse problem.Another critical issue with crosshole georadar waveform inversion schemes which clearly needs to be investigated is the consequence of the common assumption of frequency- independent electromagnetic constitutive parameters. This is crucial since in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behaviour. In particular, in the presence of water, there is a wide body of evidence showing that the dielectric permittivity can be significantly frequency dependent over the GPR frequency range, due to a variety of relaxation processes. The second part of my thesis is therefore dedicated to the evaluation of the reconstruction limits of a non-dispersive crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. I show that the inversion algorithm, combined with the iterative deconvolution-based source wavelet estimation procedure that is partially able to account for the frequency-dependent effects through an "effective" wavelet, performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.

Analysis of an iterative deconvolution approach for estimating the source wavelet during waveform inversion of crosshole ground-penetrating radar data

A major issue in the application of waveform inversion methods to crosshole georadar data is the accurate estimation of the source wavelet. Here, we explore the viability and robustness of incorporating this step into a time-domain waveform inversion procedure through an iterative deconvolution approach. Our results indicate that, at least in non-dispersive electrical environments, such an approach provides remarkably accurate and robust estimates of the source wavelet even in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity. Our results also indicate that the proposed source wavelet estimation approach is relatively insensitive to ambient noise and to the phase characteristics of the starting wavelet. Finally, there appears to be little-to-no trade-off between the wavelet estimation and the tomographic imaging procedures.

A wearable system to time gate crossing during alpine skiing slalom

INTRODUCTION: In alpine skiing, chronometry analysis is currently the most common tool to assess performance. It is widely used to rank competitors during races, as well as to manage athletes training and to evaluate material. Usually, this measurement is accurately realized using timing cells. Nevertheless, these devices are too complex and expensive to allow chronometry of every gates crossing. On the other side, differential GPS can be used for measuring gate crossing time (Waegli et al). However, this is complex (e.g. recording gate position with GPS) and mainly used in research applications. The aim of the study was to propose a wearable system to time gates crossing during alpine skiing slalom (SL), which is suitable for routine uses. METHODS: The proposed system was composed of a 3D accelerometer (ADXL320®, Analog Device, USA) placed at the sacrum of the athlete, a matrix of force sensors (Flexiforce®, Tekscan, USA) fixed on the right shin guard and a data logger (Physilog®, BioAGM, Switzerland). The sensors were sampled at 500 Hz. The crossing time were calculated in two phases. First, the accelerometer was used to detect the curves by considering the maximum of the mediolateral peak acceleration. Then, the force sensors were used to detect the impacts with the gates by considering maximum force variation. In case of non impact, the detection was realized based on the acceleration and features measured at the other gates. In order to assess the efficiency of the system, two different SL were monitored twice for two world cup level skiers, a male SL expert and a female downhill expert. RESULTS AND DISCUSSION: The combination of the accelerometer and force sensors allowed to clearly identify the gate crossing times. When comparing the runs of the SL expert and the downhill expert, we noticed that the SL expert was faster. For example for the first SL, the overall difference between the best run of each athlete was of 5.47s. At each gate, the SL expert increased the time difference slower at the beginning (0.27s/gate) than at the end (0.34s/gate). Furthermore, when comparing the runs of the SL expert, a maximum time difference of 20ms at each gate was noticed. This showed high repeatability skills of the SL expert. In opposite, the downhill expert with a maximum difference time of 1s at each gate was clearly less repeatable. Both skiers were not disturbed by the system. CONCLUSION: This study proposed a new wearable system to automatically time gates crossing during alpine skiing slalom combining force and accelerometer sensors. The system was evaluated with two professional world cup skiers and showed a high potential. This system could be extended to time other parameters. REFERENCES Waegli A, Skaloud J (2007). Inside GNSS, Spring, 24-34.

Which prior knowledge? Quantification of in vivo brain 13C MR spectra following 13C glucose infusion using AMARES.

The recent developments in high magnetic field 13C magnetic resonance spectroscopy with improved localization and shimming techniques have led to important gains in sensitivity and spectral resolution of 13C in vivo spectra in the rodent brain, enabling the separation of several 13C isotopomers of glutamate and glutamine. In this context, the assumptions used in spectral quantification might have a significant impact on the determination of the 13C concentrations and the related metabolic fluxes. In this study, the time domain spectral quantification algorithm AMARES (advanced method for accurate, robust and efficient spectral fitting) was applied to 13 C magnetic resonance spectroscopy spectra acquired in the rat brain at 9.4 T, following infusion of [1,6-(13)C2 ] glucose. Using both Monte Carlo simulations and in vivo data, the goal of this work was: (1) to validate the quantification of in vivo 13C isotopomers using AMARES; (2) to assess the impact of the prior knowledge on the quantification of in vivo 13C isotopomers using AMARES; (3) to compare AMARES and LCModel (linear combination of model spectra) for the quantification of in vivo 13C spectra. AMARES led to accurate and reliable 13C spectral quantification similar to those obtained using LCModel, when the frequency shifts, J-coupling constants and phase patterns of the different 13C isotopomers were included as prior knowledge in the analysis.

Quantification of brain glycogen concentration and turnover through localized 13C NMR of both the C1 and C6 resonances.

We have recently shown that at isotopic steady state (13)C NMR can provide a direct measurement of glycogen concentration changes, but that the turnover of glycogen was not accessible with this protocol. The aim of the present study was to design, implement and apply a novel dual-tracer infusion protocol to simultaneously measure glycogen concentration and turnover. After reaching isotopic steady state for glycogen C1 using [1-(13)C] glucose administration, [1,6-(13)C(2)] glucose was infused such that isotopic steady state was maintained at the C1 position, but the C6 position reflected (13)C label incorporation. To overcome the large chemical shift displacement error between the C1 and C6 resonances of glycogen, we implemented 2D gradient based localization using the Fourier series window approach, in conjunction with time-domain analysis of the resulting FIDs using jMRUI. The glycogen concentration of 5.1 +/- 1.6 mM measured from the C1 position was in excellent agreement with concomitant biochemical determinations. Glycogen turnover measured from the rate of label incorporation into the C6 position of glycogen in the alpha-chloralose anesthetized rat was 0.7 micromol/g/h.

Heated Optical Fiber for Distributed Soil-Moisture Measurements: A Lysimeter Experiment

An Actively Heated Fiber Optics (AHFO) method to estimate soil moisture is tested and the analysis technique improved on. The measurements were performed in a lysimeter uniformly packed with loam soil with variable water content profiles. In the first meter of the soil profi le, 30 m of fiber optic cable were installed in a 12 loops coil. The metal sheath armoring the fiber cable was used as an electrical resistance heater to generate a heat pulse, and the soil response was monitored with a Distributed Temperature Sensing (DTS) system. We study the cooling following three continuous heat pulses of 120 s at 36 W m(-1) by means of long-time approximation of radial heat conduction. The soil volumetric water contents were then inferred from the estimated thermal conductivities through a specifically calibrated model relating thermal conductivity and volumetric water content. To use the pre-asymptotic data we employed a time correction that allowed the volumetric water content to be estimated with a precision of 0.01-0.035 (m(3) m(-3)). A comparison of the AHFO measurements with soil-moisture measurements obtained with calibrated capacitance-based probes gave good agreement for wetter soils [discrepancy between the two methods was less than 0.04 (m(3) m(-3))]. In the shallow drier soils, the AHFO method underestimated the volumetric water content due to the longertime required for the temperature increment to become asymptotic in less thermally conductive media [discrepancy between the two methods was larger than 0.1 (m(3) m(-3))]. The present work suggests that future applications of the AHFO method should include longer heat pulses, that longer heating and cooling events are analyzed, and, temperature increments ideally be measured with higher frequency.

Source-wavelet estimation during full-waveform inversion of ground-penetrating radar data

A major issue in the application of waveform inversion methods to crosshole ground-penetrating radar (GPR) data is the accurate estimation of the source wavelet. Here, we explore the viability and robustness of incorporating this step into a recently published time-domain inversion procedure through an iterative deconvolution approach. Our results indicate that, at least in non-dispersive electrical environments, such an approach provides remarkably accurate and robust estimates of the source wavelet even in the presence of strong heterogeneity of both the dielectric permittivity and electrical conductivity. Our results also indicate that the proposed source wavelet estimation approach is relatively insensitive to ambient noise and to the phase characteristics of the starting wavelet. Finally, there appears to be little to no trade-off between the wavelet estimation and the tomographic imaging procedures.

Estimating vadose zone hydraulic properties using ground penetrating radar: The impact of prior information

A number of geophysical methods, such as ground-penetrating radar (GPR), have the potential to provide valuable information on hydrological properties in the unsaturated zone. In particular, the stochastic inversion of such data within a coupled geophysical-hydrological framework may allow for the effective estimation of vadose zone hydraulic parameters and their corresponding uncertainties. A critical issue in stochastic inversion is choosing prior parameter probability distributions from which potential model configurations are drawn and tested against observed data. A well chosen prior should reflect as honestly as possible the initial state of knowledge regarding the parameters and be neither overly specific nor too conservative. In a Bayesian context, combining the prior with available data yields a posterior state of knowledge about the parameters, which can then be used statistically for predictions and risk assessment. Here we investigate the influence of prior information regarding the van Genuchten-Mualem (VGM) parameters, which describe vadose zone hydraulic properties, on the stochastic inversion of crosshole GPR data collected under steady state, natural-loading conditions. We do this using a Bayesian Markov chain Monte Carlo (MCMC) inversion approach, considering first noninformative uniform prior distributions and then more informative priors derived from soil property databases. For the informative priors, we further explore the effect of including information regarding parameter correlation. Analysis of both synthetic and field data indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when we combine these data with a realistic, informative prior.

Full-waveform inversion of cross-hole ground-penetrating radar data to characterize a gravel aquifer close to the Thur River, Switzerland

Cross-hole radar tomography is a useful tool for mapping shallow subsurface electrical properties viz. dielectric permittivity and electrical conductivity. Common practice is to invert cross-hole radar data with ray-based tomographic algorithms using first arrival traveltimes and first cycle amplitudes. However, the resolution of conventional standard ray-based inversion schemes for cross-hole ground-penetrating radar (GPR) is limited because only a fraction of the information contained in the radar data is used. The resolution can be improved significantly by using a full-waveform inversion that considers the entire waveform, or significant parts thereof. A recently developed 2D time-domain vectorial full-waveform crosshole radar inversion code has been modified in the present study by allowing optimized acquisition setups that reduce the acquisition time and computational costs significantly. This is achieved by minimizing the number of transmitter points and maximizing the number of receiver positions. The improved algorithm was employed to invert cross-hole GPR data acquired within a gravel aquifer (4-10 m depth) in the Thur valley, Switzerland. The simulated traces of the final model obtained by the full-waveform inversion fit the observed traces very well in the lower part of the section and reasonably well in the upper part of the section. Compared to the ray-based inversion, the results from the full-waveform inversion show significantly higher resolution images. At either side, 2.5 m distance away from the cross-hole plane, borehole logs were acquired. There is a good correspondence between the conductivity tomograms and the natural gamma logs at the boundary of the gravel layer and the underlying lacustrine clay deposits. Using existing petrophysical models, the inversion results and neutron-neutron logs are converted to porosity. Without any additional calibration, the values obtained for the converted neutron-neutron logs and permittivity results are very close and similar vertical variations can be observed. The full-waveform inversion provides in both cases additional information about the subsurface. Due to the presence of the water table and associated refracted/reflected waves, the upper traces are not well fitted and the upper 2 m in the permittivity and conductivity tomograms are not reliably reconstructed because the unsaturated zone is not incorporated into the inversion domain.

On the diurnal soil water content dynamics during evaporation using dielectric methods

The water content dynamics in the upper soil surface during evaporation is a key element in land-atmosphere exchanges. Previous experimental studies have suggested that the soil water content increases at the depth of 5 to 15 cm below the soil surface during evapo- ration, while the layer in the immediate vicinity of the soil surface is drying. In this study, the dynamics of water content profiles exposed to solar radiative forcing was monitored at a high temporal resolution using dielectric methods both in the presence and absence of evaporation. A 4-d comparison of reported moisture content in coarse sand in covered and uncovered buckets using a commercial dielectric-based probe (70 MHz ECH2O-5TE, Decagon Devices, Pullman, WA) and the standard 1-GHz time domain reflectometry method. Both sensors reported a positive correlation between temperature and water content in the 5- to 10-cm depth, most pronounced in the morning during heating and in the afternoon during cooling. Such positive correlation might have a physical origin induced by evaporation at the surface and redistribution due to liquid water fluxes resulting from the temperature- gradient dynamics within the sand profile at those depths. Our experimental data suggest that the combined effect of surface evaporation and temperature-gradient dynamics should be considered to analyze experimental soil water profiles. Additional effects related to the frequency of operation and to protocols for temperature compensation of the dielectric sensors may also affect the probes' response during large temperature changes.

Monitoring Fatigue Status with HRV Measures in Elite Athletes : An Avenue Beyond RMSSD?

Among the tools proposed to assess the athlete's "fatigue," the analysis of heart rate variability (HRV) provides an indirect evaluation of the settings of autonomic control of heart activity. HRV analysis is performed through assessment of time-domain indices, the square root of the mean of the sum of the squares of differences between adjacent normal R-R intervals (RMSSD) measured during short (5 min) recordings in supine position upon awakening in the morning and particularly the logarithm of RMSSD (LnRMSSD) has been proposed as the most useful resting HRV indicator. However, if RMSSD can help the practitioner to identify a global "fatigue" level, it does not allow discriminating different types of fatigue. Recent results using spectral HRV analysis highlighted firstly that HRV profiles assessed in supine and standing positions are independent and complementary; and secondly that using these postural profiles allows the clustering of distinct sub-categories of "fatigue." Since, cardiovascular control settings are different in standing and lying posture, using the HRV figures of both postures to cluster fatigue state embeds information on the dynamics of control responses. Such, HRV spectral analysis appears more sensitive and enlightening than time-domain HRV indices. The wealthier information provided by this spectral analysis should improve the monitoring of the adaptive training-recovery process in athletes.

Monitoring Fatigue Status with HRV Measures in Elite Athletes: An Avenue Beyond RMSSD?

Among the tools proposed to assess the athlete's "fatigue," the analysis of heart rate variability (HRV) provides an indirect evaluation of the settings of autonomic control of heart activity. HRV analysis is performed through assessment of time-domain indices, the square root of the mean of the sum of the squares of differences between adjacent normal R-R intervals (RMSSD) measured during short (5 min) recordings in supine position upon awakening in the morning and particularly the logarithm of RMSSD (LnRMSSD) has been proposed as the most useful resting HRV indicator. However, if RMSSD can help the practitioner to identify a global "fatigue" level, it does not allow discriminating different types of fatigue. Recent results using spectral HRV analysis highlighted firstly that HRV profiles assessed in supine and standing positions are independent and complementary; and secondly that using these postural profiles allows the clustering of distinct sub-categories of "fatigue." Since, cardiovascular control settings are different in standing and lying posture, using the HRV figures of both postures to cluster fatigue state embeds information on the dynamics of control responses. Such, HRV spectral analysis appears more sensitive and enlightening than time-domain HRV indices. The wealthier information provided by this spectral analysis should improve the monitoring of the adaptive training-recovery process in athletes.

Increased risks of early sexual initiators: time makes a difference.

BACKGROUND: Years since onset of sexual intercourse (YSSI) is a rarely used variable when studying adolescents- sexual outcomes. The aim of this study is to evaluate the influence of YSSI on the adverse sexual outcomes of early sexual initiators.METHODS: Data were drawn from the 2002 Swiss Multicenter Adolescent Survey on Health database, a nationally representative cross-sectional survey including 7429 adolescents in post mandatory school aged 16-20 years. Only adolescents reporting sexual intercourse (SI) were included (N=4388; 45% females) and divided by age of onset of SI (early initiators, age<16: N=1469, 44% females; and late initiators, age?16: N=2919, 46% females). Analyses were done separately by gender. Groups were compared for personal characteristics at the bivariate level. We analyzed three sexual outcomes (?4 sexual partners, pregnancy and non-use of condom at last SI) controlling for all significant personal variables with two logistic regressions first using age, then YSSI as one of the confounding variables. Results are given as adjusted odds ratios (aOR) using lSI as the reference category.RESULTS: After adjusting for YSSI instead of age, negative sexual outcomes among early initiators were no longer significant, except for multiple sexual partners among females, although at a much lower level. Early initiators were less likely to report non-use of condom at last SI when adjusting for YSSI (females: aOR=0.59 [0.44-0.79]; p<0.001; males aOR=0.71 [0.50-1.00]; p=0.053).CONCLUSION: YSSI is an important explanatory variable when studying adolescents- sexuality and needs to be included in future research on adolescents- sexual health.