Biomagnetic inverse problems:
|
|
- Roderick Bryan
- 5 years ago
- Views:
Transcription
1 Biomagnetic inverse problems: Magnetic resonance electrical property tomography (MREPT) and magnetoencephalography (MEG) 2018 Aug. 16 The University of Tokyo Takaaki Nara 1
2 Contents Measurement What is reconstructed physically mathematically MREPT Magnetic field inside the body Electrical conductivity and permittivity Coefficients of the timeharmonic Maxwell equation MEG Type of inverse problem Active Passive Magnetic field outside the head Current in the brain Source of the Poisson equation 2
3 I. Magnetic resonance electrical property tomography (MREPT) Joint work with M. Fushimi (The University of Tokyo) 3
4 Background Electrical conductivity and permittivity imaging inside the human body pathological information: breast carcinoma (tumor) exhibits conductivity changes up to a factor of 10. [Katscher, 2009] Small tumor/cancer tissues can be detected 4
5 Magnetic resonance electrical property tomography (MREPT) `Birdcage coil generates the rotating magnetic field with a frequency of e.g. 128MHz MRI scanner What is measured: magnetic field inside the body What should be reconstructed: electric conductivity permittivity 5
6 How to measure the magnetic field using MRI 1) Actively flip the magnetic moment 2) Fourier transform of the xy comp. of magnetic moment is measured. z y The flip angle is affected by the magnetic field. amplitude x phase 8
7 Governing equation: time harmonic Maxwell s equations Faraday s law Ampere s law Silver Muller radiation condition z conductivity : admittivity permittivity MREPT inverse problem: given inside the body, reconstruct. 9
8 Conventional method (1) `standard method in U.Katscher, IEEE Tr. Medical Imaging, By eliminating, a non linear PDE for is obtained: By neglecting, a reconstruction formula is obtained: 10
9 Significant errors occur where J.K.Seo, IEEE Tr. Medical Imaging, Conductivity, S/m true reconstructed 11
10 Conventional methods (2) for inohomogeneous Y.Song and J.K.Seo, SIAM J. Appl. Math, 2013: iterative method for semilinear elliptic PDE for 2D case H.Ammari et al., Inverse Problems, 2015: iterative method for semilinear elliptic PDE for 3D case F.S.Hafalir, et al., IEEE Trans. Med. Imaging, 2014: FEM for a linear PDE for 2 nd order spatial derivatives are required E.Balidemaj, et al., IEEE Trans. Med. Imaging, 2015: integral equation based approach Iterative algorithm E and H without the human body should be computed. 12
11 Objective Derivation of an explicit reconstruction formula for inhomogeneous admittivity from the measured magnetic field 1. Derivation of a direct reconstruction formula for the 2D problem 2. Stability and regularization 3. Extension to the 3D problem 13
12 1. Derivation of a direct reconstruction formula Assumptions: 1) 2) Usual assumptions from the structure of `birdcage coils does not change w.r.t the z axis z First, we consider the 2D problem and derive a direct reconstruction formula. Second, we remove Assumption 2) and propose an iterative algorithm. 14
13 Dbar equation forms of Maxwell s equations Define the compelx derivatives: Let Then, under the two assumptions, Maxwell s equations in are written: Dbar equation of Ampere s law Dbar equation of Faraday s law 15
14 Ampere s law Complex derivative of the measured magnetic field z component of the electric field Faraday s law We use Faraday s equation to represent in terms of 16
15 Generalized Cauchy formula for Dbar equations, : domain bounded by a Jordan contour satisfies a Dbar equation: (inhomogeneous Cauchy Riemann system) Generalized Cauchy formula Cf. (Cauchy Riemann equation) Cauchy formula 17
16 Generalized Cauchy formula for Dbar equations, : domain bounded by a Jordan contour satisfies a Dbar equation: (inhomogeneous Cauchy Riemann system) Generalized Cauchy formula Faraday s law The z comp. of the electric field can be expressed as 18
17 By substituting to Ampere s law, we have Dirichlet boundary value of 19
18 Direct reconstruction formula of the admittivity Given in and on, then, the adtmivity at an arbitrary point is reconstructed by 20
19 How to give the boundary value of Scenario 1: Take as a whole body and measure near the body surface Scenario 2:Take as a local domain inside which a small cancer may exist and give a normal value of on :body body 21
20 Numerical simulations Three cylinders (conductivity = 1S/m, diameters = 10mm, 20mm, 30mm) in a cylinder (conductivity = 0.5 S/m). Forward solution was computed by an FEM software (COMSOL). f=123mhz (corresponding to 2.89T) : 82mm x 82mm 82mm Pixel: 1.4mm x 1.4mm 22
21 Singular integral For computation of and, Gaussian filter + Savitzky Golay Filter were used. (2 nd degree polynomial approximation) 23
22 Reconstruction (noiseless) Proposed Conventional (standard) The three domains were well reconstructed. The large errors were observed at the boundary. 24
23 Reconstruction (1%noise) Proposed Conventional (standard) 25
24 Phantom experiments ROI Abnormal 1S/m Normal 0.5 S/m 3T MRI Phantom 26
25 Proposed Conventional (standard) discontinuous jump of the conductivity robustness against noise 27
26 2. Stability and regularization The relative error of : Measurement noise is enhanced by division by 28
27 Zero point of the electric field: Computed in the phantom exp. Reconstructed conductivity When using the birdcage coil, becomes zero at around the center of body. Reconstruction error becomes large due to division by the small 29
28 Regularization Setting sample positions, the point by point reconstruction formula is expressed as TV regularization can be easily incorporated: 30
29 `Physical regularization :control of the location where Putting an object with high permittivity The zero of Ez can be moved. 31
30 true explicit with regularization 32
31 3. Extension to the 3D problem Assumptions: z 1) 2) Usual assumptions from the structure of `birdcage coils does not change w.r.t the z axis Ampere s law (z compnents) + i Ampere s law (xy compnent) Faraday s law (xy compnents) 33
32 Extension without the assumption The Dbar equation is changed to: Iterative Algorithm: 1) Use the direct reconstruction formula to obtain an initial estimate. 2) Compute with the reconstructed. 3) Go to 2). 34
33 Numerical simulations 64MHz, plane wave input Z axis sampling: 2mm : central difference Gaussian filter: 3x3x3, sd=1pixel Convergence criteria: z=14 mm z=6 mm z 20 mm 35
34 Noise 1% 36
35 Summary of part I For the Magnetic Resonance Electrical Property Tomography (MREPT) inverse problem, the direct reconstruction formula was derived where the admittivity was explicitly reconstructed in terms of the magnetic field. Ill posedness can be physically controlled. TV regularization is easily incorporated. An iterative algorithm for the 3D problem was proposed based on the direct reconstruction formula. 37
36 II. Magnetoencephalography (MEG) inverse problem Joint work with K. Kabashima (The University of Tokyo) K. Watanabe (The University of Tokyo) K. Amano (CiNet, NICT) 38
37 II. Magnetoencephalography (MEG) inverse problem Neural currents in the brain forward problem inverse problem magnetic field Applications analysis of the brain function diagnosis of epilepsy T.Takeda, Univ. Tokyo 39
38 MEG inverse problem Spherically head model Given at outside the head, determine. 40
39 Conventional approaches to MEG inverse problems Parametric approach Equivalent current dipoles Imaging approach Current distribution on gird points The size and shape of a source domain cannot be determined. L2 reg.: over smoothed source L1 reg.: patch is not reconstructed 41
40 L1 norm regularization Reg. parameter was chosen by the generalized Cross Validation The reconstructed source is scattered around the true patch source 42
41 Objective Propose a method that determines patch sources, i.e., domains with explicit boundaries. Difficulty How do we express focal domains on the cortical surface? How do we move them in optimization constrained on the cortical surface? 43
42 Mapping from a sphere to the cortical surface mapping We construct a mapping from a unit sphere to the cortical surface Domains on the cortical surface are mapped from domains on the sphere 44 Optimization of the shapes of the domain on the sphere is much easier.
43 Conformal mapping from a cortical surface to a sphere X.GU, Y.Wang, T.F.Chan, P.M.Tompson, and S T.Yau, Genus zero surface conformal mapping and its application to brain surface mapping, IEEE Trans. Biomedical Eng., 23, :points on the cortical surface :points on the unit sphere Since the genus of is zero, : harmonic mapping conformal mapping : harmonic mapping 45
44 Discrete harmonic mapping 1) Initial points on the unit sphere are generated by the Gauss map: 2) Renew the points on the unit sphere by the Ricci flow: Mesh generated by MRI s.t 46
45 Parametric expression of domains mapping center radius We assume a domain on the cortical surface mapped from the circle on the sphere. The shape is limited, but the spatial extent of the domain is well represented with a few parameters. In order to express more complex domain, the radius of a domain on the sphere is expressed with the truncated Fourier series. 47
46 Parametric expression of domains mapping center radius Cost function: Unknowns: Optimization in the cartesian products Global minimum can be efficiently obtained by a Lipschitzian optimization algorithm. 48
47 Global optimization for parameters in a hyper rectangular space D. R. Jones, C. D. Perttunen, and B. E. Stuckman, Lipschitzian Optimization Without the Lipschitz Constant, J. Optim. Theory Appl., vol. 79, no. 1, pp , Global minimum in a Cartesian product can be efficiently searched. 49
48 Parametric expression of domains mapping center radius Cost function: Unknowns: the center and radius of the circle on the sphere: Optimization in the cartesian products Global minimum can be efficiently obtained by a Lipschitzian optimization algorithm. 50
49 Numerical simulations Sensors: 306ch (Elekta Neuromag) Noise: 10% Gaussian 51
50 True Proposed L1 52
51 radius 3deg MEG for visual stimulus Joint work with K. Amano (CiNet) and K. Watanabe (U. Tokyo) fmri result: Activated domains get larger as the size of stimulation increases. 6deg 12 deg
52 7.2 e e e 1 radius on the 54 sphere [rad]
53 Summary of part II Domains on the cortical surface are expressed as images mapped from the unit sphere. The center and radius of a circle on the unit sphere are obtained by a Lipschitzian optimization algorithm. A focal and connected source domain can be identified by using real visual stimulation / epileptic MEG data. 55
Advanced 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 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 informationDiffuse Optical Tomography, Inverse Problems, and Optimization. Mary Katherine Huffman. Undergraduate Research Fall 2011 Spring 2012
Diffuse Optical Tomography, Inverse Problems, and Optimization Mary Katherine Huffman Undergraduate Research Fall 11 Spring 12 1. Introduction. This paper discusses research conducted in order to investigate
More informationA Model-Independent, Multi-Image Approach to MR Inhomogeneity Correction
Tina Memo No. 2007-003 Published in Proc. MIUA 2007 A Model-Independent, Multi-Image Approach to MR Inhomogeneity Correction P. A. Bromiley and N.A. Thacker Last updated 13 / 4 / 2007 Imaging Science and
More informationAnalysis of Functional MRI Timeseries Data Using Signal Processing Techniques
Analysis of Functional MRI Timeseries Data Using Signal Processing Techniques Sea Chen Department of Biomedical Engineering Advisors: Dr. Charles A. Bouman and Dr. Mark J. Lowe S. Chen Final Exam October
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 informationTheory: modeling, localization and imaging
Electromagnetic Brain Mapping with MEG/EEG Theory: modeling, localization and imaging Sylvain Baillet Imaging Group Cognitive Neuroscience & Brain Imaging Lab. Hôpital de la Salpêtrière CNRS UPR640 - LENA
More informationBiomedical Image Processing
Biomedical Image Processing Jason Thong Gabriel Grant 1 2 Motivation from the Medical Perspective MRI, CT and other biomedical imaging devices were designed to assist doctors in their diagnosis and treatment
More informationBrain Surface Conformal Spherical Mapping
Brain Surface Conformal Spherical Mapping Min Zhang Department of Industrial Engineering, Arizona State University mzhang33@asu.edu Abstract It is well known and proved that any genus zero surface can
More informationShape-based Diffeomorphic Registration on Hippocampal Surfaces Using Beltrami Holomorphic Flow
Shape-based Diffeomorphic Registration on Hippocampal Surfaces Using Beltrami Holomorphic Flow Abstract. Finding meaningful 1-1 correspondences between hippocampal (HP) surfaces is an important but difficult
More informationImage Segmentation. Shengnan Wang
Image Segmentation Shengnan Wang shengnan@cs.wisc.edu Contents I. Introduction to Segmentation II. Mean Shift Theory 1. What is Mean Shift? 2. Density Estimation Methods 3. Deriving the Mean Shift 4. Mean
More informationSPECIAL TECHNIQUES-II
SPECIAL TECHNIQUES-II Lecture 19: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay Method of Images for a spherical conductor Example :A dipole near aconducting sphere The
More informationBioEST ver. 1.5beta User Manual Chang-Hwan Im (Ph.D.)
BioEST ver. 1.5beta User Manual Chang-Hwan Im (Ph.D.) http://www.bioest.com E-mail: ichism@elecmech.snu.ac.kr imxxx010@umn.edu 1. What is BioEST? BioEST is a special edition of SNUEEG(http://www.bioinverse.com),
More informationUnfitted Discontinuous Galerkin in Brain Research and the Topology Preserving Marching Cubes DUNE User Meeting 2015
MÜNSTER Unfitted Discontinuous Galerkin in Brain Research and the Topology Preserving Marching Cubes DUNE User Meeting 2015 Andreas Nüßing September 28th, 2015 MÜNSTER UDG in Brain Research and the TPMC
More informationSpherical Microphone Arrays
Spherical Microphone Arrays Acoustic Wave Equation Helmholtz Equation Assuming the solutions of wave equation are time harmonic waves of frequency ω satisfies the homogeneous Helmholtz equation: Boundary
More informationResearch Proposal: Computational Geometry with Applications on Medical Images
Research Proposal: Computational Geometry with Applications on Medical Images MEI-HENG YUEH yueh@nctu.edu.tw National Chiao Tung University 1 Introduction My research mainly focuses on the issues of computational
More informationDigital Geometry Processing Parameterization I
Problem Definition Given a surface (mesh) S in R 3 and a domain find a bective F: S Typical Domains Cutting to a Disk disk = genus zero + boundary sphere = closed genus zero Creates artificial boundary
More informationIndex. aliasing artifacts and noise in CT images, 200 measurement of projection data, nondiffracting
Index Algebraic equations solution by Kaczmarz method, 278 Algebraic reconstruction techniques, 283-84 sequential, 289, 293 simultaneous, 285-92 Algebraic techniques reconstruction algorithms, 275-96 Algorithms
More informationImplicit Surfaces or Level Sets of Functions *
Implicit Surfaces or Level Sets of Functions * Surfaces in R 3 are either described as parametrized images F : D 2 R 3 or as implicit surfaces, i.e., as levels of functions f : R 3 R, as the set of points
More informationFoundations of Analytical and Numerical Field Computation
Foundations of Analytical and Numerical Field Computation Stephan Russenschuck, CERN-AT-MEL Stephan Russenschuck CERN, TE-MCS, 1211 Geneva, Switzerland 1 Permanent Magnet Circuits 2 Rogowski profiles Pole
More informationAspects of RF Simulation and Analysis Software Methods. David Carpenter. Remcom. B = t. D t. Remcom (Europe)
Remcom (Europe) Central Boulevard Blythe Valley Park Solihull West Midlands England, B90 8AG www.remcom.com +44 870 351 7640 +44 870 351 7641 (fax) Aspects of RF Simulation and Analysis Software Methods
More informationCombination of Parallel Imaging and Compressed Sensing for high acceleration factor at 7T
Combination of Parallel Imaging and Compressed Sensing for high acceleration factor at 7T DEDALE Workshop Nice Loubna EL GUEDDARI (NeuroSPin) Joint work with: Carole LAZARUS, Alexandre VIGNAUD and Philippe
More informationHFSS Ansys ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary
HFSS 12.0 Ansys 2009 ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary Comparison of HFSS 11 and HFSS 12 for JSF Antenna Model UHF blade antenna on Joint Strike Fighter Inherent improvements in
More informationMODELING MIXED BOUNDARY PROBLEMS WITH THE COMPLEX VARIABLE BOUNDARY ELEMENT METHOD (CVBEM) USING MATLAB AND MATHEMATICA
A. N. Johnson et al., Int. J. Comp. Meth. and Exp. Meas., Vol. 3, No. 3 (2015) 269 278 MODELING MIXED BOUNDARY PROBLEMS WITH THE COMPLEX VARIABLE BOUNDARY ELEMENT METHOD (CVBEM) USING MATLAB AND MATHEMATICA
More informationBasics of treatment planning II
Basics of treatment planning II Sastry Vedam PhD DABR Introduction to Medical Physics III: Therapy Spring 2015 Dose calculation algorithms! Correction based! Model based 1 Dose calculation algorithms!
More informationModeling the Acoustic Scattering from Axially Symmetric Fluid, Elastic, and Poroelastic Objects due to Nonsymmetric Forcing Using COMSOL Multiphysics
Modeling the Acoustic Scattering from Axially Symmetric Fluid, Elastic, and Poroelastic Objects due to Nonsymmetric Forcing Using COMSOL Multiphysics Anthony L. Bonomo *1 and Marcia J. Isakson 1 1 Applied
More informationCoordinate Transformations in Advanced Calculus
Coordinate Transformations in Advanced Calculus by Sacha Nandlall T.A. for MATH 264, McGill University Email: sacha.nandlall@mail.mcgill.ca Website: http://www.resanova.com/teaching/calculus/ Fall 2006,
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 informationRecovery of Piecewise Smooth Images from Few Fourier Samples
Recovery of Piecewise Smooth Images from Few Fourier Samples Greg Ongie*, Mathews Jacob Computational Biomedical Imaging Group (CBIG) University of Iowa SampTA 2015 Washington, D.C. 1. Introduction 2.
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationConformal Mapping and Fluid Mechanics
Conformal Mapping and Fluid Mechanics Monday, December 02, 2013 1:56 PM Homework 4 due (hard deadline) Wednesday, December 11 at 5 PM. more problems available rapid grading window: Friday, December 5 -
More informationExercise 16: Magnetostatics
Exercise 16: Magnetostatics Magnetostatics is part of the huge field of electrodynamics, founding on the well-known Maxwell-equations. Time-dependent terms are completely neglected in the computation of
More informationImage Reconstruction from Multiple Projections ECE 6258 Class project
Image Reconstruction from Multiple Projections ECE 658 Class project Introduction: The ability to reconstruct an object from multiple angular projections is a powerful tool. What this procedure gives people
More informationFinite Element Convergence for Time-Dependent PDEs with a Point Source in COMSOL 4.2
Finite Element Convergence for Time-Dependent PDEs with a Point Source in COMSOL 4.2 David W. Trott and Matthias K. Gobbert Department of Mathematics and Statistics, University of Maryland, Baltimore County,
More informationDeconvolution with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC
Deconvolution with curvelet-domain sparsity Vishal Kumar, EOS-UBC and Felix J. Herrmann, EOS-UBC SUMMARY We use the recently introduced multiscale and multidirectional curvelet transform to exploit the
More information3D Surface Reconstruction of the Brain based on Level Set Method
3D Surface Reconstruction of the Brain based on Level Set Method Shijun Tang, Bill P. Buckles, and Kamesh Namuduri Department of Computer Science & Engineering Department of Electrical Engineering University
More informationLecture 3.4 Differential Equation Based Schemes
Lecture 3.4 Differential Equation Based Schemes 1 Differential Equation Based Schemes As stated in the previous lecture, another important and most widel used method of structured mesh generation is based
More informationDual modality electrical impedance and ultrasound reflection tomography to improve image quality
Dual modality electrical impedance and ultrasound reflection tomography to improve image quality K. Ain 1,4, D. Kurniadi 2, S. Suprijanto 2, O. Santoso 3 1. Physics Department - Airlangga University, Surabaya,
More informationParallel hp-finite Element Simulations of 3D Resistivity Logging Instruments
Parallel hp-finite Element Simulations of 3D Resistivity Logging Instruments M. Paszyński 1,3, D. Pardo 1,2, L. Demkowicz 1, C. Torres-Verdin 2 1 Institute for Computational Engineering and Sciences 2
More informationA New Approach to 3D Shape Recovery of Local Planar Surface Patches from Shift-Variant Blurred Images
A New Approach to 3D Shape Recovery of Local Planar Surface Patches from Shift-Variant Blurred Images Xue Tu, Murali Subbarao, Youn-Sik Kang {tuxue, murali, yskang}@ecesunysbedu, Dept of Electrical and
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 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 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 informationModeling Focused Beam Propagation in scattering media. Janaka Ranasinghesagara
Modeling Focused Beam Propagation in scattering media Janaka Ranasinghesagara Teaching Objectives Understand the need of focused beam computational models. Understand the concepts, equations and principles
More informationAn introduction to mesh generation Part IV : elliptic meshing
Elliptic An introduction to mesh generation Part IV : elliptic meshing Department of Civil Engineering, Université catholique de Louvain, Belgium Elliptic Curvilinear Meshes Basic concept A curvilinear
More information2½D Finite Element Method for Electrical Impedance Tomography Considering the Complete Electrode Model
2½D Finite Element Method for Electrical Impedance Tomography Considering the Complete Electrode Model by Navid Bahrani A Thesis submitted to the Faculty of Graduate and Postdoctoral Affairs in partial
More informationFast marching methods
1 Fast marching methods Lecture 3 Alexander & Michael Bronstein tosca.cs.technion.ac.il/book Numerical geometry of non-rigid shapes Stanford University, Winter 2009 Metric discretization 2 Approach I:
More informationGraphs of Equations. MATH 160, Precalculus. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Graphs of Equations
Graphs of Equations MATH 160, Precalculus J. Robert Buchanan Department of Mathematics Fall 2011 Objectives In this lesson we will learn to: sketch the graphs of equations, find the x- and y-intercepts
More information= f (a, b) + (hf x + kf y ) (a,b) +
Chapter 14 Multiple Integrals 1 Double Integrals, Iterated Integrals, Cross-sections 2 Double Integrals over more general regions, Definition, Evaluation of Double Integrals, Properties of Double Integrals
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 informationBiophysical Techniques (BPHS 4090/PHYS 5800)
Biophysical Techniques (BPHS 4090/PHYS 5800) Instructors: Prof. Christopher Bergevin (cberge@yorku.ca) Schedule: MWF 1:30-2:30 (CB 122) Website: http://www.yorku.ca/cberge/4090w2017.html York University
More informationAdvanced phase retrieval: maximum likelihood technique with sparse regularization of phase and amplitude
Advanced phase retrieval: maximum likelihood technique with sparse regularization of phase and amplitude A. Migukin *, V. atkovnik and J. Astola Department of Signal Processing, Tampere University of Technology,
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 informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
More informationLecture 2: Introduction
Lecture 2: Introduction v2015.0 Release ANSYS HFSS for Antenna Design 1 2015 ANSYS, Inc. Multiple Advanced Techniques Allow HFSS to Excel at a Wide Variety of Applications Platform Integration and RCS
More informationA comparison of minimum norm and MUSIC for a combined MEG/EEG sensor array
Adv. Radio Sci., 10, 99 104, 2012 doi:10.5194/ars-10-99-2012 Author(s) 2012. CC Attribution 3.0 License. Advances in Radio Science A comparison of minimum norm and MUSIC for a combined MEG/EEG sensor array
More informationNumerical methods in plasmonics. The method of finite elements
Numerical methods in plasmonics The method of finite elements Outline Why numerical methods The method of finite elements FDTD Method Examples How do we understand the world? We understand the world through
More informationFully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon
Fully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon Tamal Pramanick 1,a) 1 Department of Mathematics, Indian Institute of Technology Guwahati, Guwahati
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 informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More informationG009 Scale and Direction-guided Interpolation of Aliased Seismic Data in the Curvelet Domain
G009 Scale and Direction-guided Interpolation of Aliased Seismic Data in the Curvelet Domain M. Naghizadeh* (University of Alberta) & M. Sacchi (University of Alberta) SUMMARY We propose a robust interpolation
More informationA Diagonal Split-cell Model for the High-order Symplectic FDTD Scheme
PIERS ONLINE, VOL. 2, NO. 6, 2006 715 A Diagonal Split-cell Model for the High-order Symplectic FDTD Scheme Wei Sha, Xianliang Wu, and Mingsheng Chen Key Laboratory of Intelligent Computing & Signal Processing
More informationMATLAB. Advanced Mathematics and Mechanics Applications Using. Third Edition. David Halpern University of Alabama CHAPMAN & HALL/CRC
Advanced Mathematics and Mechanics Applications Using MATLAB Third Edition Howard B. Wilson University of Alabama Louis H. Turcotte Rose-Hulman Institute of Technology David Halpern University of Alabama
More informationAn Automated Image-based Method for Multi-Leaf Collimator Positioning Verification in Intensity Modulated Radiation Therapy
An Automated Image-based Method for Multi-Leaf Collimator Positioning Verification in Intensity Modulated Radiation Therapy Chenyang Xu 1, Siemens Corporate Research, Inc., Princeton, NJ, USA Xiaolei Huang,
More informationResearch on the New Image De-Noising Methodology Based on Neural Network and HMM-Hidden Markov Models
Research on the New Image De-Noising Methodology Based on Neural Network and HMM-Hidden Markov Models Wenzhun Huang 1, a and Xinxin Xie 1, b 1 School of Information Engineering, Xijing University, Xi an
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 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 informationSparse sampling in MRI: From basic theory to clinical application. R. Marc Lebel, PhD Department of Electrical Engineering Department of Radiology
Sparse sampling in MRI: From basic theory to clinical application R. Marc Lebel, PhD Department of Electrical Engineering Department of Radiology Objective Provide an intuitive overview of compressed sensing
More informationMath 231E, Lecture 34. Polar Coordinates and Polar Parametric Equations
Math 231E, Lecture 34. Polar Coordinates and Polar Parametric Equations 1 Definition of polar coordinates Let us first recall the definition of Cartesian coordinates: to each point in the plane we can
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 informationA1:Orthogonal Coordinate Systems
A1:Orthogonal Coordinate Systems A1.1 General Change of Variables Suppose that we express x and y as a function of two other variables u and by the equations We say that these equations are defining a
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,00,000 M Open access books available International authors and editors Downloads Our authors
More informationA Transmission Line Matrix Model for Shielding Effects in Stents
A Transmission Line Matrix Model for hielding Effects in tents Razvan Ciocan (1), Nathan Ida (2) (1) Clemson University Physics and Astronomy Department Clemson University, C 29634-0978 ciocan@clemon.edu
More informationFundamentals of 3D. Lecture 3: Debriefing: Lecture 2 Rigid transformations Quaternions Iterative Closest Point (+Kd-trees)
INF555 Fundamentals of 3D Lecture 3: Debriefing: Lecture 2 Rigid transformations Quaternions Iterative Closest Point (+Kd-trees) Frank Nielsen nielsen@lix.polytechnique.fr Harris-Stephens' combined corner/edge
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 informationRECONSTRUCTION AND ENHANCEMENT OF CURRENT DISTRIBUTION ON CURVED SURFACES FROM BIOMAGNETIC FIELDS USING POCS
DRAFT: October, 4: File: ramon-et-al pp.4 Page 4 Sheet of 8 CANADIAN APPLIED MATHEMATICS QUARTERLY Volume, Number, Summer RECONSTRUCTION AND ENHANCEMENT OF CURRENT DISTRIBUTION ON CURVED SURFACES FROM
More informationEnhao Gong, PhD Candidate, Electrical Engineering, Stanford University Dr. John Pauly, Professor in Electrical Engineering, Stanford University Dr.
Enhao Gong, PhD Candidate, Electrical Engineering, Stanford University Dr. John Pauly, Professor in Electrical Engineering, Stanford University Dr. Greg Zaharchuk, Associate Professor in Radiology, Stanford
More informationInformation Processing for Remote Sensing Ed. Chen CH, World Scientific, New Jersey (1999)
Information Processing for Remote Sensing Ed. Chen CH, World Scientific, New Jersey (1999) DISCRIMINATION OF BURIED PLASTIC AND METAL OBJECTS IN SUBSURFACE SOIL D. C. CHIN, DR. R. SRINIVASAN, AND ROBERT
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 informationMöbius Transformations in Scientific Computing. David Eppstein
Möbius Transformations in Scientific Computing David Eppstein Univ. of California, Irvine School of Information and Computer Science (including joint work with Marshall Bern from WADS 01 and SODA 03) Outline
More informationAn Investigation of a Model of Percentage Depth Dose for Irregularly Shaped Fields
Int. J. Cancer (Radiat. Oncol. Invest): 96, 140 145 (2001) 2001 Wiley-Liss, Inc. Publication of the International Union Against Cancer An Investigation of a Model of Percentage Depth Dose for Irregularly
More informationWEINER FILTER AND SUB-BLOCK DECOMPOSITION BASED IMAGE RESTORATION FOR MEDICAL APPLICATIONS
WEINER FILTER AND SUB-BLOCK DECOMPOSITION BASED IMAGE RESTORATION FOR MEDICAL APPLICATIONS ARIFA SULTANA 1 & KANDARPA KUMAR SARMA 2 1,2 Department of Electronics and Communication Engineering, Gauhati
More informationRole of Parallel Imaging in High Field Functional MRI
Role of Parallel Imaging in High Field Functional MRI Douglas C. Noll & Bradley P. Sutton Department of Biomedical Engineering, University of Michigan Supported by NIH Grant DA15410 & The Whitaker Foundation
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 informationLet s review the four equations we now call Maxwell s equations. (Gauss s law for magnetism) (Faraday s law)
Electromagnetic Waves Let s review the four equations we now call Maxwell s equations. E da= B d A= Q encl ε E B d l = ( ic + ε ) encl (Gauss s law) (Gauss s law for magnetism) dφ µ (Ampere s law) dt dφ
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 informationRectangular Coordinates in Space
Rectangular Coordinates in Space Philippe B. Laval KSU Today Philippe B. Laval (KSU) Rectangular Coordinates in Space Today 1 / 11 Introduction We quickly review one and two-dimensional spaces and then
More informationSolved with COMSOL Multiphysics 4.3a
Magnetic Lens Introduction Scanning electron microscopes image samples by scanning with a high-energy beam of electrons. The subsequent electron interactions produce signals such as secondary and back-scattered
More informationInterior Reconstruction Using the Truncated Hilbert Transform via Singular Value Decomposition
Interior Reconstruction Using the Truncated Hilbert Transform via Singular Value Decomposition Hengyong Yu 1, Yangbo Ye 2 and Ge Wang 1 1 CT Laboratory, Biomedical Imaging Division, VT-WFU School of Biomedical
More informationMRI Brain Image Segmentation Using an AM-FM Model
MRI Brain Image Segmentation Using an AM-FM Model Marios S. Pattichis', Helen Petropoulos2, and William M. Brooks2 1 Department of Electrical and Computer Engineering, The University of New Mexico, Albuquerque,
More informationIndex. C m (Ω), 141 L 2 (Ω) space, 143 p-th order, 17
Bibliography [1] J. Adams, P. Swarztrauber, and R. Sweet. Fishpack: Efficient Fortran subprograms for the solution of separable elliptic partial differential equations. http://www.netlib.org/fishpack/.
More informationElectromagnetic Field Modeling in Human Tissue
Electromagnetic Field Modeling in Human Tissue Iliana Marinova, Valentin Mateev Abstract For investigations of electromagnetic field distributions in biological structures by Finite Element Method (FEM),
More informationAn Iterative Approach for Reconstruction of Arbitrary Sparsely Sampled Magnetic Resonance Images
An Iterative Approach for Reconstruction of Arbitrary Sparsely Sampled Magnetic Resonance Images Hamed Pirsiavash¹, Mohammad Soleymani², Gholam-Ali Hossein-Zadeh³ ¹Department of electrical engineering,
More informationTumor Detection and classification of Medical MRI UsingAdvance ROIPropANN Algorithm
International Journal of Engineering Research and Advanced Technology (IJERAT) DOI:http://dx.doi.org/10.31695/IJERAT.2018.3273 E-ISSN : 2454-6135 Volume.4, Issue 6 June -2018 Tumor Detection and classification
More informationCreating an Automated Blood Vessel. Diameter Tracking Tool
Medical Biophysics 3970Z 6 Week Project: Creating an Automated Blood Vessel Diameter Tracking Tool Peter McLachlan - 250068036 April 2, 2013 Introduction In order to meet the demands of tissues the body
More informationRay Optics I. Last time, finished EM theory Looked at complex boundary problems TIR: Snell s law complex Metal mirrors: index complex
Phys 531 Lecture 8 20 September 2005 Ray Optics I Last time, finished EM theory Looked at complex boundary problems TIR: Snell s law complex Metal mirrors: index complex Today shift gears, start applying
More informationDevelopment of a Maxwell Equation Solver for Application to Two Fluid Plasma Models. C. Aberle, A. Hakim, and U. Shumlak
Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models C. Aberle, A. Hakim, and U. Shumlak Aerospace and Astronautics University of Washington, Seattle American Physical Society
More informationA fast breast nonlinear elastography reconstruction technique using the Veronda-Westman model
A fast breast nonlinear elastography reconstruction technique using the Veronda-Westman model Mohammadhosein Amooshahi a and Abbas Samani abc a Department of Electrical & Computer Engineering, University
More informationNew Approaches for EEG Source Localization and Dipole Moment Estimation. Shun Chi Wu, Yuchen Yao, A. Lee Swindlehurst University of California Irvine
New Approaches for EEG Source Localization and Dipole Moment Estimation Shun Chi Wu, Yuchen Yao, A. Lee Swindlehurst University of California Irvine Outline Motivation why EEG? Mathematical Model equivalent
More informationOff-the-Grid Compressive Imaging: Recovery of Piecewise Constant Images from Few Fourier Samples
Off-the-Grid Compressive Imaging: Recovery of Piecewise Constant Images from Few Fourier Samples Greg Ongie PhD Candidate Department of Applied Math and Computational Sciences University of Iowa April
More information