Improved 3D image plane parallel magnetic resonance imaging (pmri) method
|
|
- Marvin Nichols
- 5 years ago
- Views:
Transcription
1 B. Wu, R. P. Millane, R. Watts, P. J. Bones, Improved 3D Image Plane Parallel Magnetic Resonance Imaging (pmri) Method, Proceedings of Image and Vision Computing New Zealand 2007, pp , Hamilton, New Zealand, December Improved 3D image plane parallel magnetic resonance imaging (pmri) method B. Wu 1, R. P. Millane 1, R. Watts 2 and P. J. Bones 1 1 Computational Imaging Group, Dept. Electrical & Computer Engineering, University of Canterbury, Christchurch, New Zealand. 2 Dept. Physics & Astronomy, University of Canterbury, Christchurch, New Zealand. bing.wu@elec.canterbury.ac.nz Abstract A new image plane parallel magnetic resonance imaging (pmri) method, Generalized Unaliasing Incorporating object Support constraint and sensitivity Encoding (GUISE), was recently introduced. It allows images to be recovered from arbitrary cartesian k-space (spatial frequency space) sampling patterns, restricted for efficiency to be those consisting of repeated blocks. Sensitivity encoding (SENSE) is seen as a special case of GUISE where a minimal repeating block size is used. An efficient sampling pattern design method for GUISE is proposed in order to reduce the overall noise level in the reconstructed image. Experimental results show that GUISE with this sampling pattern method provides better image reconstruction than SENSE. The improvement is attributed to better exploitation of the differences in spatial sensitivity across the set of receiver coils and the object support constraint. Keywords: 3D MRI, pmri, object support 1 INTRODUCTION 3D volumetric magnetic resonance imaging (MRI) allows images to be viewed in different planes independent of the order used for the data acquisition, and is widely used for breast imaging [1], brain imaging [2] and angiogram imaging [3]. However, clinical usage is limited by the long scan time. Recent development of parallel MRI (pmri) techniques have significantly reduced the necessary data acquisition time of 3D MRI examinations [3]. In general, pmri methods employ multiple receiver coils during data acquisition, and the different data sets gained from the coils, each of which has different spatial sensitivity variation, are exploited in post-processing to compensate for the skipped phase encoding (PE) steps. Among various pmri methods proposed to date, SENSitivity Encoding (SENSE) [4] has received much attention and has been implemented on commercial MR scanners because of its simplistic nature and high signal-tonoise ratio (SNR) performance. We have recently introduced a new pmri method, Generalized Unaliasing Incorporating object Support constraint and sensitivity Encoding (GUISE) [5], which generalizes the image plane pmri methods. In both SENSE and GUISE, a series of matrix inversions are to be performed to separate the folded voxels using knowledge of the coil sensitivity weighting for direct image recovery. The key difference is that SENSE results in fixed aliasing pattern whereas the aliasing pattern in GUISE is determined by the sampling pattern adopted. We postulate that this freedom in the aliasing pattern design allows more efficient exploitation of coil sensitivity and the object support region, which leads to improved SNR in the reconstructed image. In this paper, we present an overview of GUISE and propose a sampling pattern design method. We then compare the performance of GUISE with SENSE based on experimental data. 2 THEORY The term k-space is widely used in the MRI community to refer to the spatial frequency space in which data is acquired. In 3D MRI, the measured k-space data represent a 3D Fourier encoding of the magnetization signal. However, inverse discrete Fourier transformation along the fully sampled frequency encoding (FE) direction allows the 3D k-space to be decomposed to its 2D multi-slice equivalent. Each resulting image slice is thereby defined in terms of the two PE directions and reconstructions can be performed on individual slices to restore the entire 3D image. For ease of illustration and without loss of generality, we there- 311
2 fore present our work based on a single 2D plane and considerations for 3D space are given where needed. 2.1 Formulation of GUISE Consider the case of M receiver coils. The k-space data from the mth coil represents the 2D Fourier transform of the coil sensitivity weighted object magnetization. In the case of a full set of samples at Nyquist spacing, the system can be represented conveniently as F m = W diag(c m ) f, (1) where f, c m, and F m represent the stacked column vectors of the 2D object magnetization, the mth coil sensitivity map, and the acquired samples from the mth receiver coil, respectively. Diag(x) refers to a diagonal matrix with vector x on the main diagonal and zero elsewhere. W is the stacked Fourier matrix, left multiplication by which performs the 2D DFT. Under-sampling in k-space is equivalent to deleting the entries corresponding to skipped sample locations, which can be expressed as diag(h) F m = diag(h) W diag(c m ) f, (2) where H is a 2D k-space binary sampling mask stacked into a column vector. Inverse DFT of the under-sampled data set gives an image vector, f m, that contains aliasing artifacts: fm = W 1 diag(h) F m = W 1 diag(h) W diag(c m ) f = A diag(c m ) f, (3) where A is a matrix that maps the coil sensitivity weighted object magnetization to the aliased image, which we call a Fourier encoding matrix. The aliasing artifacts cannot be removed by inverting A directly because A is rank deficient. For multiple receiver coils, we can create a new hybrid encoding system, B, that maps the object magnetization to the aliased versions in different coils. Mathematically, this is achieved by stacking the aliased images from different coils and the corresponding encoding functions, i.e., B f = f, B = f = A diag(c 1 ) A diag(c 2 ) A diag(c M ) f1 f2 fm,. (4) B is invertible as long as the number of independent rows is greater than or equal to the number of independent columns. Due to distinct coil spatial sensitivity weightings, the simultaneously acquired samples in different coils are independent and the row rank increases by the number of phase encoding steps for every coil used. Thus an acceleration factor up to M is possible. Knowledge of the object support allows us to eliminate the regions outside the object from reconstruction, as they are known to generate zero signal. Only those entries in f that lie within the support region and the corresponding columns in B are retained. In the interests of brevity, we henceforth denote the reduced versions of f and B as f[ρ] and B[ρ], respectively, thus B[ρ] f[ρ] = f. (5) Denoting the ratio of the object support area to the total FOV area as φ, where φ 1, the maximum possible acceleration factor for recovering the 2D image is improved by the use of the support constraint to α max = M/φ. 2.2 System decomposition The system in (5) poses a problem for efficient matrix inversion as a result of its size, but this can be relieved by employing a periodic k-space sampling pattern that decomposes the overall system [6]. In a FOV of size N 1 N 2, if the k-space sampling pattern corresponding to H consists of a repeating pattern of size C 1 C 2 whose stacked column vector is G, then its point spread function has a sparse structure with nonzero values separated by L 1 = N 1 /C 1 and L 2 = N 2 /C 2 voxels in each respective direction. As a result, aliasing takes place in a well controlled manner, as depicted in Figure 1: only voxels that are separated by integer multiples of L 1 and L 2 in each direction respectively are aliased together, forming what we denote as a subsequence. Voxels within each such aliased subsequence are independent from voxels in all other sub-sequences [6]. Thus the overall system can be decomposed into L = L 1 L 2 sub-sequences in the above manner, and a sub-system matrix can be set up for the lth subsequence: B l [ρ l ] f l [ρ l ] = f l, fl = B l [ρ l ] = fl,1 fl,2 fl,m Â[ρ l ] diag(c l,1 [ρ l ]) Â[ρ l ] diag(c l,2 [ρ l ]) Â[ρ l ] diag(c l,m [ρ l ]),, (6) 312
3 2.3.1 Choosing the block size Figure 1: A periodic k-space sampling pattern leads to decomposition of the image plane system. On the left, the central region of the k-space sampling pattern, consisting of a 4 4 repeating block, is shown (bright and dark parts represent acquired and skipped samples, respectively). Aliasing in the resulting image plane is limited to occur within sub-sequences of voxels, one of which is depicted on the right as a regular set of white dots. For each such sub-sequence only those voxels which fall within the support region (outlined in white) need to be recovered. The the final 2D image is reconstructed by assembling all the recovered sub-sequences. Each image plane has its own support region. where  l [ρ l ] = (W 1 diag(g)w)[ρ l ], l = 0, 1, 2,..., L 1 L 2 1, and l is used to denote quantities related to the lth sub-sequence only. For the lth sub-sequence we describe the ratio of the number of voxels within the object support to the total number of sub-sequence voxels as the subsupport ratio, denoted φ l. The entire image can be recovered by inverting a series of sub-matrices followed by appropriate reordering. 2.3 Noise propagation We have shown in [5] that the overall noise amplification level in the reconstruction can be written as: P 1 = trace(b[ρ]b[ρ] H ) +, (7) σ 2 i=1 i where σ denotes the singular value of B and P is the total number of voxels to be recovered. We refer this quantity as inverse trace metric (IT metric) and use it as a measurement tool to judge the noise immunity of different sampling strategies. In the case of using periodic sampling pattern, this quantity can be written as the sum of the IT metrics of individual sub-systems: trace(b[ρ]b[ρ] H ) + = L trace(b l [ρ l ]B l [ρ l ] H ) +. l=1 (8) Recall that the block has dimensions C 1 C 2. To achieve an acceleration factor of α, appropriate block sizes are restricted to those which make C 1 C 2 an integer multiple of α. This still leaves considerable freedom in choosing the block size. We now consider the factors in choosing the block parameters, first in terms of shape, then in terms of overall size. As can be seen in Figure 2, the values of C 1 and C 2 are directly reflected in the number of voxels being aliased in the corresponding direction. There is an advantage therefore in choosing the larger of C 1 or C 2 to correspond to the direction in which there is a greater variation in coil sensitivity. There is also an advantage to be gained by choosing the larger of C 1 or C 2 to correspond to the direction for which the out-of-support region is larger (i.e., for which the distance from the edge of the frame to the support is greater.) Assuming that the constraints above have been satisfied, we now consider the effects of choosing an overall large block size. As illustrated in Figure 2, larger k-space block sizes lead to a reduced number of sub-sequences in the spatial domain. Likewise the spatial extent of each sub-sequence is increased and sub-sequence elements are more closely packed. Based on the knowledge that coil sensitivities are smoothly varying spatial functions and object support areas are generally contiguous, we postulate that the use of large block sizes has the following influence: 1. The range of sensitivity variation within each sub-sequence is enhanced and finer. This can be exploited by using appropriate sampling patterns for improved reconstruction, however it also means that bad sampling patterns would lead to worse results compared to the case of using smaller block sizes. 2. Voxels within the object support are more evenly distributed among the resulting subsequences. Applying the support constraint therefore increases the degree of determinacy of the sub-systems, and the determinacy for each individual sub-system is determined by the sub-support ratio, φ l. Thus the distribution of φ l has a large impact on the noise profile and an uneven distribution of φ l values, as illustrated in Figure 2(a), leads to a nonuniform noise profile in the reconstruction. Furthermore, because the noise amplification level (IT metric) of a particular sub-system grows rapidly with increasing φ l, a narrow distribution of the φ l values is also a favorable condition for a lower overall noise level. 313
4 Figure 2: In an image plane of size , sub-systems are shown which correspond to k-space block sizes of (a) 2 4, (b) The histogram of all the φ l values indicates that using larger block size leads to an even distribution of the voxels within the support region (the central shaded circular region) and also a lower peak φ l. Thus expanding the k-space block size promotes the maximum acceleration factor achievable Sample selection The actual repeating pattern determines how the voxels are aliased and is a critical factor in the conditioning. Although periodic sampling reduces the problem size very significantly, an exhaustive search for the optimal pattern is still computationally intractable. Rather we propose to use the sequential forward selection (SFS) algorithm originally proposed by Gao et al [7], which sequentially selects sample positions based on minimizing a specified cost function in each step. However, the IT metric cost function involves a computationally expensive matrix inverse. We have shown previously [5] that the sum of singular values raised to the fourth power, i.e., P σi 4 = trace(b[ρ]b[ρ] H ) 2, (9) i=1 which we refer to as the squared trace metric (ST metric), exhibits similar general behavior to the IT metric under the process of minimization. Although the ST metric appears computationally expensive, very efficient calculation can be achieved in a sequential algorithm. 3 METHODS To compare the performance of GUISE and SENSE, a 3D brain scan of a healthy adult volunteer was performed on a GE 1.5T scanner using an 8-channel head coil. Informed consent was obtained prior to the study. The two PE directions were set to leftright (LR) and anterior-posterior (AP), which correspond to the two directions with significant coil sensitivity variation, and FE is set to the superiorinferior (SI) direction. T 1 weighted 3D Spoiled Gradient Recalled (SPGR) was used with the following sequence parameters: T R = 23ms, T E = 10ms, flip angle = 15, matrix = , resolution = 1.5mm 0.75mm 1.5mm. A full coronal plane data set was acquired to allow reconstructions with different sampling patterns to be performed in post-processing and the inverse DFT was performed along the FE direction prior to further reconstruction in individual axial slices. Image reconstructions using GUISE and SENSE methods were implemented in Matlab R (The Mathworks, Natick, MA, USA). 4 RESULTS Reconstructions of axial plane images using SENSE and GUISE at a sampling acceleration factor of 8 are shown in Figure 3. In SENSE, asymmetrical block sizes with a greater LR dimension than the AP dimension are chosen since there is more nonobject margin in the LR direction. In GUISE, a repeating block size of 8 8 is used and samples have been sequentially selected based on minimizing the overall ST metric. It is seen that in both cases GUISE method results in comparatively more uniform noise profile. The overall noise amplification levels are compared using the normalized root mean squared distance (NRMS), defined by 314
5 Figure 3: Comparing the reconstructed axial plane image slices using SENSE and GUISE at an overall acceleration factor of 8. Reconstructed axial slice of the brain region: (a) reconstruction from full data acquisition (b) using GUISE (c) using SENSE. Reconstructed axial slice through the eyeball: (d) reconstruction from full data acquisition (e) using GUISE (f) using SENSE. The NRMS for each accelerated reconstruction is shown in the bottom corner. NRMS = P i=1 ( Ireference i I reconstructed i ) 2 P, i=1 Ireference i 2 where I i is the ith voxel and P denotes the number of voxels recovered within the support region. The reconstructions gained from full sample acquisition were used as reference images (Figure 3(a) and Figure 3(d)). The NRMS values in the lower right of Figure 3(b), (c), (e) and (f) indicate that GUISE also out-performed SENSE in providing a lower overall noise level. Reconstructed brain slice images had better SNR compared to reconstructed eyeball slice images as the former features a smaller object support region and therefore lower φ value than the latter. 5 DISCUSSION GUISE treats the image plane aliasing due to undersampling in k-space as convolving the image with the point spread function of the sampling pattern. In general for pmri methods, the knowledge of image plane coil sensitivity weighting is used to remove the aliasing artifacts. GUISE explicitly incorporates the knowledge of the object s spatial support constraint in the image recovery process, which improves the determinacy of the linear system. By introducing periodicity in the k-space sampling pattern, the overall image plane system can be decomposed into independent sub-systems so that image recovery efficiency is greatly improved. The size of the repeating block in k-space sampling pattern determines the spatial extent of the resulting sub-systems and consequently their coil sensitivity weighting and sub-support ratios. Sufficient exploitation of the coil sensitivity distinctness and the finite object spatial extent is restricted by the small repeating block size used in SENSE. However, the repeating sampling block needs to welldesigned to keep the keep the sub-systems wellconditioned, and it is a difficult task to choose a suitable pattern for the block. The SFS method based on an approximated cost function achieves a satisfactory balance between efficiency and optimality. Since decreasing marginal gain is expected as the block size is further expanded in return of dramatically increased sample selection complexity, a moderate block size should be used in practice. 315
6 6 CONCLUSION GUISE is a new parallel MRI method with high flexibility in sampling pattern design. GUISE with appropriate sampling patterns can achieve better image reconstruction than the currently favoured method in clinical practice, SENSE, at high acceleration factors. References [1] K. Akazawa, Y. Tamaki, T. Taguchi, Y. Tanji, Y. Miyoshi, and S. Kim, Preoperative evaluation of residual tumor extent by threedimensional magnetic resonance imaging in breast cancer patients treated with neoadjuvant chemotherapy, The Breast Journal, vol. 12(2), pp , [2] C. Helmchen, H. Rambold, A. Sprenger, C. Erdmann, and F. Binkofski, Cerebellar activation in opsoclonus: an fmri study, Neurology, vol. 61(3), pp , [3] R. Muthupillai, G. Vick, S. Flamm, and T. Chung, Time-resolved contrast-enhanced magnetic resonance angiography in pediatric patients using sensitivity encoding, Journal of Magnetic Resonance Imaging, vol. 17, pp , [4] M. Weiger, K. Pruessmann, and P. Boesiger, 2D SENSE for faster 3D MRI, MAGMA, vol. 14, pp , [5] P. Bones, B. Wu, R. Millane, and R. Watts, Support constraint in 3-D magnetic resonance imaging, in Signal Recovery and Synthesis (Postdeadline Papers), The Optical Society of America, Washington DC, United States of America, April 2007, p. SMD6. [6] P. Bones, N. Alwesh, T. Connolly, and N. Blakeley, Recovery of limited-extent images aliased because of spectral undersampling, JOSA A, vol. 18, pp , [7] Y. Gao and S. Reeves, Optimal k-space sampling in MRSI for images with a limited region of support, IEEE Trans. Med. Imaging, vol. 12, pp ,
Parallel Imaging. Marcin.
Parallel Imaging Marcin m.jankiewicz@gmail.com Parallel Imaging initial thoughts Over the last 15 years, great progress in the development of pmri methods has taken place, thereby producing a multitude
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 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 informationCompressed Sensing Reconstructions for Dynamic Contrast Enhanced MRI
1 Compressed Sensing Reconstructions for Dynamic Contrast Enhanced MRI Kevin T. Looby klooby@stanford.edu ABSTRACT The temporal resolution necessary for dynamic contrast enhanced (DCE) magnetic resonance
More informationZigzag Sampling for Improved Parallel Imaging
Magnetic Resonance in Medicine 60:474 478 (2008) Zigzag Sampling for Improved Parallel Imaging Felix A. Breuer, 1 * Hisamoto Moriguchi, 2 Nicole Seiberlich, 3 Martin Blaimer, 1 Peter M. Jakob, 1,3 Jeffrey
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 informationClassification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging
1 CS 9 Final Project Classification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging Feiyu Chen Department of Electrical Engineering ABSTRACT Subject motion is a significant
More informationControlled Aliasing in Parallel Imaging Results in Higher Acceleration (CAIPIRINHA)
www.siemens.com/magnetom-world Controlled Aliasing in Parallel Imaging Results in Higher Acceleration (CAIPIRINHA) Felix Breuer; Martin Blaimer; Mark Griswold; Peter Jakob Answers for life. Controlled
More informationDynamic Autocalibrated Parallel Imaging Using Temporal GRAPPA (TGRAPPA)
Magnetic Resonance in Medicine 53:981 985 (2005) Dynamic Autocalibrated Parallel Imaging Using Temporal GRAPPA (TGRAPPA) Felix A. Breuer, 1 * Peter Kellman, 2 Mark A. Griswold, 1 and Peter M. Jakob 1 Current
More informationLab Location: MRI, B2, Cardinal Carter Wing, St. Michael s Hospital, 30 Bond Street
Lab Location: MRI, B2, Cardinal Carter Wing, St. Michael s Hospital, 30 Bond Street MRI is located in the sub basement of CC wing. From Queen or Victoria, follow the baby blue arrows and ride the CC south
More informationRedundancy Encoding for Fast Dynamic MR Imaging using Structured Sparsity
Redundancy Encoding for Fast Dynamic MR Imaging using Structured Sparsity Vimal Singh and Ahmed H. Tewfik Electrical and Computer Engineering Dept., The University of Texas at Austin, USA Abstract. For
More informationTwo-Parameter Selection Techniques for Projection-based Regularization Methods: Application to Partial-Fourier pmri
Two-Parameter Selection Techniques for Projection-based Regularization Methods: Application to Partial-Fourier pmri Misha E. Kilmer 1 Scott Hoge 1 Dept. of Mathematics, Tufts University Medford, MA Dept.
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 informationWeighted-CS for reconstruction of highly under-sampled dynamic MRI sequences
Weighted- for reconstruction of highly under-sampled dynamic MRI sequences Dornoosh Zonoobi and Ashraf A. Kassim Dept. Electrical and Computer Engineering National University of Singapore, Singapore E-mail:
More informationA Study of Medical Image Analysis System
Indian Journal of Science and Technology, Vol 8(25), DOI: 10.17485/ijst/2015/v8i25/80492, October 2015 ISSN (Print) : 0974-6846 ISSN (Online) : 0974-5645 A Study of Medical Image Analysis System Kim Tae-Eun
More informationControlled Aliasing in Volumetric Parallel Imaging (2D CAIPIRINHA)
Magnetic Resonance in Medicine 55:549 556 (2006) Controlled Aliasing in Volumetric Parallel Imaging (2D CAIPIRINHA) Felix A. Breuer,* Martin Blaimer, Matthias F. Mueller, Nicole Seiberlich, Robin M. Heidemann,
More informationNew Technology Allows Multiple Image Contrasts in a Single Scan
These images were acquired with an investigational device. PD T2 T2 FLAIR T1 MAP T1 FLAIR PSIR T1 New Technology Allows Multiple Image Contrasts in a Single Scan MR exams can be time consuming. A typical
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 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 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 informationSparse MRI: The Application of Compressed Sensing for Rapid MR Imaging
Sparse MRI: The Application of Compressed Sensing for Rapid MR Imaging Michael Lustig, David Donoho and John M. Pauly Magnetic Resonance Systems Research Laboratory, Department of Electrical Engineering,
More informationCT NOISE POWER SPECTRUM FOR FILTERED BACKPROJECTION AND ITERATIVE RECONSTRUCTION
CT NOISE POWER SPECTRUM FOR FILTERED BACKPROJECTION AND ITERATIVE RECONSTRUCTION Frank Dong, PhD, DABR Diagnostic Physicist, Imaging Institute Cleveland Clinic Foundation and Associate Professor of Radiology
More informationEE369C: Assignment 4 Solutions
EE369C Fall 2017-18 Medical Image Reconstruction 1 EE369C: Assignment 4 Solutions Due Wednesday, Oct. 25 The problems this week will be concerned with 2DFT SENSE reconstruction. The data is an axial brain
More informationSampling, Ordering, Interleaving
Sampling, Ordering, Interleaving Sampling patterns and PSFs View ordering Modulation due to transients Temporal modulations Slice interleaving Sequential, Odd/even, bit-reversed Arbitrary Other considerations:
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 informationMODEL-BASED FREE-BREATHING CARDIAC MRI RECONSTRUCTION USING DEEP LEARNED & STORM PRIORS: MODL-STORM
MODEL-BASED FREE-BREATHING CARDIAC MRI RECONSTRUCTION USING DEEP LEARNED & STORM PRIORS: MODL-STORM Sampurna Biswas, Hemant K. Aggarwal, Sunrita Poddar, and Mathews Jacob Department of Electrical and Computer
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 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 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 informationSupplementary methods
Supplementary methods This section provides additional technical details on the sample, the applied imaging and analysis steps and methods. Structural imaging Trained radiographers placed all participants
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 informationarxiv: v2 [physics.med-ph] 22 Jul 2014
Multichannel Compressive Sensing MRI Using Noiselet Encoding arxiv:1407.5536v2 [physics.med-ph] 22 Jul 2014 Kamlesh Pawar 1,2,3, Gary Egan 4, and Jingxin Zhang 1,5,* 1 Department of Electrical and Computer
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 informationM R I Physics Course
M R I Physics Course Multichannel Technology & Parallel Imaging Nathan Yanasak, Ph.D. Jerry Allison Ph.D. Tom Lavin, B.S. Department of Radiology Medical College of Georgia References: 1) The Physics of
More informationSparse MRI: The Application of Compressed Sensing for Rapid MR Imaging
Magnetic Resonance in Medicine 58:1182 1195 (2007) Sparse MRI: The Application of Compressed Sensing for Rapid MR Imaging Michael Lustig, 1 David Donoho, 2 and John M. Pauly 1 The sparsity which is implicit
More informationMulti-azimuth velocity estimation
Stanford Exploration Project, Report 84, May 9, 2001, pages 1 87 Multi-azimuth velocity estimation Robert G. Clapp and Biondo Biondi 1 ABSTRACT It is well known that the inverse problem of estimating interval
More informationImproved Spatial Localization in 3D MRSI with a Sequence Combining PSF-Choice, EPSI and a Resolution Enhancement Algorithm
Improved Spatial Localization in 3D MRSI with a Sequence Combining PSF-Choice, EPSI and a Resolution Enhancement Algorithm L.P. Panych 1,3, B. Madore 1,3, W.S. Hoge 1,3, R.V. Mulkern 2,3 1 Brigham and
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 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 informationRecovery of Piecewise Smooth Images from Few Fourier Samples
Recovery of Piecewise Smooth Images from Few Fourier Samples Greg Ongie*, Mathews Jacob Computational Biomedical Imaging Group (CBIG) University of Iowa SampTA 2015 Washington, D.C. 1. Introduction 2.
More 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 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 informationA Bayesian Approach to SENSE Image Reconstruction in FMRI
A Bayesian Approach to SENSE Image Reconstruction in FMRI Daniel B. Rowe Program in Computational Sciences Department of Mathematics, Statistics, and Computer Science Marquette University July 31, 017
More informationA Study of Nonlinear Approaches to Parallel Magnetic Resonance Imaging
University of Wisconsin Milwaukee UWM Digital Commons Theses and Dissertations December 2012 A Study of Nonlinear Approaches to Parallel Magnetic Resonance Imaging Yuchou Chang University of Wisconsin-Milwaukee
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 informationA low rank based seismic data interpolation via frequencypatches transform and low rank space projection
A low rank based seismic data interpolation via frequencypatches transform and low rank space projection Zhengsheng Yao, Mike Galbraith and Randy Kolesar Schlumberger Summary We propose a new algorithm
More informationAdaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited
Adaptive Waveform Inversion: Theory Mike Warner*, Imperial College London, and Lluís Guasch, Sub Salt Solutions Limited Summary We present a new method for performing full-waveform inversion that appears
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 informationMidterm Review
Midterm Review - 2017 EE369B Concepts Noise Simulations with Bloch Matrices, EPG Gradient Echo Imaging 1 About the Midterm Monday Oct 30, 2017. CCSR 4107 Up to end of C2 1. Write your name legibly on this
More informationHigh dynamic range magnetic resonance flow imaging in the abdomen
High dynamic range magnetic resonance flow imaging in the abdomen Christopher M. Sandino EE 367 Project Proposal 1 Motivation Time-resolved, volumetric phase-contrast magnetic resonance imaging (also known
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 informationHead motion in diffusion MRI
Head motion in diffusion MRI Anastasia Yendiki HMS/MGH/MIT Athinoula A. Martinos Center for Biomedical Imaging 11/06/13 Head motion in diffusion MRI 0/33 Diffusion contrast Basic principle of diffusion
More informationChapter 3 Set Redundancy in Magnetic Resonance Brain Images
16 Chapter 3 Set Redundancy in Magnetic Resonance Brain Images 3.1 MRI (magnetic resonance imaging) MRI is a technique of measuring physical structure within the human anatomy. Our proposed research focuses
More informationMorphological Pyramids in Multiresolution MIP Rendering of. Large Volume Data: Survey and New Results
Morphological Pyramids in Multiresolution MIP Rendering of Large Volume Data: Survey and New Results Jos B.T.M. Roerdink Institute for Mathematics and Computing Science University of Groningen P.O. Box
More informationCalculating the Distance Map for Binary Sampled Data
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Calculating the Distance Map for Binary Sampled Data Sarah F. Frisken Gibson TR99-6 December 999 Abstract High quality rendering and physics-based
More informationsurface Image reconstruction: 2D Fourier Transform
2/1/217 Chapter 2-3 K-space Intro to k-space sampling (chap 3) Frequenc encoding and Discrete sampling (chap 2) Point Spread Function K-space properties K-space sampling principles (chap 3) Basic Contrast
More informationCHAPTER 3 DIFFERENT DOMAINS OF WATERMARKING. domain. In spatial domain the watermark bits directly added to the pixels of the cover
38 CHAPTER 3 DIFFERENT DOMAINS OF WATERMARKING Digital image watermarking can be done in both spatial domain and transform domain. In spatial domain the watermark bits directly added to the pixels of the
More informationWEINER FILTER AND SUB-BLOCK DECOMPOSITION BASED IMAGE RESTORATION FOR MEDICAL APPLICATIONS
WEINER FILTER AND SUB-BLOCK DECOMPOSITION BASED IMAGE RESTORATION FOR MEDICAL APPLICATIONS ARIFA SULTANA 1 & KANDARPA KUMAR SARMA 2 1,2 Department of Electronics and Communication Engineering, Gauhati
More 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 information(a Scrhon5 R2iwd b. P)jc%z 5. ivcr3. 1. I. ZOms Xn,s. 1E IDrAS boms. EE225E/BIOE265 Spring 2013 Principles of MRI. Assignment 8 Solutions
EE225E/BIOE265 Spring 2013 Principles of MRI Miki Lustig Assignment 8 Solutions 1. Nishimura 7.1 P)jc%z 5 ivcr3. 1. I Due Wednesday April 10th, 2013 (a Scrhon5 R2iwd b 0 ZOms Xn,s r cx > qs 4-4 8ni6 4
More informationImage Enhancement Techniques for Fingerprint Identification
March 2013 1 Image Enhancement Techniques for Fingerprint Identification Pankaj Deshmukh, Siraj Pathan, Riyaz Pathan Abstract The aim of this paper is to propose a new method in fingerprint enhancement
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 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 informationA Nuclear Norm Minimization Algorithm with Application to Five Dimensional (5D) Seismic Data Recovery
A Nuclear Norm Minimization Algorithm with Application to Five Dimensional (5D) Seismic Data Recovery Summary N. Kreimer, A. Stanton and M. D. Sacchi, University of Alberta, Edmonton, Canada kreimer@ualberta.ca
More informationParallel Magnetic Resonance Imaging (pmri): How Does it Work, and What is it Good For?
Parallel Magnetic Resonance Imaging (pmri): How Does it Work, and What is it Good For? Nathan Yanasak, Ph.D. Chair, AAPM TG118 Department of Radiology Georgia Regents University Overview Phased-array coils
More informationCompressed Sensing Algorithm for Real-Time Doppler Ultrasound Image Reconstruction
Mathematical Modelling and Applications 2017; 2(6): 75-80 http://www.sciencepublishinggroup.com/j/mma doi: 10.11648/j.mma.20170206.14 ISSN: 2575-1786 (Print); ISSN: 2575-1794 (Online) Compressed Sensing
More 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 informationA Transmission Line Matrix Model for Shielding Effects in Stents
A Transmission Line Matrix Model for hielding Effects in tents Razvan Ciocan (1), Nathan Ida (2) (1) Clemson University Physics and Astronomy Department Clemson University, C 29634-0978 ciocan@clemon.edu
More informationA MORPHOLOGY-BASED FILTER STRUCTURE FOR EDGE-ENHANCING SMOOTHING
Proceedings of the 1994 IEEE International Conference on Image Processing (ICIP-94), pp. 530-534. (Austin, Texas, 13-16 November 1994.) A MORPHOLOGY-BASED FILTER STRUCTURE FOR EDGE-ENHANCING SMOOTHING
More informationIterative CT Reconstruction Using Curvelet-Based Regularization
Iterative CT Reconstruction Using Curvelet-Based Regularization Haibo Wu 1,2, Andreas Maier 1, Joachim Hornegger 1,2 1 Pattern Recognition Lab (LME), Department of Computer Science, 2 Graduate School in
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 informationUniqueness in Bilinear Inverse Problems with Applications to Subspace and Sparsity Models
Uniqueness in Bilinear Inverse Problems with Applications to Subspace and Sparsity Models Yanjun Li Joint work with Kiryung Lee and Yoram Bresler Coordinated Science Laboratory Department of Electrical
More informationSampling, Ordering, Interleaving
Sampling, Ordering, Interleaving Sampling patterns and PSFs View ordering Modulation due to transients Temporal modulations Timing: cine, gating, triggering Slice interleaving Sequential, Odd/even, bit-reversed
More informationClinical Importance. Aortic Stenosis. Aortic Regurgitation. Ultrasound vs. MRI. Carotid Artery Stenosis
Clinical Importance Rapid cardiovascular flow quantitation using sliceselective Fourier velocity encoding with spiral readouts Valve disease affects 10% of patients with heart disease in the U.S. Most
More informationSkåne University Hospital Lund, Lund, Sweden 2 Deparment of Numerical Analysis, Centre for Mathematical Sciences, Lund University, Lund, Sweden
Volume Tracking: A New Method for Visualization of Intracardiac Blood Flow from Three-Dimensional, Time-Resolved, Three-Component Magnetic Resonance Velocity Mapping Appendix: Theory and Numerical Implementation
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 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 informationAn Intuitive Explanation of Fourier Theory
An Intuitive Explanation of Fourier Theory Steven Lehar slehar@cns.bu.edu Fourier theory is pretty complicated mathematically. But there are some beautifully simple holistic concepts behind Fourier theory
More informationImproving the quality of compressed sensing MRI that exploits adjacent slice similarity
Improving the quality of compressed sensing MRI that exploits adjacent slice similarity Norihito Inamuro and Akira Hirabayashi College of Information Science and Engineering, Ritsumeikan University, Kusatsu,
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 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 informationStereo Vision. MAN-522 Computer Vision
Stereo Vision MAN-522 Computer Vision What is the goal of stereo vision? The recovery of the 3D structure of a scene using two or more images of the 3D scene, each acquired from a different viewpoint in
More informationComputational Medical Imaging Analysis Chapter 4: Image Visualization
Computational Medical Imaging Analysis Chapter 4: Image Visualization Jun Zhang Laboratory for Computational Medical Imaging & Data Analysis Department of Computer Science University of Kentucky Lexington,
More informationEE795: Computer Vision and Intelligent Systems
EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 WRI C225 Lecture 04 130131 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Review Histogram Equalization Image Filtering Linear
More information29 th NATIONAL RADIO SCIENCE CONFERENCE (NRSC 2012) April 10 12, 2012, Faculty of Engineering/Cairo University, Egypt
K1. High Performance Compressed Sensing MRI Image Reconstruction Ahmed Abdel Salam, Fadwa Fawzy, Norhan Shaker, Yasser M.Kadah Biomedical Engineering, Cairo University, Cairo, Egypt, ymk@k-space.org Computer
More informationACQUIRING AND PROCESSING SUSCEPTIBILITY WEIGHTED IMAGING (SWI) DATA ON GE 3.0T
ACQUIRING AND PROCESSING SUSCEPTIBILITY WEIGHTED IMAGING (SWI) DATA ON GE 3.0T Revision date: 12/13/2010 Overview Susceptibility Weighted Imaging (SWI) is a relatively new data acquisition and processing
More informationVD-AUTO-SMASH Imaging
Magnetic Resonance in Medicine 45:1066 1074 (2001) VD-AUTO-SMASH Imaging Robin M. Heidemann, Mark A. Griswold, Axel Haase, and Peter M. Jakob* Recently a self-calibrating SMASH technique, AUTO-SMASH, was
More informationIndex. aliasing artifacts and noise in CT images, 200 measurement of projection data, nondiffracting
Index Algebraic equations solution by Kaczmarz method, 278 Algebraic reconstruction techniques, 283-84 sequential, 289, 293 simultaneous, 285-92 Algebraic techniques reconstruction algorithms, 275-96 Algorithms
More informationSpatial Enhancement Definition
Spatial Enhancement Nickolas Faust The Electro- Optics, Environment, and Materials Laboratory Georgia Tech Research Institute Georgia Institute of Technology Definition Spectral enhancement relies on changing
More informationECE 8201: Low-dimensional Signal Models for High-dimensional Data Analysis
ECE 8201: Low-dimensional Signal Models for High-dimensional Data Analysis Yuejie Chi Departments of ECE and BMI The Ohio State University September 24, 2015 Time, location, and office hours Time: Tue/Thu
More informationSegmentation and Modeling of the Spinal Cord for Reality-based Surgical Simulator
Segmentation and Modeling of the Spinal Cord for Reality-based Surgical Simulator Li X.C.,, Chui C. K.,, and Ong S. H.,* Dept. of Electrical and Computer Engineering Dept. of Mechanical Engineering, National
More informationDiffuse Optical Tomography, Inverse Problems, and Optimization. Mary Katherine Huffman. Undergraduate Research Fall 2011 Spring 2012
Diffuse Optical Tomography, Inverse Problems, and Optimization Mary Katherine Huffman Undergraduate Research Fall 11 Spring 12 1. Introduction. This paper discusses research conducted in order to investigate
More 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 informationFast methods for magnetic resonance angiography (MRA)
Fast methods for magnetic resonance angiography (MRA) Bahareh Vafadar Department of Electrical and Computer Engineering A thesis presented for the degree of Doctor of Philosophy University of Canterbury,
More informationModule 5: Dynamic Imaging and Phase Sharing. (true-fisp, TRICKS, CAPR, DISTAL, DISCO, HYPR) Review. Improving Temporal Resolution.
MRES 7005 - Fast Imaging Techniques Module 5: Dynamic Imaging and Phase Sharing (true-fisp, TRICKS, CAPR, DISTAL, DISCO, HYPR) Review Improving Temporal Resolution True-FISP (I) True-FISP (II) Keyhole
More informationFeature Selection for fmri Classification
Feature Selection for fmri Classification Chuang Wu Program of Computational Biology Carnegie Mellon University Pittsburgh, PA 15213 chuangw@andrew.cmu.edu Abstract The functional Magnetic Resonance Imaging
More informationVolume 2, Issue 9, September 2014 ISSN
Fingerprint Verification of the Digital Images by Using the Discrete Cosine Transformation, Run length Encoding, Fourier transformation and Correlation. Palvee Sharma 1, Dr. Rajeev Mahajan 2 1M.Tech Student
More informationImage Compression Algorithm and JPEG Standard
International Journal of Scientific and Research Publications, Volume 7, Issue 12, December 2017 150 Image Compression Algorithm and JPEG Standard Suman Kunwar sumn2u@gmail.com Summary. The interest in
More informationThe Benefit of Tree Sparsity in Accelerated MRI
The Benefit of Tree Sparsity in Accelerated MRI Chen Chen and Junzhou Huang Department of Computer Science and Engineering, The University of Texas at Arlington, TX, USA 76019 Abstract. The wavelet coefficients
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 information