Advanced methods for image reconstruction in fmri
|
|
- Jeffery Williamson
- 6 years ago
- Views:
Transcription
1 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,
2 Outline MR image reconstruction problem description Overview of image reconstruction methods Model-based image reconstruction Iterations and computation (NUFFT etc.) Regularization Myths about iterative reconstruction Field inhomogeneity correction Parallel (sensitivity encoded) imaging Image reconstruction toolbox: fessler 1
3 Why Iterative Image Reconstruction? Statistical modeling may reduce noise Incorporate prior information, e.g.: support constraints (piecewise) smoothness phase constraints No density compensation needed Non-Fourier physical effects such as field inhomogeneity Incorporation of coil sensitivity maps Improved results for under-sampled trajectories (?)... ( Avoiding k-space interpolation is not a compelling reason!) 2
4 Primary drawbacks of Iterative Methods Choosing regularization parameter(s) Algorithm speed 3
5 Example: Iterative Reconstruction under B 4
6 Introduction to Reconstruction 5
7 Standard MR Image Reconstruction MR k space data Reconstructed Image Cartesian sampling in k-space. An inverse FFT. End of story. Commercial MR system quotes 4 FFTs (256 2 ) per second. 6
8 Non-Cartesian MR Image Reconstruction k-space data y = (y 1,...,y M ) image f( r) k y k x k-space trajectory: κ(t) = (k x (t),k y (t)) = spatial coordinates: d r R 7
9 Textbook MRI Measurement Model Ignoring lots of things, the standard measurement model is: y i = s(t i )+noise i, s(t) = Z i = 1,...,M f( r)e ı2π κ(t) r d r = F( κ(t)). r: spatial coordinates κ(t): k-space trajectory of the MR pulse sequence f( r): object s unknown transverse magnetization F( κ): Fourier transform of f( r). We get noisy samples of this! e ı2π κ(t) r provides spatial information = Nobel Prize Goal of image reconstruction: find f( r) from measurements {y i } M i=1. The unknown object f( r) is a continuous-space function, but the recorded measurements y = (y 1,...,y M ) are finite. Under-determined (ill posed) problem = no canonical solution. All MR scans provide only partial k-space data. 8
10 Image Reconstruction Strategies Continuous-continuous formulation Pretend that a continuum of measurements are available: F( κ) = Z f( r)e ı2π κ r d r. The solution is an inverse Fourier transform: f( r) = Now discretize the integral solution: ˆf( r) = Z F( κ)e ı2π κ r d κ. M F( κ i )e ı2π κ i r w i i=1 M y i w i e ı2π κ i r, i=1 where w i values are sampling density compensation factors. Numerous methods for choosing w i values in the literature. For Cartesian sampling, using w i = 1/N suffices, and the summation is an inverse FFT. For non-cartesian sampling, replace summation with gridding. 9
11 Continuous-discrete formulation Use many-to-one linear model: y = A f +ε, where A : L 2 (R d ) C M. Minimum norm solution (cf. natural pixels ): min ˆf subject to y=a ˆf ˆf ˆf = A (AA ) 1 y = M i=1 c i e ı2π κ i r, where AA c = y. Discrete-discrete formulation Assume parametric model for object: f( r) = N j=1 f j p j ( r). Estimate parameter vector f = ( f 1,..., f N ) from data vector y.
12 Model-Based Image Reconstruction: Overview 11
13 Model-Based Image Reconstruction MR signal equation with more complete physics: s(t) = Z f( r)s coil ( r)e ıω( r)t e R 2 ( r)t e ı2π κ(t) r d r y i = s(t i )+noise i, i = 1,...,M s coil ( r) Receive-coil sensitivity pattern(s) (for SENSE) ω( r) Off-resonance frequency map (due to field inhomogeneity / magnetic susceptibility) R 2 ( r) Relaxation map Other physical factors (?) Eddy current effects; in κ(t) Concomitant gradient terms Chemical shift Motion Goal? (it depends) 12
14 Field Inhomogeneity-Corrected Reconstruction s(t) = Z f( r)s coil ( r)e ıω( r)t e R 2 ( r)t e ı2π κ(t) r d r Goal: reconstruct f( r) given field map ω( r). (Assume all other terms are known or unimportant.) Motivation Essential for functional MRI of brain regions near sinus cavities! (Sutton et al., ISMRM 21; T-MI 23) 13
15 Sensitivity-Encoded (SENSE) Reconstruction s(t) = Z f( r)s coil ( r)e ıω( r)t e R 2 ( r)t e ı2π κ(t) r d r Goal: reconstruct f( r) given sensitivity maps s coil ( r). (Assume all other terms are known or unimportant.) Can combine SENSE with field inhomogeneity correction easily. (Sutton et al., ISMRM 21, Olafsson et al., ISBI 26) 14
16 Joint Estimation of Image and Field-Map s(t) = Z f( r)s coil ( r)e ıω( r)t e R 2 ( r)t e ı2π κ(t) r d r Goal: estimate both the image f( r) and the field map ω( r) (Assume all other terms are known or unimportant.) Analogy: joint estimation of emission image and attenuation map in PET. (Sutton et al., ISMRM Workshop, 21; ISBI 22; ISMRM 22; ISMRM 23; MRM 24) 15
17 s(t) = Z The Kitchen Sink f( r)s coil ( r)e ıω( r)t e R 2 ( r)t e ı2π κ(t) r d r Goal: estimate image f( r), field map ω( r), and relaxation map R 2 ( r) Requires suitable k-space trajectory. (Sutton et al., ISMRM 22; Twieg, MRM, 23) 16
18 Estimation of Dynamic Maps s(t) = Z f( r)s coil ( r)e ıω( r)t e R 2 ( r)t e ı2π κ(t) r d r Goal: estimate dynamic field map ω( r) and BOLD effect R 2 ( r) given baseline image f( r) in fmri. Motion... 17
19 Model-Based Image Reconstruction: Details 18
20 Basic Signal Model y i = s(t i )+ε i, E[y i ] = s(t i ) = Goal: reconstruct f( r) from y = (y 1,...,y M ). Series expansion of unknown object: f( r) N j=1 Z f( r)e ı2π κ i r d r f j p( r r j ) usually 2D rect functions. Substituting into signal model yields Z [ ] N E[y i ] = f j p( r r j ) e ı2π κ i r d r = j=1 = N [ Z p( r r j )e ı2π κ i d r] r f j j=1 N a i j f j, a i j = P( κ i )e ı2π κ i r j, p( r) FT P( κ). j=1 Discrete-discrete measurement model with system matrix A = {a i j }: y = A f + ε. Goal: estimate coefficients (pixel values) f = ( f 1,..., f N ) from y. 19
21 Small Pixel Size Need Not Matter 127 x true N= N= N= N= N=
22 Profiles true N=32 N=64 N=128 N=256 N=512 f(x,) horizontal position [mm] 21
23 Regularized Least-Squares Estimation Estimate object by minimizing a cost function: ˆf = argmin f C N Ψ( f), Ψ( f) = y A f 2 + αr( f) data fit term y A f 2 corresponds to negative log-likelihood of Gaussian distribution regularizing term R( f) controls noise by penalizing roughness, e.g. : R( f) Z f 2 d r regularization parameter α > controls tradeoff between spatial resolution and noise αr( f) Equivalent to Bayesian MAP estimation with prior e Issues: choosing R( f) choosing α computing minimizer rapidly. 22
24 Quadratic regularization 1D example: squared differences between neighboring pixel values: R( f) = N 1 j=2 2 f j f j 1 2. In matrix-vector notation, R( f) = 1 2 C f 2 where C =, so C f = f 2 f 1. f N f N 1 For 2D and higher-order differences, modify differencing matrix C. Leads to closed-form solution: ˆf = argmin y A f 2 + α C f 2 f = [ A A+αC C ] 1 A y.. (a formula of limited practical use for computing ˆf ) 23
25 Choosing the Regularization Parameter Spatial resolution analysis (Fessler & Rogers, IEEE T-IP, 1996): ˆf = [ A A+αC C ] 1 A y ] E[ ˆf = [ A A+αC C ] 1 A E[y] ] E[ ˆf = [ A A+αC C ] 1 A A f }{{} blur A A and C C are Toeplitz = blur is approximately shift-invariant. Frequency response of blur: L(ω) = H(ω k ) = FFT(A Ae j ) (lowpass) R(ω k ) = FFT(C C e j ) (highpass) H(ω) H(ω)+αR(ω) Adjust α to achieve desired spatial resolution. 24
26 Spatial Resolution Example 1 A A e j 1 α C C e j 1 PSF π H(ω) π R(ω) π L=H/(H+R) ω Y ω Y ω Y π π π ω X π π π ω X π π π Radial k-space trajectory, FWHM of PSF is 1.2 pixels ω X 25
27 Spatial Resolution Example: Profiles 1 x 15 H(ω) R(ω) L(ω).8.6 π π ω 26
28 Resolution/noise tradeoffs Noise analysis: } Cov{ ˆf = [ A A+αC C ] 1 A Cov{y} A [ A A+αC C ] 1 Using circulant approximations to A A and C C yields: Var { ˆf j } σ 2 ε k H(ω k ) = FFT(A Ae j ) (lowpass) R(ω k ) = FFT(C C e j ) (highpass) H(ω k ) (H(ω k )+αr(ω k )) 2 = Predicting reconstructed image noise requires just 2 FFTs. (cf. gridding approach?) Adjust α to achieve desired spatial resolution / noise tradeoff. 27
29 Resolution/Noise Tradeoff Example Relative standard deviation α Under sampled radial Nyquist sampled radial Cartesian α PSF FWHM [pixels] In short: one can choose α rapidly and predictably for quadratic regularization. 28
30 Iterative Minimization by Conjugate Gradients Choose initial guess f () (e.g., fast conjugate phase / gridding). Iteration (unregularized): g (n) = Ψ ( f (n)) = A (A f (n) y) gradient p (n) = Pg (n) precondition, n = γ n = g (n), p (n) g (n 1), p (n 1), n > d (n) = p (n) + γ n d (n 1) search direction v (n) = Ad (n) α n = d (n), g (n) / v (n) 2 step size f (n+1) = f (n) + α n d (n) update Bottlenecks: computing A f (n) and A r. A is too large to store explicitly (not sparse) Even if A were stored, directly computing A f is O(MN) per iteration, whereas FFT is only O(M logm). 29
31 [A f] i = Computing A f Rapidly N a i j f j = P( κ i ) j=1 N e ı2π κ i r j f j, j=1 i = 1,...,M Pixel locations { r j } are uniformly spaced k-space locations { κ i } are unequally spaced = needs nonuniform fast Fourier transform (NUFFT) 3
32 NUFFT (Type 2) Compute over-sampled FFT of equally-spaced signal samples Interpolate onto desired unequally-spaced frequency locations Dutt & Rokhlin, SIAM JSC, 1993, Gaussian bell interpolator Fessler & Sutton, IEEE T-SP, 23, min-max interpolator and min-max optimized Kaiser-Bessel interpolator. NUFFT toolbox: fessler/code 1 X(ω) 5 π π/2 ω? π/2 π 31
33 Worst-Case NUFFT Interpolation Error Maximum error for K/N= E max Min Max (uniform) Gaussian (best) Min Max (best L=2) Kaiser Bessel (best) Min Max (L=13, β=1 fit) J 32
34 Further Acceleration using Toeplitz Matrices Cost-function gradient: g (n) = A (A f (n) y) = T f (n) b, where T A A, b A y. In the absence of field inhomogeneity, the Gram matrix T is Toeplitz: [A A] jk = M P( κ i ) 2 e ı2π κ i ( r j r k ). i=1 Computing T f (n) requires an ordinary (2 over-sampled) FFT. (Chan & Ng, SIAM Review, 1996) In 2D: block Toeplitz with Toeplitz blocks (BTTB). Precomputing the first column of T and b requires a couple NUFFTs. (Wajer, ISMRM 21, Eggers ISMRM 22, Liu ISMRM 25) This formulation seems ideal for hardware FFT systems. 33
35 Toeplitz Acceleration Example: image. radial trajectory, 2 angular under-sampling. True CP 1 True CP.6 CG NUFFT CG Toeplitz CG NUFFT CG Toeplitz (Iterative provides reduced aliasing energy.) 34
36 Toeplitz Acceleration Method A Dy b = A y T 2 iter Total Time NRMS (5dB) Conj. Phase % CG-NUFFT % CG-Toeplitz % Toeplitz aproach reduces CPU time by more than 2 on conventional workstation (Xeon 3.4GHz) Eliminates k-space interpolations = ideal for FFT hardware No SNR compromise CG reduces NRMS error relative to CP, but 15 slower... (More dramatic improvements seen in fmri when correcting for field inhomogeneity.) 35
37 NUFFT with Field Inhomogeneity? Combine NUFFT with min-max temporal interpolator (Sutton et al., IEEE T-MI, 23) (forward version of time segmentation, Noll, T-MI, 1991) Recall signal model including field inhomogeneity: s(t) = Z f( r)e ıω( r)t e ı2π κ(t) r d r. Temporal interpolation approximation (aka time segmentation ): e ıω( r)t L a l (t)e ıω( r)τ l l=1 for min-max optimized temporal interpolation functions {a l ( )} L l=1. L Z [ ] s(t) a l (t) f( r)e ıω( r)τ l e ı2π κ(t) r d r l=1 Linear combination of L NUFFT calls. 36
38 Field Corrected Reconstruction Example Simulation using known field map ω( r). 37
39 Simulation Quantitative Comparison Computation time? NRMSE between ˆf and f true? Reconstruction Method Time (s) NRMSE NRMSE complex magnitude No Correction Full Conjugate Phase Fast Conjugate Phase Fast Iterative (1 iters) Exact Iterative (1 iters)
40 Human Data: Field Correction 39
41 Acceleration using Toeplitz Approximations In the presence of field inhomogeneity, the system matrix is: a i j = P( κ i )e ıω( r j)t i e ı2π κ i r j The Gram matrix T = A A is not Toeplitz: [A A] jk = Approximation ( time segmentation ): M P( κ i ) 2 e ı2π κ i ( r j r k ) e ı (ω( r j ) ω( r k ))t i. i=1 e ı L (ω( r j ) ω( r k ))t i b il e ı (ω( r j ) ω( r k ))τ l l=1 T = A L A D D l diag { } e ıω( r j)τ l lt l D l, [T l=1 l ] jk M i=1 P( κ i ) 2 b il e ı2π κ i ( r j r k ). Each T l is Toeplitz = T f using L pairs of FFTs. (Fessler et al., IEEE T-SP, Sep. 25, brain imaging special issue) 4
42 Toeplitz Results 1 Image and support 2.5 Uncorrected Conj. Phase, L= Fieldmap: Brain Hz CG NUFFT L=6 CG Toeplitz L=8 41
43 Toeplitz Acceleration Precomputation NRMS % vs SNR Method L B,C A Dy b = A y T l 15 iter Total Time 5 db 4 db 3 db 2 db Conj. Phase CG-NUFFT CG-Toeplitz Reduces CPU time by 2 on conventional workstation (Mac G5) No SNR compromise Eliminates k-space interpolations = ideal for FFT hardware 42
44 Joint Field-Map / Image Reconstruction Signal model: After discretization: y i = s(t i )+ε i, s(t) = Z f( r)e ıω( r)t e ı2π κ(t) r d r. y = A(ω) f + ε, a i j (ω) = P( κ i )e ıω jt i e ı2π κ i r j. Joint estimation via regularized (nonlinear) least-squares: ( ˆf, ˆω) = argmin f C N,ω R N y A(ω) f 2 + β 1 R 1 ( f)+β 2 R 2 (ω). Alternating minimization: Using current estimate of fieldmap ˆω, update ˆf using CG algorithm. Using current estimate ˆf of image, update fieldmap ˆω using gradient descent. Use spiral-in / spiral-out sequence or racetrack EPI. (Sutton et al., MRM, 24) 43
45 Joint Estimation Example (a) uncorr., (b) std. map, (c) joint map, (d) T1 ref, (e) using std, (f) using joint. 44
46 Activation Results: Static vs Dynamic Field Maps 45
47 Functional results for the two reconstructions for 3 human subjects. Reconstruction using the standard field map for (a) subject 1, (b) subject 2, and (c) subject 3. Reconstruction using the jointly estimated field map for (d) subject 1, (e) subject 2, and (f) subject 3. Number of pixels with correlation coefficients higher than thresholds for (g) subject 1, (h) subject 2, and (i) subject 3. Take home message: dynamic field mapping is possible, using iterative reconstruction as an essential tool. (Standard field maps based on echo-time differences work poorly for spiral-in / spiral-out sequences due to phase discrepancies.)
48 Tracking Respiration-Induced Field Changes 47
49 Choosing α is difficult Myths Sample density weighting is desirable 48
50 Sampling density weighted LS Some researchers recommend using a weighted LS cost function: Ψ( f) = y A f 2 W where the weighting matrix W is related to the k-space sample density compensation factors (DCF). Purported benefits: Faster convergence Better conditioning But, Gauss-Markov theorem from statistical estimation theory states that lowest estimator variance is realized by choosing W = σ 2 ε I, the inverse of the data noise covariance. 49
51 Resolution/Noise Tradeoff: Example with Weighting Relative standard deviation Nyquist sampled radial, DCF weights Nyquist sampled radial Cartesian cond(a WA) e+6 cond(a A) e PSF FWHM [pixels] 5
52 Don t just take it from me... Fig. 5 of Gabr et al., MRM, Dec
53 Acceleration via Weighting? Fig. 8 of Gabr et al., MRM, Dec. 26. Zero initialization! 52
54 Acceleration via Initialization 1 8 Unweighted, initialized Weighted, initialized Unweighted, CP initialized NRMS Error Iteration 53
55 Parallel Imaging 54
56 Sensitivity encoded (SENSE) imaging Use multiple receive coils (requires multiple RF channels). Exploit spatial localization of sensitivity pattern of each coil. Note: at 1.5T, RF is about 6MHz. = RF wavelength is about m/s/6 1 6 Hz = 5 meters RF coil sensitivity patterns 1 Array coil images 1 Regularized sensitivity maps Pruessmann et al., MRM,
57 SENSE Model Multiple coil data: y li = s l (t i )+ε li, s l (t) = Z f( r)s coil l ( r)e ı2π κ(t) r d r, l = 1,...,L = N coil Goal: reconstruct f( r) from coil data y 1,..., y L given sensitivity maps { s coil l ( r) } L l=1. Benefit: reduced scan time. Left: sum of squares; right: SENSE. 56
58 SENSE Reconstruction Signal model: s l (t) = Z f( r)s coil l ( r)e ı2π κ(t) r d r Discretized form: y l = AD l f + ε l, l = 1,...,L, where A is the usual frequency/phase encoding matrix and D l contains the sensitivity pattern of the lth coil: D l = diag { s coil l ( r j ) }. Regularized least-squares estimation: ˆf = argmin f L y l AD l f 2 + βr( f). l=1 Can generalize to account for noise correlation due to coil coupling. Easy to apply CG algorithm, including Toeplitz/NUFFT acceleration. For Cartesian SENSE, iterations are not needed. (Solve small system of linear equations for each voxel.) 57
59 Summary Iterative reconstruction: much potential in MRI Quadratic regularization parameter selection is tractable Computation: reduced by tools like NUFFT / Toeplitz But optimization algorithm design remains important (cf. Shepp and Vardi, 1982, PET) Some current challenges Sensitivity pattern mapping for SENSE Through-voxel field inhomogeneity gradients Motion / dynamics / partial k-space data Establishing diagnostic efficacy with clinical data... Image reconstruction toolbox: fessler 58
Tradeoffs and complexities in model-based MR image reconstruction
Tradeoffs and complexities in model-based MR image reconstruction Jeffrey A. Fessler EECS Dept., The University of Michigan fessler@umich.edu December 6, 2007 Abstract Model-based methods for image reconstruction
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 informationSteen Moeller Center for Magnetic Resonance research University of Minnesota
Steen Moeller Center for Magnetic Resonance research University of Minnesota moeller@cmrr.umn.edu Lot of material is from a talk by Douglas C. Noll Department of Biomedical Engineering Functional MRI Laboratory
More informationFast, Iterative Image Reconstruction for MRI in the Presence of Field Inhomogeneities
178 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 22, NO. 2, FEBRUARY 2003 Fast, Iterative Image Reconstruction for MRI in the Presence of Field Inhomogeneities Bradley P. Sutton*, Student Member, IEEE, Douglas
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 informationFOV. ] are the gradient waveforms. The reconstruction of this signal proceeds by an inverse Fourier Transform as:. [2] ( ) ( )
Gridding Procedures for Non-Cartesian K-space Trajectories Douglas C. Noll and Bradley P. Sutton Dept. of Biomedical Engineering, University of Michigan, Ann Arbor, MI, USA 1. Introduction The data collected
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 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 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 informationNUFFT for Medical and Subsurface Image Reconstruction
NUFFT for Medical and Subsurface Image Reconstruction Qing H. Liu Department of Electrical and Computer Engineering Duke University Duke Frontiers 2006 May 16, 2006 Acknowledgment Jiayu Song main contributor
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 informationA Novel Iterative Thresholding Algorithm for Compressed Sensing Reconstruction of Quantitative MRI Parameters from Insufficient Data
A Novel Iterative Thresholding Algorithm for Compressed Sensing Reconstruction of Quantitative MRI Parameters from Insufficient Data Alexey Samsonov, Julia Velikina Departments of Radiology and Medical
More informationNuts & Bolts of Advanced Imaging. Image Reconstruction Parallel Imaging
Nuts & Bolts of Advanced Imaging Image Reconstruction Parallel Imaging Michael S. Hansen, PhD Magnetic Resonance Technology Program National Institutes of Health, NHLBI Declaration of Financial Interests
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 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 informationSPM8 for Basic and Clinical Investigators. 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 informationG Practical Magnetic Resonance Imaging II Sackler Institute of Biomedical Sciences New York University School of Medicine. Compressed Sensing
G16.4428 Practical Magnetic Resonance Imaging II Sackler Institute of Biomedical Sciences New York University School of Medicine Compressed Sensing Ricardo Otazo, PhD ricardo.otazo@nyumc.org Compressed
More informationMotion Correction in fmri by Mapping Slice-to-Volume with Concurrent Field-Inhomogeneity Correction
Motion Correction in fmri by Mapping Slice-to-Volume with Concurrent Field-Inhomogeneity Correction Desmond T.B. Yeo 1,2, Jeffery A. Fessler 2, and Boklye Kim 1 1 Department of Radiology, University of
More informationDynamic Image and Fieldmap Joint Estimation Methods for MRI Using Single-Shot Trajectories
Dynamic Image and Fieldmap Joint Estimation Methods for MRI Using Single-Shot Trajectories by Antonis Matakos A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor
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 informationIMAGE reconstruction in conventional magnetic resonance
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 24, NO. 3, MARCH 2005 325 Conjugate Phase MRI Reconstruction With Spatially Variant Sample Density Correction Douglas C. Noll*, Member, IEEE, Jeffrey A. Fessler,
More informationFast and Motion Robust Dynamic R 2 Reconstruction for Functional MRI
Fast and Motion Robust Dynamic R 2 Reconstruction for Functional MRI by Valur Thor Olafsson A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy (Electrical
More informationMotion Robust Magnetic Susceptibility and Field Inhomogeneity Estimation Using Regularized Image Restoration Techniques for fmri
Motion Robust Magnetic Susceptibility and Field Inhomogeneity Estimation Using Regularized Image Restoration Techniques for fmri Desmond Tec Beng Yeo 1,, Jeffrey A. Fessler 1,, and Bolye Kim 1 1 Department
More informationRegularized Estimation of Main and RF Field Inhomogeneity and Longitudinal Relaxation Rate in Magnetic Resonance Imaging
Regularized Estimation of Main and RF Field Inhomogeneity and Longitudinal Relaxation Rate in Magnetic Resonance Imaging by Amanda K. Funai A dissertation submitted in partial fulfillment of the requirements
More information, Member, IEEE, Douglas C. Noll, Senior Member, IEEE, and Jeffrey A. Fessler, Fellow, IEEE
604 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 37, NO., FEBRUARY 018 Fast Spatial Resolution Analysis of Quadratic Penalized Least-Squares Image Reconstruction With Separate Real and Imaginary Roughness
More informationTITLE: Regularized Reconstruction of Dynamic Contrast-Enhanced MR Images for Evaluation of Breast Lesions
AD Award Number: W81XWH-08-1-0273 TITLE: Regularized Reconstruction of Dynamic Contrast-Enhanced MR Images for Evaluation of Breast Lesions PRINCIPAL INVESTIGATOR: Kimberly A. Khalsa CONTRACTING ORGANIZATION:
More informationBias-Variance Tradeos Analysis Using Uniform CR Bound. Mohammad Usman, Alfred O. Hero, Jerey A. Fessler and W. L. Rogers. University of Michigan
Bias-Variance Tradeos Analysis Using Uniform CR Bound Mohammad Usman, Alfred O. Hero, Jerey A. Fessler and W. L. Rogers University of Michigan ABSTRACT We quantify fundamental bias-variance tradeos for
More informationEstimating 3D Respiratory Motion from Orbiting Views
Estimating 3D Respiratory Motion from Orbiting Views Rongping Zeng, Jeffrey A. Fessler, James M. Balter The University of Michigan Oct. 2005 Funding provided by NIH Grant P01 CA59827 Motivation Free-breathing
More informationAdvanced Imaging Trajectories
Advanced Imaging Trajectories Cartesian EPI Spiral Radial Projection 1 Radial and Projection Imaging Sample spokes Radial out : from k=0 to kmax Projection: from -kmax to kmax Trajectory design considerations
More information6 credits. BMSC-GA Practical Magnetic Resonance Imaging II
BMSC-GA 4428 - Practical Magnetic Resonance Imaging II 6 credits Course director: Ricardo Otazo, PhD Course description: This course is a practical introduction to image reconstruction, image analysis
More informationAccelerated MRI Techniques: Basics of Parallel Imaging and Compressed Sensing
Accelerated MRI Techniques: Basics of Parallel Imaging and Compressed Sensing Peng Hu, Ph.D. Associate Professor Department of Radiological Sciences PengHu@mednet.ucla.edu 310-267-6838 MRI... MRI has low
More informationFast Imaging Trajectories: Non-Cartesian Sampling (1)
Fast Imaging Trajectories: Non-Cartesian Sampling (1) M229 Advanced Topics in MRI Holden H. Wu, Ph.D. 2018.05.03 Department of Radiological Sciences David Geffen School of Medicine at UCLA Class Business
More informationImaging Notes, Part IV
BME 483 MRI Notes 34 page 1 Imaging Notes, Part IV Slice Selective Excitation The most common approach for dealing with the 3 rd (z) dimension is to use slice selective excitation. This is done by applying
More informationAN ALGORITHM FOR BLIND RESTORATION OF BLURRED AND NOISY IMAGES
AN ALGORITHM FOR BLIND RESTORATION OF BLURRED AND NOISY IMAGES Nader Moayeri and Konstantinos Konstantinides Hewlett-Packard Laboratories 1501 Page Mill Road Palo Alto, CA 94304-1120 moayeri,konstant@hpl.hp.com
More informationComputer Vision I. Announcements. Fourier Tansform. Efficient Implementation. Edge and Corner Detection. CSE252A Lecture 13.
Announcements Edge and Corner Detection HW3 assigned CSE252A Lecture 13 Efficient Implementation Both, the Box filter and the Gaussian filter are separable: First convolve each row of input image I with
More informationEfficient MR Image Reconstruction for Compressed MR Imaging
Efficient MR Image Reconstruction for Compressed MR Imaging Junzhou Huang, Shaoting Zhang, and Dimitris Metaxas Division of Computer and Information Sciences, Rutgers University, NJ, USA 08854 Abstract.
More informationSome Blind Deconvolution Techniques in Image Processing
Some Blind Deconvolution Techniques in Image Processing Tony Chan Math Dept., UCLA Joint work with Frederick Park and Andy M. Yip IPAM Workshop on Mathematical Challenges in Astronomical Imaging July 26-30,
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 informationField Maps. 1 Field Map Acquisition. John Pauly. October 5, 2005
Field Maps John Pauly October 5, 25 The acquisition and reconstruction of frequency, or field, maps is important for both the acquisition of MRI data, and for its reconstruction. Many of the imaging methods
More informationMRI Physics II: Gradients, Imaging
MRI Physics II: Gradients, Imaging Douglas C., Ph.D. Dept. of Biomedical Engineering University of Michigan, Ann Arbor Magnetic Fields in MRI B 0 The main magnetic field. Always on (0.5-7 T) Magnetizes
More informationComputational Aspects of MRI
David Atkinson Philip Batchelor David Larkman Programme 09:30 11:00 Fourier, sampling, gridding, interpolation. Matrices and Linear Algebra 11:30 13:00 MRI Lunch (not provided) 14:00 15:30 SVD, eigenvalues.
More informationMarquette University Iain P. Bruce Marquette University, M. Muge Karaman Marquette University
Marquette University e-publications@marquette Mathematics, Statistics and omputer Science Faculty Research and Publications Mathematics, Statistics and omputer Science, Department of 10-1-2012 The SENSE-Isomorphism
More informationEvaluations of k-space Trajectories for Fast MR Imaging for project of the course EE591, Fall 2004
Evaluations of k-space Trajectories for Fast MR Imaging for project of the course EE591, Fall 24 1 Alec Chi-Wah Wong Department of Electrical Engineering University of Southern California 374 McClintock
More informationFmri Spatial Processing
Educational Course: Fmri Spatial Processing Ray Razlighi Jun. 8, 2014 Spatial Processing Spatial Re-alignment Geometric distortion correction Spatial Normalization Smoothing Why, When, How, Which Why is
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 informationProjection Space Maximum A Posterior Method for Low Photon Counts PET Image Reconstruction
Proection Space Maximum A Posterior Method for Low Photon Counts PET Image Reconstruction Liu Zhen Computer Department / Zhe Jiang Wanli University / Ningbo ABSTRACT In this paper, we proposed a new MAP
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 informationGPR Migration Imaging Algorithm Based on NUFFT
PIERS ONLINE, VOL. 6, NO. 1, 010 16 GPR Migration Imaging Algorithm Based on NUFFT Hao Chen, Renbiao Wu, Jiaxue Liu, and Zhiyong Han Tianjin Key Laboratory for Advanced Signal Processing, Civil Aviation
More informationK-Space Trajectories and Spiral Scan
K-Space and Spiral Scan Presented by: Novena Rangwala nrangw2@uic.edu 1 Outline K-space Gridding Reconstruction Features of Spiral Sampling Pulse Sequences Mathematical Basis of Spiral Scanning Variations
More informationImage Registration: Warping Without Folding
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
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 informationNonuniform Fast Fourier Transforms Using Min-Max Interpolation
560 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 51, NO. 2, FEBRUARY 2003 Nonuniform Fast Fourier Transforms Using Min-Max Interpolation Jeffrey A. Fessler, Senior Member, IEEE, and Bradley P. Sutton,
More informationHST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008
MIT OpenCourseWare http://ocw.mit.edu HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationRefractive Index Map Reconstruction in Optical Deflectometry Using Total-Variation Regularization
Refractive Index Map Reconstruction in Optical Deflectometry Using Total-Variation Regularization L. Jacques 1, A. González 2, E. Foumouo 2 and P. Antoine 2 1 Institute of Information and Communication
More informationChallenges in Improved Sensitivity of Quantification of PET Data for Alzheimer s Disease Studies
Challenges in Improved Sensitivity of Quantification of PET Data for Alzheimer s Disease Studies Rosemary Renaut, Hongbin Guo, Kewei Chen and Wolfgang Stefan Supported by Arizona Alzheimer s Research Center
More informationThis Time. fmri Data analysis
This Time Reslice example Spatial Normalization Noise in fmri Methods for estimating and correcting for physiologic noise SPM Example Spatial Normalization: Remind ourselves what a typical functional image
More informationarxiv: v3 [math.na] 14 Jun 2018
LAP: A LINEARIZE AND PROJECT METHOD FOR SOLVING INVERSE PROBLEMS WITH COUPLED VARIABLES JAMES L. HERRING, JAMES G. NAGY, AND LARS RUTHOTTO arxiv:1705.09992v3 [math.na] 14 Jun 2018 Abstract. Many inverse
More informationSeparate Magnitude and Phase Regularization via Compressed Sensing
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 31, NO. 9, SEPTEMBER 2012 1713 Separate Magnitude and Phase Regularization via Compressed Sensing Feng Zhao*, Douglas C. Noll, Senior Member, IEEE, Jon-Fredrik
More informationAugmented Lagrangian with Variable Splitting for Faster Non-Cartesian L 1 -SPIRiT MR Image Reconstruction: Supplementary Material
IEEE TRANSACTIONS ON MEDICAL IMAGING: SUPPLEMENTARY MATERIAL 1 Augmented Lagrangian with Variable Splitting for Faster Non-Cartesian L 1 -SPIRiT MR Image Reconstruction: Supplementary Material Daniel S.
More informationEE290T: Advanced Reconstruction Methods for Magnetic Resonance Imaging. Martin Uecker
EE290T: Advanced Reconstruction Methods for Magnetic Resonance Imaging Martin Uecker Tentative Syllabus 01: Jan 27 Introduction 02: Feb 03 Parallel Imaging as Inverse Problem 03: Feb 10 Iterative Reconstruction
More informationConstrained Reconstruction of Sparse Cardiac MR DTI Data
Constrained Reconstruction of Sparse Cardiac MR DTI Data Ganesh Adluru 1,3, Edward Hsu, and Edward V.R. DiBella,3 1 Electrical and Computer Engineering department, 50 S. Central Campus Dr., MEB, University
More informationScanner Parameter Estimation Using Bilevel Scans of Star Charts
ICDAR, Seattle WA September Scanner Parameter Estimation Using Bilevel Scans of Star Charts Elisa H. Barney Smith Electrical and Computer Engineering Department Boise State University, Boise, Idaho 8375
More informationMachine Learning / Jan 27, 2010
Revisiting Logistic Regression & Naïve Bayes Aarti Singh Machine Learning 10-701/15-781 Jan 27, 2010 Generative and Discriminative Classifiers Training classifiers involves learning a mapping f: X -> Y,
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 informationXI Signal-to-Noise (SNR)
XI Signal-to-Noise (SNR) Lecture notes by Assaf Tal n(t) t. Noise. Characterizing Noise Noise is a random signal that gets added to all of our measurements. In D it looks like this: while in D
More informationHigher Degree Total Variation for 3-D Image Recovery
Higher Degree Total Variation for 3-D Image Recovery Greg Ongie*, Yue Hu, Mathews Jacob Computational Biomedical Imaging Group (CBIG) University of Iowa ISBI 2014 Beijing, China Motivation: Compressed
More informationQuantitative Statistical Methods for Image Quality Assessment
Quantitative Statistical Methods for Image Quality Assessment The Harvard community has made this article openly available. Please share how this access benefits you. Your story matters Citation Dutta,
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 informationAn Introduction to Image Reconstruction, Processing, and their Effects in FMRI
An Introduction to Image Reconstruction, Processing, and their Effects in FMRI Daniel B. Rowe Program in Computational Sciences Department of Mathematics, Statistics, and Computer Science Marquette University
More informationPET Image Reconstruction using Anatomical Information through Mutual Information Based Priors
2005 IEEE Nuclear Science Symposium Conference Record M11-354 PET Image Reconstruction using Anatomical Information through Mutual Information Based Priors Sangeetha Somayajula, Evren Asma, and Richard
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 informationSlide 1. Technical Aspects of Quality Control in Magnetic Resonance Imaging. Slide 2. Annual Compliance Testing. of MRI Systems.
Slide 1 Technical Aspects of Quality Control in Magnetic Resonance Imaging Slide 2 Compliance Testing of MRI Systems, Ph.D. Department of Radiology Henry Ford Hospital, Detroit, MI Slide 3 Compliance Testing
More informationAccelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method
Accelerating the Hessian-free Gauss-Newton Full-waveform Inversion via Preconditioned Conjugate Gradient Method Wenyong Pan 1, Kris Innanen 1 and Wenyuan Liao 2 1. CREWES Project, Department of Geoscience,
More informationStatistical Analysis of Image Reconstructed Fully-Sampled and Sub-Sampled fmri Data
Statistical Analysis of Image Reconstructed Fully-Sampled and Sub-Sampled fmri Data Daniel B. Rowe Program in Computational Sciences Department of Mathematics, Statistics, and Computer Science Marquette
More informationIterative reconstruction method for threedimensional non-cartesian parallel MRI
University of Iowa Iowa Research Online Theses and Dissertations Spring 2011 Iterative reconstruction method for threedimensional non-cartesian parallel MRI Xuguang Jiang University of Iowa Copyright 2011
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 informationSTATISTICAL image reconstruction methods for X-ray
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 28, NO. 5, MAY 2009 645 Quadratic Regularization Design for 2-D CT Hugo R. Shi*, Student Member, IEEE, and Jeffrey A. Fessler, Fellow, IEEE Abstract Statistical
More informationModule 4. K-Space Symmetry. Review. K-Space Review. K-Space Symmetry. Partial or Fractional Echo. Half or Partial Fourier HASTE
MRES 7005 - Fast Imaging Techniques Module 4 K-Space Symmetry Review K-Space Review K-Space Symmetry Partial or Fractional Echo Half or Partial Fourier HASTE Conditions for successful reconstruction Interpolation
More informationGuided Image Super-Resolution: A New Technique for Photogeometric Super-Resolution in Hybrid 3-D Range Imaging
Guided Image Super-Resolution: A New Technique for Photogeometric Super-Resolution in Hybrid 3-D Range Imaging Florin C. Ghesu 1, Thomas Köhler 1,2, Sven Haase 1, Joachim Hornegger 1,2 04.09.2014 1 Pattern
More informationThe Fast Multipole Method (FMM)
The Fast Multipole Method (FMM) Motivation for FMM Computational Physics Problems involving mutual interactions of N particles Gravitational or Electrostatic forces Collective (but weak) long-range forces
More informationDigital Image Processing Laboratory: MAP Image Restoration
Purdue University: Digital Image Processing Laboratories 1 Digital Image Processing Laboratory: MAP Image Restoration October, 015 1 Introduction This laboratory explores the use of maximum a posteriori
More informationOptimal Sampling Geometries for TV-Norm Reconstruction of fmri Data
Optimal Sampling Geometries for TV-Norm Reconstruction of fmri Data Oliver M. Jeromin, Student Member, IEEE, Vince D. Calhoun, Senior Member, IEEE, and Marios S. Pattichis, Senior Member, IEEE Abstract
More informationAdvanced Image Reconstruction Methods for Photoacoustic Tomography
Advanced Image Reconstruction Methods for Photoacoustic Tomography Mark A. Anastasio, Kun Wang, and Robert Schoonover Department of Biomedical Engineering Washington University in St. Louis 1 Outline Photoacoustic/thermoacoustic
More informationA Novel Image Super-resolution Reconstruction Algorithm based on Modified Sparse Representation
, pp.162-167 http://dx.doi.org/10.14257/astl.2016.138.33 A Novel Image Super-resolution Reconstruction Algorithm based on Modified Sparse Representation Liqiang Hu, Chaofeng He Shijiazhuang Tiedao University,
More informationSpatial Regularization of Functional Connectivity Using High-Dimensional Markov Random Fields
Spatial Regularization of Functional Connectivity Using High-Dimensional Markov Random Fields Wei Liu 1, Peihong Zhu 1, Jeffrey S. Anderson 2, Deborah Yurgelun-Todd 3, and P. Thomas Fletcher 1 1 Scientific
More informationPhase Difference Reconstruction. Outline
Advanced MRI Phase Difference Reconstruction Faik Can MERAL Outline Introduction Quantitative Description Arctangent operation ATAN2 Phased-Array Multiple Coil Data Correction of Predictable Phase Errors
More informationCollaborative Sparsity and Compressive MRI
Modeling and Computation Seminar February 14, 2013 Table of Contents 1 T2 Estimation 2 Undersampling in MRI 3 Compressed Sensing 4 Model-Based Approach 5 From L1 to L0 6 Spatially Adaptive Sparsity MRI
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 informationA Weighted Least Squares PET Image Reconstruction Method Using Iterative Coordinate Descent Algorithms
A Weighted Least Squares PET Image Reconstruction Method Using Iterative Coordinate Descent Algorithms Hongqing Zhu, Huazhong Shu, Jian Zhou and Limin Luo Department of Biological Science and Medical Engineering,
More informationIEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 23, NO. 4, APRIL
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 23, NO. 4, APRIL 2004 401 Iterative Tomographic Image Reconstruction Using Fourier-Based Forward and Back-Projectors Samuel Matej*, Senior Member, IEEE, Jeffrey
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 informationA Nonparametric Bayesian Approach to Detecting Spatial Activation Patterns in fmri Data
A Nonparametric Bayesian Approach to Detecting Spatial Activation Patterns in fmri Data Seyoung Kim, Padhraic Smyth, and Hal Stern Bren School of Information and Computer Sciences University of California,
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 informationHST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008
MIT OpenCourseWare http://ocw.mit.edu HST.583 Functional Magnetic Resonance Imaging: Data Acquisition and Analysis Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationLaser sensors. Transmitter. Receiver. Basilio Bona ROBOTICA 03CFIOR
Mobile & Service Robotics Sensors for Robotics 3 Laser sensors Rays are transmitted and received coaxially The target is illuminated by collimated rays The receiver measures the time of flight (back and
More informationPassive Differential Matched-field Depth Estimation of Moving Acoustic Sources
Lincoln Laboratory ASAP-2001 Workshop Passive Differential Matched-field Depth Estimation of Moving Acoustic Sources Shawn Kraut and Jeffrey Krolik Duke University Department of Electrical and Computer
More informationCompressed Sensing for Rapid MR Imaging
Compressed Sensing for Rapid Imaging Michael Lustig1, Juan Santos1, David Donoho2 and John Pauly1 1 Electrical Engineering Department, Stanford University 2 Statistics Department, Stanford University rapid
More informationParametric Texture Model based on Joint Statistics
Parametric Texture Model based on Joint Statistics Gowtham Bellala, Kumar Sricharan, Jayanth Srinivasa Department of Electrical Engineering, University of Michigan, Ann Arbor 1. INTRODUCTION Texture images
More information3-D Tomographic Reconstruction
Mitglied der Helmholtz-Gemeinschaft 3-D Tomographic Reconstruction of Atmospheric Trace Gas Concentrations for Infrared Limb-Imagers February 21, 2011, Nanocenter, USC Jörn Ungermann Tomographic trace
More informationPicture quality requirements and NUT proposals for JPEG AIC
Mar. 2006, Cupertino Picture quality requirements and NUT proposals for JPEG AIC Jae-Jeong Hwang, Young Huh, Dai-Gyoung Kim Kunsan National Univ., KERI, Hanyang Univ. hwang@kunsan.ac.kr Contents 1. Picture
More information