skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Hydraulic Inverse Modeling Using Total-Variation Regularization with Relaxed Variable-Splitting

ORCiD logo [1]; ORCiD logo [1]; ORCiD logo [1]; ORCiD logo [1]
  1. Los Alamos National Laboratory
Publication Date:
Research Org.:
Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
Sponsoring Org.:
OSTI Identifier:
Report Number(s):
DOE Contract Number:
Resource Type:
Resource Relation:
Conference: SIAM Computational Science and Engineering ; 2017-02-27 - 2017-02-27 ; Atlanta, Georgia, United States
Country of Publication:
United States
Earth Sciences

Citation Formats

Lin, Youzuo, Vesselinov, Velimir Valentinov, O'Malley, Daniel, and Wohlberg, Brendt Egon. Hydraulic Inverse Modeling Using Total-Variation Regularization with Relaxed Variable-Splitting. United States: N. p., 2017. Web.
Lin, Youzuo, Vesselinov, Velimir Valentinov, O'Malley, Daniel, & Wohlberg, Brendt Egon. Hydraulic Inverse Modeling Using Total-Variation Regularization with Relaxed Variable-Splitting. United States.
Lin, Youzuo, Vesselinov, Velimir Valentinov, O'Malley, Daniel, and Wohlberg, Brendt Egon. Wed . "Hydraulic Inverse Modeling Using Total-Variation Regularization with Relaxed Variable-Splitting". United States. doi:.
title = {Hydraulic Inverse Modeling Using Total-Variation Regularization with Relaxed Variable-Splitting},
author = {Lin, Youzuo and Vesselinov, Velimir Valentinov and O'Malley, Daniel and Wohlberg, Brendt Egon},
abstractNote = {},
doi = {},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Wed Mar 01 00:00:00 EST 2017},
month = {Wed Mar 01 00:00:00 EST 2017}

Other availability
Please see Document Availability for additional information on obtaining the full-text document. Library patrons may search WorldCat to identify libraries that hold this conference proceeding.

