Diffuse Optical Tomography, Inverse Problems, and Optimization. Mary Katherine Huffman. Undergraduate Research Fall 2011 Spring 2012
|
|
- Jacob Parsons
- 6 years ago
- Views:
Transcription
1 Diffuse Optical Tomography, Inverse Problems, and Optimization Mary Katherine Huffman Undergraduate Research Fall 11 Spring Introduction. This paper discusses research conducted in order to investigate improvements to an optimization problem regarding diffuse optical tomography, an imaging technique in which near-infrared light illuminates tissue. The scattering and absorption of the light is recorded and studied, and this light allows for a picture of the inner characteristics of the tissue to be formed [2]. Specifically, the optimization problem arises after a parametric level sets-based regularization method adequately models an anomaly and its background. This research examines possible improvements regarding decreasing the amount of parameters, generating a good initial guess as to the location of the anomaly, and introducing a second optimization phase after a first has been completed. 2. Background. Before describing the analysis and results of the research, this paper will discuss three background concepts. Section 2.1 will address diffuse optical tomography and its beginnings. Section 2.2 will then explain the nature of inverse problems and illustrate the parametric level sets method as a representation of an anomaly and a medium using a modest number of parameters. Finally, 2.3 will discuss the resulting optimization problem. 2.1 Diffuse Optical Tomography. As previously stated, diffuse optical tomography (DOT) is a technique used to illuminate tissue and develop a model to reconstruct a picture of the inner characteristics. From [2], sources and detectors are placed around the tissue, and this model is created based on the measurements that are calculated. Currently, diffuse optical tomography is being applied to the detection of tumors in the brain and the breast. Before DOT was used, photon migration imaging was explored, and to calculate photon measurements, either illumination by pulses of light, continuous wave illumination, or radio-frequency amplitude modulated illumination was employed. Other pre-dot applications include pulse oximetry and near-infrared spectroscopy. Different choices of layouts for the sources and detectors have been proposed. Light sources, times, and frequencies can deviate; also, the use of virtual sources can be considered. Sources and detectors may be located on opposite sides of the tissue, gathered on the same side of the tissue, or placed in rings around the tissue. The illumination is typically either continuous wave or radio-frequency amplitude modulated illumination. The concept of virtual sources addresses the idea of an optical source inside the tissue; however, since the implementation of sources inside human beings is not always practical, frequencies are used instead. Ultimately, to retrieve the wanted information from DOT, an inverse problem must be solved to recover the desired parameters. Recent research explores parametric level set methods for the solving of inverse problems. 2.2 Parametric Level Sets. Inverse problems are abundant in many fields of science and engineering, including DOT [1]. Generally, a measure is taken from some medium, and the inverse problem serves as a means of finding a physical model to explain the observations.
2 Many times in real-world applications, the inverse problem is ill-conditioned, meaning that it is sensitive to noise; thus, retrieving true, accurate information is often extremely difficult. To deal with the ill-conditioned inverse problem, regularization is used, which stabilizes the reconstruction. The goal of solving the inverse problem is to identify characteristics about a particular region in a medium. As a way of approximating the region of interest along with the adjacent regions in the medium, the parametric level sets (PaLS) method has been introduced. For inverse problems regarding the retrieving of the geometry of a region, the foreground and background are viewed as separate cases, and the boundary between the two is represented as the level set of a function. For images involving DOT, we are interested in the background tissue and the anomaly located inside of the region. The PaLS method assumes that the characteristics for the scattering and absorption of light are significantly different in the anomaly and the background. The level set function is defined as where c ϵ R and is a parameter vector. If the function value is greater than c, we consider that region inside the anomaly. If the function value is equal to c, the function is on the boundary between the anomaly and the background, and if the function value is less than c, the function is outside of the anomaly. We assume the unknown property p is defined as Above, denotes the property values inside of the domain, and denotes the values outside of the domain. H is the Heaviside function, defined as H is approximated by a smooth Heaviside function,. The PaLS method represents the level set function as a summation of radial basis functions, defined as.. where α is the parameter vector of expansion coefficients, β is the vector of dilation coefficients, χ denotes the vector of the centers of the radial basis functions, and ψ denotes a smooth radial basis function. The norm utilized is an approximated Euclidean norm, which we define as where v is a small number. This parameterization regularizes the problem, and now the reconstruction of the image is an optimization problem.
3 2.3 Optimization. Specifically, the solution to an inverse problem is based on a measurement u of the region of interest and is defined by a model M, a function of the property p of the region, and noise η [3]. We represent u as where u is the actual measurement, and p, we solve for the residual operator, defined as is the assumed model. In order to approximate For ease, we will refer to any function as. The function then to minimize is To solve this minimization problem, iterative methods such as the Gauss-Newton (GN) method, Levenberg-Marquardt (LM) method, or TREGS method can be used [3]. The GN method approximates the function as.. where is the current approximate solution and is the ill-conditioned Jacobian of the residual operator. The LM method replaces in the GN method with where λ > 0 and controls the step size. TREGS uses the ideas of singular value decomposition and a trust region approach to approximate more accurate updates. For further detail regarding these methods, refer to [3]. 3. Analysis Overview. Current MATLAB code utilizes the PaLS and TREGS methods to reconstruct an image for a DOT model. The region of interest is constructed in a grid with both the x and y axes ranging from 0 to. The anomaly is within this region. The code assumes distinct property values for both the anomaly and the background and begins the process using an initialized guess of the anomaly location based on the radial basis functions. Then, µ is updated until the norm of the function is below an inputted tolerance. My research aimed to improve the optimization method in several ways. Section 3.1 will discuss holding certain parameter values constant and thus inhibiting the optimization to change them. Section 3.2 examines the results of using only a single radial basis function to define the parametric level sets function to locate the anomaly. Once the method produced a result using this single function, based on the reconstruction I strategically placed radial basis functions around the specified region and began a second optimization. I refer to this process as reseeding. Finally, section 3.3 will discuss three other reseeding options using more than one radial basis function for the first phase of optimization. 3.1 Constant Parameters. Recall that the parameters of the parametric level sets function include an α vector, β vector, and coordinates of the centers of the radial basis functions. By removing one of these vectors from the function, the optimization updates µ with respect to the remaining parameters. Thus, instead of optimizing with minor updates to both the dilation and expansion coefficients, the method makes larger updates to only one vector. Running the
4 optimization method once with a constant parameter vector did not always produce a better reconstruction, as seen in Figure 1 below. Figure 1. The actual image (A), the reconstruction after the original optimization (B), the reconstruction after the optimization with a constant β vector A. B. C. 3.2 A Single Radial Basis Locating Approach. In order to cheaply locate the anomaly in the region, I experimented with reconstructing the image in a first optimization phase with only one radial basis function. The theory was that the optimization would update this single function until the general area of the anomaly was located; then, in a second optimization phase, I would reseed the region with more strategically placed radial basis functions based on the first optimization of the single function. I began the optimization with a single function in the middle of the region, giving the radial basis function a center of (2., 2.). Figure 2 provides an example. The single radial basis function was able to locate the anomaly when only one anomaly existed. After this first reconstruction, the code updated five different radial basis functions. These functions were given centers corresponding to the final update of the first reconstruction. The x and y coordinates of the centers were a combination of the x and y coordinates of this final update and those values plus or minus 0.. Let (x*, y*) denote the center of the final update of the first function. Then, the second phase included five functions with centers (x*, y*+0.), (x*, y* 0.), (x*,y*), (x*+0., y*), and (x* 0., y*). Figure 2. The actual image (A), the single radial basis function (B), the first reconstruction with the single function (C), the second reconstruction with five functions (D) A. B. C. D. While the single function was able to locate the anomaly and approximate its size, the second optimization was able to capture the shape to a greater extent. However, if more than one anomaly was present in the region, the single radial basis function was not able to find the distinct locations. Instead, the optimization updated the function to ultimately locate only one anomaly or match the background. Figures 3 and 4 demonstrate these findings respectively.
5 Figure 3. The actual image (A), the reconstructed image (B). 4 x A. B. Figure 4. The actual image (A), the reconstructed image (B). x A. B. 3.3 Other Reseeding Methods. Additionally, I examined three options regarding reseeding the method with radial basis functions after a first optimization phase using more than one function. All examples begin with the same original image, and the first phase of the optimization utilizes fifteen radial basis functions distributed uniformly throughout the grid region. One option included using the centers of the final update of the radial basis functions after the first reconstruction as the centers of the functions for the second reconstruction. As shown in Figure, using radial basis functions already known to be close to the anomaly provided a more accurate reconstruction. The second reconstruction results in less background noise. Figure. The actual image (A), the reconstruction after a first optimization phase (B), the reconstruction after using the centers of the final update of the first phase functions as the centers for the radial basis functions of the second phase A. B. C.
6 Another option used the centers of the final update of the radial basis functions after the first reconstruction that also had positive expansion coefficients as the centers of the functions for the second reconstruction. These results also showed much improvement, as demonstrated in Figure 6. The higher accuracy could be due to the fact that the functions with positive expansion coefficients are the functions that are seen in the reconstructed image. Figure 6. The actual image (A), the reconstruction after a first optimization phase (B), the reconstruction after using the centers of the final update of the first phase functions with positive expansion coefficients as the centers for the radial basis functions of the second phase A. B. C. The third option used the centers of the final update of the radial basis functions after the first reconstruction with dilation coefficients greater than 0.9 as the centers of the functions for the second reconstruction. For the image used, this option appeared to be the most accurate. A possible explanation could be that since the functions with a dilation coefficient of a certain value or higher result from the optimization of the first reconstruction, these functions were clearly those most utilized in the first reimaging. Figure 7 provides the graphic results. Figure 7. The actual image (A), the reconstruction after a first optimization phase (B), the reconstruction after using the centers of the final update of the first phase functions with dilation coefficients greater than 0.9 as the centers for the radial basis functions of the second phase A. B. C. 4. Conclusions. After reviewing the results of the experiments, this research suggests that in general, a second optimization phase with reseeded radial basis functions produces a better reconstruction than a single phase of optimization. The three options for the reseeding procedure examined were (1) using the centers of the final updates of the functions from the first phase of optimization as the beginning centers of the functions for the second phase, (2) using the centers
7 of the final updates of the functions from the first phase of optimization with positive expansion coefficients as the beginning centers of the functions for the second phase, and (3) using the centers of the final updates of the functions from the first phase of optimization with dilation coefficients greater than 0.9 as the beginning centers of the functions for the second phase. In all three instances, a more accurate reconstruction does in fact appear to make sense. In (1), the first optimization phase produces an accurate reconstruction with radial basis functions spread uniformly throughout the region; therefore, the second phase should naturally produce accurate results when the radial basis functions are located more specifically around the anomaly. In (2), the reseeding process utilizes the updated functions of the first phase with positive expansion coefficients. This accurate image results from the fact that the functions with positive expansion coefficients are the functions that are above a specific level and thus contribute much to what is seen of the image. In (3), I postulate that using the functions with dilation coefficients greater than 0.9 created a more accurate image because this value would indicate that the corresponding function was utilized to a greater extent in the first optimization phase. The dilation coefficients are all initialized at 0.8. Therefore, with a final value of 0.9, these functions were beneficial to the reconstructed image and would be of greater use than the others. The method of using a single radial basis function is advantageous when only one anomaly exists in the region of interest. The single function accurately captures the location and relative size of the anomaly; however, when more than one anomaly exists, the single function is optimized to match the background or locates only one anomaly. I would not recommend this method of location in human applications unless it is known that only one anomaly is in existence. Regarding keeping a parameter vector constant, it appears that this method is not as accurate as utilizing all parameters. If only one parameter vector is available, this procedure can provide an approximate location for the anomalies. Unfortunately, the sizes and shapes of the anomalies are not as accurately displayed as with other methods. Finally, further research aims to investigate improving the optimization by reducing the Jacobian matrix to exclude certain columns. This method would eliminate the columns with norms below a tolerance determined by the matrix norm. The reduced matrix would in turn produce a more useful singular value decomposition for solving ill-conditioned matrix equations.
8 Works Cited [1] Aghasi, A., Kilmer, M., & Miller, E. L. (11). Parametric Level Set Methods for Inverse Problems. SIAM J. Imaging Sciences, [2] Boas, D. A., Brooks, D. H., Miller, E. L., DiMarzio, C. A., Kilmer, M., Gaudette, R. J., et al. (01). Imaging the Body with Diffuse Optical Tomography. IEEE Signal Processing Magazine, [3] de Sturler, E. & Kilmer, M. E. (11). A Regularized Gauss-Newton Trust Region Approach to Imaging in Diffuse Optical Tomography. SIAM J. Scientific Computing, 7-86.
REDUCED ORDER MODELING IN MULTISPECTRAL PHOTOACOUSTIC TOMOGRAPHY
REDUCED ORDER MODELING IN MULTISPECTRAL PHOTOACOUSTIC TOMOGRAPHY Arvind Saibaba Sarah Vallélian Statistical and Applied Mathematical Sciences Institute & North Carolina State University May 26, 2016 OUTLINE
More informationParametric Level Set Methods for Inverse Problems
arxiv:1007.2467v2 [math.ap] 6 Oct 2011 Parametric Level Set Methods for Inverse Problems Alireza Aghasi, Misha Kilmer, Eric L. Miller October 7, 2011 Abstract In this paper, a parametric level set method
More informationCharacterization of a Time-of-Flight Camera for Use in Diffuse Optical Tomography
Characterization of a Time-of-Flight Camera for Use in Diffuse Optical Tomography A thesis submitted in partial fulfillment of the requirements for the degree of Bachelor of Science degree in Physics from
More informationUltrasonic Multi-Skip Tomography for Pipe Inspection
18 th World Conference on Non destructive Testing, 16-2 April 212, Durban, South Africa Ultrasonic Multi-Skip Tomography for Pipe Inspection Arno VOLKER 1, Rik VOS 1 Alan HUNTER 1 1 TNO, Stieltjesweg 1,
More informationALIREZA AGHASI, MISHA KILMER, ERIC L. MILLER
PARAMETRIC LEVEL SET METHODS FOR INVERSE PROBLEMS ALIREZA AGHASI, MISHA KILMER, ERIC L. MILLER AAGHAS1@ECE.TUFTS.EDU, MISHA.KILMER@TUFTS.EDU, ELMILLER@ECE.TUTS.EDU Abstract. In this paper, a parametric
More informationAssessment of Diffuse Optical Tomography Image Reconstruction Methods Using a Photon Transport Model
Assessment of Diffuse Optical Tomography Image Reconstruction Methods Using a Photon Transport Model Murad M. Althobaiti 1 *, Hassan S. Salehi 2, and Quing Zhu 1,2 1 Department of Biomedical Engineering,
More informationModel parametrization strategies for Newton-based acoustic full waveform
Model parametrization strategies for Newton-based acoustic full waveform inversion Amsalu Y. Anagaw, University of Alberta, Edmonton, Canada, aanagaw@ualberta.ca Summary This paper studies the effects
More informationParameter Estimation in Differential Equations: A Numerical Study of Shooting Methods
Parameter Estimation in Differential Equations: A Numerical Study of Shooting Methods Franz Hamilton Faculty Advisor: Dr Timothy Sauer January 5, 2011 Abstract Differential equation modeling is central
More informationSUPPLEMENTARY INFORMATION
Supplementary Information Compact spectrometer based on a disordered photonic chip Brandon Redding, Seng Fatt Liew, Raktim Sarma, Hui Cao* Department of Applied Physics, Yale University, New Haven, CT
More informationBayesian estimation of optical properties of the human head via 3D structural MRI p.1
Bayesian estimation of optical properties of the human head via 3D structural MRI June 23, 2003 at ECBO 2003 Alex Barnett Courant Institute, New York University Collaborators (NMR Center, Mass. Gen. Hosp.,
More informationBiomagnetic inverse problems:
Biomagnetic inverse problems: Magnetic resonance electrical property tomography (MREPT) and magnetoencephalography (MEG) 2018 Aug. 16 The University of Tokyo Takaaki Nara 1 Contents Measurement What is
More informationOutline. Level Set Methods. For Inverse Obstacle Problems 4. Introduction. Introduction. Martin Burger
For Inverse Obstacle Problems Martin Burger Outline Introduction Optimal Geometries Inverse Obstacle Problems & Shape Optimization Sensitivity Analysis based on Gradient Flows Numerical Methods University
More informationSource Array. anomaly. Tissue. Receiver Array
A Shape-based Reconstruction Technique for DPDW Data Misha E. Kilmer Department of Mathematics, Tufts University, Medford MA 1 mkilme1@tufts.edu Eric L. Miller Department of Electrical and Computer Engineering,
More informationWEIGHTED L1 AND L2 NORMS FOR IMAGE RECONSTRUCTION: FIRST CLINICAL RESULTS OF EIT LUNG VENTILATION DATA
WEIGHTED L1 AND L2 NORMS FOR IMAGE RECONSTRUCTION: FIRST CLINICAL RESULTS OF EIT LUNG VENTILATION DATA PRESENTER: PEYMAN RAHMATI, PHD CANDIDATE, Biomedical Eng., Dept. of Systems and Computer Eng., Carleton
More informationDiffuse light tomography to detect blood vessels using Tikhonov regularization Huseyin Ozgur Kazanci* a, Steven L. Jacques b a
Diffuse light tomography to detect blood vessels using Tikhonov regularization Huseyin Ozgur Kazanci* a, Steven L. Jacques b a Biomedical Engineering, Faculty of Engineering, Akdeniz University, 07058
More informationImproving Forward Matrix Generation and Utilization for Time Domain Diffuse Optical Tomography. Damon Hyde
Improving Forward Matrix Generation and Utilization for Time Domain Diffuse Optical Tomography by Damon Hyde B.S., Worcester Polytechnic Institute, A thesis submitted to the Department of Electrical and
More informationParametric estimation of 3D tubular structures for diffuse optical tomography
Parametric estimation of 3D tubular structures for diffuse optical tomography Fridrik Larusson, 1,4, Pamela G. Anderson, 2 Elizabeth Rosenberg, 2 Misha E. Kilmer, 3 Angelo Sassaroli, 2 Sergio Fantini,
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More information(x, y, z) m 2. (x, y, z) ...] T. m 2. m = [m 1. m 3. Φ = r T V 1 r + λ 1. m T Wm. m T L T Lm + λ 2. m T Hm + λ 3. t(x, y, z) = m 1
Class 1: Joint Geophysical Inversions Wed, December 1, 29 Invert multiple types of data residuals simultaneously Apply soft mutual constraints: empirical, physical, statistical Deal with data in the same
More informationParameterization. Michael S. Floater. November 10, 2011
Parameterization Michael S. Floater November 10, 2011 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to generate from point
More informationInverse Scattering. Brad Nelson 3/9/2012 Math 126 Final Project
Inverse Scattering Brad Nelson /9/ Math 6 Final Project This paper investigates the ability to determine the boundary of an unknown object scattering waves using three different scenarios. In the first,
More informationAccelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method
Accelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method Wenyong Pan 1, Kris Innanen 1 and Wenyuan Liao 2 1. CREWES Project, Department of Geoscience,
More informationElectromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah
Electromagnetic migration of marine CSEM data in areas with rough bathymetry Michael S. Zhdanov and Martin Čuma*, University of Utah Summary In this paper we present a new approach to the interpretation
More informationImage Reconstruction in Optical Tomography : Utilizing Large Data Sets
Image Reconstruction in Optical Tomography : Utilizing Large Data Sets Vadim A. Markel, Ph.D. Department of Radiology John C. Schotland, M.D., Ph.D. Department of Bioengineering http://whale.seas.upenn.edu
More information3-D Tomographic Reconstruction
Mitglied der Helmholtz-Gemeinschaft 3-D Tomographic Reconstruction of Atmospheric Trace Gas Concentrations for Infrared Limb-Imagers February 21, 2011, Nanocenter, USC Jörn Ungermann Tomographic trace
More informationCritical computational aspects of near infrared circular tomographic imaging: Analysis of measurement number, mesh resolution and reconstruction basis
Critical computational aspects of near infrared circular tomographic imaging: Analysis of measurement number, mesh resolution and reconstruction basis Phaneendra K. Yalavarthy, Hamid Dehghani,, Brian W.
More informationAdaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited
Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited Summary We present a new method for performing full-waveform inversion that appears
More informationTomographic Algorithm for Industrial Plasmas
Tomographic Algorithm for Industrial Plasmas More info about this article: http://www.ndt.net/?id=22342 1 Sudhir K. Chaudhary, 1 Kavita Rathore, 2 Sudeep Bhattacharjee, 1 Prabhat Munshi 1 Nuclear Engineering
More informationCalibrating an Overhead Video Camera
Calibrating an Overhead Video Camera Raul Rojas Freie Universität Berlin, Takustraße 9, 495 Berlin, Germany http://www.fu-fighters.de Abstract. In this section we discuss how to calibrate an overhead video
More informationStatistical Analysis and Evaluation of Near Infrared Tomographic Imaging System
Statistical Analysis and Evaluation of Near Infrared Tomographic Imaging System A Thesis Submitted to the Faculty In partial fulfillment of the requirements for the Degree of Doctor of Philosophy by Xiaomei
More informationVisual Tracking (1) Pixel-intensity-based methods
Intelligent Control Systems Visual Tracking (1) Pixel-intensity-based methods Shingo Kagami Graduate School of Information Sciences, Tohoku University swk(at)ic.is.tohoku.ac.jp http://www.ic.is.tohoku.ac.jp/ja/swk/
More informationThe computer training program on tomography
The computer training program on tomography A.V. Bulba, L.A. Luizova, A.D. Khakhaev Petrozavodsk State University, 33 Lenin St., Petrozavodsk, 185640, Russia alim@karelia.ru Abstract: Computer training
More informationTwo Dimensional Microwave Imaging Using a Divide and Unite Algorithm
Two Dimensional Microwave Imaging Using a Divide and Unite Algorithm Disha Shur 1, K. Yaswanth 2, and Uday K. Khankhoje 2 1 Indian Institute of Engineering Science and Technology, Shibpur, India 2 Indian
More information26257 Nonlinear Inverse Modeling of Magnetic Anomalies due to Thin Sheets and Cylinders Using Occam s Method
26257 Nonlinear Inverse Modeling of Anomalies due to Thin Sheets and Cylinders Using Occam s Method R. Ghanati* (University of Tehran, Insitute of Geophysics), H.A. Ghari (University of Tehran, Insitute
More informationIterative resolution estimation in Kirchhoff imaging
Stanford Exploration Project, Report SERGEY, November 9, 2000, pages 369?? Iterative resolution estimation in Kirchhoff imaging Robert G. Clapp, Sergey Fomel, and Marie Prucha 1 ABSTRACT We apply iterative
More informationSIMULATION OF ARTIFICIAL SYSTEMS BEHAVIOR IN PARAMETRIC EIGHT-DIMENSIONAL SPACE
78 Proceedings of the 4 th International Conference on Informatics and Information Technology SIMULATION OF ARTIFICIAL SYSTEMS BEHAVIOR IN PARAMETRIC EIGHT-DIMENSIONAL SPACE D. Ulbikiene, J. Ulbikas, K.
More informationLinear DOT Reconstruction Comparison Abstract. We compare, through simulations, the performance of four linear algorithms for diffuse optical tomograp
A comparison study of linear reconstruction techniques for diffuse optical tomographic imaging of absorption coefficient Richard J. Gaudette, Dana H. Brooks, Charles A. DiMarzio, Misha E. Kilmer, Eric
More information13. Learning Ballistic Movementsof a Robot Arm 212
13. Learning Ballistic Movementsof a Robot Arm 212 13. LEARNING BALLISTIC MOVEMENTS OF A ROBOT ARM 13.1 Problem and Model Approach After a sufficiently long training phase, the network described in the
More informationJournal of Physics: Conference Series. Recent citations. To cite this article: Inga Berre et al 2008 J. Phys.: Conf. Ser.
Journal of Physics: Conference Series A multi-level strategy for regularized parameter identification using nonlinear reparameterization with sample application for reservoir characterization Recent citations
More informationAbstract. Three-dimensional diffuse optical tomography (DOT) attempts to map in vivo tissue blood oxygenation
Abstract Three-dimensional diffuse optical tomography (DOT) attempts to map in vivo tissue blood oxygenation and blood volume levels by reconstructing the spatial distribution of the optical absorption
More information2.710 Optics Spring 09 Solutions to Problem Set #1 Posted Wednesday, Feb. 18, 2009
MASSACHUSETTS INSTITUTE OF TECHNOLOGY.70 Optics Spring 09 Solutions to Problem Set # Posted Wednesday, Feb. 8, 009 Problem : Spherical waves and energy conservation In class we mentioned that the radiation
More informationPlane wave solutions for the radiative transport equation
Plane wave solutions for the radiative transport equation Arnold D. Kim adkim@ucmerced.edu School of Natural Sciences University of California, Merced Overview Introduction to plane wave solutions Green
More informationFourier Transformation Methods in the Field of Gamma Spectrometry
International Journal of Pure and Applied Physics ISSN 0973-1776 Volume 3 Number 1 (2007) pp. 132 141 Research India Publications http://www.ripublication.com/ijpap.htm Fourier Transformation Methods in
More informationAdvanced Image Reconstruction Methods for Photoacoustic Tomography
Advanced Image Reconstruction Methods for Photoacoustic Tomography Mark A. Anastasio, Kun Wang, and Robert Schoonover Department of Biomedical Engineering Washington University in St. Louis 1 Outline Photoacoustic/thermoacoustic
More informationHartley - Zisserman reading club. Part I: Hartley and Zisserman Appendix 6: Part II: Zhengyou Zhang: Presented by Daniel Fontijne
Hartley - Zisserman reading club Part I: Hartley and Zisserman Appendix 6: Iterative estimation methods Part II: Zhengyou Zhang: A Flexible New Technique for Camera Calibration Presented by Daniel Fontijne
More informationSparse wavelet expansions for seismic tomography: Methods and algorithms
Sparse wavelet expansions for seismic tomography: Methods and algorithms Ignace Loris Université Libre de Bruxelles International symposium on geophysical imaging with localized waves 24 28 July 2011 (Joint
More informationVisual Tracking (1) Tracking of Feature Points and Planar Rigid Objects
Intelligent Control Systems Visual Tracking (1) Tracking of Feature Points and Planar Rigid Objects Shingo Kagami Graduate School of Information Sciences, Tohoku University swk(at)ic.is.tohoku.ac.jp http://www.ic.is.tohoku.ac.jp/ja/swk/
More informationCamera calibration. Robotic vision. Ville Kyrki
Camera calibration Robotic vision 19.1.2017 Where are we? Images, imaging Image enhancement Feature extraction and matching Image-based tracking Camera models and calibration Pose estimation Motion analysis
More informationForward Problem Solution as the Operator of Filtered and Back Projection Matrix to Reconstruct the Various Method of Collecting Data and the Object Element in Electrical Impedance Tomography K. Ain 1,3,
More informationCS231A Course Notes 4: Stereo Systems and Structure from Motion
CS231A Course Notes 4: Stereo Systems and Structure from Motion Kenji Hata and Silvio Savarese 1 Introduction In the previous notes, we covered how adding additional viewpoints of a scene can greatly enhance
More informationFundamentals of Digital Image Processing
\L\.6 Gw.i Fundamentals of Digital Image Processing A Practical Approach with Examples in Matlab Chris Solomon School of Physical Sciences, University of Kent, Canterbury, UK Toby Breckon School of Engineering,
More informationCoE4TN4 Image Processing. Chapter 5 Image Restoration and Reconstruction
CoE4TN4 Image Processing Chapter 5 Image Restoration and Reconstruction Image Restoration Similar to image enhancement, the ultimate goal of restoration techniques is to improve an image Restoration: a
More informationMonte Carlo Method for Solving Inverse Problems of Radiation Transfer
INVERSE AND ILL-POSED PROBLEMS SERIES Monte Carlo Method for Solving Inverse Problems of Radiation Transfer V.S.Antyufeev. ///VSP/// UTRECHT BOSTON KÖLN TOKYO 2000 Contents Chapter 1. Monte Carlo modifications
More informationLecture 8 Object Descriptors
Lecture 8 Object Descriptors Azadeh Fakhrzadeh Centre for Image Analysis Swedish University of Agricultural Sciences Uppsala University 2 Reading instructions Chapter 11.1 11.4 in G-W Azadeh Fakhrzadeh
More informationEngineered Diffusers Intensity vs Irradiance
Engineered Diffusers Intensity vs Irradiance Engineered Diffusers are specified by their divergence angle and intensity profile. The divergence angle usually is given as the width of the intensity distribution
More informationTotal variation tomographic inversion via the Alternating Direction Method of Multipliers
Total variation tomographic inversion via the Alternating Direction Method of Multipliers Landon Safron and Mauricio D. Sacchi Department of Physics, University of Alberta Summary Geophysical inverse problems
More informationVisual Tracking (1) Feature Point Tracking and Block Matching
Intelligent Control Systems Visual Tracking (1) Feature Point Tracking and Block Matching Shingo Kagami Graduate School of Information Sciences, Tohoku University swk(at)ic.is.tohoku.ac.jp http://www.ic.is.tohoku.ac.jp/ja/swk/
More informationIMAGE DE-NOISING IN WAVELET DOMAIN
IMAGE DE-NOISING IN WAVELET DOMAIN Aaditya Verma a, Shrey Agarwal a a Department of Civil Engineering, Indian Institute of Technology, Kanpur, India - (aaditya, ashrey)@iitk.ac.in KEY WORDS: Wavelets,
More informationAmplitude and kinematic corrections of migrated images for non-unitary imaging operators
Stanford Exploration Project, Report 113, July 8, 2003, pages 349 363 Amplitude and kinematic corrections of migrated images for non-unitary imaging operators Antoine Guitton 1 ABSTRACT Obtaining true-amplitude
More informationSYSTEMS OF NONLINEAR EQUATIONS
SYSTEMS OF NONLINEAR EQUATIONS Widely used in the mathematical modeling of real world phenomena. We introduce some numerical methods for their solution. For better intuition, we examine systems of two
More informationElectrical Impedance Tomography: advances in Image Reconstruction
Electrical Impedance Tomography: advances in Image Reconstruction Andy Adler Systems and Computer Engineering, Carleton U, Ottawa, Canada Outline Electrical Impedance Tomography Imaging the lungs Measurement
More informationICRA 2016 Tutorial on SLAM. Graph-Based SLAM and Sparsity. Cyrill Stachniss
ICRA 2016 Tutorial on SLAM Graph-Based SLAM and Sparsity Cyrill Stachniss 1 Graph-Based SLAM?? 2 Graph-Based SLAM?? SLAM = simultaneous localization and mapping 3 Graph-Based SLAM?? SLAM = simultaneous
More informationMorphological Image Processing
Morphological Image Processing Ranga Rodrigo October 9, 29 Outline Contents Preliminaries 2 Dilation and Erosion 3 2. Dilation.............................................. 3 2.2 Erosion..............................................
More informationIndex. 3D reconstruction, point algorithm, point algorithm, point algorithm, point algorithm, 253
Index 3D reconstruction, 123 5+1-point algorithm, 274 5-point algorithm, 260 7-point algorithm, 255 8-point algorithm, 253 affine point, 43 affine transformation, 55 affine transformation group, 55 affine
More informationOverview. Spectral Processing of Point- Sampled Geometry. Introduction. Introduction. Fourier Transform. Fourier Transform
Overview Spectral Processing of Point- Sampled Geometry Introduction Fourier transform Spectral processing pipeline Spectral filtering Adaptive subsampling Summary Point-Based Computer Graphics Markus
More informationPlane Wave Imaging Using Phased Array Arno Volker 1
11th European Conference on Non-Destructive Testing (ECNDT 2014), October 6-10, 2014, Prague, Czech Republic More Info at Open Access Database www.ndt.net/?id=16409 Plane Wave Imaging Using Phased Array
More informationApplications of the generalized norm solver
Applications of the generalized norm solver Mandy Wong, Nader Moussa, and Mohammad Maysami ABSTRACT The application of a L1/L2 regression solver, termed the generalized norm solver, to two test cases,
More informationVision par ordinateur
Epipolar geometry π Vision par ordinateur Underlying structure in set of matches for rigid scenes l T 1 l 2 C1 m1 l1 e1 M L2 L1 e2 Géométrie épipolaire Fundamental matrix (x rank 2 matrix) m2 C2 l2 Frédéric
More informationGeometric theory of inversion and seismic imaging II: INVERSION + DATUMING + STATIC + ENHANCEMENT. August Lau and Chuan Yin.
Geometric theory of inversion and seismic imaging II: INVERSION + DATUMING + STATIC + ENHANCEMENT August Lau and Chuan Yin January 6, 2017 Abstract The goal of seismic processing is to convert input data
More informationCHAPTER 6 DETECTION OF MASS USING NOVEL SEGMENTATION, GLCM AND NEURAL NETWORKS
130 CHAPTER 6 DETECTION OF MASS USING NOVEL SEGMENTATION, GLCM AND NEURAL NETWORKS A mass is defined as a space-occupying lesion seen in more than one projection and it is described by its shapes and margin
More informationLecture 04. Fundamentals of Lidar Remote Sensing (2)
Lecture 04. Fundamentals of Lidar Remote Sensing (2) Lidar Equation Introduction Physical Picture of Lidar Equation Fundamental Lidar Equation Different Forms of Lidar Equation Illustration of Lidar Equation
More informationAbstract. Keywords. Computer Vision, Geometric and Morphologic Analysis, Stereo Vision, 3D and Range Data Analysis.
Morphological Corner Detection. Application to Camera Calibration L. Alvarez, C. Cuenca and L. Mazorra Departamento de Informática y Sistemas Universidad de Las Palmas de Gran Canaria. Campus de Tafira,
More informationWavefront-based models for inverse electrocardiography. Alireza Ghodrati (Draeger Medical) Dana Brooks, Gilead Tadmor (Northeastern University)
Wavefront-based models for inverse electrocardiography Alireza Ghodrati (Draeger Medical) Dana Brooks, Gilead Tadmor (Northeastern University) Rob MacLeod (University of Utah) Inverse ECG Basics Problem
More informationContents. I Basics 1. Copyright by SIAM. Unauthorized reproduction of this article is prohibited.
page v Preface xiii I Basics 1 1 Optimization Models 3 1.1 Introduction... 3 1.2 Optimization: An Informal Introduction... 4 1.3 Linear Equations... 7 1.4 Linear Optimization... 10 Exercises... 12 1.5
More information2 T. x + 2 T. , T( x, y = 0) = T 1
LAB 2: Conduction with Finite Difference Method Objective: The objective of this laboratory is to introduce the basic steps needed to numerically solve a steady state two-dimensional conduction problem
More informationFull waveform inversion by deconvolution gradient method
Full waveform inversion by deconvolution gradient method Fuchun Gao*, Paul Williamson, Henri Houllevigue, Total), 2012 Lei Fu Rice University November 14, 2012 Outline Introduction Method Implementation
More informationDescriptive Statistics, Standard Deviation and Standard Error
AP Biology Calculations: Descriptive Statistics, Standard Deviation and Standard Error SBI4UP The Scientific Method & Experimental Design Scientific method is used to explore observations and answer questions.
More informationXI Conference "Medical Informatics & Technologies" VALIDITY OF MRI BRAIN PERFUSION IMAGING METHOD
XI Conference "Medical Informatics & Technologies" - 2006 medical imaging, MRI, brain perfusion Bartosz KARCZEWSKI 1, Jacek RUMIŃSKI 1 VALIDITY OF MRI BRAIN PERFUSION IMAGING METHOD Brain perfusion imaging
More informationEE795: Computer Vision and Intelligent Systems
EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 FDH 204 Lecture 14 130307 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Review Stereo Dense Motion Estimation Translational
More informationLecture 24. Lidar Simulation
Lecture 24. Lidar Simulation q Introduction q Lidar Modeling via Lidar Simulation & Error Analysis q Functions of Lidar Simulation and Error Analysis q How to Build up Lidar Simulation? q Range-resolved
More informationAngle-gather time migration a
Angle-gather time migration a a Published in SEP report, 1, 141-15 (1999) Sergey Fomel and Marie Prucha 1 ABSTRACT Angle-gather migration creates seismic images for different reflection angles at the reflector.
More informationTraining of Neural Networks. Q.J. Zhang, Carleton University
Training of Neural Networks Notation: x: input of the original modeling problem or the neural network y: output of the original modeling problem or the neural network w: internal weights/parameters of
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Volume 19, 2013 http://acousticalsociety.org/ ICA 2013 Montreal Montreal, Canada 2-7 June 2013 Biomedical Acoustics Session 1aBA: Ultrasound Tomography 1aBA3. A contrast
More informationCCSSM Curriculum Analysis Project Tool 1 Interpreting Functions in Grades 9-12
Tool 1: Standards for Mathematical ent: Interpreting Functions CCSSM Curriculum Analysis Project Tool 1 Interpreting Functions in Grades 9-12 Name of Reviewer School/District Date Name of Curriculum Materials:
More informationSUMMARY. method to synthetic datasets is discussed in the present paper.
Geophysical modeling through simultaneous Joint Inversion of Seismic, Gravity and Magnetotelluric data Michele De Stefano (1), Daniele Colombo (1) WesternGeco EM - Geosystem, via Clericetti 42/A, 20133
More informationJournal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp
Journal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp. 59 74 SOLID MECHANICS WAVE PROPAGATION DUE TO AN EMBEDDED SEISMIC SOURCE IN A GRADED HALF-PLANE WITH RELIEF PECULIARITIES.
More informationMultiple View Geometry in Computer Vision
Multiple View Geometry in Computer Vision Prasanna Sahoo Department of Mathematics University of Louisville 1 Structure Computation Lecture 18 March 22, 2005 2 3D Reconstruction The goal of 3D reconstruction
More informationSimulation of Multipoint Ultrasonic Flowmeter
Simulation of Multipoint Ultrasonic Flowmeter Jakub Filipský 1,*, Jiří Nožička 2 1,2 CTU in Prague, Faculty of Mechanical Engineering, Department of Fluid Mechanics and Thermodynamics, Technická 4, 166
More informationcalibrated coordinates Linear transformation pixel coordinates
1 calibrated coordinates Linear transformation pixel coordinates 2 Calibration with a rig Uncalibrated epipolar geometry Ambiguities in image formation Stratified reconstruction Autocalibration with partial
More informationEE795: Computer Vision and Intelligent Systems
EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 FDH 204 Lecture 10 130221 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Review Canny Edge Detector Hough Transform Feature-Based
More informationA projected Hessian matrix for full waveform inversion Yong Ma and Dave Hale, Center for Wave Phenomena, Colorado School of Mines
A projected Hessian matrix for full waveform inversion Yong Ma and Dave Hale, Center for Wave Phenomena, Colorado School of Mines SUMMARY A Hessian matrix in full waveform inversion (FWI) is difficult
More informationIdentifying Car Model from Photographs
Identifying Car Model from Photographs Fine grained Classification using 3D Reconstruction and 3D Shape Registration Xinheng Li davidxli@stanford.edu Abstract Fine grained classification from photographs
More informationPlanes Intersecting Cones: Static Hypertext Version
Page 1 of 12 Planes Intersecting Cones: Static Hypertext Version On this page, we develop some of the details of the plane-slicing-cone picture discussed in the introduction. The relationship between the
More informationStabilization of the Euler deconvolution algorithm by means of a two steps regularization approach
Stabilization of the Euler deconvolution algorithm by means of a two steps regularization approach R. Pašteka ( 1,2 ), D. Kušnirák ( 1 ), H.-J. Götze ( 2 ) ( 1 )Department of Applied Geophysics, Comenius
More informationMULTI-RESOLUTION STATISTICAL ANALYSIS ON GRAPH STRUCTURED DATA IN NEUROIMAGING
MULTI-RESOLUTION STATISTICAL ANALYSIS ON GRAPH STRUCTURED DATA IN NEUROIMAGING, Vikas Singh, Moo Chung, Nagesh Adluru, Barbara B. Bendlin, Sterling C. Johnson University of Wisconsin Madison Apr. 19, 2015
More informationA Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods
A Random Variable Shape Parameter Strategy for Radial Basis Function Approximation Methods Scott A. Sarra, Derek Sturgill Marshall University, Department of Mathematics, One John Marshall Drive, Huntington
More informationRobust PDF Table Locator
Robust PDF Table Locator December 17, 2016 1 Introduction Data scientists rely on an abundance of tabular data stored in easy-to-machine-read formats like.csv files. Unfortunately, most government records
More informationTomographic Reconstruction
Tomographic Reconstruction 3D Image Processing Torsten Möller Reading Gonzales + Woods, Chapter 5.11 2 Overview Physics History Reconstruction basic idea Radon transform Fourier-Slice theorem (Parallel-beam)
More informationAPPLICATION OF ALGORITHMS FOR AUTOMATIC GENERATION OF HEXAHEDRAL FINITE ELEMENT MESHES
MESTRADO EM ENGENHARIA MECÂNICA November 2014 APPLICATION OF ALGORITHMS FOR AUTOMATIC GENERATION OF HEXAHEDRAL FINITE ELEMENT MESHES Luís Miguel Rodrigues Reis Abstract. The accuracy of a finite element
More informationStereo Vision. MAN-522 Computer Vision
Stereo Vision MAN-522 Computer Vision What is the goal of stereo vision? The recovery of the 3D structure of a scene using two or more images of the 3D scene, each acquired from a different viewpoint in
More information