Image Registration: Warping Without Folding
|
|
- Silvia Parks
- 5 years ago
- Views:
Transcription
1 Image Registration: Warping Without Folding Jeffrey A. Fessler EECS Department BME Department, Dept. of Radiology The University of Michigan UM BME 500 Oct 7, 2009 Acknowledgments: Se Young Chun, Matt Jacobson James Balter, Marc Kessler 1
2 Introduction Outline Image registration Enforcing / encouraging local invertibility (diffeomorphism) IEEE J. Selected Topics in Signal Processing, Feb Motion-compensated image reconstruction Conventional Model based Temporal regularization Image reconstruction (and registration) toolbox: fessler 2
3 Image Reconstruction Overview object f System L noise measured data y Reconstruction ˆf X-ray CT Sinogram X-ray CT Image 3
4 Image Reconstruction Formulations Static f( r) Dynamic f( r,t) contrast changes object motion (synergy with image registration) Image Reconstruction Research Topics Models for imaging physics Models for measurement statistics Models for object Algorithm formulation Algorithm acceleration Algorithm analysis / evaluation... 4
5 Why FBP (thin slice) MBIR (120 kvp, ma with modulation, 0.4 s rotation, pitch 1.375, slice thickness mm) 5
6 Part 1 Nonrigid image registration ensuring local invertibility 6
7 Image Registration Many applications, e.g., forensics, remote sensing, medicine... rigid transformations nonrigid transformations (warps) Example: Respiratory motion (application: radiotherapy planning) Target Source Inhale Exhale 7
8 Image registration: Overview Given two images (or image volumes): f( r) and g( r), r = (x,y,z), find a spatial transformation T( r), where T : R 3 R 3, such that f( r) is similar to the warped image g( T( r)) Usual steps: parameterize by α the spatial transformation: T( r; α) ( ) choose a similarity measure Ψ f( ),g( T( )) find optimal deformation parameters α numerically: ( ) ˆα = argmax Ψ f( ),g( T( ; α)) α Challenge: want estimated transformation T( r; ˆα) to be plausible. Typically we want it to be diffeomorphic, or topology preserving, or invertible, or at least locally invertible. 8
9 Image registration: Similarity measures sum of squared differences Ψ( f 1, f 2 ) = f 1 (r i ) f 2 (r i ) 2 i correlation ( )( ) Ψ( f 1, f 2 ) = f1 (r i ) f 1 f2 (r i ) f 2 i mutual information... Ψ( f 1, f 2 ) = MI( f 1, f 2 ) = H( f 1 )+H( f 2 ) H( f 1, f 2 ) 9
10 Image registration: B-spline deformations Nonrigid spatial transformation: T( r; α) = r+(d x ( r; α x ),d y ( r; α y ),d z ( r; α z )), }{{} deformation where α = (α x, α y, α z ) denotes unknown deformation coefficients. Tensor-product B-spline deformation model: d x ( r; α x ) = α x i jk β(x/m x i)β(y/m y j)β(z/m z k) i, j,k α y i jk β(x/m x i)β(y/m y j)β(z/m z k) d y ( r; α y ) = i, j,k d z ( r; α z ) = i, j,k α z i jk β(x/m x i)β(y/m y j)β(z/m z k) m x,m y,m z denote the knot spacing in each dimension. These spacings determine the spatial scale of the deformation. 10
11 Cubic B-spline Kernel 2/3 1/
12 B-spline deformations: Benefits differentiable (smooth) local support recursive filters for computations piecewise polynomial hierarchical Nonrigid image registration similarity measures usually have many local maximizers. To help find a good local maximum, one usually uses coarse-to-fine search. This is easy with B-spline deformations. (Thevenaz & Unser, IEEE T-IP, 2000) 12
13 B-spline deformations illustrated knot locations, mx=8 my=16 Knot locations Local support 13
14 B-spline deformations illustrated Invertible Not invertible 14
15 Invertibility When is r T( r) = r + d( r) a locally invertible transformation? By the inverse function theorem, it is sufficient for T to be continuously differentiable, and { } have positive Jacobian determinant: det T( r) > 0 for all r. Jacobian of transformation/deformation ( ) T( r) = r+ d( r) = x dx y dx z dx x dy y dy z dy x dz y dz z dz (mathematical theory vs practice) 15
16 Ensuring local invertibility in image registration We need to estimate B-spline deformation coefficients α subject to some local invertiblity constraint α C: ˆα = argmax Ψ(α). α C Ideal sufficient condition ensuring local invertibility for parametric deformation model: { { } } α C 0 = α : det T( r; α) > 0, r R 3. This condition is very difficult to implement. Conventional relaxed local invertibility condition: C 0 C 1 { { } } α C 1 = α : det T( r; α) > 0, r grid points. This condition does not ensure local invertibility everywhere. It is also computationally demanding. We seek simpler sufficient conditions for local invertibility: C C 0. 16
17 Unconstrained vs constrained optimization Source No constraint Jacobian constraint Target No constraint Jacobian penalty Image registration is an ill-posed problem. Jacobian constraint on grid required > 3 computation as unconstrained case. Nevertheless, some negative Jacobians remain (between grid points) because C 0 C 1. We need a simpler constraint that ensures positive Jacobian determinants everywhere. 17
18 A sufficient condition: Box constraints Simple lower/upper bounds on B-spline coefficients: { α C 3 = α : α x i jk m x K, α y i jk m y K, α z i jk m } z K, i, j,k, where K 2.05 in 2D and K 2.48 in 3D. Choi et al., 2000; Rueckert et al., MICCAI 2006 Fact: C 3 C 0. So constraining α C 3 ensures local invertibility everywhere. Box constraints are particularly simple for optimization. However, C 3 is a very restrictive set of deformations. Maximum displacement is only about half the knot spacing. Precludes even simple (large) global translations. 18
19 Proposed sufficient condition for invertibility Theorem (IEEE JSTSP, Feb. 2009): Suppose 0 k q < 1 for q {x,y,z}. Define the set: 2 C 4 = {α : m x k x α x i+1, j,k α x i, j,k m x K x, m y k y α y i, j+1,k αy i, j,k m yk y, If α C 4, then r R 3 : { } 1 (k x + k y + k z ) det T( r; α) m z k z α z i, j,k+1 α z i, j,k m z K z, α q i+1, j,k αq i, j,k m qk q for q = y,z, α q i, j+1,k αq i, j,k m qk q for q = x,z, α q i, j,k+1 αq i, j,k m qk q for q = x,y, i, j,k}. (1+K x )(1+K y )(1+K z )+(1+K x )k y k z + k x (1+K y )k z + k x k y (1+K z ). Corollary: { } Choosing k x = k y = k z = 1/3 ε ensures that 0 < det T( r; α), r. 19
20 Comparing sufficient conditions 1D example with two coefficients: α 1, α 2, for n = 2 (quadratic B-splines) α 2 2 C 0 = C 4 C α 1 20
21 Limitations of sufficient conditions 2D simulations using augmented Lagrange multiplier approach to enforce the constraint α C 3 or α C 4. Source Box constraint Proposed constraint Target Box constraint Proposed constraint Clearly C 3 C 0 and C 4 C 0. 21
22 Solution: composition of transformations Composing multiple transformations can overcome the limitations of sufficient conditions, e.g., T T α3 T α2 T α1 where α 1, α 2, α 3 C 4. Proposed constraint Proposed constraint Proposed constraint 22
23 Composition for box constraints Requires many more compositions: Rueckert et al., MICCAI 2006 Box constraint Box constraint Box constraint Box constraint Box constraint Box constraint Each of the 30 warps used many augmented Lagrangian iterations. Tradeoff: simplicity of constraint and its flexibility. 23
24 Simplyfing further via regularization Idea: replace constrained optimization ˆα = argmax α C 4 Ψ(α) with simpler unconstrained, but regularized, optimization: ˆα = argmax Ψ(α) γr(α) α where R(α) is zero if α C 4 but is large otherwise. This encourages local invertibility, but does not enforce it strictly. The regularization parameter γ controls the tradeoff between image similarity regularity of the deformation (local invertibility). 24
25 Proposed regularizer proposed quadratic α x i+1,j,k αx i,j,k α y i+1,j,k αy i,j,k α z i+1,j,k αz i,j,k α x i,j+1,k αx i,j,k α y i,j+1,k αy i,j,k α z i,j+1,k αz i,j,k α x i,j,k+1 αx i,j,k α y i,j,k+1 αy i,j,k α z i,j,k+1 αz i,j,k Interval constraints in C 4 replaced by piecewise quadratic penalty function of differences of neighboring B-spline coefficients. cf. conventional quadratic roughness regularization m x = m y = m z = 1,k x = k y = k z = 1/3, and K x = K y = K z = 4/3. 25
26 Regularization tradeoffs As regularization parameter γ # of negative Jacobian determinants so RMS difference between images log of (1 + # of negative Jacobian determinants) Proposed quadratic log of RMS difference between target and deformed images Proposed regularizer: good image similarity, few negative Jacobian determinants. 26
27 Regularization tradeoffs: Details Proportion of negative Jacobian determinants 15 x no constraint Jacobian proposed quadratic kim RMS error of deformed image 27
28 3D registration of CT inhale/exhale scans 3D CT scans of a cancer patient at exhale and inhale, for radiation treatment planning voxels. Source: Coronal Target: Coronal 28
29 CT inhale/exhale scans: Sagittal Source: Sagittal Target: Sagittal 29
30 3D registration results: Unconstrained No constraint: Sagittal No constraint: Sagittal 500 No constraint: Sagittal Relatively small difference image values, but many negative Jacobian determinants 30
31 3D registration results: Jacobian Jacobian penalty: Sagittal Jacobian penalty: Sagittal 500 Jacobian penalty: Sagittal Jacobian penalty based on C 1 (grid) Better behaved warp and reasonable difference image. But slow. 31
32 3D registration results: Proposed regularizer Proposed method: Sagittal Proposed method: Sagittal 500 Proposed method: Sagittal Proposed regularizer based on C 4. Tailored design of constraint/penalty: k x = k y = 1/4, k z = 1/2. Similar warp and difference image, but faster. 32
33 Quantification Method CPU time RMS difference # negative (seconds) (HU) Jacobians Unconstrained Jacobian penalty Proposed penalty Computation time per iteration (in seconds) at the finest level Regularization parameter adjusted empirically in both penalized cases to be the smallest value that yields no negative Jacobian determinants on the voxel grid. 3 multiresolution levels: knot spacings 8 pixels with downsampled images, 8 pixels with original images, 4 pixels with original images. 120 iterations of CG at each level. Work in progress to compose a couple coarse-scale deformations before refining to fine scale to reduce RMS differences. Code on web site: fessler 33
34 Diffeomorphisms: To be or not to be... Source Target Invertibility penalty Sliding at diaphragm / rib cage interface. Enforcing smoothness leads to bone warping. 34
35 Work in progress... Requiring the warp to be smooth everywhere seems suboptimal. One possible solution presented at SPIE Stay tuned... 35
36 Summary Simple condition for ensuring local invertibility everywhere Admits more deformations than conventional box constraints Simple regularizer requires comparable computation as unregularized image registration and much less computation than Jacobian determinant constraints / penalties Rigid structures (bones) Sliding tissue interfaces Parameter selection Computation (GPUs?) Performance characterization Open problems 36
37 Part 2 Model-based image reconstruction with motion compensation 37
38 Motion in image reconstruction Object being scanned: f( r,t) Measured data vector: y = (y 1,...,y M ) Static image reconstruction: Assume f( r,t) = f( r,t 0 ) = f( r) during scan. Estimate f( r) from measurements y. (Ill posed.) Dynamic image reconstruction ( List-mode data model) Assume each data point y i is acquired instantaneously at a corresponding time instant t i Relate y i to object at time t i, e.g., y i = Z a i ( r) }{{} physics f( r,t i )d r + }{{} ε i. statistics More generally: p(y i f(,t i )). Estimate f( r,t) from measurements y. (Even more ill posed!) 38
39 Gated data model Group data into K vectors, e.g., K respiratory phases: y 1,..., y K Assume f( r,t) is stationary during kth phase of data acquisition Relate y k to f k ( r) f( r,t k ) using physics and statistics From K data vectors y 1,..., y K, reconstruct object:? 39
40 Gated data model: Illustration f 1 y 1 f 2 y 2 Imaging System Reconstruction Method ˆf? f K y K 40
41 Gated data: Image reconstruction options Pool all data and ignore motion fast low noise variance motion induced blur Frame-wise: reconstruct each gate/frame separately: y k ˆf k simple high noise variance no motion blur (except within-gate motion) Frame-wise with post-reconstruction averaging (FW-PRA) map each reconstructed frame onto 1st frame, then average averaging should reduce noise should avoid motion blur if registration is accurate registration accuracy limited by noise in the individual gates FW-PRA with motion estimates from a separate modality use another modality (e.g., PET-CT) to estimate motion performance depends on consistency of motion between modalities Thorndyke et al., Med Phys,
42 Frame-wise with post-reconstruction averaging y 1 common target image 1 K ˆα 2 ˆf y 2 1 K ˆα K y K SSD registration 1 K Image Reconstruction Post Registration Consolidation 42
43 Forward model with motion Post-reconstruction averaging assumes an implicit model that relates the frames f 2,..., f K to the first frame f 1. We now make the (motion) model explicit: f k = W(α k ) f 1, k = 2,...,K. W(α k ) is the linear (!) transformation of the image values corresponding to motion α k. This model suggests additional image reconstruction approaches. 43
44 Linear interpolation and Nonrigid deformations B-spline interpolation model for continuous-space image: f( r) = f(x,y,z) = c nml β(x n)β(y m)β(z l) n,m,l Find coefficients c = {c nml } by prefiltering digital image f[n,m,l]. Nonrigid deformation of f : ) g( r) = f( T( r; α) = c nml β(t x ( r; α) n)β(t y ( r; α) m)β(t z ( r; α) l). n,m,l Resample warped image on grid (of target image): g( r j ) = c nml W r j ;n,m,l(α), j = 1,...,N n,m,l W r;n,m,l (α) β(t x ( r; α) n)β(t y ( r; α) m)β(t z ( r; α) l) In matrix-vector form, where g = {g( r j )} and f = { f[n,m,l]}: g = W(α)c, c = W 1 (0) f = g = W(α)W 1 (0) f. 44
45 Forward model with motion f 1 y 1 A ε 1 W(α 2 ) f 1 y 2 A ε 2 W(α K ) f 1 y K A ε K Motion Model Acquisition Process Measured Data 45
46 Gated data: More image reconstruction options Model-based image reconstruction with motion compensation given motion estimates, from FW-PRA or from separate modality, compensate for motion in reconstruction process. Qiao et al., PMB 2006; Taguchi et al., SPIE Model-based image reconstruction jointly with registration Jacobson & Fessler, IEEE NSS-MIC 2003, IEEE SSP 2003, ISBI 2006 Odille et al., MRM 2008 estimate jointly the first frame and the motion from all data Model-based image reconstruction with temporal regularization Mair et al., IEEE T-MI, 2006 estimate all frames and the motion between frames from all data 46
47 Model-based image reconstruction with motion compensation Given: motion estimates ˆα k for k = 2,...,K, from FW approach or from a separate modality, model for system physics / statistics: p(y k f k ) = p(y k W( ˆα k ) f 1 ). Perform penalized-likelihood (aka MAP) estimation of one image: ˆf 1 = argmax f 1 Ψ( f 1 ; ˆα 2,..., ˆα K ; y 1,..., y K ) Ψ( f 1 ; ˆα 2,..., ˆα K ; y 1,..., y K ) K logp(y k W( ˆα k ) f 1 ) βr( f 1 ). k=1 R( f) is optional regularization to control noise in ill-posed image reconstruction problems. For linear model with additive gaussian noise y k = A k f k +ε k : ˆf 1 = argmin f 1 K k=1 y k A k W( ˆα k ) f βr( f 1 ). 47
48 Motion compensated image reconstruction y 1 common source image argmax f1 Ψ( f 1 ; ˆα 2,..., ˆα K ; y 1,..., y K ) ˆα 2 ˆf 1 y 2 ˆα K y K Image Reconstruction SSD registration Post Registration Consolidation 48
49 Joint image reconstruction / registration Previous approach used possibly suboptimal motion estimates: ˆf 1 = argmax f 1 Ψ( f 1 ; ˆα 2,..., ˆα K ; y 1,..., y K ) Alternative: jointly estimate one image and K 1 deformation parameters: ( ˆf 1, ˆα 2,..., ˆα K ) = argmax f 1,α 2,...,α K Ψ( f 1 ; α 2,..., α K ; y 1,..., y K ) K Ψ( f 1 ; α 2,..., α K ; y 1,..., y K ) = logp(y k W(α k ) f 1 ) βr( f 1 ) k=1 Natural optimization strategy is to alternate between: updating image estimate ˆf 1 using current motion parameters, updating motion estimates { ˆα k } using current image estimate. Can initialize motion parameters using frame-wise method. 49
50 Joint estimation illustrated y 1 ˆα y 2 argmax Ψ( f 1 ; α 2,..., α K ; y 1,..., y K ) ˆf 1 y K Goal: find image estimate and motion parameters that best fit all measured data. 50
51 Motion-compensated temporal regularization Previous joint estimation approach: ( ˆf 1, ˆα 2,..., ˆα K ) = argmax f 1,α 2,...,α K Ψ( f 1 ; α 2,..., α K ; y 1,..., y K ) Alternative approach based on temporal regularization: ( ˆf 1,..., ˆf K ; ˆα 2,..., ˆα K ) = argmax f 1,..., ˆf K ;α 2,...,α K Ψ( f 1,..., f K ; α 2,..., α K ; y 1,..., y K ) Ψ( f 1,..., f K ; α 2,..., α K ; y 1,..., y K ) = K logp(y k f k ) βr( f k ) k=1 K γ f k+1 W(α k ) f k 2 k=2 } {{ } temporal regularization with motion effects Pro: no warp in log-likelihood. Con: more unknowns; γ choice? 51
52 Ten-Gate 3D PET Simulation (a) True Image (b) Mean Pre Image of Ungated Data 80K total counts/axial mm and 30% randoms (ECAT HR+), divided across 10 gates. Derived from 17 slices of real thorax anatomy. B-spline deformations (11x14x5x3 control grid), derived from helical CT scans at multiple inspirations (Matt Jacobson, 2006 thesis) 52
53 Sample Reconstructed Images (a) True Activity Image (b) PL with Known Motion (c) Fully Joint Estimation (JEDM) (d) Frame Wise Semi Statistical (FWPR PLC) (e) Frame Wise Post Averaging (FWPR PA) (f) Ungated PL 53
54 Mean Reconstructed Images (a) True Activity Image (b) PL with Known Motion (c) Fully Joint Estimation (JEDM) (d) Frame Wise Semi Statistical (FWPR PLC) (e) Frame Wise Post Averaging (FWPR PA) (f) Ungated PL 54
55 Lesion Recovery Comparison Study Parameters: 80K total counts/axial mm, 10 gates, 30% randoms Known Motion JEDM FWPR PA FWPR PLC Ungated 60 Cumulatively Averaged %Error No. of Realizations 55
56 Motion Tracking Performance Front Back motion Left Right motion Axial motion True Motion JEDM FWPR PLC FWPR PA True Motion JEDM FWPR PLC FWPR PA 20 True Motion JEDM FWPR PLC FWPR PA Lesion Position (mm) Lesion Position (mm) Lesion Position (mm) Gate No Gate No Gate No. 56
57 Regularization Using PET-CT Side Info. Relax regularization strength in neighborhood of lesion. (a) (b) 57
58 Regularization Using PET-CT Side Info. (cont d) (a) Uniformly Penalized Ungated PL (b) Uniformly Penalized JEDM (c) JEDM with Unpenalized ROI Uptake Error = 54.65% Uptake Error = 36.89% Uptake Error = 13.85% Very weak use of boundary side information = robust to mis-registration. 58
59 Summary Several possible methods for motion-compensated image reconstruction Model-based approaches such as joint estimation have potential Repeated motion estimation steps necessitate simple invertibility regularizers More work needed on algorithms, acceleration, evaluation,... 59
Estimating 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 informationImage Registration using Constrained Optimization
Image Registration using Constrained Optimization Jeongtae Kim and Jeffrey A. Fessler Information Electronics Dept., Ewha Womans University, Korea EECS Dept., The University of Michigan, Ann Arbor SIAM
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 informationIMAGE registration is a core tool in many medical imaging
JOURNAL OF SELECTED TOPICS IN SIGNAL PROCESSING, TO APPEAR. 1 A simple regularizer for B-spline nonrigid image registration that encourages local invertibility Se Young Chun, Student Member, IEEE, and
More informationMOTION ASPECTS IN JOINT IMAGE RECONSTRUCTION AND NONRIGID MOTION ESTIMATION
MOTION ASPECTS IN JOINT IMAGE RECONSTRUCTION AND NONRIGID MOTION ESTIMATION by Se Young Chun A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy (Electrical
More informationIEEE JOURNAL OF SELECTED TOPICS IN SIGNAL PROCESSING, VOL. 3, NO. 1, FEBRUARY
IEEE JOURNAL OF SELECTED TOPICS IN SIGNAL PROCESSING, VOL. 3, NO. 1, FEBRUARY 2009 159 A Simple Regularizer B-spline Nonrigid Image Registration That Encourages Local Invertibility Se Young Chun, Student
More informationNonrigid Registration using Free-Form Deformations
Nonrigid Registration using Free-Form Deformations Hongchang Peng April 20th Paper Presented: Rueckert et al., TMI 1999: Nonrigid registration using freeform deformations: Application to breast MR images
More informationGood Morning! Thank you for joining us
Good Morning! Thank you for joining us Deformable Registration, Contour Propagation and Dose Mapping: 101 and 201 Marc Kessler, PhD, FAAPM The University of Michigan Conflict of Interest I receive direct
More informationDeformable Image Registration, Contour Propagation and Dose Mapping: 101 and 201. Please do not (re)redistribute
Deformable Registration, Contour Deformable Registration, Contour Propagation and Dose Mapping: 101 and 201 Marc Kessler, PhD The University of Michigan Jean Pouliot, PhD University of California Learning
More informationAxial block coordinate descent (ABCD) algorithm for X-ray CT image reconstruction
Axial block coordinate descent (ABCD) algorithm for X-ray CT image reconstruction Jeffrey A. Fessler and Donghwan Kim EECS Department University of Michigan Fully 3D Image Reconstruction Conference July
More information2005 IEEE Nuclear Science Symposium Conference Record M10-2
25 IEEE Nuclear Science Symposium Conference Record M1-2 Estimating 3D Respiratory Motion from Orbiting Views Rongping Zeng 1, Jeffrey A. Fessler 1, and James M. Balter 2 rzeng@eecs.umich.edu, fessler@eecs.umich.edu,
More informationMOTION of a patient induces artifacts such as streaks
Motion-compensated Image Reconstruction for Cardiac CT with Sinogram-Based Motion Estimation JangHwanChoMember IEEE and Jeffrey A. Fessler Fellow IEEE Abstract Motion-compensated image reconstruction ()
More informationMethodological progress in image registration for ventilation estimation, segmentation propagation and multi-modal fusion
Methodological progress in image registration for ventilation estimation, segmentation propagation and multi-modal fusion Mattias P. Heinrich Julia A. Schnabel, Mark Jenkinson, Sir Michael Brady 2 Clinical
More informationImplementation of Advanced Image Guided Radiation Therapy
Image Acquisition Course Outline Principles, characteristics& applications of the available modalities Image Processing in the T x room Image guided treatment delivery What can / can t we do in the room
More informationDeviceless respiratory motion correction in PET imaging exploring the potential of novel data driven strategies
g Deviceless respiratory motion correction in PET imaging exploring the potential of novel data driven strategies Presented by Adam Kesner, Ph.D., DABR Assistant Professor, Division of Radiological Sciences,
More informationNon-Rigid Multimodal Medical Image Registration using Optical Flow and Gradient Orientation
M. HEINRICH et al.: MULTIMODAL REGISTRATION USING GRADIENT ORIENTATION 1 Non-Rigid Multimodal Medical Image Registration using Optical Flow and Gradient Orientation Mattias P. Heinrich 1 mattias.heinrich@eng.ox.ac.uk
More informationImage Segmentation and Registration
Image Segmentation and Registration Dr. Christine Tanner (tanner@vision.ee.ethz.ch) Computer Vision Laboratory, ETH Zürich Dr. Verena Kaynig, Machine Learning Laboratory, ETH Zürich Outline Segmentation
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 informationIs deformable image registration a solved problem?
Is deformable image registration a solved problem? Marcel van Herk On behalf of the imaging group of the RT department of NKI/AVL Amsterdam, the Netherlands DIR 1 Image registration Find translation.deformation
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 informationRobust Lung Ventilation Assessment
Fifth International Workshop on Pulmonary Image Analysis -75- Robust Lung Ventilation Assessment Sven Kabus 1, Tobias Klinder 1, Tokihiro Yamamoto 2, Paul J. Keall 3, Billy W. Loo, Jr. 4, and Cristian
More informationNonrigid Registration Using a Rigidity Constraint
Nonrigid Registration Using a Rigidity Constraint Marius Staring, Stefan Klein and Josien P.W. Pluim Image Sciences Institute, University Medical Center Utrecht, P.O. Box 85500, 3508 GA, Room Q0S.459,
More informationDense Image-based Motion Estimation Algorithms & Optical Flow
Dense mage-based Motion Estimation Algorithms & Optical Flow Video A video is a sequence of frames captured at different times The video data is a function of v time (t) v space (x,y) ntroduction to motion
More informationBasic fmri Design and Analysis. Preprocessing
Basic fmri Design and Analysis Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial filtering
More informationRespiratory Motion Compensation for Simultaneous PET/MR Based on Strongly Undersampled Radial MR Data
Respiratory Motion Compensation for Simultaneous PET/MR Based on Strongly Undersampled Radial MR Data Christopher M Rank 1, Thorsten Heußer 1, Andreas Wetscherek 1, and Marc Kachelrieß 1 1 German Cancer
More information3/27/2012 WHY SPECT / CT? SPECT / CT Basic Principles. Advantages of SPECT. Advantages of CT. Dr John C. Dickson, Principal Physicist UCLH
3/27/212 Advantages of SPECT SPECT / CT Basic Principles Dr John C. Dickson, Principal Physicist UCLH Institute of Nuclear Medicine, University College London Hospitals and University College London john.dickson@uclh.nhs.uk
More informationIntroduction to Image Super-resolution. Presenter: Kevin Su
Introduction to Image Super-resolution Presenter: Kevin Su References 1. S.C. Park, M.K. Park, and M.G. KANG, Super-Resolution Image Reconstruction: A Technical Overview, IEEE Signal Processing Magazine,
More informationEECS 556 Image Processing W 09. Interpolation. Interpolation techniques B splines
EECS 556 Image Processing W 09 Interpolation Interpolation techniques B splines What is image processing? Image processing is the application of 2D signal processing methods to images Image representation
More informationImage Registration with Local Rigidity Constraints
Image Registration with Local Rigidity Constraints Jan Modersitzki Institute of Mathematics, University of Lübeck, Wallstraße 40, D-23560 Lübeck Email: modersitzki@math.uni-luebeck.de Abstract. Registration
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 informationBasic principles of MR image analysis. Basic principles of MR image analysis. Basic principles of MR image analysis
Basic principles of MR image analysis Basic principles of MR image analysis Julien Milles Leiden University Medical Center Terminology of fmri Brain extraction Registration Linear registration Non-linear
More informationUnion of Learned Sparsifying Transforms Based Low-Dose 3D CT Image Reconstruction
Union of Learned Sparsifying Transforms Based Low-Dose 3D CT Image Reconstruction Xuehang Zheng 1, Saiprasad Ravishankar 2, Yong Long 1, Jeff Fessler 2 1 University of Michigan - Shanghai Jiao Tong University
More informationSPM8 for Basic and Clinical Investigators. Preprocessing. fmri Preprocessing
SPM8 for Basic and Clinical Investigators Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization Spatial
More information1 Trajectories. Class Notes, Trajectory Planning, COMS4733. Figure 1: Robot control system.
Class Notes, Trajectory Planning, COMS4733 Figure 1: Robot control system. 1 Trajectories Trajectories are characterized by a path which is a space curve of the end effector. We can parameterize this curve
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 informationOptimization of CT Simulation Imaging. Ingrid Reiser Dept. of Radiology The University of Chicago
Optimization of CT Simulation Imaging Ingrid Reiser Dept. of Radiology The University of Chicago Optimization of CT imaging Goal: Achieve image quality that allows to perform the task at hand (diagnostic
More informationImage Registration. Prof. Dr. Lucas Ferrari de Oliveira UFPR Informatics Department
Image Registration Prof. Dr. Lucas Ferrari de Oliveira UFPR Informatics Department Introduction Visualize objects inside the human body Advances in CS methods to diagnosis, treatment planning and medical
More informationLearning-based Neuroimage Registration
Learning-based Neuroimage Registration Leonid Teverovskiy and Yanxi Liu 1 October 2004 CMU-CALD-04-108, CMU-RI-TR-04-59 School of Computer Science Carnegie Mellon University Pittsburgh, PA 15213 Abstract
More informationGPU acceleration of 3D forward and backward projection using separable footprints for X-ray CT image reconstruction
GPU acceleration of 3D forward and backward projection using separable footprints for X-ray CT image reconstruction Meng Wu and Jeffrey A. Fessler EECS Department University of Michigan Fully 3D Image
More informationNon-Rigid Image Registration III
Non-Rigid Image Registration III CS6240 Multimedia Analysis Leow Wee Kheng Department of Computer Science School of Computing National University of Singapore Leow Wee Kheng (CS6240) Non-Rigid Image Registration
More informationSuper-resolution Reconstruction of Fetal Brain MRI
Super-resolution Reconstruction of Fetal Brain MRI Ali Gholipour and Simon K. Warfield Computational Radiology Laboratory Children s Hospital Boston, Harvard Medical School Worshop on Image Analysis for
More informationVariational Lung Registration With Explicit Boundary Alignment
Variational Lung Registration With Explicit Boundary Alignment Jan Rühaak and Stefan Heldmann Fraunhofer MEVIS, Project Group Image Registration Maria-Goeppert-Str. 1a, 23562 Luebeck, Germany {jan.ruehaak,stefan.heldmann}@mevis.fraunhofer.de
More informationBiomedical Imaging Registration Trends and Applications. Francisco P. M. Oliveira, João Manuel R. S. Tavares
Biomedical Imaging Registration Trends and Applications Francisco P. M. Oliveira, João Manuel R. S. Tavares tavares@fe.up.pt, www.fe.up.pt/~tavares Outline 1. Introduction 2. Spatial Registration of (2D
More informationFunctional MRI in Clinical Research and Practice Preprocessing
Functional MRI in Clinical Research and Practice Preprocessing fmri Preprocessing Slice timing correction Geometric distortion correction Head motion correction Temporal filtering Intensity normalization
More informationUvA-DARE (Digital Academic Repository) Motion compensation for 4D PET/CT Kruis, M.F. Link to publication
UvA-DARE (Digital Academic Repository) Motion compensation for 4D PET/CT Kruis, M.F. Link to publication Citation for published version (APA): Kruis, M. F. (2014). Motion compensation for 4D PET/CT General
More informationHybrid Spline-based Multimodal Registration using a Local Measure for Mutual Information
Hybrid Spline-based Multimodal Registration using a Local Measure for Mutual Information Andreas Biesdorf 1, Stefan Wörz 1, Hans-Jürgen Kaiser 2, Karl Rohr 1 1 University of Heidelberg, BIOQUANT, IPMB,
More informationSpatio-Temporal Registration of Biomedical Images by Computational Methods
Spatio-Temporal Registration of Biomedical Images by Computational Methods Francisco P. M. Oliveira, João Manuel R. S. Tavares tavares@fe.up.pt, www.fe.up.pt/~tavares Outline 1. Introduction 2. Spatial
More informationWhite Pixel Artifact. Caused by a noise spike during acquisition Spike in K-space <--> sinusoid in image space
White Pixel Artifact Caused by a noise spike during acquisition Spike in K-space sinusoid in image space Susceptibility Artifacts Off-resonance artifacts caused by adjacent regions with different
More informationArtefakt-resistente Bewegungsschätzung für die bewegungskompensierte CT
Artefakt-resistente Bewegungsschätzung für die bewegungskompensierte CT Marcus Brehm 1,2, Thorsten Heußer 1, Pascal Paysan 3, Markus Oehlhafen 3, and Marc Kachelrieß 1,2 1 German Cancer Research Center
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 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 informationAdvanced methods for image reconstruction in fmri
Advanced methods for image reconstruction in fmri Jeffrey A. Fessler EECS Department The University of Michigan Regional Symposium on MRI Sep. 28, 27 Acknowledgements: Doug Noll, Brad Sutton, Outline MR
More informationEPI Data Are Acquired Serially. EPI Data Are Acquired Serially 10/23/2011. Functional Connectivity Preprocessing. fmri Preprocessing
Functional Connectivity Preprocessing Geometric distortion Head motion Geometric distortion Head motion EPI Data Are Acquired Serially EPI Data Are Acquired Serially descending 1 EPI Data Are Acquired
More informationFMA901F: Machine Learning Lecture 3: Linear Models for Regression. Cristian Sminchisescu
FMA901F: Machine Learning Lecture 3: Linear Models for Regression Cristian Sminchisescu Machine Learning: Frequentist vs. Bayesian In the frequentist setting, we seek a fixed parameter (vector), with value(s)
More informationAn introduction to interpolation and splines
An introduction to interpolation and splines Kenneth H. Carpenter, EECE KSU November 22, 1999 revised November 20, 2001, April 24, 2002, April 14, 2004 1 Introduction Suppose one wishes to draw a curve
More informationNonrigid Registration with Adaptive, Content-Based Filtering of the Deformation Field
Nonrigid Registration with Adaptive, Content-Based Filtering of the Deformation Field Marius Staring*, Stefan Klein and Josien P.W. Pluim Image Sciences Institute, University Medical Center Utrecht, P.O.
More informationConstrained and Unconstrained Optimization
Constrained and Unconstrained Optimization Carlos Hurtado Department of Economics University of Illinois at Urbana-Champaign hrtdmrt2@illinois.edu Oct 10th, 2017 C. Hurtado (UIUC - Economics) Numerical
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 informationSimplified statistical image reconstruction algorithm for polyenergetic X-ray CT. y i Poisson I i (E)e R } where, b i
implified statistical image reconstruction algorithm for polyenergetic X-ray C omesh rivastava, tudent member, IEEE, and Jeffrey A. Fessler, enior member, IEEE Abstract In X-ray computed tomography (C),
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 informationRespiratory Motion Estimation using a 3D Diaphragm Model
Respiratory Motion Estimation using a 3D Diaphragm Model Marco Bögel 1,2, Christian Riess 1,2, Andreas Maier 1, Joachim Hornegger 1, Rebecca Fahrig 2 1 Pattern Recognition Lab, FAU Erlangen-Nürnberg 2
More informationCluster of Workstation based Nonrigid Image Registration Using Free-Form Deformation
Cluster of Workstation based Nonrigid Image Registration Using Free-Form Deformation Xiaofen Zheng, Jayaram K. Udupa, and Xinjian Chen Medical Image Processing Group, Department of Radiology 423 Guardian
More informationPenalized-Likelihood Reconstruction for Sparse Data Acquisitions with Unregistered Prior Images and Compressed Sensing Penalties
Penalized-Likelihood Reconstruction for Sparse Data Acquisitions with Unregistered Prior Images and Compressed Sensing Penalties J. W. Stayman* a, W. Zbijewski a, Y. Otake b, A. Uneri b, S. Schafer a,
More informationNonrigid Surface Modelling. and Fast Recovery. Department of Computer Science and Engineering. Committee: Prof. Leo J. Jia and Prof. K. H.
Nonrigid Surface Modelling and Fast Recovery Zhu Jianke Supervisor: Prof. Michael R. Lyu Committee: Prof. Leo J. Jia and Prof. K. H. Wong Department of Computer Science and Engineering May 11, 2007 1 2
More informationImage Registration + Other Stuff
Image Registration + Other Stuff John Ashburner Pre-processing Overview fmri time-series Motion Correct Anatomical MRI Coregister m11 m 21 m 31 m12 m13 m14 m 22 m 23 m 24 m 32 m 33 m 34 1 Template Estimate
More informationMotion. 1 Introduction. 2 Optical Flow. Sohaib A Khan. 2.1 Brightness Constancy Equation
Motion Sohaib A Khan 1 Introduction So far, we have dealing with single images of a static scene taken by a fixed camera. Here we will deal with sequence of images taken at different time intervals. Motion
More informationTG 132: Use of Image Registration and Fusion in RT
TG 132: Use of Image Registration and Fusion in RT Kristy K Brock, PhD, DABR, FAAPM Associate Professor Department of Radiation Oncology, University of Michigan Chair, AAPM TG 132: Image Registration and
More informationREAL-TIME ADAPTIVITY IN HEAD-AND-NECK AND LUNG CANCER RADIOTHERAPY IN A GPU ENVIRONMENT
REAL-TIME ADAPTIVITY IN HEAD-AND-NECK AND LUNG CANCER RADIOTHERAPY IN A GPU ENVIRONMENT Anand P Santhanam Assistant Professor, Department of Radiation Oncology OUTLINE Adaptive radiotherapy for head and
More informationAlgorithms for Markov Random Fields in Computer Vision
Algorithms for Markov Random Fields in Computer Vision Dan Huttenlocher November, 2003 (Joint work with Pedro Felzenszwalb) Random Field Broadly applicable stochastic model Collection of n sites S Hidden
More informationInterpolating and approximating scattered 3D-data with hierarchical tensor product B-splines
Interpolating and approximating scattered 3D-data with hierarchical tensor product B-splines Günther Greiner Kai Hormann Abstract In this note we describe surface reconstruction algorithms based on optimization
More informationWilliam Yang Group 14 Mentor: Dr. Rogerio Richa Visual Tracking of Surgical Tools in Retinal Surgery using Particle Filtering
Mutual Information Computation and Maximization Using GPU Yuping Lin and Gérard Medioni Computer Vision and Pattern Recognition Workshops (CVPR) Anchorage, AK, pp. 1-6, June 2008 Project Summary and Paper
More informationRespiratory Motion Compensation for C-arm CT Liver Imaging
Respiratory Motion Compensation for C-arm CT Liver Imaging Aline Sindel 1, Marco Bögel 1,2, Andreas Maier 1,2, Rebecca Fahrig 3, Joachim Hornegger 1,2, Arnd Dörfler 4 1 Pattern Recognition Lab, FAU Erlangen-Nürnberg
More informationBiomedical Image Analysis based on Computational Registration Methods. João Manuel R. S. Tavares
Biomedical Image Analysis based on Computational Registration Methods João Manuel R. S. Tavares tavares@fe.up.pt, www.fe.up.pt/~tavares Outline 1. Introduction 2. Methods a) Spatial Registration of (2D
More information1. Suppose that the equation F (x, y, z) = 0 implicitly defines each of the three variables x, y, and z as functions of the other two:
Final Solutions. Suppose that the equation F (x, y, z) implicitly defines each of the three variables x, y, and z as functions of the other two: z f(x, y), y g(x, z), x h(y, z). If F is differentiable
More informationRIGID IMAGE REGISTRATION
RIGID IMAGE REGISTRATION Duygu Tosun-Turgut, Ph.D. Center for Imaging of Neurodegenerative Diseases Department of Radiology and Biomedical Imaging duygu.tosun@ucsf.edu What is registration? Image registration
More information8/2/2016. Measures the degradation/distortion of the acquired image (relative to an ideal image) using a quantitative figure-of-merit
Ke Li Assistant Professor Department of Medical Physics and Department of Radiology School of Medicine and Public Health, University of Wisconsin-Madison This work is partially supported by an NIH Grant
More informationRegularization and Markov Random Fields (MRF) CS 664 Spring 2008
Regularization and Markov Random Fields (MRF) CS 664 Spring 2008 Regularization in Low Level Vision Low level vision problems concerned with estimating some quantity at each pixel Visual motion (u(x,y),v(x,y))
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 information3-D Compounding of B-Scan Ultrasound Images
3-D Compounding of B-Scan Ultrasound Images Jochen F. Krücker, Charles R. Meyer, Theresa A. Tuthill, Gerald L. LeCarpentier, J. Brian Fowlkes, Paul L. Carson University of Michigan, Dept. of Radiology,
More informationImage Registration Lecture 4: First Examples
Image Registration Lecture 4: First Examples Prof. Charlene Tsai Outline Example Intensity-based registration SSD error function Image mapping Function minimization: Gradient descent Derivative calculation
More informationMultimodal Elastic Image Matching
Research results based on my diploma thesis supervised by Prof. Witsch 2 and in cooperation with Prof. Mai 3. 1 February 22 nd 2011 1 Karlsruhe Institute of Technology (KIT) 2 Applied Mathematics Department,
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 1 3.1 Linearization and Optimization of Functions of Vectors 1 Problem Notation 2 Outline 3.1.1 Linearization 3.1.2 Optimization of Objective Functions 3.1.3 Constrained
More informationThe Near Future in Cardiac CT Image Reconstruction
SCCT 2010 The Near Future in Cardiac CT Image Reconstruction Marc Kachelrieß Institute of Medical Physics (IMP) Friedrich-Alexander Alexander-University Erlangen-Nürnberg rnberg www.imp.uni-erlangen.de
More informationDigital Image Processing
Digital Image Processing SPECIAL TOPICS CT IMAGES Hamid R. Rabiee Fall 2015 What is an image? 2 Are images only about visual concepts? We ve already seen that there are other kinds of image. In this lecture
More informationSuper-Resolution. Many slides from Miki Elad Technion Yosi Rubner RTC and more
Super-Resolution Many slides from Mii Elad Technion Yosi Rubner RTC and more 1 Example - Video 53 images, ratio 1:4 2 Example Surveillance 40 images ratio 1:4 3 Example Enhance Mosaics 4 5 Super-Resolution
More informationImage 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 informationBildverarbeitung für die Medizin 2007
Bildverarbeitung für die Medizin 2007 Image Registration with Local Rigidity Constraints Jan Modersitzki Institute of Mathematics, University of Lübeck, Wallstraße 40, D-23560 Lübeck 1 Summary Registration
More informationSIFT: SCALE INVARIANT FEATURE TRANSFORM SURF: SPEEDED UP ROBUST FEATURES BASHAR ALSADIK EOS DEPT. TOPMAP M13 3D GEOINFORMATION FROM IMAGES 2014
SIFT: SCALE INVARIANT FEATURE TRANSFORM SURF: SPEEDED UP ROBUST FEATURES BASHAR ALSADIK EOS DEPT. TOPMAP M13 3D GEOINFORMATION FROM IMAGES 2014 SIFT SIFT: Scale Invariant Feature Transform; transform image
More informationMR IMAGE SEGMENTATION
MR IMAGE SEGMENTATION Prepared by : Monil Shah What is Segmentation? Partitioning a region or regions of interest in images such that each region corresponds to one or more anatomic structures Classification
More informationMotion artifact detection in four-dimensional computed tomography images
Motion artifact detection in four-dimensional computed tomography images G Bouilhol 1,, M Ayadi, R Pinho, S Rit 1, and D Sarrut 1, 1 University of Lyon, CREATIS; CNRS UMR 5; Inserm U144; INSA-Lyon; University
More informationA Fast GPU-Based Approach to Branchless Distance-Driven Projection and Back-Projection in Cone Beam CT
A Fast GPU-Based Approach to Branchless Distance-Driven Projection and Back-Projection in Cone Beam CT Daniel Schlifske ab and Henry Medeiros a a Marquette University, 1250 W Wisconsin Ave, Milwaukee,
More informationDigital Image Processing COSC 6380/4393
Digital Image Processing COSC 6380/4393 Lecture 21 Nov 16 th, 2017 Pranav Mantini Ack: Shah. M Image Processing Geometric Transformation Point Operations Filtering (spatial, Frequency) Input Restoration/
More informationOverview of Proposed TG-132 Recommendations
Overview of Proposed TG-132 Recommendations Kristy K Brock, Ph.D., DABR Associate Professor Department of Radiation Oncology, University of Michigan Chair, AAPM TG 132: Image Registration and Fusion Conflict
More informationThe goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a
The goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a coordinate system and then the measuring of the point with
More informationHYBRID MULTISCALE LANDMARK AND DEFORMABLE IMAGE REGISTRATION. Dana Paquin. Doron Levy. Lei Xing. (Communicated by Yang Kuang)
MATHEMATICAL BIOSCIENCES http://www.mbejournal.org/ AND ENGINEERING Volume 4, Number 4, October 2007 pp. 711 737 HYBRID MULTISCALE LANDMARK AND DEFORMABLE IMAGE REGISTRATION Dana Paquin Department of Mathematics,
More informationCS 565 Computer Vision. Nazar Khan PUCIT Lectures 15 and 16: Optic Flow
CS 565 Computer Vision Nazar Khan PUCIT Lectures 15 and 16: Optic Flow Introduction Basic Problem given: image sequence f(x, y, z), where (x, y) specifies the location and z denotes time wanted: displacement
More informationComputational Neuroanatomy
Computational Neuroanatomy John Ashburner john@fil.ion.ucl.ac.uk Smoothing Motion Correction Between Modality Co-registration Spatial Normalisation Segmentation Morphometry Overview fmri time-series kernel
More informationCentral Slice Theorem
Central Slice Theorem Incident X-rays y f(x,y) R x r x Detected p(, x ) The thick line is described by xcos +ysin =R Properties of Fourier Transform F [ f ( x a)] F [ f ( x)] e j 2 a Spatial Domain Spatial
More informationLecture VIII. Global Approximation Methods: I
Lecture VIII Global Approximation Methods: I Gianluca Violante New York University Quantitative Macroeconomics G. Violante, Global Methods p. 1 /29 Global function approximation Global methods: function
More informationImage registration for motion estimation in cardiac CT
Image registration for motion estimation in cardiac CT Bibo Shi a, Gene Katsevich b, Be-Shan Chiang c, Alexander Katsevich d, and Alexander Zamyatin c a School of Elec. Engi. and Comp. Sci, Ohio University,
More information