A High-Resolution 2D/3D Seismic Study of a Thrust Fault Zone in Lake Geneva, Switzerland

P130 A HIGH-RESOLUTION 2D/3D SEISMIC STUDY OF A THRUST FAULT ZONE IN LAKE GENEVA SWITZERLAND M. SCHEIDHAUER M. BERES D. DUPUY and F. MARILLIER Institute of Geophysics University of Lausanne 1015 Lausanne, Switzerland Summary A high-resolution three-dimensional (3D) seismic reflection survey has been conducted in Lake Geneva near the city of Lausanne Switzerland where the faulted molasse basement (Tertiary sandstones) is overlain by complex Quaternary sedimentary structures. Using a single 48-channel streamer an area of 1200 m x 600 m was surveyed in 10 days. With a 5-m shot spacing and a receiver spacing of 2.5 m in the inline direction and 7.5 m in the crossline direction, a 12-fold data coverage was achieved. A maximum penetration depth of ~150 m was achieved with a 15 cu. in. water gun operated at 140 bars. The multi-channel data allow the determination of an accurate velocity field for 3D processing, and they show particularly clean images of the fault zone and the overlying sediments in horizontal and vertical sections. In order to compare different sources, inline 55 was repeated with a 30/30 and a 15/15 cu. in. double-chamber air gun (Mini GI) operated at 100 and 80 bars, respectively. A maximum penetration depth of ~450 m was achieved with this source.

Model-based Segmentation and Fusion of 3D Computed Tomography and 3D Ultrasound of the Eye for Radiotherapy Planning

Computed Tomography (CT) represents the standard imaging modality for tumor volume delineation for radiotherapy treatment planning of retinoblastoma despite some inherent limitations. CT scan is very useful in providing information on physical density for dose calculation and morphological volumetric information but presents a low sensitivity in assessing the tumor viability. On the other hand, 3D ultrasound (US) allows a highly accurate definition of the tumor volume thanks to its high spatial resolution but it is not currently integrated in the treatment planning but used only for diagnosis and follow-up. Our ultimate goal is an automatic segmentation of gross tumor volume (GTV) in the 3D US, the segmentation of the organs at risk (OAR) in the CT and the registration of both modalities. In this paper, we present some preliminary results in this direction. We present 3D active contour-based segmentation of the eye ball and the lens in CT images; the presented approach incorporates the prior knowledge of the anatomy by using a 3D geometrical eye model. The automated segmentation results are validated by comparing with manual segmentations. Then, we present two approaches for the fusion of 3D CT and US images: (i) landmark-based transformation, and (ii) object-based transformation that makes use of eye ball contour information on CT and US images.

Quantitative analysis of crystal/grain sizes and their distributions in 2D and 3D

We review methods to estimate the average crystal (grain) size and the crystal (grain) size distribution in solid rocks. Average grain sizes often provide the base for stress estimates or rheological calculations requiring the quantification of grain sizes in a rock's microstructure. The primary data for grain size data are either 1D (i.e. line intercept methods), 2D (area analysis) or 3D (e.g., computed tomography, serial sectioning). These data have been used for different data treatments over the years, whereas several studies assume a certain probability function (e.g., logarithm, square root) to calculate statistical parameters as the mean, median, mode or the skewness of a crystal size distribution. The finally calculated average grain sizes have to be compatible between the different grain size estimation approaches in order to be properly applied, for example, in paleo-piezometers or grain size sensitive flow laws. Such compatibility is tested for different data treatments using one- and two-dimensional measurements. We propose an empirical conversion matrix for different datasets. These conversion factors provide the option to make different datasets compatible with each other, although the primary calculations were obtained in different ways. In order to present an average grain size, we propose to use the area-weighted and volume-weighted mean in the case of unimodal grain size distributions, respectively, for 2D and 3D measurements. The shape of the crystal size distribution is important for studies of nucleation and growth of minerals. The shape of the crystal size distribution of garnet populations is compared between different 2D and 3D measurements, which are serial sectioning and computed tomography. The comparison of different direct measured 3D data; stereological data and direct presented 20 data show the problems of the quality of the smallest grain sizes and the overestimation of small grain sizes in stereological tools, depending on the type of CSD. (C) 2011 Published by Elsevier Ltd.

Direct Fourier Tomographic Reconstruction Image-To-Image Filter

