497 research outputs found

    The D-bar Method for Diffuse Optical Tomography: a computational study

    Get PDF
    The D-bar method at negative energy is numerically implemented. Using the method we are able to numerically reconstruct potentials and investigate exceptional points at negative energy. Subsequently, applying the method to Diffusive Optical Tomography, a new way of reconstructing the diffusion coefficient from the associated Complex Geometrics Optics solution is suggested and numerically validated

    The D-Bar Method for Diffuse Optical Tomography : A Computational Study

    Get PDF
    The D-bar method at negative energy is numerically implemented. Using the method, we are able to numerically reconstruct potentials and investigate exceptional points at negative energy. Subsequently, applying the method to diffuse optical tomography, a new way of reconstructing the diffusion coefficient from the associated Complex Geometrics Optics solution is suggested and numerically validated.Peer reviewe

    Image reconstruction in quantitative photoacoustic tomography using adaptive optical Monte Carlo

    Get PDF
    In quantitative photoacoustic tomography (QPAT), distributions of optical parameters inside the target are reconstructed from photoacoustic images. In this work, we utilize the Monte Carlo (MC) method for light transport in the image reconstruction of QPAT. Modeling light transport accurately with the MC requires simulating a large number of photon packets, which can be computationally expensive. On the other hand, too low number of photon packets results in a high level of stochastic noise, which can lead to significant errors in reconstructed images. In this work, we use an adaptive approach, where the number of simulated photon packets is adjusted during an iterative image reconstruction. It is based on a norm test where the expected relative error of the minimization direction is controlled. The adaptive approach automatically determines the number of simulated photon packets to provide sufficiently accurate light transport modeling without unnecessary computational burden. The presented approach is studied with two-dimensional simulations

    Adaptive stochastic Gauss-Newton method with optical Monte Carlo for quantitative photoacoustic tomography

    Get PDF
    SIGNIFICANCE: The image reconstruction problem in quantitative photoacoustic tomography (QPAT) is an ill-posed inverse problem. Monte Carlo method for light transport can be utilized in solving this image reconstruction problem. AIM: The aim was to develop an adaptive image reconstruction method where the number of photon packets in Monte Carlo simulation is varied to achieve a sufficient accuracy with reduced computational burden. APPROACH: The image reconstruction problem was formulated as a minimization problem. An adaptive stochastic Gauss-Newton (A-SGN) method combined with Monte Carlo method for light transport was developed. In the algorithm, the number of photon packets used on Gauss-Newton (GN) iteration was varied utilizing a so-called norm test. RESULTS: The approach was evaluated with numerical simulations. With the proposed approach, the number of photon packets needed for solving the inverse problem was significantly smaller than in a conventional approach where the number of photon packets was fixed for each GN iteration. CONCLUSIONS: The A-SGN method with a norm test can be utilized in QPAT to provide accurate and computationally efficient solutions

    On Learned Operator Correction in Inverse Problems

    Get PDF
    We discuss the possibility of learning a data-driven explicit model correction for inverse problems and whether such a model correction can be used within a variational framework to obtain regularized reconstructions. This paper discusses the conceptual difficulty of learning such a forward model correction and proceeds to present a possible solution as a forward-adjoint correction that explicitly corrects in both data and solution spaces. We then derive conditions under which solutions to the variational problem with a learned correction converge to solutions obtained with the correct operator. The proposed approach is evaluated on an application to limited view photoacoustic tomography and compared to the established framework of the Bayesian approximation error method

    Why we can no longer ignore consecutive disasters

    Get PDF
    In recent decades, a striking number of countries have suffered from consecutive disasters: events whose impacts overlap both spatially and temporally, while recovery is still under way. The risk of consecutive disasters will increase due to growing exposure, the interconnectedness of human society, and the increased frequency and intensity of nontectonic hazard. This paper provides an overview of the different types of consecutive disasters, their causes, and impacts. The impacts can be distinctly different from disasters occurring in isolation (both spatially and temporally) from other disasters, noting that full isolation never occurs. We use existing empirical disaster databases to show the global probabilistic occurrence for selected hazard types. Current state‐of‐the art risk assessment models and their outputs do not allow for a thorough representation and analysis of consecutive disasters. This is mainly due to the many challenges that are introduced by addressing and combining hazards of different nature, and accounting for their interactions and dynamics. Disaster risk management needs to be more holistic and codesigned between researchers, policy makers, first responders, and companies
    corecore