Multilevel Optimization for Multi-Modal X-ray Data Analysis
|
|
- Sydney Hopkins
- 6 years ago
- Views:
Transcription
1 Multilevel Optimization for Multi-Modal X-ray Data Analysis Zichao (Wendy) Di Mathematics & Computer Science Division Argonne National Laboratory May 25, 2016
2 2 / 35 Outline Multi-Modality Imaging Example: X-ray Tomography Forward Model Optimization Algorithms Numerical Results Multilevel-Based Acceleration Introduction of MG/OPT Application of Multilevel in Tomographic Reconstruction Numerical Results Summary and Outlook
3 3 / 35 Outline Multi-Modality Imaging Example: X-ray Tomography Forward Model Optimization Algorithms Numerical Results Multilevel-Based Acceleration Introduction of MG/OPT Application of Multilevel in Tomographic Reconstruction Numerical Results Summary and Outlook
4 More is Better 1 CT:hard-tissue; T2(MRI):soft tissue; PET: functional characteristics. 1 Boss, A. and et al. (2010). Hybrid PET/MRI of Intracranial Masses: Initial Experiences and Comparison to PET/CT, JNM 4 / 35
5 5 / 35 Outline Multi-Modality Imaging Example: X-ray Tomography Forward Model Optimization Algorithms Numerical Results Multilevel-Based Acceleration Summary and Outlook
6 6 / 35 Tomographic Image Reconstruction Non-invasive imaging technique to visualize internal structures of object. Tomography applications: physics, chemistry, astronomy, geophysics, medicine, etc. Tomographic imaging modalities: X-ray transmission, ultrasound, magnetic resonance, X-ray fluorescence, etc. Task: estimate distribution of physical quantities in sample from measurements. Limited angle tomography reconstruction is naturally ill-conditioned.
7 Schematic Experimental Setup with Two Modalities 7 / 35
8 8 / 35 Reconstruction Approaches Traditional: filtered backprojection Restricted to simple tomographic model Requires a large number of projections Alternatively, iterative reconstruction from single data modality Requires much less data acquisition, results in higher accuracy Our Goal: Formulate a joint inversion integrating XRF and XRT data to improve the reconstruciton quality of elemental map.
9 9 / 35 Outline Multi-Modality Imaging Example: X-ray Tomography Forward Model Optimization Algorithms Numerical Results Multilevel-Based Acceleration Summary and Outlook
10 10 / 35 X-ray Transmission (XRT) Traditionally, the XRT projection of the object from beam line (θ, τ) is modeled as { F T θ,τ ( µ E ) = I 0 exp } L v µ E v. v to directly solve the linear attenuation coefficient µ E v for each voxel v, In our approach,notice µ E v = e W v,e µ E e, { F T θ,τ (W) = I 0 exp } L v µ E e W v,e v,e I 0 : incident photon flux µ E e : mass attenuation coefficient of element e E at beam incident energy E W = W v,e: tensor denoting how much of element e is in voxel v L = [L v ]: tensor of intersection length of beam line (θ, τ) with the voxel v
11 X-ray Fluorescence (XRF) 11 / 35
12 12 / 35 Mathematical Model of XRF First, we obtain the unit fluorescence spectrum ( { })) 1 x M e = F (F(I 1 2 e ) F σ 2π exp 2σ 2 Then, the XRF spectrum, F R θ,τ = L v W v,e M e exp n d ( W v,e µ E e L v I v U v + µee v,e,d v,e e P v,v,d ) P = [P v,v,d ]: tensor of intersection length of fluorescence detectorlet path d with the voxel v
13 13 / 35 Outline Multi-Modality Imaging Example: X-ray Tomography Forward Model Optimization Algorithms Numerical Results Multilevel-Based Acceleration Summary and Outlook
14 14 / 35 Resulting Optimization: Joint Reconstruction (JRT) Goal: Find W so that F R θ,τ (W) = DR θ,τ and FT θ,τ (W) = DT θ,τ min φ(w) = ( 1 F R θ,τ (W) D R θ,τ 2 + β F T θ,τ (W) D T θ,τ 2) W θ,τ where D R θ,τ Rn E : measurement data of XRF signal detected at angle θ from light beam τ R: the measurement data of XRT signal detected at angle θ from light beam τ D T θ,τ β > 0: a scaling parameter to balance the two modalities Note: Optimization differs on how F R θ,τ and FT θ,τ are combined
15 How Multimodality Can Help Consider an overdetermined (more equations than the number of unknowns) but rank deficient system which has an infinitude of solutions: 1 2 [ ] 2 4 x1 = x Another such rank-deficient (not full-rank) system: [ ] [ ] [ ] 2 3 x1 1.5 = However, combine these two can form a full rank and consistent system with a unique solution x [ ] x1 2 = 2 3 x x 2 15 / 35
16 16 / 35 Optimization Solver In experiments, we use truncated-newton (TN) method with preconditioned conjugate gradient (PCG) to provide a search direction: To satisfy the bound constraints, the projected PCG is applied to the reduced Newton system One TN iteration typically requires κ(o(10)) PCG iterations Main expense of each outer iteration is κ + 2 function-gradient evaluations.
17 17 / 35 Outline Multi-Modality Imaging Example: X-ray Tomography Forward Model Optimization Algorithms Numerical Results Multilevel-Based Acceleration Summary and Outlook
18 JRT versus Single XRF Reconstruction (Synthetic) 1 Element Unit: g/µm 2 1 Di, Leyffer, and Wild (2016). An Optimization-Based Approach for Tomographic Inversion from Multiple Data Modalities, SIAMIS 18 / 35
19 Convergence 19 / 35
20 20 / 35 Addressing Self-Absorption Effect Sample: Solid glassrod (Silicon) with 2 wires (Tungsten, Gold).
21 21 / 35 Outline Multi-Modality Imaging Example: X-ray Tomography Forward Model Optimization Algorithms Numerical Results Multilevel-Based Acceleration Introduction of MG/OPT Application of Multilevel in Tomographic Reconstruction Numerical Results Summary and Outlook
22 22 / 35 Motivation N = 9 2, flops = N = 5 2, flops = N = 3 2, flops = N = 2 2, flops =
23 23 / 35 Outline Multi-Modality Imaging Multilevel-Based Acceleration Introduction of MG/OPT Application of Multilevel in Tomographic Reconstruction Numerical Results Summary and Outlook
24 Introduction of MG/OPT Multigrid optimization algorithm (MG/OPT) is a general framework to accelerate a traditional optimization algorithm 1. Recursively use coarse problems to generate search directions for fine problems. MG/OPT can deal with more general problems in an optimization perspective, in particular, it is able to handle inequality constraints in a natural way. Multiple options to design the hierarchy of the problem: through image space. through data space. 1 Nash, S. G. (2000). A Multigrid Approach to Discretized Optimization Problems, OMS 24 / 35
25 MG/OPT Given: High-resolution model f h (z h ), easier-to-solve low-resolution model f H (z H ) z + OPT(f (z), z, k) A restriction operator Ih H and an interpolation operator IH h An initial estimate z 0 h of the solution z h on the fine level Integers k 1 and k 2 satisfying k 1 + k 2 > 0 25 / 35
26 26 / 35 Outline Multi-Modality Imaging Multilevel-Based Acceleration Introduction of MG/OPT Application of Multilevel in Tomographic Reconstruction Numerical Results Summary and Outlook
27 27 / 35 Interpolation/Restriction Operators Restriction on 2D parameter W h R 2 to produce W H using full weighted matrix: Restriction on gradient ÎH h = C IH h where C balances the order difference between φ h ( W h ) and φ H (Ih H W h ) Interpolation operator: I h H = 4(IH h )T.
28 28 / 35 Coarse Grid Surrogate Model Shift experimental data for coarse level: D D D θ,τ = D θ,τ (F,h θ,τ ( W h ) F,H θ,τ (IH h W h )) The surrogate model: φ H (W H ) = θ,τ ( 1 2 F,H θ,τ (W H ) D D D θ,τ 2 ) ( φ H (Ih H W h ) ÎH h φ h( W T h )) WH
29 29 / 35 Outline Multi-Modality Imaging Multilevel-Based Acceleration Introduction of MG/OPT Application of Multilevel in Tomographic Reconstruction Numerical Results Summary and Outlook
30 30 / 35 MG/OPT Search Directions versus Error Problem Size: [33 17] Fine Level Error Interpolated Coarse Level Error
31 31 / 35 2-level MG/OPT Reconstruction 10 2 Single Level Multigrid Reconstruction Error Reduction Equivalent number of f g evaluations
32 32 / 35 Outline Multi-Modality Imaging Example: X-ray Tomography Forward Model Optimization Algorithms Numerical Results Multilevel-Based Acceleration Introduction of MG/OPT Application of Multilevel in Tomographic Reconstruction Numerical Results Summary and Outlook
33 33 / 35 Summary Established a link between X-ray transmission and X-ray fluorescence datasets by reformulating their corresponding physical models. Developed a simultaneous optimization approach for the joint inversion, and achieved a dramatic improvement of reconstruction quality with no extra computational cost. Proposed a multigrid-based optimization framework to further reduce the computational cost of the reconstruction problem. Preliminary results show that coarsening in voxel space improves accuracy/speed.
34 34 / 35 Making More of More Extend our multimodal analysis tool to data from different instruments, involving varying spatial resolution and contrast mechanisms. Guided by the hierarchical nature of our multilevel algorithm, we will investigate new data acquisition strategies and allow for flexible and adaptive sampling approaches. Enable a true real-time feedback.
35 35 / 35 Acknowledgement This work was funded under the project The Tao of Fusion: Pathways for Big-data Analysis of Energy Materials at Work and Base optimization (DOE-ASCR DE-AC02-06CH11357).
Image Acquisition Systems
Image Acquisition Systems Goals and Terminology Conventional Radiography Axial Tomography Computer Axial Tomography (CAT) Magnetic Resonance Imaging (MRI) PET, SPECT Ultrasound Microscopy Imaging ITCS
More informationSIMULTANEOUS SENSING ERROR RECOVERY AND TOMOGRAPHIC INVERSION USING AN OPTIMIZATION-BASED APPROACH
SIMULTANEOUS SENSING ERROR RECOVERY AND TOMOGRAPHIC INVERSION USING AN OPTIMIZATION-BASED APPROACH ANTHONY P. AUSTIN, ZICHAO (WENDY) DI, SVEN LEYFFER, AND STEFAN M. WILD. Abstract. Tomography can be used
More informationUNIVERSITY OF SOUTHAMPTON
UNIVERSITY OF SOUTHAMPTON PHYS2007W1 SEMESTER 2 EXAMINATION 2014-2015 MEDICAL PHYSICS Duration: 120 MINS (2 hours) This paper contains 10 questions. Answer all questions in Section A and only two questions
More informationEmission Computed Tomography Notes
Noll (24) ECT Notes: Page 1 Emission Computed Tomography Notes Introduction Emission computed tomography (ECT) is the CT applied to nuclear medicine. There are two varieties of ECT: 1. SPECT single-photon
More informationImaging of flow in porous media - from optimal transport to prediction
Imaging of flow in porous media - from optimal transport to prediction Eldad Haber Dept of EOS and Mathematics, UBC October 15, 2013 With Rowan Lars Jenn Cocket Ruthotto Fohring Outline Prediction is very
More informationsmooth coefficients H. Köstler, U. Rüde
A robust multigrid solver for the optical flow problem with non- smooth coefficients H. Köstler, U. Rüde Overview Optical Flow Problem Data term and various regularizers A Robust Multigrid Solver Galerkin
More informationSplitting-Based Statistical X-Ray CT Image Reconstruction with Blind Gain Correction
Splitting-Based Statistical X-Ray CT Image Reconstruction with Blind Gain Correction Hung Nien and Jeffrey A. Fessler Department of Electrical Engineering and Computer Science University of Michigan, Ann
More informationREDUCED 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 informationLow-Dose Dual-Energy CT for PET Attenuation Correction with Statistical Sinogram Restoration
Low-Dose Dual-Energy CT for PET Attenuation Correction with Statistical Sinogram Restoration Joonki Noh, Jeffrey A. Fessler EECS Department, The University of Michigan Paul E. Kinahan Radiology Department,
More informationReconstruction in CT and relation to other imaging modalities
Reconstruction in CT and relation to other imaging modalities Jørgen Arendt Jensen November 1, 2017 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound
More informationComputer Graphics. - Volume Rendering - Philipp Slusallek
Computer Graphics - Volume Rendering - Philipp Slusallek Overview Motivation Volume Representation Indirect Volume Rendering Volume Classification Direct Volume Rendering Applications: Bioinformatics Image
More informationRadon Transform and Filtered Backprojection
Radon Transform and Filtered Backprojection Jørgen Arendt Jensen October 13, 2016 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound Imaging Department
More informationAlgebraic Iterative Methods for Computed Tomography
Algebraic Iterative Methods for Computed Tomography Per Christian Hansen DTU Compute Department of Applied Mathematics and Computer Science Technical University of Denmark Per Christian Hansen Algebraic
More informationChapter6 Image Reconstruction
Chapter6 Image Reconstruction Preview 61I 6.1 Introduction 6.2 Reconstruction by Fourier Inversion 6.3 Reconstruction by convolution and backprojection 6.4 Finite series-expansion 1 Preview Reconstruction
More informationMultigrid Pattern. I. Problem. II. Driving Forces. III. Solution
Multigrid Pattern I. Problem Problem domain is decomposed into a set of geometric grids, where each element participates in a local computation followed by data exchanges with adjacent neighbors. The grids
More informationMoscow-Bavarian Joint Advanced Student School 2006 / Medical Imaging Principles of Computerized Tomographic Imaging and Cone-Beam Reconstruction
Line Integrals Line integrals represent the integral of some parameter of the object along the line (e.g. attenuation of x-rays) Object: f(x,y) Line: x cosθ + y sinθ = t Line integral / Radon transform:
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 informationMulti-slice CT Image Reconstruction Jiang Hsieh, Ph.D.
Multi-slice CT Image Reconstruction Jiang Hsieh, Ph.D. Applied Science Laboratory, GE Healthcare Technologies 1 Image Generation Reconstruction of images from projections. textbook reconstruction advanced
More informationarxiv: v2 [cond-mat.mtrl-sci] 5 Jan 2010
Gamma scattering scanning of concrete block for detection of voids. Shivaramu 1, Arijit Bose 2 and M. Margret 1 1 Radiological Safety Division, Safety Group, IGCAR, Kalpakaam - 63 12 (India) 2 Chennai
More informationl ealgorithms for Image Registration
FAIR: exib Image Registration l F l ealgorithms for Jan Modersitzki Computing And Software, McMaster University 1280 Main Street West, Hamilton On, L8S 4K1, Canada modersit@cas.mcmaster.ca August 13, 2008
More informationContinuous and Discrete Image Reconstruction
25 th SSIP Summer School on Image Processing 17 July 2017, Novi Sad, Serbia Continuous and Discrete Image Reconstruction Péter Balázs Department of Image Processing and Computer Graphics University of
More informationMéthodes d imagerie pour les écoulements et le CND
Méthodes d imagerie pour les écoulements et le CND Journée scientifique FED3G CEA LIST/Lab Imagerie Tomographie et Traitement Samuel Legoupil 15 juin 2012 2D/3D imaging tomography Example Petrochemical
More informationGradient Free Design of Microfluidic Structures on a GPU Cluster
Gradient Free Design of Microfluidic Structures on a GPU Cluster Austen Duffy - Florida State University SIAM Conference on Computational Science and Engineering March 2, 2011 Acknowledgements This work
More informationLINE SEARCH MULTILEVEL OPTIMIZATION AS COMPUTATIONAL METHODS FOR DENSE OPTICAL FLOW
LINE SEARCH MULTILEVEL OPTIMIZATION AS COMPUTATIONAL METHODS FOR DENSE OPTICAL FLOW EL MOSTAFA KALMOUN, LUIS GARRIDO, AND VICENT CASELLES Abstract. We evaluate the performance of different optimization
More informationReconstruction in CT and relation to other imaging modalities
Reconstruction in CT and relation to other imaging modalities Jørgen Arendt Jensen November 16, 2015 Center for Fast Ultrasound Imaging, Build 349 Department of Electrical Engineering Center for Fast Ultrasound
More informationIEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 14, NO. 1, JANUARY
IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL 14, NO 1, JANUARY 2005 125 A General Framework for Nonlinear Multigrid Inversion Seungseok Oh, Student Member, IEEE, Adam B Milstein, Student Member, IEEE, Charles
More informationRipplet: a New Transform for Feature Extraction and Image Representation
Ripplet: a New Transform for Feature Extraction and Image Representation Dr. Dapeng Oliver Wu Joint work with Jun Xu Department of Electrical and Computer Engineering University of Florida Outline Motivation
More informationWE consider the gate-sizing problem, that is, the problem
2760 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: REGULAR PAPERS, VOL 55, NO 9, OCTOBER 2008 An Efficient Method for Large-Scale Gate Sizing Siddharth Joshi and Stephen Boyd, Fellow, IEEE Abstract We consider
More informationX-ray tomography. X-ray tomography. Applications in Science. X-Rays. Notes. Notes. Notes. Notes
X-ray tomography Important application of the Fast Fourier transform: X-ray tomography. Also referred to as CAT scan (Computerized Axial Tomography) X-ray tomography This has revolutionized medical diagnosis.
More informationMEDICAL IMAGE ANALYSIS
SECOND EDITION MEDICAL IMAGE ANALYSIS ATAM P. DHAWAN g, A B IEEE Engineering in Medicine and Biology Society, Sponsor IEEE Press Series in Biomedical Engineering Metin Akay, Series Editor +IEEE IEEE PRESS
More informationTomographic reconstruction of shock layer flows
Tomographic reconstruction of shock layer flows A thesis submitted for the degree of Doctor of Philosophy at The Australian National University May 2005 Rado Faletič Department of Physics Faculty of Science
More informationA new calibration-free beam hardening reduction method for industrial CT
A new calibration-free beam hardening reduction method for industrial CT ECC 2 for industrial CT Tobias Würfl 1, Nicole Maaß 2, Frank Dennerlein 2, Andreas K. Maier 1 1 Pattern Recognition Lab, FAU Erlangen-Nürnberg;
More informationτ-extrapolation on 3D semi-structured finite element meshes
τ-extrapolation on 3D semi-structured finite element meshes European Multi-Grid Conference EMG 2010 Björn Gmeiner Joint work with: Tobias Gradl, Ulrich Rüde September, 2010 Contents The HHG Framework τ-extrapolation
More informationA Projection Access Scheme for Iterative Reconstruction Based on the Golden Section
A Projection Access Scheme for Iterative Reconstruction Based on the Golden Section Thomas Köhler Philips Research Laboratories Roentgenstrasse - Hamburg Germany Abstract A new access scheme for projections
More informationDIGITAL TERRAIN MODELLING. Endre Katona University of Szeged Department of Informatics
DIGITAL TERRAIN MODELLING Endre Katona University of Szeged Department of Informatics katona@inf.u-szeged.hu The problem: data sources data structures algorithms DTM = Digital Terrain Model Terrain function:
More informationEstimating 3D Respiratory Motion from Orbiting Views
Estimating 3D Respiratory Motion from Orbiting Views Rongping Zeng, Jeffrey A. Fessler, James M. Balter The University of Michigan Oct. 2005 Funding provided by NIH Grant P01 CA59827 Motivation Free-breathing
More informationThe Spherical Harmonics Discrete Ordinate Method for Atmospheric Radiative Transfer
The Spherical Harmonics Discrete Ordinate Method for Atmospheric Radiative Transfer K. Franklin Evans Program in Atmospheric and Oceanic Sciences University of Colorado, Boulder Computational Methods in
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 informationAlgebraic Iterative Methods for Computed Tomography
Algebraic Iterative Methods for Computed Tomography Per Christian Hansen DTU Compute Department of Applied Mathematics and Computer Science Technical University of Denmark Per Christian Hansen Algebraic
More informationMedical Image Analysis
Computer assisted Image Analysis VT04 29 april 2004 Medical Image Analysis Lecture 10 (part 1) Xavier Tizon Medical Image Processing Medical imaging modalities XRay,, CT Ultrasound MRI PET, SPECT Generic
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 informationDEVELOPMENT OF CONE BEAM TOMOGRAPHIC RECONSTRUCTION SOFTWARE MODULE
Rajesh et al. : Proceedings of the National Seminar & Exhibition on Non-Destructive Evaluation DEVELOPMENT OF CONE BEAM TOMOGRAPHIC RECONSTRUCTION SOFTWARE MODULE Rajesh V Acharya, Umesh Kumar, Gursharan
More informationMedical Image Reconstruction Term II 2012 Topic 6: Tomography
Medical Image Reconstruction Term II 2012 Topic 6: Tomography Professor Yasser Mostafa Kadah Tomography The Greek word tomos means a section, a slice, or a cut. Tomography is the process of imaging a cross
More informationData mining with sparse grids
Data mining with sparse grids Jochen Garcke and Michael Griebel Institut für Angewandte Mathematik Universität Bonn Data mining with sparse grids p.1/40 Overview What is Data mining? Regularization networks
More informationME964 High Performance Computing for Engineering Applications
ME964 High Performance Computing for Engineering Applications Outlining Midterm Projects Topic 3: GPU-based FEA Topic 4: GPU Direct Solver for Sparse Linear Algebra March 01, 2011 Dan Negrut, 2011 ME964
More informationCompressed Sensing for Electron Tomography
University of Maryland, College Park Department of Mathematics February 10, 2015 1/33 Outline I Introduction 1 Introduction 2 3 4 2/33 1 Introduction 2 3 4 3/33 Tomography Introduction Tomography - Producing
More informationEvaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT
Evaluation of Spectrum Mismatching using Spectrum Binning Approach for Statistical Polychromatic Reconstruction in CT Qiao Yang 1,4, Meng Wu 2, Andreas Maier 1,3,4, Joachim Hornegger 1,3,4, Rebecca Fahrig
More informationA Non-Linear Image Registration Scheme for Real-Time Liver Ultrasound Tracking using Normalized Gradient Fields
A Non-Linear Image Registration Scheme for Real-Time Liver Ultrasound Tracking using Normalized Gradient Fields Lars König, Till Kipshagen and Jan Rühaak Fraunhofer MEVIS Project Group Image Registration,
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 informationMultigrid Methods for Markov Chains
Multigrid Methods for Markov Chains Hans De Sterck Department of Applied Mathematics, University of Waterloo collaborators Killian Miller Department of Applied Mathematics, University of Waterloo, Canada
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 informationRADIOMICS: potential role in the clinics and challenges
27 giugno 2018 Dipartimento di Fisica Università degli Studi di Milano RADIOMICS: potential role in the clinics and challenges Dr. Francesca Botta Medical Physicist Istituto Europeo di Oncologia (Milano)
More informationEfficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers
Efficient Beam Velocity Model Building with Tomography Designed to Accept 3d Residuals Aligning Depth Offset Gathers J.W.C. Sherwood* (PGS), K. Sherwood (PGS), H. Tieman (PGS), R. Mager (PGS) & C. Zhou
More informationGPU implementation for rapid iterative image reconstruction algorithm
GPU implementation for rapid iterative image reconstruction algorithm and its applications in nuclear medicine Jakub Pietrzak Krzysztof Kacperski Department of Medical Physics, Maria Skłodowska-Curie Memorial
More informationRecent Developments in Model-based Derivative-free Optimization
Recent Developments in Model-based Derivative-free Optimization Seppo Pulkkinen April 23, 2010 Introduction Problem definition The problem we are considering is a nonlinear optimization problem with constraints:
More informationNONLINEAR MULTIGRID INVERSION ALGORITHMS WITH APPLICATIONS TO STATISTICAL IMAGE RECONSTRUCTION. A Thesis. Submitted to the Faculty.
NONLINEAR MULTIGRID INVERSION ALGORITHMS WITH APPLICATIONS TO STATISTICAL IMAGE RECONSTRUCTION A Thesis Submitted to the Faculty of Purdue University by Seungseok Oh In Partial Fulfillment of the Requirements
More informationBME I5000: Biomedical Imaging
1 Lucas Parra, CCNY BME I5000: Biomedical Imaging Lecture 4 Computed Tomography Lucas C. Parra, parra@ccny.cuny.edu some slides inspired by lecture notes of Andreas H. Hilscher at Columbia University.
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationEfficient Multigrid based solvers for Isogeometric Analysis
Efficient Multigrid based solvers for Isogeometric Analysis R. Tielen, M. Möller and C. Vuik Delft Institute of Applied Mathematics (DIAM) Numerical Analysis / 22 Isogeometric Analysis (IgA) Extension
More informationMemory Efficient Adaptive Mesh Generation and Implementation of Multigrid Algorithms Using Sierpinski Curves
Memory Efficient Adaptive Mesh Generation and Implementation of Multigrid Algorithms Using Sierpinski Curves Michael Bader TU München Stefanie Schraufstetter TU München Jörn Behrens AWI Bremerhaven Abstract
More informationFast Reconstruction of CFRP X-ray Images based on a Neural Network Filtered Backprojection Approach
Fast Reconstruction of CFRP X-ray Images based on a Neural Network Filtered Backprojection Approach More info about this article: http://www.ndt.net/?id=20852 Eline Janssens 1, Sascha Senck 2, Christoph
More informationCh. 4 Physical Principles of CT
Ch. 4 Physical Principles of CT CLRS 408: Intro to CT Department of Radiation Sciences Review: Why CT? Solution for radiography/tomography limitations Superimposition of structures Distinguishing between
More informationImplementation of a backprojection algorithm on CELL
Implementation of a backprojection algorithm on CELL Mario Koerner March 17, 2006 1 Introduction X-ray imaging is one of the most important imaging technologies in medical applications. It allows to look
More informationAn Introduc+on to Mathema+cal Image Processing IAS, Park City Mathema2cs Ins2tute, Utah Undergraduate Summer School 2010
An Introduc+on to Mathema+cal Image Processing IAS, Park City Mathema2cs Ins2tute, Utah Undergraduate Summer School 2010 Luminita Vese Todd WiCman Department of Mathema2cs, UCLA lvese@math.ucla.edu wicman@math.ucla.edu
More informationIterative methods for use with the Fast Multipole Method
Iterative methods for use with the Fast Multipole Method Ramani Duraiswami Perceptual Interfaces and Reality Lab. Computer Science & UMIACS University of Maryland, College Park, MD Joint work with Nail
More informationProgramming, numerics and optimization
Programming, numerics and optimization Lecture C-4: Constrained optimization Łukasz Jankowski ljank@ippt.pan.pl Institute of Fundamental Technological Research Room 4.32, Phone +22.8261281 ext. 428 June
More informationPhase problem and the Radon transform
Phase problem and the Radon transform Andrei V. Bronnikov Bronnikov Algorithms The Netherlands The Radon transform and applications Inverse problem of phase-contrast CT Fundamental theorem Image reconstruction
More informationDiscrete Estimation of Data Completeness for 3D Scan Trajectories with Detector Offset
Discrete Estimation of Data Completeness for 3D Scan Trajectories with Detector Offset Andreas Maier 1, Patrick Kugler 2, Günter Lauritsch 2, Joachim Hornegger 1 1 Pattern Recognition Lab and SAOT Erlangen,
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 informationMatrix-free IPM with GPU acceleration
Matrix-free IPM with GPU acceleration Julian Hall, Edmund Smith and Jacek Gondzio School of Mathematics University of Edinburgh jajhall@ed.ac.uk 29th June 2011 Linear programming theory Primal-dual pair
More informationRefractive Index Map Reconstruction in Optical Deflectometry Using Total-Variation Regularization
Refractive Index Map Reconstruction in Optical Deflectometry Using Total-Variation Regularization L. Jacques 1, A. González 2, E. Foumouo 2 and P. Antoine 2 1 Institute of Information and Communication
More informationMedical Images Analysis and Processing
Medical Images Analysis and Processing - 25642 Emad Course Introduction Course Information: Type: Graduated Credits: 3 Prerequisites: Digital Image Processing Course Introduction Reference(s): Insight
More informationAnisotropy-preserving 5D interpolation by hybrid Fourier transform
Anisotropy-preserving 5D interpolation by hybrid Fourier transform Juefu Wang and Shaowu Wang, CGG Summary We present an anisotropy-preserving interpolation method based on a hybrid 5D Fourier transform,
More informationDiscontinuous Galerkin Sparse Grid method for Maxwell s equations
Discontinuous Galerkin Sparse Grid method for Maxwell s equations Student: Tianyang Wang Mentor: Dr. Lin Mu, Dr. David L.Green, Dr. Ed D Azevedo, Dr. Kwai Wong Motivation u We consider the Maxwell s equations,
More informationGE s Revolution CT MATLAB III: CT. Kathleen Chen March 20, 2018
GE s Revolution CT MATLAB III: CT Kathleen Chen chens18@rpi.edu March 20, 2018 https://www.zmescience.com/medicine/inside-human-body-real-time-gifs-demo-power-ct-scan/ Reminders Make sure you have MATLAB
More informationMedical Imaging BMEN Spring 2016
Name Medical Imaging BMEN 420-501 Spring 2016 Homework #4 and Nuclear Medicine Notes All questions are from the introductory Powerpoint (based on Chapter 7) and text Medical Imaging Signals and Systems,
More informationAn Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems
An Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems Long Chen University of California, Irvine chenlong@math.uci.edu Joint work with: Huayi Wei (Xiangtan University),
More informationRefraction Corrected Transmission Ultrasound Computed Tomography for Application in Breast Imaging
Refraction Corrected Transmission Ultrasound Computed Tomography for Application in Breast Imaging Joint Research With Trond Varslot Marcel Jackowski Shengying Li and Klaus Mueller Ultrasound Detection
More informationGraphics 2009/2010, period 1. Lecture 8: ray tracing
Graphics 2009/2010, period 1 Lecture 8 Ray tracing Original plan Finish the basic stuff (book ch 1-9) today: We will not look into signal processing (chapter 4) We will look at triangle rasterization (cf.
More informationDipartimento di Ingegneria Informatica, Automatica e Gestionale A. Ruberti, SAPIENZA, Università di Roma, via Ariosto, Roma, Italy.
Data article Title: Data and performance profiles applying an adaptive truncation criterion, within linesearchbased truncated Newton methods, in large scale nonconvex optimization. Authors: Andrea Caliciotti
More informationEstimation of Tikhonov Regularization Parameter for Image Reconstruction in Electromagnetic Geotomography
Rafał Zdunek Andrze Prałat Estimation of ikhonov Regularization Parameter for Image Reconstruction in Electromagnetic Geotomograph Abstract: he aim of this research was to devise an efficient wa of finding
More informationTomographic Image Reconstruction in Noisy and Limited Data Settings.
Tomographic Image Reconstruction in Noisy and Limited Data Settings. Syed Tabish Abbas International Institute of Information Technology, Hyderabad syed.abbas@research.iiit.ac.in July 1, 2016 Tabish (IIIT-H)
More informationStudies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization
Studies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization Siva Nadarajah Antony Jameson Stanford University 15th AIAA Computational Fluid Dynamics Conference
More informationComputational Medical Imaging Analysis
Computational Medical Imaging Analysis Chapter 1: Introduction to Imaging Science Jun Zhang Laboratory for Computational Medical Imaging & Data Analysis Department of Computer Science University of Kentucky
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 informationModel Reduction for Variable-Fidelity Optimization Frameworks
Model Reduction for Variable-Fidelity Optimization Frameworks Karen Willcox Aerospace Computational Design Laboratory Department of Aeronautics and Astronautics Massachusetts Institute of Technology Workshop
More informationTheoretical Investigations of Tomographic Methods used for Determination of the Integrity of Spent BWR Nuclear Fuel
a UPPSALA UNIVERSITY Department of Radiation Sciences Box 535, S-751 1 Uppsala, Sweden http://www.tsl.uu.se/ Internal report ISV-6/97 August 1996 Theoretical Investigations of Tomographic Methods used
More informationData mining with sparse grids using simplicial basis functions
Data mining with sparse grids using simplicial basis functions Jochen Garcke and Michael Griebel Institut für Angewandte Mathematik Universität Bonn Part of the work was supported within the project 03GRM6BN
More informationBackground. Outline. Radiographic Tomosynthesis: Image Quality and Artifacts Reduction 1 / GE /
Radiographic Tomosynthesis: Image Quality and Artifacts Reduction Baojun Li, Ph.D Department of Radiology Boston University Medical Center 2012 AAPM Annual Meeting Background Linear Trajectory Tomosynthesis
More informationCompressed Sensing Algorithm for Real-Time Doppler Ultrasound Image Reconstruction
Mathematical Modelling and Applications 2017; 2(6): 75-80 http://www.sciencepublishinggroup.com/j/mma doi: 10.11648/j.mma.20170206.14 ISSN: 2575-1786 (Print); ISSN: 2575-1794 (Online) Compressed Sensing
More informationEnhanced material contrast by dual-energy microct imaging
Enhanced material contrast by dual-energy microct imaging Method note Page 1 of 12 2 Method note: Dual-energy microct analysis 1. Introduction 1.1. The basis for dual energy imaging Micro-computed tomography
More informationMultidimensional scaling
Multidimensional scaling Lecture 5 Alexander & Michael Bronstein tosca.cs.technion.ac.il/book Numerical geometry of non-rigid shapes Stanford University, Winter 2009 Cinderella 2.0 2 If it doesn t fit,
More informationRegularization by multigrid-type algorithms
Regularization by multigrid-type algorithms Marco Donatelli Department of Physics and Mathematics University of Insubria Joint work with S. Serra-Capizzano Outline 1 Restoration of blurred and noisy images
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 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 informationLecture 15: More Iterative Ideas
Lecture 15: More Iterative Ideas David Bindel 15 Mar 2010 Logistics HW 2 due! Some notes on HW 2. Where we are / where we re going More iterative ideas. Intro to HW 3. More HW 2 notes See solution code!
More informationGeometric Partitioning for Tomography
Geometric Partitioning for Tomography Jan-Willem Buurlage, CWI Amsterdam Rob Bisseling, Utrecht University Joost Batenburg, CWI Amsterdam 2018-09-03, Tokyo, Japan SIAM Conference on Parallel Processing
More informationNew Mexico Tech Hyd 510
Numerics Motivation Modeling process (JLW) To construct a model we assemble and synthesize data and other information to formulate a conceptual model of the situation. The model is conditioned on the science
More informationAdvanced Graphics. Path Tracing and Photon Mapping Part 2. Path Tracing and Photon Mapping
Advanced Graphics Path Tracing and Photon Mapping Part 2 Path Tracing and Photon Mapping Importance Sampling Combine importance sampling techniques Reflectance function (diffuse + specular) Light source
More informationMassively Parallel Phase Field Simulations using HPC Framework walberla
Massively Parallel Phase Field Simulations using HPC Framework walberla SIAM CSE 2015, March 15 th 2015 Martin Bauer, Florian Schornbaum, Christian Godenschwager, Johannes Hötzer, Harald Köstler and Ulrich
More information