We present an open-source ITK implementation of a directFourier method for tomographic reconstruction, applicableto parallel-beam x-ray images. Direct Fourierreconstruction makes use of the central-slice theorem tobuild a polar 2D Fourier space from the 1D transformedprojections of the scanned object, that is resampled intoa Cartesian grid. Inverse 2D Fourier transform eventuallyyields the reconstructed image. Additionally, we providea complex wrapper to the BSplineInterpolateImageFunctionto overcome ITKâeuro?s current lack for image interpolatorsdealing with complex data types. A sample application ispresented and extensively illustrated on the Shepp-Loganhead phantom. We show that appropriate input zeropaddingand 2D-DFT oversampling rates together with radial cubicb-spline interpolation improve 2D-DFT interpolationquality and are efficient remedies to reducereconstruction artifacts.

A fast Bayesian reconstruction algorithm for emission tomography with entropy prior converging to feasible images

The development and tests of an iterative reconstruction algorithm for emission tomography based on Bayesian statistical concepts are described. The algorithm uses the entropy of the generated image as a prior distribution, can be accelerated by the choice of an exponent, and converges uniformly to feasible images by the choice of one adjustable parameter. A feasible image has been defined as one that is consistent with the initial data (i.e. it is an image that, if truly a source of radiation in a patient, could have generated the initial data by the Poisson process that governs radioactive disintegration). The fundamental ideas of Bayesian reconstruction are discussed, along with the use of an entropy prior with an adjustable contrast parameter, the use of likelihood with data increment parameters as conditional probability, and the development of the new fast maximum a posteriori with entropy (FMAPE) Algorithm by the successive substitution method. It is shown that in the maximum likelihood estimator (MLE) and FMAPE algorithms, the only correct choice of initial image for the iterative procedure in the absence of a priori knowledge about the image configuration is a uniform field.

Bayesian image reconstruction with space-variant noise suppression

In this paper we present a Bayesian image reconstruction algorithm with entropy prior (FMAPE) that uses a space-variant hyperparameter. The spatial variation of the hyperparameter allows different degrees of resolution in areas of different statistical characteristics, thus avoiding the large residuals resulting from algorithms that use a constant hyperparameter. In the first implementation of the algorithm, we begin by segmenting a Maximum Likelihood Estimator (MLE) reconstruction. The segmentation method is based on using a wavelet decomposition and a self-organizing neural network. The result is a predetermined number of extended regions plus a small region for each star or bright object. To assign a different value of the hyperparameter to each extended region and star, we use either feasibility tests or cross-validation methods. Once the set of hyperparameters is obtained, we carried out the final Bayesian reconstruction, leading to a reconstruction with decreased bias and excellent visual characteristics. The method has been applied to data from the non-refurbished Hubble Space Telescope. The method can be also applied to ground-based images.

Super-resolution of remotely sensed images with variable-pixel linear reconstruction

This paper describes the development and applications of a super-resolution method, known as Super-Resolution Variable-Pixel Linear Reconstruction. The algorithm works combining different lower resolution images in order to obtain, as a result, a higher resolution image. We show that it can make significant spatial resolution improvements to satellite images of the Earth¿s surface allowing recognition of objects with size approaching the limiting spatial resolution of the lower resolution images. The algorithm is based on the Variable-Pixel Linear Reconstruction algorithm developed by Fruchter and Hook, a well-known method in astronomy but never used for Earth remote sensing purposes. The algorithm preserves photometry, can weight input images according to the statistical significance of each pixel, and removes the effect of geometric distortion on both image shape and photometry. In this paper, we describe its development for remote sensing purposes, show the usefulness of the algorithm working with images as different to the astronomical images as the remote sensing ones, and show applications to: 1) a set of simulated multispectral images obtained from a real Quickbird image; and 2) a set of multispectral real Landsat Enhanced Thematic Mapper Plus (ETM+) images. These examples show that the algorithm provides a substantial improvement in limiting spatial resolution for both simulated and real data sets without significantly altering the multispectral content of the input low-resolution images, without amplifying the noise, and with very few artifacts.

Three-dimensional ordered-subset expectation maximization iterative protocol for evaluation of left ventricular volumes and function by quantitative gated SPECT: a dynamic phantom study.

