Statistical computing on manifolds and data assimilation: from medical images to anatomical and physiological models
Centrale OMA/MA33 - MVA 2016-2017
X. Pennec Introduction to medical image acquisition and processing
Asclepios team 2004, route des Lucioles B.P. 93 06902 Sophia Antipolis Cedex
http://www-sop.inria.fr/asclepios 1
•4 INRIA’s Research Centers
INRIA Lille Nord Europe
INRIA Paris Rocquencourt •INRIA Nancy •Grand Est INRIA •INRIA Rennes Saclay •Bretagne Atlantique Île-de-France
2,800 scientists in INRIA Grenoble computer science Rhône-Alpes INRIA Bordeaux 1300 researchers and teaching- Sud-Ouest researchers •INRIA Sophia Antipolis 1000 Ph.D. students •Méditerranée 500 postdoctoral fellows
Statistical computing on manifolds and data assimilation: from medical images to anatomical and physiological models Medical Imaging – Module 2 - MVA 2016-2017
Fri Jan 6: Introduction to Medical Image Analysis and Processing [XP]
Fri Jan 13: Statistics on Riemannian manifolds and Lie groups [XP]
Fri Jan 20: Biomechanical Modelling and Simulation [HD]
Fri Jan 27: Diffeomorphic transformations for image registration [XP]
Fri Feb 17: Statistics on deformations for computational anatomy [XP]
Fri Feb 24: Manifold-valued image processing: the tensor example [XP]
Fri Mar 3: Data assimilation methods for physiological modeling [HD]
Fri Mar 10: Surgery Simulation [HD]
Fri Mar 17: Exam [Hervé Delingette, Xavier Pennec]
3
1 Course overview
The data of Computational Anatomy & Physiology
Image acquisition
Medical Image Analysis processing
4
1895
Roentgen
First Nobel prize in Physics in 1901
5
Todays’s Medical Imaging modalities
CT Scan MRI PET Ultrasound
•Source :T. Peters
2 Dynamic Images (4-D)
CT Scan MRI
Bio-signals
ECG Pressure Sensor
Temperature
Medical Imaging in clinical practice
TEACHING
Image DataBase
PREVENTION DIAGNOSIS THERAPY Pre-operative Patient Patient Patient
Post-operative robot
•CT •US •X-Ray •MRI •X-Ray •US •US •Video •NM •iMRI
Per-operative
3 Course overview
The data of Computational Anatomy & Physiology
Image acquisition
Introduction
Tomography
Nuclear medicine
MRI
Medical Image Analysis processing
10
Today MRI X-Scan
Density and X-ray structure of absorption protons density
PET / SPECT Ultrasounds
Density of Variations of Radioactive acoustic isotopes impedance
11
Different imaging modalities
X-ray
Magnetic resonance imaging
anatomic, functional, angiographic, diffusion, spectroscopic, tagged
Transmission Tomography (X Scan)
Nuclear Medicine :
Positron emission tomography (PET)
Single photon emission tomography (SPECT)
Ultrasonography
Histological Imaging, confocal in-vivo microscopy, molecular imaging,…
12
4 Characteristics of medical images
Intensity values are related to physical tissue characteristics which in turn may relate to a physiological phenomenon
Anatomy
Physics
Physiology
Volumetric medical images
Very often medical images are volumetric Voxel Representation
I (x,y,z)
Example of volumetric images : CT-scan (Scanner)
Size: 512 x 512 x 128 Resolution: 0.5 x 0.5 x 1 mm
5 Volumetric images
T2 MRI
16
Course overview
The data of Computational Anatomy & Physiology
Image acquisition
Introduction
Tomography
Nuclear medicine
MRI
Medical Image Analysis processing
17
Principle of CT Imaging (1)
X-Ray X-Ray
18
6 Principle of CT Imaging (2)
Input X-Ray intensity : Ni
Measured Output X-Ray intensity : No
Ni N0
Exponential attenuation : Objective :
measure (x) = absorption coefficient of X-ray
19
Computed Tomography
Principle : • Reconstruct n dimensional function (image) from projected data of (n–1) dimension
Radon Transform (1917) • “Two dimension and three dimension object can be reconstructed from the infinite set of projection data”.
Fourier Slice Theorem 1D Fourier Transform of Projected slice of 2D field (x,y)
is equal to 1D slice in the same direction of the 2D Fourier Transform of (x,y)
21
7 Basic principle of CT -Reconstruction of 2 dimensional image-
Projection Data curvilinear integral of absorption coefficient regarding Y
y y
object x x
X-ray tube Data Acquisition field Reconstruction field Simple Backprojection
Reconstruction Principle Backprojection based on inverse Radon Transform