Fast Implementation of Iterative Image Reconstruction
|
|
- Arnold Dickerson
- 5 years ago
- Views:
Transcription
1 Fast Implementation of Iterative Image Reconstruction Tobias Schaeffter 1, Michael Schacht Hansen 2, Thomas Sangild Sørensen 3 1 Division of Imaging Sciences and Biomedical Engineering, Kings College London 2 National Heart, Lung, and Blood Institute, NIH, Bethesda, MD, USA 3 Department of Computer Science and Institute of Clinical Medicine, Aarhus University Background Iterative reconstruction is becoming an increasingly important topic in medical imaging. In computed tomography (CT) iterative reconstruction algorithms are becoming standard and they are implemented on current CT-scanners to reduce radiation dose [1]. This is not a new trend; the very first CT machines used an iterative technique, the algebraic reconstruction technique (ART) [2]. In emission tomography in nuclear medicine iterative reconstructions were used in the early days to integrate prior knowledge about the noise statistics into the reconstruction [3] and more advanced iterative algorithms have become a standard in nuclear medicine [4]. Also in MRI there is a trend for advanced reconstruction algorithms to integrate knowledge about coil sensitivities in parallel MRI [5], about motion [6-9] and other prior information into the reconstruction process. All these approaches result in large set of linear equations Ax=b, which need to be solved. A direct inversion of this equation system is impractical, except in the most trivial cases, since it requires a large amount of memory and a prohibitive number of operations. However, this complexity can be reduced by using iterative reconstruction algorithms, where direct inversion of the system matrix is not needed. These algorithms start with an initial guess for the solution and then iteratively progress towards an approximate solution. A variety of such techniques exist for the treatment of large linear systems. For many cases, the so-called conjugate gradient (CG) method [11] particularly suited. It converges to a minimum in N 2 steps (where N is the rank of the matrix) if the matrix A is symmetric positive definite. The rate of convergence can be improved if preconditioning of the matrix A is performed. However, also in well-conditioned systems still around iterations need to be performed. Therefore one of the main limitations of introducing iterative algorithms in image reconstruction has been the long computation time so far hindering the clinical translations of many new reconstruction algorithms. Modern algorithms can make use of acceleration techniques and parallelization to speed up the reconstruction. These techniques make use of the on-going developments in computer processors (CPU) and especially in commodity graphics hardware (GPU). General purpose computation on the GPU is still an emerging research area, but many algorithms have been accelerated significantly compared to their CPU counterparts [12, 13]. One of the most recent GPUs, the Geforce 560 (Nvidia, US), has 352 stream processors for these computations allowing significant acceleration if the computational problem can be solved in parallel. More importantly, both major hardware manufacturers (AMD and Nvidia, US) have released high-level programming environments requiring less knowledge of graphics programming in general. The GPU implementation of MR-reconstruction algorithms has been demonstrated for a number of applications [14-20]. In this presentation the GPU-implementation of an iterative reconstruction algorithm is presented for the example of non-cartesian parallel imaging. However, the basic concepts can be applied to other applications.
2 Parallel Imaging using radial trajectories Non-uniform sampling schemes are more and more used in MR-scanning since they are often more robust against motion and flow. In particular, radial imaging results in more acceptable artifact due to residual motion and changes of the imaging plane during interactive scanning [22]. Furthermore, for radial acquisitions the central part of k-space is fully sampled, which can be used directly for estimating the coil sensitivity and regularization [9]. Non-Cartesian sampling can be combined with parallel imaging techniques such as SENSE to reconstruct images from highly undersampled data [23]. In parallel imaging the data are acquired simultaneously by multiple receiver coils with unique spatial sensitivity profiles S c (r). These profiles encode the unknown image v(r) such that a measured complex value m c (k) sampled by the c-th coil in k-space by: m c (k) = v(r) S c (r) e ikr dr (1) For the discrete case this equation can be written as the linear system of equations: m = Ev (2) using the encoding matrix E: E (γ,κ ),ρ = S γ, (r ρ )e ik κr ρ (3) After noise-decorrelation of the data m one solution of this system is given in ref. 23 by: ( E H E)v = E H m (4) where H is the conjugate transpose of the matrix. In order to solve this equation system (4) by an iterative conjugate gradient optimization, in each step of the iteration multiplications of the encoding matrices E and E H with some temporary vectors have to be performed. These operations are given by equation (3), which are multiplications with the coil sensitivity and Fourier transforms from k-space to image space and vice versa. However, since a non-uniform sampling is used a gridding procedure (including density correction) in combination with a Fast Fourier Transform (FFT) can be used. In Figure 1 it is shown that in each iteration of the conjugate gradient algorithm the residuum image of each coil element undergoes subsequent operations in the image domain and in k-space: re-gridding in k-space from non-uniform to uniform sampling Fast Fourier Transform multiplication with the coil sensitivity (and its conjugate) (inverse) Fast Fourier Transform gridding from uniform to non-uniform sampling
3 Figure 1: Iterative reconstruction using a conjugate gradient algorithm: S coil sensitivity; FFT Fast Fourier Transform; GRID regridding from non-cartesian to cartesian grid or vice versa; Figure adapted from ref. (5) From a computational load perspective most of the operations are simple multiplications and thus the gridding operations and the FFTs are the more demanding. A disadvantage of the proposed algorithm is that the convergence speed depends on how the encoding matrix is conditioned. Recently, a reconstruction algorithm based on Tikhonov regularization has been proposed that reduces the SNR loss due to geometric correlations in the spatial information [24]. The Tikhonov regularization provides a framework for stabilizing the solution of ill-conditioned linear equations by solving an equation system with a least squared estimate: 2 v λ = argmin Ev m 2 + λ 2 v 1 2 { regv 2 } (5) With v reg the inverse of a reference image for regularization and λ the regularization parameter. This approach results in the following equation system [16]: D E H ϒE + λ 2 2 ( v reg )D(D 1 v) = DE H ϒm (6) This equation system is similar to (4) but includes a reference image as a priori information for regularization. Furthermore, the matrix D is a diagonal preconditioning matrix aiming to equalize the diagonal elements of the system matrix on the left side. ϒ is another diagonal matrix that performs density compensation to account the non-uniform sampling in k-space. It should be noted that a similar system needs to be solved for k-t SENSE using a training image data set for regularization. Again in order to solve the equation system by a conjugate gradient optimization, for each iteration different vector matrix multiplications have to be performed. In particular the multiplication of the current search vector with the encoding matrix includes the multiplication the coil sensitivity matrix S followed by a non-equispaced Fourier transform from image space to k-space for each coil, resulting in distinct signals. The conjugate transpose of the encoding matrix E H, consists of a nonequispaced Fourier transform from k-space to image space followed by a
4 multiplication of the complex conjugate of the coil matrix for each coil. Since these computational expensive operations have to be performed for each iteration step, they need be optimized by an efficient implementation. Implementation of Non- equidistant Fourier transform on GPU A description of the implementation and an evaluation of the GPU-based nonequidistant fast Fourier transform (NFFT or gridding) can be found in references [14,16]. Here we summarize the fundamental considerations. The NFFT consists of three steps, namely 1) convolution, 2) FFT, and 3) deapodization. The deapodization step is merely element-wise multiplication of two vectors and thus trivial to implement in parallel. An efficient GPU implementation of the FFT is typically provided by the programming API (e.g. CUDA) and therefore also straightforward to apply. The implementation of the convolution operation on the other hand needs careful consideration for a data-parallel implementation. First, let us consider how a typical sequential CPU implementation would proceed: For both the NFFT and NFFT H implementation we iterate over the acquired samples one by one. For every sample we then iterative over all grid cells within the range of the convolution kernel and compute each convolution contribution. For the NFFT H this contribution is added to the resulting grid cells one by one within the inner loop. For the NFFT on the other hand the contributions of the inner loop are accumulated within the loop and written to the corresponding index in the resulting samples array. For a data-parallel approach these operations must take place in parallel. For the forwards convolution (the NFFT) this is straightforward; we can start a thread for each sample and have each thread perform the corresponding inner loop of the sequential approach. The NFFT H convolution is unfortunately more problematic. If we again parallelized by starting a thread per sample and each sample performs the inner loop of the sequential approach we encounter write conflicts (and undefined behaviour) as multiple threads attempt to update the same Cartesian grid cell in parallel. The solution, in sequential terms, is to interchange the order of the two loops. Correspondingly, in parallel terms, we invoke a thread for each grid cell. Each thread then accumulates the contribution from the samples within the radius of the convolution kernel. This strategy imposes a new challenge however; how to determine (in constant time) which samples that fall within the convolution radius of each cell? We refer the reader to reference [16] for a description of an appropriate pre-processing step. To maximize the performance of the GPU implementation of both the NFFT and NFFT H it is crucial to reduce the obvious bottleneck: the required memory bandwidth. Several strategies can be combined: 1) Use online computation of the convolution kernel instead of predefined lookup tables. Much of the cost of the computation can be hidden by the memory lookup of the corresponding sample or cell. 2) Use on-chip memory to hold accumulation variables and avoid unnecessary writes to main memory. References [2xTMI] describe a strategy in which each thread accumulates the contribution from multiple threads or cells. This reduces memory bandwidth requirements at the expense of parallelism as fewer threads are launched. For small image sizes the reduction in the number of concurrent threads launched is undesired. For parallel MRI one can consider to instead use the sparse on-chip memory to hold the convolution results from the multiple coils.
5 Applications There is wide range of applications that benefit from fast reconstruction. However, reconstructions using the GPU is not always possible due to memory constrains. In the presentation two different applications are demonstrated that benefit from fast reconstruction: 1) Real-time 2D imaging with real-time reconstruction. For this reconstruction computer is connected directly to a streaming output with of a scanner. 2) Offline reconstruction of 3D data that can be separated into a number of independent 2D iterative reconstruction problems. For the real-time cardiac MR-application data was acquired on a Siemens Avanto 1.5 T system using 4 coil-elements. A radial steady state free precession (SSFP) pulse sequence (FOV=300mm, PBW=1500Hz, Flip-angle=60, TR=2.5ms TE=1.2ms) was used to sample projections (each 256 samples) with a golden angular spacing of Figure 2 shows a real-time reconstruction of cardiac function obtained in a clinically referred pediatric patient during free breathing. For each reconstructed frame either 34 or 50 projections were used and a moderate regularization (λ 2 =0.05) was applied. A fixed number of iterations (15 and 20) in the iterative SENSE algorithm ensured that the GPU-reconstruction of all channels was always faster than data acquisition. The temporal fidelity of the reconstruction is demonstrated in a separated plot of one line spatial line over time. Figure 2: Real-time frame of radial SSFP acquisition. The GPU implementation of iterative SENSE reconstruction allows real-time reconstruction, i.e. faster than the data acquisition. a) Sliding window reconstruction using 34 projections and 15 iterations to update an image frame takes 83ms (< 85ms acquisition time). b) Sliding window reconstruction using 50 projections and 20 iterations to update an image frame takes 120ms (125ms acquisition time). For the 3D reconstruction example a radial phase encoding data reconstruction was tested. This encoding scheme combines Cartesian sampling along the readout direction with a radial sampling scheme in the phase encoding plane was suggested [25]. In particular, to reconstruct the 3D volume, initially the fast Fourier Transform (FFT) is applied in the fully sampled readout direction, converting the remaining reconstruction problem into a series (i.e. umber of samples along the readout directions) of two-dimensional radial reconstructions. Cardiac-triggered SSFP imaging (90 flip angle, TR/TE=5.0 ms/2.28 ms) with radial phase encoding was performed during free breathing (6 mm navigator window) to acquire a 3D dataset in end diastole using a resolution of mm 3 (in 4 minutes assuming 100% navigator efficiency) on a 1.5T system (Philips Achieva). The dataset contains 64*128=8192 lines of k-space data using Cartesian sampling in the readout (feet-head) direction (32 receive channels). In the phase encoding plane these were
6 laid out as 64 radial profiles of 256 samples (undersampling 2) each. Data were obtained on Iterative SENSE is required to resolve undersampling along the radial (Rr=2) and angular direction (Ra=4). Usually, this reconstruction is very time consuming; and reconstruction time of 3 minutes/slice using 32 coils was reported for a Matlab implementation resulting in total reconstruction time more than 12 hours [25]. Therefore the performance of GPU-implementation of the iterative SENSE reconstruction was tested. For this each 2D radial iterative SENSE reconstruction was performed on the GPU adapting parts of the implementation using a target matrix size of 256x256 for 256 slices. Furthermore, the coil sensitivity maps were acquired in a separate reference scan and no explicit regularization term was included in the reconstructions. We used 10 iterations of the conjugate gradient solver. The data was reconstructed on an Nvidia Geforce 280 GTX GPU with 1GB of memory. Figure 3 shows three reformatted slices (transversal, coronal, and sagittal view) from the three reconstructed, isotropic 3D volumes. The reconstruction time was reduced to 1.2s per slice. These numbers include data transfers from disc taking roughly 0.5s. In total, the full 3D datasets (256 slices) were each reconstructed in 5 minutes. Figure 3: Multiplanar reformatted slices of a 3D whole heart dataset using an iterative SENSE reconstruction. The GPU-implementation reduced the reconstruction per slice to 1.2 seconds. Conclusion Iterative reconstruction algorithms are becoming increasingly important in advanced MR-acquisition schemes. However, the clinical application is often hindered by long reconstruction times. In particular, the performance growth of individual CPU cores has stagnated and is unlikely to improve notably in the future. The iterative SENSE algorithm presented in this presentation is good candidate for implementation on GPU. This concept can be applied to other iterative algorithms making the application of advanced reconstructions feasible in clinical applications. References [1] Fleischmann D, Boas FE. Computed tomography-old ideas and new technology. Eur Radiol. 2011; 21(3): [2] Macovski A, Herman GT. Principles of reconstruction algorithms. In Newton TH, Potts DG (eds) Radiology of the skull and brain: technical aspects of computed tomography. Mosby, Saint Louis, MO, 1981; [3] Rockmore AJ, Macovski A. A maximum likelihood approach to emission image reconstruction from projections. IEEE Trans Nucl Sci : [4] Vandenberghe S, D'Asseler Y, Van de Walle R, Kauppinen T, Koole M, Bouwens L, Van Laere K,
7 Lemahieua I, Dierckx RA. Iterative reconstruction algorithms in nuclear medicine. Computerized medical imaging and graphics, (2), [5] Pruessmann KP. Encoding and reconstruction in parallel MRI. NMR Biomed. 2006;19(3): [6] Tsao J, Boesiger P, Pruessmann KP, k-t BLAST and k-t SENSE: dynamic MRI with high frame rate exploiting spatiotemporal correlations, Magnetic Resonance in Medicine 2003 vol. 50, [7] Hansen M, Baltes C, Tsao J, Kozerke S, Pruessmann KP, Eggers H, kt BLAST reconstruction from non-cartesian k-t SPACE sampling, Magn. Reson. Med 2006, 55, 1, [8] Batchelor PG, Atkinson D, Irarrazaval P, Hill DLG, Hajnal J, Larkman D, Matrix description of general motion correction applied to multishot images, Magnetic Resonance in Medicine 2005 vol. 54, no. 5, p [9] Odille F, Vuissoz PA, Marie PY, Felblinger J. Generalized reconstruction by inversion of coupled systems (GRICS) applied to free-breathing MRI. Magn Reson Med Jul;60(1): [10] Prieto C, Uribe S, Razavi R, Atkinson D, Schaeffter T. 3D undersampled golden-radial phase encoding for DCE-MRA using inherently regularized iterative SENSE. Magn Reson Med Aug;64(2): [11] Shewchuk JR, School of Computer Science, An Introduction to the Conjugate Gradient Method Without the Agonizing Pain, Carnegie Mellon University Pittsburgh, 1994, [12] Owens JD, Luebke D, Govindaraju N, Harris M, Krüger J, Lefohn AE, Purcell TJ, A survey of general-purpose computation on graphics hardware, Computer Graphics Forum 2007 vol. 26, no. 1, pp , [13] Sorensen TS, Mosegaard J, An Introduction to GPU Accelerated Surgical Simulation. Berlin, Germany: Springer, 2006, vol. 4072, pp [14] Sorensen TS, Schaeffter T, Noe KO, Hansen MS.Accelerating the nonequispaced fast Fourier transform on commodity graphics hardware. IEEE Trans Med Imaging. 2008;27(4): [15] Hansen MS, Atkinson D, Sorensen TS. Cartesian SENSE and k-t SENSE reconstruction using commodity graphics hardware. Magn Reson Med Mar;59(3):463-8 [16] Sørensen, T. S., Atkinson, D., Schaeffter, T., & Hansen, M. S. (2009). Real-time reconstruction of sensitivity encoded radial magnetic resonance imaging using a graphics processing unit IEEE transactions on medical imaging 2009, 28(12), [17] Knoll F, Unger M, Diwoky C, Clason C, Pock T, Stollberger R. Fast reduction of undersampling artifacts in radial MR angiography with 3D total variation on graphics hardware. MAGMA Apr;23(2): [18] Roujol S, de Senneville BD, Vahala E, Sørensen TS, Moonen C, Ries M. Online real-time reconstruction of adaptive TSENSE with commodity CPU/GPU hardware. Magn Reson Med. 2009; 62(6): [19] Kim D, Trzasko J, Smelyanskiy M, Haider C, Dubey P, Manduca A.High-Performance 3D Compressive Sensing MRI Reconstruction Using Many-Core Architectures. Int J Biomed Imaging. 2011:473128; [20] Thomas Sangild Sørensen, Claudia Prieto, David Atkinson, Michael Schacht Hansen, Tobias Schaeffter. GPU accelerated iterative SENSE reconstruction of radial phase encoded whole-heart MRI Proc. ISMRM 2010 [21] Eggers H, Mazurkewitz P, Boesiger P. Noise amplification in parallel imaging. Proc. ISMRM Workshop Non-Cartesian Imaging, [22] Schaeffter T, Weiss S, Eggers H, Rasche V. Projection reconstruction balanced fast field echo for interactive real-time cardiac imaging. Magn. Reson. Med., 2001; 46, 6, [23] Pruessmann KP, Weiger M, Börnert P, Boesiger P. Advances in sensitivity encoding with arbitrary k-space trajectories. Magnetic Resonance in Medicine 2001, 46(4), [24] Lin, F.-H., Kwong, K. K., Belliveau, J. W., Wald, L. L. Parallel imaging reconstruction using automatic regularization. Magnetic Resonance in Medicine 2004, 51(3), [25] Boubertakh R, Prieto C, Batchelor PG, Uribe S, Atkinson D, Eggers H, Sørensen TS, Hansen MS, Razavi RS, Schaeffter T. Whole-heart imaging using undersampled radial phase encoding (RPE) and iterative sensitivity encoding (SENSE) reconstruction. Magn Reson Med. 2009;62(5):
Role 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 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 informationInstitute of Cardiovascular Science, UCL Centre for Cardiovascular Imaging, London, United Kingdom, 2
Grzegorz Tomasz Kowalik 1, Jennifer Anne Steeden 1, Bejal Pandya 1, David Atkinson 2, Andrew Taylor 1, and Vivek Muthurangu 1 1 Institute of Cardiovascular Science, UCL Centre for Cardiovascular Imaging,
More informationDevelopment of fast imaging techniques in MRI From the principle to the recent development
980-8575 2-1 2012 10 13 Development of fast imaging techniques in MRI From the principle to the recent development Yoshio MACHIDA and Issei MORI Health Sciences, Tohoku University Graduate School of Medicine
More informationk-t BLAST Reconstruction From Non-Cartesian k-t Space Sampling
Magnetic Resonance in Medicine 55:85 91 (006) k-t BLAST Reconstruction From Non-Cartesian k-t Space Sampling Michael S. Hansen, 1, * Christof Baltes, 1 Jeffrey Tsao, 3 Sebastian Kozerke, 1 Klaas P. Pruessmann,
More informationTHIS PAPER describes a parallel implementation of the
538 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 27, NO. 4, APRIL 2008 Accelerating the Nonequispaced Fast Fourier Transform on Commodity Graphics Hardware Thomas Sangild Sørensen*, Tobias Schaeffter, Karsten
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 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 informationFree-Breathing Whole-Heart Coronary MRA: Motion Compensation Integrated into 3D Cartesian Compressed Sensing Reconstruction
Free-Breathing Whole-Heart Coronary MRA: Motion Compensation Integrated into 3D Cartesian Compressed Sensing Reconstruction Christoph Forman 1,2, Robert Grimm 1, Jana Hutter 1,2, Andreas Maier 1,2, Joachim
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 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 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 informationSpiral keyhole imaging for MR fingerprinting
Spiral keyhole imaging for MR fingerprinting Guido Buonincontri 1, Laura Biagi 1,2, Pedro A Gómez 3,4, Rolf F Schulte 4, Michela Tosetti 1,2 1 IMAGO7 Research Center, Pisa, Italy 2 IRCCS Stella Maris,
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 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 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 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 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 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 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 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 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 informationGadgetron: An Open Source Framework for Medical Image Reconstruction
COMPUTER PROCESSING AND MODELING - Full Papers Magnetic Resonance in Medicine 69:1768 1776 (2013) Gadgetron: An Open Source Framework for Medical Image Reconstruction Michael Schacht Hansen 1 * and Thomas
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 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 informationTHE STEP BY STEP INTERACTIVE GUIDE
COMSATS Institute of Information Technology, Islamabad PAKISTAN A MATLAB BASED INTERACTIVE GRAPHICAL USER INTERFACE FOR ADVANCE IMAGE RECONSTRUCTION ALGORITHMS IN MRI Medical Image Processing Research
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 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 informationSingle Breath-hold Abdominal T 1 Mapping using 3-D Cartesian Sampling and Spatiotemporally Constrained Reconstruction
Single Breath-hold Abdominal T 1 Mapping using 3-D Cartesian Sampling and Spatiotemporally Constrained Reconstruction Felix Lugauer 1,3, Jens Wetzl 1, Christoph Forman 2, Manuel Schneider 1, Berthold Kiefer
More informationRespiratory Motion Compensation for Simultaneous PET/MR Based on Strongly Undersampled Radial MR Data
Respiratory Motion Compensation for Simultaneous PET/MR Based on Strongly Undersampled Radial MR Data Christopher M Rank 1, Thorsten Heußer 1, Andreas Wetscherek 1, and Marc Kachelrieß 1 1 German Cancer
More 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 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 informationMotion Artifacts and Suppression in MRI At a Glance
Motion Artifacts and Suppression in MRI At a Glance Xiaodong Zhong, PhD MR R&D Collaborations Siemens Healthcare MRI Motion Artifacts and Suppression At a Glance Outline Background Physics Common Motion
More informationGRAPPA Operator for Wider Radial Bands (GROWL) with Optimally Regularized Self-Calibration
GRAPPA Operator for Wider Radial Bands (GROWL) with Optimally Regularized Self-Calibration Wei Lin,* Feng Huang, Yu Li, and Arne Reykowski Magnetic Resonance in Medicine 64:757 766 (2010) A self-calibrated
More informationFrom Image to Video: Real-time Medical Imaging with MRI
From Image to Video: Real-time Medical Imaging with MRI Sebastian Schaetz, Martin Uecker BiomedNMR Forschungs GmbH at the MPI for biophysical Chemistry, Goettingen, Germany Electrical Engineering and Computer
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 informationImage reconstruction using compressed sensing for individual and collective coil methods.
Biomedical Research 2016; Special Issue: S287-S292 ISSN 0970-938X www.biomedres.info Image reconstruction using compressed sensing for individual and collective coil methods. Mahmood Qureshi *, Muhammad
More informationGADGETRON. Michael S. Hansen, PhD. Magnetic Resonance Technology Program National Institutes of Health - NHLBI
GADGETRON Michael S. Hansen, PhD Magnetic Resonance Technology Program National Institutes of Health - NHLBI QUESTIONS/COMMENTS EMAIL: michael.hansen@nih.gov Twi8er: @ReconstructThis Outline Gadgetron
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 informationNIH Public Access Author Manuscript Med Phys. Author manuscript; available in PMC 2009 March 13.
NIH Public Access Author Manuscript Published in final edited form as: Med Phys. 2008 February ; 35(2): 660 663. Prior image constrained compressed sensing (PICCS): A method to accurately reconstruct dynamic
More informationGeneral and Efficient Super-Resolution Method for Multi-slice MRI
General and Efficient Super-Resolution Method for Multi-slice MRI D.H.J. Poot 1,2,V.VanMeir 2, and J. Sijbers 2 1 BIGR, Erasmus Medical Center, Rotterdam 2 Visionlab, University of Antwerp, Antwerp Abstract.
More informationGraz, Austria 2 Institute of Computer Graphics and Vision, Graz University of Technology, Inffeldgasse
Magnetic Resonance Materials in Physics, Biology and Medicine Fast Reduction of Undersampling Artifacts in Radial MR Angiography with 3D Total Variation on Graphics Hardware Florian Knoll 1, Markus Unger
More informationMotion compensated reconstruction
Motion compensated reconstruction Educational course: Image acquisition and Reconstruction 25th Annual meeting of the ISMRM, Honolulu, 2017 Sajan Goud Lingala Siemens Healthineers, Princeton, USA Declaration
More informationSpread Spectrum Using Chirp Modulated RF Pulses for Incoherent Sampling Compressive Sensing MRI
Spread Spectrum Using Chirp Modulated RF Pulses for Incoherent Sampling Compressive Sensing MRI Sulaiman A. AL Hasani Department of ECSE, Monash University, Melbourne, Australia Email: sulaiman.alhasani@monash.edu
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 informationFOREWORD TO THE SPECIAL ISSUE ON MOTION DETECTION AND COMPENSATION
Philips J. Res. 51 (1998) 197-201 FOREWORD TO THE SPECIAL ISSUE ON MOTION DETECTION AND COMPENSATION This special issue of Philips Journalof Research includes a number of papers presented at a Philips
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 information3D MAGNETIC RESONANCE IMAGING OF THE HUMAN BRAIN NOVEL RADIAL SAMPLING, FILTERING AND RECONSTRUCTION
3D MAGNETIC RESONANCE IMAGING OF THE HUMAN BRAIN NOVEL RADIAL SAMPLING, FILTERING AND RECONSTRUCTION Maria Magnusson 1,2,4, Olof Dahlqvist Leinhard 2,4, Patrik Brynolfsson 2,4, Per Thyr 2,4 and Peter Lundberg
More informationhigh performance medical reconstruction using stream programming paradigms
high performance medical reconstruction using stream programming paradigms This Paper describes the implementation and results of CT reconstruction using Filtered Back Projection on various stream programming
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 informationParallel 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 informationQualitative Comparison of Conventional and Oblique MRI for Detection of Herniated Spinal Discs
Qualitative Comparison of Conventional and Oblique MRI for Detection of Herniated Spinal Discs Doug Dean Final Project Presentation ENGN 2500: Medical Image Analysis May 16, 2011 Outline Review of the
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 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 informationSuper-resolution Reconstruction of Fetal Brain MRI
Super-resolution Reconstruction of Fetal Brain MRI Ali Gholipour and Simon K. Warfield Computational Radiology Laboratory Children s Hospital Boston, Harvard Medical School Worshop on Image Analysis for
More informationOnline Real-Time Reconstruction of Adaptive TSENSE With Commodity CPU/GPU Hardware
Online Real-Time Reconstruction of Adaptive TSENSE With Commodity CPU/GPU Hardware Sébastien Roujol, 1 Baudouin Denis de Senneville, 1 Erkki Vahala, 2 Thomas Sangild Sørensen, 3 Chrit Moonen, 1 and Mario
More informationA practical acceleration algorithm for real-time imaging
A practical acceleration algorithm for real-time imaging The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation As Published Publisher
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 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 informationGPUs Open New Avenues in Medical MRI
GPUs Open New Avenues in Medical MRI Chris A. Cocosco D. Gallichan, F. Testud, M. Zaitsev, and J. Hennig Dept. of Radiology, Medical Physics, UNIVERSITY MEDICAL CENTER FREIBURG 1 Our research group: Biomedical
More informationDYNAMIC magnetic resonance imaging (MRI), which captures
IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL 26, NO 7, JULY 2007 917 Accelerating Dynamic Spiral MRI by Algebraic Reconstruction From Undersampled k t Space Taehoon Shin*, Student Member, IEEE, Jon-Fredrik
More informationGeorgia Institute of Technology, August 17, Justin W. L. Wan. Canada Research Chair in Scientific Computing
Real-Time Rigid id 2D-3D Medical Image Registration ti Using RapidMind Multi-Core Platform Georgia Tech/AFRL Workshop on Computational Science Challenge Using Emerging & Massively Parallel Computer Architectures
More informationWhole-Heart Cine MRI Using Real-Time Respiratory Self-Gating
Magnetic Resonance in Medicine 57:606 613 (2007) Whole-Heart Cine MRI Using Real-Time Respiratory Self-Gating Sergio Uribe, 1,2 Vivek Muthurangu, 2 Redha Boubertakh, 1,2 Tobias Schaeffter, 2 Reza Razavi,
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 informationTomographic Reconstruction
Tomographic Reconstruction 3D Image Processing Torsten Möller Reading Gonzales + Woods, Chapter 5.11 2 Overview Physics History Reconstruction basic idea Radon transform Fourier-Slice theorem (Parallel-beam)
More informationcuibm A GPU Accelerated Immersed Boundary Method
cuibm A GPU Accelerated Immersed Boundary Method S. K. Layton, A. Krishnan and L. A. Barba Corresponding author: labarba@bu.edu Department of Mechanical Engineering, Boston University, Boston, MA, 225,
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 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 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 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 informationEnhao Gong, PhD Candidate, Electrical Engineering, Stanford University Dr. John Pauly, Professor in Electrical Engineering, Stanford University Dr.
Enhao Gong, PhD Candidate, Electrical Engineering, Stanford University Dr. John Pauly, Professor in Electrical Engineering, Stanford University Dr. Greg Zaharchuk, Associate Professor in Radiology, Stanford
More informationConstructing System Matrices for SPECT Simulations and Reconstructions
Constructing System Matrices for SPECT Simulations and Reconstructions Nirantha Balagopal April 28th, 2017 M.S. Report The University of Arizona College of Optical Sciences 1 Acknowledgement I would like
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 informationHighly Efficient Respiratory Motion Compensated Free-Breathing Coronary MRA Using Golden-Step Cartesian Acquisition
JOURNAL OF MAGNETIC RESONANCE IMAGING 41:738 746 (2015) Technical Development Highly Efficient Respiratory Motion Compensated Free-Breathing Coronary MRA Using Golden-Step Cartesian Acquisition Claudia
More informationDeep Learning for Fast and Spatially- Constrained Tissue Quantification from Highly-Undersampled Data in Magnetic Resonance Fingerprinting (MRF)
Deep Learning for Fast and Spatially- Constrained Tissue Quantification from Highly-Undersampled Data in Magnetic Resonance Fingerprinting (MRF) Zhenghan Fang 1, Yong Chen 1, Mingxia Liu 1, Yiqiang Zhan
More informationAlgebraic Iterative Methods for Computed Tomography
Algebraic Iterative Methods for Computed Tomography Per Christian Hansen DTU Compute Department of Applied Mathematics and Computer Science Technical University of Denmark Per Christian Hansen Algebraic
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 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 informationImage Registration. Prof. Dr. Lucas Ferrari de Oliveira UFPR Informatics Department
Image Registration Prof. Dr. Lucas Ferrari de Oliveira UFPR Informatics Department Introduction Visualize objects inside the human body Advances in CS methods to diagnosis, treatment planning and medical
More 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 informationHYPERDRIVE IMPLEMENTATION AND ANALYSIS OF A PARALLEL, CONJUGATE GRADIENT LINEAR SOLVER PROF. BRYANT PROF. KAYVON 15618: PARALLEL COMPUTER ARCHITECTURE
HYPERDRIVE IMPLEMENTATION AND ANALYSIS OF A PARALLEL, CONJUGATE GRADIENT LINEAR SOLVER AVISHA DHISLE PRERIT RODNEY ADHISLE PRODNEY 15618: PARALLEL COMPUTER ARCHITECTURE PROF. BRYANT PROF. KAYVON LET S
More information3D Registration based on Normalized Mutual Information
3D Registration based on Normalized Mutual Information Performance of CPU vs. GPU Implementation Florian Jung, Stefan Wesarg Interactive Graphics Systems Group (GRIS), TU Darmstadt, Germany stefan.wesarg@gris.tu-darmstadt.de
More informationCompressed Sensing MRI with Multichannel Data Using Multicore Processors
Compressed Sensing MRI with Multichannel Data Using Multicore Processors Ching-Hua Chang and Jim Ji* Magnetic Resonance in Medicine 64:1135 1139 (2010) Compressed sensing (CS) is a promising method to
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 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 informationTradeoffs 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 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 informationImproved 3D image plane parallel magnetic resonance imaging (pmri) method
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. 311 316, Hamilton,
More informationTotal Variation Regularization Method for 3D Rotational Coronary Angiography
Total Variation Regularization Method for 3D Rotational Coronary Angiography Haibo Wu 1,2, Christopher Rohkohl 1,3, Joachim Hornegger 1,2 1 Pattern Recognition Lab (LME), Department of Computer Science,
More informationarxiv: v1 [cs.dc] 7 Jan 2013
A multi-gpu Programming Library for Real-Time Applications Sebastian Schaetz 1 and Martin Uecker 2 arxiv:1301.1215v1 [cs.dc] 7 Jan 2013 1 BiomedNMR Forschungs GmbH at the Max Planck Institute for biophysical
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 informationDepth-Layer-Based Patient Motion Compensation for the Overlay of 3D Volumes onto X-Ray Sequences
Depth-Layer-Based Patient Motion Compensation for the Overlay of 3D Volumes onto X-Ray Sequences Jian Wang 1,2, Anja Borsdorf 2, Joachim Hornegger 1,3 1 Pattern Recognition Lab, Friedrich-Alexander-Universität
More informationFiber Selection from Diffusion Tensor Data based on Boolean Operators
Fiber Selection from Diffusion Tensor Data based on Boolean Operators D. Merhof 1, G. Greiner 2, M. Buchfelder 3, C. Nimsky 4 1 Visual Computing, University of Konstanz, Konstanz, Germany 2 Computer Graphics
More informationBlind Sparse Motion MRI with Linear Subpixel Interpolation
Blind Sparse Motion MRI with Linear Subpixel Interpolation Anita Möller, Marco Maaß, Alfred Mertins Institute for Signal Processing, University of Lübeck moeller@isip.uni-luebeck.de Abstract. Vital and
More informationIterative Methods for Solving Linear Problems
Iterative Methods for Solving Linear Problems When problems become too large (too many data points, too many model parameters), SVD and related approaches become impractical. Iterative Methods for Solving
More informationGadgetron implementation for PET/MR motion correction
implementation for motion correction Thomas Küstner School of Biomedical Engineering & Imaging Sciences, King s College London, St. Thomas Hospital, London, UK Dept. of Diagnostic and Interventional Radiology,
More informationFast Cardiac CINE MRI by Iterative Truncation of Small Transformed Coefficients
pissn 2384-1095 eissn 2384-1109 imri 2015;19:19-30 http://dx.doi.org/10.13104/imri.2015.19.1.19 Fast Cardiac CINE MRI by Iterative Truncation of Small Transformed Coefficients Jinho Park, Hye-Jin Hong,
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 informationAccelerating a Simulation of Type I X ray Bursts from Accreting Neutron Stars Mark Mackey Professor Alexander Heger
Accelerating a Simulation of Type I X ray Bursts from Accreting Neutron Stars Mark Mackey Professor Alexander Heger The goal of my project was to develop an optimized linear system solver to shorten the
More informationCompressed Sensing Reconstruction for Whole-Heart Imaging with 3D Radial Trajectories: A Graphics Processing Unit Implementation
Magnetic Resonance in Medicine 69:91 102 (2013) Compressed Sensing Reconstruction for Whole-Heart Imaging with 3D Radial Trajectories: A Graphics Processing Unit Implementation Seunghoon Nam, 1,2 Mehmet
More information