The purposes of this study were to characterize the performance of a 3-dimensional (3D) ordered-subset expectation maximization (OSEM) algorithm in the quantification of left ventricular (LV) function with (99m)Tc-labeled agent gated SPECT (G-SPECT), the QGS program, and a beating-heart phantom and to optimize the reconstruction parameters for clinical applications. METHODS: A G-SPECT image of a dynamic heart phantom simulating the beating left ventricle was acquired. The exact volumes of the phantom were known and were as follows: end-diastolic volume (EDV) of 112 mL, end-systolic volume (ESV) of 37 mL, and stroke volume (SV) of 75 mL; these volumes produced an LV ejection fraction (LVEF) of 67%. Tomographic reconstructions were obtained after 10-20 iterations (I) with 4, 8, and 16 subsets (S) at full width at half maximum (FWHM) gaussian postprocessing filter cutoff values of 8-15 mm. The QGS program was used for quantitative measurements. RESULTS: Measured values ranged from 72 to 92 mL for EDV, from 18 to 32 mL for ESV, and from 54 to 63 mL for SV, and the calculated LVEF ranged from 65% to 76%. Overall, the combination of 10 I, 8 S, and a cutoff filter value of 10 mm produced the most accurate results. The plot of the measures with respect to the expectation maximization-equivalent iterations (I x S product) revealed a bell-shaped curve for the LV volumes and a reverse distribution for the LVEF, with the best results in the intermediate range. In particular, FWHM cutoff values exceeding 10 mm affected the estimation of the LV volumes. CONCLUSION: The QGS program is able to correctly calculate the LVEF when used in association with an optimized 3D OSEM algorithm (8 S, 10 I, and FWHM of 10 mm) but underestimates the LV volumes. However, various combinations of technical parameters, including a limited range of I and S (80-160 expectation maximization-equivalent iterations) and low cutoff values (< or =10 mm) for the gaussian postprocessing filter, produced results with similar accuracies and without clinically relevant differences in the LV volumes and the estimated LVEF.

A filtering method to correct time-lapse 3D ERT data and improve imaging of natural aquifer dynamics

We have developed a processing methodology that allows crosshole ERT (electrical resistivity tomography) monitoring data to be used to derive temporal fluctuations of groundwater electrical resistivity and thereby characterize the dynamics of groundwater in a gravel aquifer as it is infiltrated by river water. Temporal variations of the raw ERT apparent-resistivity data were mainly sensitive to the resistivity (salinity), temperature and height of the groundwater, with the relative contributions of these effects depending on the time and the electrode configuration. To resolve the changes in groundwater resistivity, we first expressed fluctuations of temperature-detrended apparent-resistivity data as linear superpositions of (i) time series of riverwater-resistivity variations convolved with suitable filter functions and (ii) linear and quadratic representations of river-water-height variations multiplied by appropriate sensitivity factors; river-water height was determined to be a reliable proxy for groundwater height. Individual filter functions and sensitivity factors were obtained for each electrode configuration via deconvolution using a one month calibration period and then the predicted contributions related to changes in water height were removed prior to inversion of the temperature-detrended apparent-resistivity data. Applications of the filter functions and sensitivity factors accurately predicted the apparent-resistivity variations (the correlation coefficient was 0.98). Furthermore, the filtered ERT monitoring data and resultant time-lapse resistivity models correlated closely with independently measured groundwater electrical resistivity monitoring data and only weakly with the groundwater-height fluctuations. The inversion results based on the filtered ERT data also showed significantly less inversion artefacts than the raw data inversions. We observed resistivity increases of up to 10% and the arrival time peaks in the time-lapse resistivity models matched those in the groundwater resistivity monitoring data.

Reproducibility of 3D free-breathing magnetic resonance coronary vessel wall imaging.

AIMS: Although the coronary artery vessel wall can be imaged non-invasively using magnetic resonance imaging (MRI), the in vivo reproducibility of wall thickness measures has not been previously investigated. Using a refined magnetization preparation scheme, we sought to assess the reproducibility of three-dimensional (3D) free-breathing black-blood coronary MRI in vivo. METHODS AND RESULTS: MRI vessel wall scans parallel to the right coronary artery (RCA) were obtained in 18 healthy individuals (age range 25-43, six women), with no known history of coronary artery disease, using a 3D dual-inversion navigator-gated black-blood spiral imaging sequence. Vessel wall scans were repeated 1 month later in eight subjects. The visible vessel wall segment and the wall thickness were quantitatively assessed using a semi-automatic tool and the intra-observer, inter-observer, and inter-scan reproducibilities were determined. The average imaged length of the RCA vessel wall was 44.5+/-7 mm. The average wall thickness was 1.6+/-0.2 mm. There was a highly significant intra-observer (r=0.97), inter-observer (r=0.94), and inter-scan (r=0.90) correlation for wall thickness (all P<0.001). There was also a significant agreement for intra-observer, inter-observer, and inter-scan measurements on Bland-Altman analysis. The intra-class correlation coefficients for intra-observer (r=0.97), inter-observer (r=0.92), and inter-scan (r=0.86) analyses were also excellent. CONCLUSION: The use of black-blood free-breathing 3D MRI in conjunction with semi-automated analysis software allows for reproducible measurements of right coronary arterial vessel-wall thickness. This technique may be well-suited for non-invasive longitudinal studies of coronary atherosclerosis.