Save / Share:
  • Purpose: Radiation therapy with high dose rate and flattening filter-free (FFF) beams has the potential advantage of greatly reduced treatment time and out-of-field dose. Current inverse planning algorithms are, however, not customized for beams with nonuniform incident profiles and the resultant IMRT plans are often inefficient in delivery. The authors propose a total-variation regularization (TVR)-based formalism by taking the inherent shapes of incident beam profiles into account. Methods: A novel TVR-based inverse planning formalism is established for IMRT with nonuniform beam profiles. The authors introduce a TVR term into the objective function, which encourages piecewise constant fluence in the nonuniformmore » FFF fluence domain. The proposed algorithm is applied to lung and prostate and head and neck cases and its performance is evaluated by comparing the resulting plans to those obtained using a conventional beamlet-based optimization (BBO). Results: For the prostate case, the authors' algorithm produces acceptable dose distributions with only 21 segments, while the conventional BBO requires 114 segments. For the lung case and the head and neck case, the proposed method generates similar coverage of target volume and sparing of the organs-at-risk as compared to BBO, but with a markedly reduced segment number. Conclusions: TVR-based optimization in nonflat beam domain provides an effective way to maximally leverage the technical capacity of radiation therapy with FFF fields. The technique can generate effective IMRT plans with improved dose delivery efficiency without significant deterioration of the dose distribution.« less
  • Double-difference waveform inversion is a promising tool for quantitative monitoring for enhanced geothermal systems (EGS). The method uses time-lapse seismic data to jointly inverts for reservoir changes. Due to the ill-posedness of waveform inversion, it is a great challenge to obtain reservoir changes accurately and efficiently, particularly when using timelapse seismic reflection data. To improve reconstruction, we develop a spatially-variant total-variation regularization scheme into double-difference waveform inversion to improve the inversion accuracy and robustness. The new regularization scheme employs different regularization parameters in different regions of the model to obtain an optimal regularization in each area. We compare the resultsmore » obtained using a spatially-variant parameter with those obtained using a constant regularization parameter. Utilizing a spatially-variant regularization scheme, the target monitoring regions are well reconstructed and the image noise is significantly reduced outside the monitoring regions. Our numerical examples demonstrate that the spatially-variant total-variation regularization scheme provides the flexibility to regularize local regions based on the a priori spatial information without increasing computational costs and the computer memory requirement.« less
  • Breast ultrasound tomography is a rapidly developing imaging modality that has the potential to impact breast cancer screening and diagnosis. A new ultrasound breast imaging device (CURE) with a ring array of transducers has been designed and built at Karmanos Cancer Institute, which acquires both reflection and transmission ultrasound signals. To extract the sound-speed information from the breast data acquired by CURE, we have developed an iterative sound-speed image reconstruction algorithm for breast ultrasound transmission tomography based on total-variation (TV) minimization. We investigate applicability of the TV tomography algorithm using in vivo ultrasound breast data from 61 patients, and comparemore » the results with those obtained using the Tikhonov regularization method. We demonstrate that, compared to the Tikhonov regularization scheme, the TV regularization method significantly improves image quality, resulting in sound-speed tomography images with sharp (preserved) edges of abnormalities and few artifacts.« less
  • We present a new algorithm for bound-constrained total-variation (TV) regularization that in comparison with its predecessors is simple, fast, and flexible. We use a splitting approach to decouple TV minimization from enforcing the constraints. Consequently, existing TV solvers can be employed with minimal alteration. This also makes the approach straightforward to generalize to any situation where TV can be applied. We consider deblurring of images with Gaussian or salt-and-pepper noise, as well as Abel inversion of radiographs with Poisson noise. We incorporate previous iterative reweighting algorithms to solve the TV portion.
  • Purpose: Dual-energy CT (DECT) is being increasingly used for its capability of material decomposition and energy-selective imaging. A generic problem of DECT, however, is that the decomposition process is unstable in the sense that the relative magnitude of decomposed signals is reduced due to signal cancellation while the image noise is accumulating from the two CT images of independent scans. Direct image decomposition, therefore, leads to severe degradation of signal-to-noise ratio on the resultant images. Existing noise suppression techniques are typically implemented in DECT with the procedures of reconstruction and decomposition performed independently, which do not explore the statistical propertiesmore » of decomposed images during the reconstruction for noise reduction. In this work, the authors propose an iterative approach that combines the reconstruction and the signal decomposition procedures to minimize the DECT image noise without noticeable loss of resolution. Methods: The proposed algorithm is formulated as an optimization problem, which balances the data fidelity and total variation of decomposed images in one framework, and the decomposition step is carried out iteratively together with reconstruction. The noise in the CT images from the proposed algorithm becomes well correlated even though the noise of the raw projections is independent on the two CT scans. Due to this feature, the proposed algorithm avoids noise accumulation during the decomposition process. The authors evaluate the method performance on noise suppression and spatial resolution using phantom studies and compare the algorithm with conventional denoising approaches as well as combined iterative reconstruction methods with different forms of regularization. Results: On the Catphan©600 phantom, the proposed method outperforms the existing denoising methods on preserving spatial resolution at the same level of noise suppression, i.e., a reduction of noise standard deviation by one order of magnitude. This improvement is mainly attributed to the high noise correlation in the CT images reconstructed by the proposed algorithm. Iterative reconstruction using different regularization, including quadratic orq-generalized Gaussian Markov random field regularization, achieves similar noise suppression from high noise correlation. However, the proposed TV regularization obtains a better edge preserving performance. Studies of electron density measurement also show that our method reduces the average estimation error from 9.5% to 7.1%. On the anthropomorphic head phantom, the proposed method suppresses the noise standard deviation of the decomposed images by a factor of ∼14 without blurring the fine structures in the sinus area. Conclusions: The authors propose a practical method for DECT imaging reconstruction, which combines the image reconstruction and material decomposition into one optimization framework. Compared to the existing approaches, our method achieves a superior performance on DECT imaging with respect to decomposition accuracy, noise reduction, and spatial resolution.« less