Revision anterior cruciate ligament reconstruction with quadriceps tendon-patellar bone autograft.

PURPOSE: To evaluate the cause of recurrent pathologic instability after anterior cruciate ligament (ACL) surgery and the effectiveness of revision reconstruction using a quadriceps tendon autograft using a 2-incision technique. TYPE OF STUDY: Retrospective follow-up study. METHODS: Between 1999 and 2001, 31 patients underwent ACL revision reconstruction because of recurrent pathologic instability during sports or daily activities. Twenty-eight patients were reviewed after a mean follow-up of 4.2 years (range, 3.3 to 5.6 years). The mean age at revision surgery was 27 years (range, 18 to 41 years). The average time from primary procedure to revision surgery was 26 months (range, 9 to 45 months). A clinical, functional, and radiographic evaluation was performed. Also magnetic resonance imaging (MRI) or computed tomography (CT) scanning was performed. The International Knee Documentation Committee (IKDC), Lysholm, and Tegner scales were used. A KT-1000 arthrometer measurement (MEDmetric, San Diego, CA) by an experienced physician was made. RESULTS: Of the failures, 79% had radiographic evidence of malposition of their tunnels. In only 6 cases (21%) was the radiologic anatomy of tunnel placement judged to be correct on both the femoral and tibial side. The MRI or CT showed, in 6 cases, a too-centrally placed femoral tunnel. After revision surgery, the position of tunnels was corrected. A significant improvement of Lachman and pivot-shift phenomenon was observed. In particular, 17 patients had a negative Lachman test, and 11 patients had a grade I Lachman with a firm end point. Preoperatively, the pivot-shift test was positive in all cases, and at last follow-up in 7 patients (25%) a grade 1+ was found. Postoperatively, KT-1000 testing showed a mean manual maximum translation of 8.6 mm (SD, 2.34) for the affected knee; 97% of patients had a maximum manual side-to-side translation <5 mm. At the final postoperative evaluation, 26 patients (93%) graded their knees as normal or nearly normal according to the IKDC score. The mean Lysholm score was 93.6 (SD, 8.77) and the mean Tegner activity score was 6.1 (SD, 1.37). No patient required further revision. Five patients (18%) complained of hypersensitive scars from the reconstructive surgery that made kneeling difficult. CONCLUSIONS: There were satisfactory results after ACL revision surgery using quadriceps tendon and a 2-incision technique at a minimum 3 years' follow-up; 93% of patients returned to sports activities. LEVEL OF EVIDENCE: Level IV, case series, no control group.

Aortic vessel wall magnetic resonance imaging at 3.0 Tesla: a reproducibility study of respiratory navigator gated free-breathing 3D black blood magnetic resonance imaging.

The purpose of this study was to evaluate a free-breathing three-dimensional (3D) dual inversion-recovery (DIR) segmented k-space gradient-echo (turbo field echo [TFE]) imaging sequence at 3T for the quantification of aortic vessel wall dimensions. The effect of respiratory motion suppression on image quality was tested. Furthermore, the reproducibility of the aortic vessel wall measurements was investigated. Seven healthy subjects underwent 3D DIR TFE imaging of the aortic vessel wall with and without respiratory navigator. Subsequently, this sequence with respiratory navigator was performed twice in 10 healthy subjects to test its reproducibility. The signal-to-noise (SNR), contrast-to-noise ratio (CNR), vessel wall sharpness, and vessel wall volume (VWV) were assessed. Data were compared using the paired t-test, and the reproducibility of VWV measurements was evaluated using intraclass correlation coefficients (ICCs). SNR, CNR, and vessel wall sharpness were superior in scans performed with respiratory navigator compared to scans performed without. The ICCs concerning intraobserver, interobserver, and interscan reproducibility were excellent (0.99, 0.94, and 0.95, respectively). In conclusion, respiratory motion suppression substantially improves image quality of 3D DIR TFE imaging of the aortic vessel wall at 3T. Furthermore, this optimized technique with respiratory motion suppression enables assessment of aortic vessel wall dimensions with high reproducibility.

Correction for heart rate variability during 3D whole heart MR coronary angiography.

PURPOSE: To evaluate the effect of a real-time adaptive trigger delay on image quality to correct for heart rate variability in 3D whole-heart coronary MR angiography (MRA). MATERIALS AND METHODS: Twelve healthy adults underwent 3D whole-heart coronary MRA with and without the use of an adaptive trigger delay. The moment of minimal coronary artery motion was visually determined on a high temporal resolution MRI. Throughout the scan performed without adaptive trigger delay, trigger delay was kept constant, whereas during the scan performed with adaptive trigger delay, trigger delay was continuously updated after each RR-interval using physiological modeling. Signal-to-noise, contrast-to-noise, vessel length, vessel sharpness, and subjective image quality were compared in a blinded manner. RESULTS: Vessel sharpness improved significantly for the middle segment of the right coronary artery (RCA) with the use of the adaptive trigger delay (52.3 +/- 7.1% versus 48.9 +/- 7.9%, P = 0.026). Subjective image quality was significantly better in the middle segments of the RCA and left anterior descending artery (LAD) when the scan was performed with adaptive trigger delay compared to constant trigger delay. CONCLUSION: Our results demonstrate that the use of an adaptive trigger delay to correct for heart rate variability improves image quality mainly in the middle segments of the RCA and LAD.

Whole-body 3D T1-weighted MR imaging in patients with prostate cancer: feasibility and evaluation in screening for metastatic disease.

PURPOSE: To develop and assess the diagnostic performance of a three-dimensional (3D) whole-body T1-weighted magnetic resonance (MR) imaging pulse sequence at 3.0 T for bone and node staging in patients with prostate cancer. MATERIALS AND METHODS This prospective study was approved by the institutional ethics committee; informed consent was obtained from all patients. Thirty patients with prostate cancer at high risk for metastases underwent whole-body 3D T1-weighted imaging in addition to the routine MR imaging protocol for node and/or bone metastasis screening, which included coronal two-dimensional (2D) whole-body T1-weighted MR imaging, sagittal proton-density fat-saturated (PDFS) imaging of the spine, and whole-body diffusion-weighted MR imaging. Two observers read the 2D and 3D images separately in a blinded manner for bone and node screening. Images were read in random order. The consensus review of MR images and the findings at prospective clinical and MR imaging follow-up at 6 months were used as the standard of reference. The interobserver agreement and diagnostic performance of each sequence were assessed on per-patient and per-lesion bases. RESULTS: The signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) were significantly higher with whole-body 3D T1-weighted imaging than with whole-body 2D T1-weighted imaging regardless of the reference region (bone or fat) and lesion location (bone or node) (P < .003 for all). For node metastasis, diagnostic performance (area under the receiver operating characteristic curve) was higher for whole-body 3D T1-weighted imaging (per-patient analysis; observer 1: P < .001 for 2D T1-weighted imaging vs 3D T1-weighted imaging, P = .006 for 2D T1-weighted imaging + PDFS imaging vs 3D T1-weighted imaging; observer 2: P = .006 for 2D T1-weighted imaging vs 3D T1-weighted imaging, P = .006 for 2D T1-weighted imaging + PDFS imaging vs 3D T1-weighted imaging), as was sensitivity (per-lesion analysis; observer 1: P < .001 for 2D T1-weighted imaging vs 3D T1-weighted imaging, P < .001 for 2D T1-weighted imaging + PDFS imaging vs 3D T1-weighted imaging; observer 2: P < .001 for 2D T1-weighted imaging vs 3D T1-weighted imaging, P < .001 for 2D T1-weighted imaging + PDFS imaging vs 3D T1-weighted imaging). CONCLUSION: Whole-body MR imaging is feasible with a 3D T1-weighted sequence and provides better SNR and CNR compared with 2D sequences, with a diagnostic performance that is as good or better for the detection of bone metastases and better for the detection of lymph node metastases.