Accelerated 4D flow MRI

Size: px
Start display at page:

Download "Accelerated 4D flow MRI"

Transcription

1 Accelerated 4D flow MRI Comparing SENSE, k-t PCA and Compressed Sensing Tomas Kaandorp Patient comfort can be increased by accelerating an MRI scanner to reduce scan time. 4D flow MRI scans can be accelerated by undersampling the k-space. Unfortunately this has a negative effect on the image quality and signal-to-noise ratio. To look at the effects of undersampling on the flow amplitude, SENSE, k-t PCA and Compressed Sensing (CS) are compared for acceleration factors (R): 6, 8, and 10 in five healthy volunteers. The results show that the maximum flow amplitude decreases for all scanning modalities if R increases. The average drop in flow amplitude for SENSE, k-t PCA and CS are respectively: 33 %, 20 %, and 7.5 %. Furthermore CS showed a better image quality for higher R compared to SENSE and k-t PCA, which leads to reliable flow curves and usable image quality up to R10. Thus CS shows promising results for reliable shortening of 4D flow MRI in clinical environments. Report Bachelor Project Physics and Astronomy, size 15 EC, conducted between and Supervisor: Dr. Ir. A. J. Nederveen Second assessor: Prof. Dr. Ir. G. J. Strijkers Mentor: L. M. Gottwald Department of Radiology and Nucleair medicine Amsterdam Medical Centre University of Amsterdam August 2017

2 Contents 1 Introduction 2 2 Theory MRI Acceleration methods Hypothesis Method Phantom scan In vivo SNR Results Phantom In vivo Prospective vs retrospective triggering SNR Discussion 19 6 Conclusion 20 7 References 21 1

3 1 Introduction The global average life expectancy has tremendously increased for decades. This means that the number of people dying from cardiovascular diseases (CV) is on the rise, and is now even the primary cause of premature deaths (Heissel, 2017). For the correct diagnosis and treatment a 4D flow MRI scan might be beneficial. For example, a 4D flow MRI scan allows physicians to look at blood flow, peak blood flow and wall shear stress. However these scans needed for the diagnosis are high resolution scans which means they take a long time. This in term means that patients have to lie inside the scanner for a prolonged period of time, which can be claustrophobic or uncomfortable. To increase patient well-being it is vital to reduce scanning time of 4D flow MRI. However accelerating an MRI scan has a negative trade-off to the image quality of the scan, this trade-off is expressed in terms of the Signal-to-Noise Ratio (SNR). It is therefor vital to finding a way to shorten the scan time without the drop in image quality. Over the years the scan time of 4D flow has been shortened i.e. by three methods. First by the implementation of multiple receiving coils (Pruessmann et. al., 1999), secondly by the use of k-t Principal Component Analyses (k-t PCA) and more recently by implementing Compressed Sensing (CS). Since the aorta is of interest to diagnose CV, this paper will compare the aorta blood flow for SENSE, k-t PCA and CS to see if the acceleration has an effect on the reliability of the blood flow quantization. 2 Theory 2.1 MRI A MRI scanner uses the 1 H spin of the Nucleus to acquire an image. Firstly a subject (i.e. volunteer or patient) is placed into a homogeneous magnetic field (the scanner). By doing so, all the nuclear spins of the atoms align in to the foot head direction. Once all the spins are aligned, a Radio Frequency (RF) pulse is given with a frequency matching the Larmor frequency of 1 H. This RF-pulse makes the nuclear spin resonate. When resonating the spins send out a RF signal which can be measured with the receiving coil. This can be done via T1 imaging of T2 imaging, depending on what sort of tissue that is scanned. The difference is determined by the time the 1 H nucleus remains resonating. To transfer the frequencies measured during acquisition a Fourier transform is used. The frequencies are measured in the Frequency space, also called k-space. In principle the entire k-space needs to be scanned to get a good image quality. If only parts of the k-space are scanned the image becomes blurry and other artifacts appear. Over the years a number of methods where devised to scan the k-space. Initially k-space was scanned in horizontal row (k x direction) progressing downwards in the vertical direction (k y ) after scanning each horizontal row. 2

4 Figure 1: The Nyquist criterion sets the required k-space coverage, which can be achieved using various sampling trajectories. Image resolution is determined by the extent of the k-space coverage. The supported field of view is determined by the sampling density. Violation of the Nyquist criterion causes artifacts in linear reconstructions, which depend on the sampling pattern. (Lustig et. al., 2008) Although there are multiple ways of scanning the k-space a Cartesian grid is preferred since this can be easily transformed from frequency domain to space domain by a Fourier transform. New techniques such as spiral trajectories based on radial coordinates are gaining popularity since these are less susceptible to motion artifacts and can be more easily under-sampled than Cartesian trajectories (Lustig et. al., 2008). 2.2 Acceleration methods Accelerating MRI scans is mainly done by collecting less data points. Data is collected in the frequency domain (k-space). By collecting less data points in the k-space, i.e. under sampling the k-space, the duration of the scan shortens. However this is where a trade off arises. A fully sampled k-space is needed for high resolution images which physicians like, but this takes time. Under sampling the k-space shortens the scan time but this also means the image quality of the scan deteriorates. This trade off is expressed in the Signal-toNoise Ratio (Formula 1) SN R = Signalintensity SD(N oiseintensity) 3 (1)

5 In short, if a scan is accelerated, the SNR decreases which means the image quality decreases. One way of accelerating MRI scans is by using a SENSE reconstruction. During a SENSE scan, data is collected using multiple receiver coils. This means there are multiple images created (i. e. 2). These images contain aliasing effects since they are based on only part of the k-space. By later combining these multiple images, one complete image is created and aliasing effects are minimized. This however requires the scanner to save all the images during the scan and to form a reconstructed image by superposition of these saved images weighted by the corresponding coil sensitivity maps. Figure 2: SENSE graphically displayed, two images with aliasing formed by different coils. Afterwards these images are combined using a superposition weighted by coil sensitivity maps (MRIquestions.com, 2017) Another way of reducing scan time is by under sampling the k-space over time, this method is called k-t PCA. Next an adaptive filter is used in the Fourier transformed domain to correct for the aliasing effects that occur due to the under sampling. However, the reconstruction is now under determined since the filter leads to more unknown values than equations. This has a negative effect on the temporal fidelity. To solve this problem k-t BLAST is revised to k-t PCA. Principle Component Analysis (PCA) reduces a multi-dimensional data set to a lower dimensional set by using a training image. PCA constrains the appearance of the object such that the frequency must be a linear combination of predefined basis functions. The basis functions are derived from the training image. The PCA reconstruction separates the overlapped images accurately. Since the temporal fidelity is also separated, there is no negative effect on the fidelity. Therefore, k-t PCA then leads to a shorter 4D flow scan time. Unfortunately this still leads to some aliasing in the image. This means the image quality decreases and it may become difficult to define a correct Region Of Interest (ROI). Thus, limiting the acceleration factor (R). Gwenaël Pagé et. al., 4

6 2017 showed that it is possible to accelerate the scan using R = 8 and retain a good image quality. However the flow amplitude decreases significantly above R6. Thus limiting 4D flow MRI. Figure 3: An example of k-t PCA. Firstly the training image is made (a), then the data is under sampled (b). By using a filter to correct for the aliasing effect, the image becomes under determined which is solved with the training image(c). (Petersen et. al., 2009) A more recent development is called Compressed Sensing (CS). CS also under samples the k-space but has some more requirements. Firstly the data has to be under sampled in a random way. Secondly, transform sparsity is needed and lastly, to acquire the best image quality, an iterative reconstruction method is needed. There are again a number of ways to scan the k-space, a Poisson distribution in Cartesian coordinates can be used. However a spiral shape with more data points in the centre than on the outer radius shows promising results in scan time reduction and image quality(l. M. Gottwald et. al., 2017). The collected data has to be incoherent/sparse to make it possible to separate noise from the signal. Via a wavelet transform, the data is made sparse which is necessary to do the further data analyses. In the sparse data set, the strongest signal is identified and abstracted from the whole signal. In the next iteration of data process, the signal then allows for the identification of a signal with a lower intensity. After multiple iterations only the noise remains and the complete spectrum of signals is identified. In Figure 2, the CS reconstruction is graphically 5

7 displayed. After the separation of the signals, the image can be reconstructed through the regular methods i.e. a Fast Fourier transform(fft). Because all signals are deduced from the under sampled k-space, the image quality remains acceptable for clinical use (Lustig et. al., 2008). Figure 4: Heuristic procedure for reconstruction from under sampled data. A sparse signal (a) is 8-fold under sampled in its 1-D k-space domain (b). Equispaced under sampling results in signal aliasing (d) preventing recovery. Pseudorandom under sampling results in incoherent interference (c). Some strong signal components stick above the interference level, are detected and recovered by thresholding (e) and (f). The interference of these components is computed (g) and subtracted (h), thus lowering the total interference level and enabling recovery of weaker components.(lustig et. al., Hypothesis Based on the paper from Gwenaël Pagé et. al., 2017, it is expected that for SENSE and k-t PCA the amplitude of the flow curves decrease if the acceleration factor (R) is higher then R6. Secondly, since the SNR decreases if the acceleration rate increases, it is to be expected that the image quality decreases. Since this effect is greater for SENSE than for k-t PCA or CS, it is expected that SENSE has the worst image quality and thus lowest SNR, followed by k-t PCA. Since the CS reconstruction is better in determining the signals with low intensity because it uses an iterative reconstruction method, CS should have a better image quality. Therefore the average SNR of CS should be higher than the average SNR of k-t PCA and SENSE. 6

8 3 Method 3.1 Phantom scan Firstly, an exploratory scan was made with a carotid bifurcation phantom (LifeTec, Eindhoven, Netherlands). A pulsating water flow with a simulated heart rate of 60 bpm and a variability of 5 bpm created a temporal mean flow of 300 ml/min. A time stamp from the pump signal was used for retrospective triggering. A schematic overview of the experimental setup and all measurement parameters are shown in Figure 5. Figure 5: Experimental setup used for the Phantom scans, (L. M. Gottwald, 2017) All experiments were conducted on a 3T MRI scanner (Philips Healthcare, Best, Netherlands) using an 8-channel neck coil. A total matrix size(fov) of 64x64x64 mm 3 and a voxel size of 1x1x1 mm 3 was acquired in three flow encoding directions using a VENC of 150 cm/s. Profile lists of different lengths were created with a nominal acceleration factor (R) of 2, 4, 6, 8 and 10 per acceleration method i.e. SENSE and k-t PCA. The acquired data was processed by ReconFrame (Gyrotools, Zurich, Zwiterland) and retrospectively absolute binned in 24 cardiac frames. For the flow analysis the return artery was used since this was the largest phantom artery and thus yields the highest signal. Next Matlab(The Math- Works Inc, r2016b) was used to draw Regions Of Interest (ROI). Based on these ROI, a Black and White (BW) mask was created. This BW mask was then used to make the flow curves since it multiplies the data in the aorta with 1 and all other areas with 0. The flow curves are compared to the pre-set maximum flow and pulsation of the phantom. The flow was quantified using the Phase image. The phase data was multiplied with the BW mask so only the Aorta was selected. Next the flow data was scaled to the voxel size and the preset VENC. This was done using Formula 2 and 3. V elocityroi = mean(bw mask Dataset) (2) AverageF lowroi = mean( V elocityroi BW Area ) (3) 100 7

9 3.2 In vivo Secondly the flow was compared in vivo. This was also done on a 3T MRI scanner (Phillips Healthcare, Best, Netherlands) using a 16-channel Torso and an 8-channel posterior coil. Five healthy volunteers, aged between years where scanned. A total FOV size of 315x275x60 mm 3 and a VENC of 150 cm/s was used, the resolution was set to 2.5x2.5x2.5 mm 3. Data was collected in three different flow encoding directions for each acceleration method(sense, k-t PCA and CS). Scanning R2 and R4 would lead to a longer scanning time then 1.5 hours. It becomes very uncomfortable for a volunteer to lay inside the scanner longer then 1.5 hours, therefor R2 and R4 are not included in this paper. The ROI was drawn on the Ascending Aorta (AAO) for all 24 cardiac frames. next this was scaled for the phase offset. To calculate the flow Equation 2 and 3 where used again. Next the mean flow is determined over all volunteers. The flow curves are compared based on the acceleration rate (R) with fixed modality and also compared per modality with fixed R. Furthermore it was investigated if retrospective cardiac synchronization differs from prospective cardiac synchronization. This is done by comparing data, for a fixed R, between retrospective and prospective gating for CS and k-t PCA. This can be interesting since it may occur that data is not binned correctly in retrospective binning than with prospective binning (triggering). 3.3 SNR Lastly, the SNR is compared. This is done by making two extra ROI. The first one is inside the body but on a location where there should be no significant change in the signal intensity. The second ROI is then drawn outside the body where only noise is expected. By using Formula 1 the SNR is calculated. In contrary to the flow, the SNR is calculated based on the magnitude image and not on the phase image. This is done since the intensity in a phase image can change around the blood vessels due to blood flow. 4 Results 4.1 Phantom Firstly, a Phantom was examined for k-t PCA and SENSE, the flow curves per acceleration rate are displayed in Figure 6. As Pagé et. al., 2017, predicted, the amplitude of the flow curves declines as the acceleration rate (R) increases. Furthermore it can be seen that the lines become less smooth for increased R. Both these observations are due to the undersampling of the data. 8

10 Figure 6: Flow comparison in the Phantom, expressed for fixed R. Notice the drop in maximum flow amplitude. 9

11 4.2 In vivo First the image quality of the magnitude images are compared to see if the data is usable and how the undersampling effects the image quality. This comparison is done by eye based on Figure 7 Figure 7: on the x-axis R 6, 8, and 10 are displayed. On the y-axis CS, k-t PCA and Sense are displayed. On the top row it is visible that SENSE is very noisy and unclear. Below kt-pca shows a good image quality for all R although some streaking appears. On the bottom CS also shows good image quality. The image is a bit blurry due to scaling. CS is scaled by multiplying with 0.2 and SENSE by multiplying with 2 10

12 During the drawing of the ROI it became clear that the aorta was very difficult to locate on the accelerated SENSE scans. Since the volunteers where scanned for multiple scanning modalities, sometimes k-t PCA or CS scans where used to draw the ROI. Firstly, the differences in the flow curves is compared within each scanning modality, in Figure 8 the flow for all R is plotted for SENSE. Figure 8: Mean flow for SENSE. R 6, 8, and 10 are compared. R6 shows no real flow. This may be due to a small and incorrect data set. The max flow amplitude of R10 is lower then R8. It can be seen that the flow curves are very rough, also R6 preforms very bad compared to R8 and R10. This might be because the mean R6 flow is determined on less than five volunteers due to some failed scans. The decrease in amplitude between R8 and R10 is neatly visible. In Figure 9, k-t PCA is compared between acceleration rate (R). The flow curves preform better and are smoother compared to SENSE. Again the drop in amplitude is visible as predicted by Pagé et. al. On the x-axis the cardiac frames are displayed. In all flow curves we can see when the heart makes a beat. First a peak is visible since the heart is pumping, next the a slow relaxation occurs. During the relaxation some low peaks are still caused by the relaxation 11

13 of the slightly expanded aorta due to the blood pressure. On the y-axis, the Mean blood flow is given in ml/s. Figure 9: Mean flow for k-t PCA. R 6, 8, and 10 are compared. The drop in amplitude is visible and for low R the graphs are smoother due to less undersampling. Figure 10 shows the mean flow curves of CS. Here it can be seen that the max amplitude is almost constant. Again as expected the amplitude does drop, the difference between the highest peak (R6) and the lowest peak (R8) is ml/s. The amplitude difference between R8 and R10 can not be calculated due to a difference in cardiac frames. However since the peak flow amplitude of R10 lies between the peak flow of R8 and R6 this is still smaller then ml/s. 12

14 Figure 10: Mean flow for CS. R 6, 8, and 10 are compared. The difference in max flow between R8 and R6 is ml/s Next the differences per modality are compared. 13

15 Figure 11: Mean flow for SENSE, k-t PCA, and CS for R6 Figure 12: Mean flow for SENSE, k-t PCA, and CS for R8 14

16 Figure 13: Mean flow for SENSE, k-t PCA, and CS for R10. Again CS preforms stable with a maximum flow around 160 ml/s SENSE preforms well for R8 but decreases signifcantly compared to CS for R10. For k-t PCA a steady drop can be seen in maximum amplitude for all R compared to SENSE and CS. Also the curve remains reasonably smooth. The maximum difference in flow amplitude between R6 and R10 for SENSE, k-t PCA and CS are respectively: 33%, 20%, and 7.5%. 15

17 4.3 Prospective vs retrospective triggering It can be seen in Figures 14, 15, 16, and 17 that retrospective binning leads to a binning error at the end of the flow curves. A ascending line can be seen which should actually be at the beginning of the flow curve. Furthermore it can be seen that the amplitude for retrospective binning matches more to the amplitude of CS and SENSE. Figure 14: Mean flow for retrospective and prospective binning in the AAo accelerated to R8. A small binning error can be seen at the end of the flow curve for retrospective binning. 16

18 Figure 15: Velocity for retrospective and prospective binning in the AAo accelerated to R8. A small binning error can be seen at the end of the flow curve for retrospective binning. Figure 16: Mean flow for retrospective and prospective binning in the DAo accelerated to R8. A small binning error can be seen at the end of the flow curve for retrospective binning. 17

19 Figure 17: Velocity for retrospective and prospective binning in the DAo accelerated to R8. A small binning error can be seen at the end of the flow curve for retrospective binning. 18

20 4.4 SNR As the acceleration rate (R) increases, the increased under sampling of the k- space leads to a decrease in signal intensity. Therefor, the SNR calculated with Formula 1 should drop. Since the image quality for CS is also better then for SENSE or k-t PCA, it is also expected that the SNR is higher for CS than for SENSE and k-t PCA. However this is not what is seen in the data. The SNR values are given below in Table 1. SENSE k-t PCA CS R R R Tabel 1: average SNR expressed per method and R As seen in the table, CS has an higher average SNR than k-t PCA or SENSE over all volunteers. Furthermore there is an unexpected peak in the average SNR for R8 in SENSE and k-t PCA. As expected the SNR value for R10 are below R6 for k-t PCA and SENSE as should be based on the fact that the signal decreases if the undersampling increases. 5 Discussion Not all scans could be performed one after the other due to a software change during the examination, resulting in a scan break. Moreover, some exams had to be spitted into two exams, when the gating efficiency of the subject had been very poor. The composition of the subject scans and masks are shown in Table 2. 19

21 Volunteer: Mask 1 SENSE + kt-pca 4 SENSE + k-t PCA SENSE k-t PCA SENSE + k-t PCA Kt-PCA Mask 2 CS CS CS + SENSE 2.8, 3.2 CS CS + SENSE Mask 3 k-t PCA: 6, 8 Table 2: Composition of scanned volunteer and corresponding BW masks The drawing of the ROI remains a human task. This means it is susceptible to flaws such as drawing the ROI to big or to small. Drawing the ROI wrong can lead to misguided flow curves. One other solution is to only look at the velocity since this is not scaled for the area of the ROI. Furthermore, in determining the Signal-to-Noise Ratio (SNR), it was assumed that the signal outside of the body would only contain noise. It may be possible that the signal outside the body is also scaled depending on the scan modality used. Thus the used baseline of noise would differ per modality and it would not be able to make a good comparison. The average SNR for CS is about 20 times higher than the SNR for SENSE and k-t PCA. This might be explained by a scaling factor. There are also strong indications that the SNR is influenced by the multiple scanning sessions of volunteers and the software update of the scanner. For example the CS SNR for volunteer two is around 17 while all the other SNR are around 50. Volunteer two was scanned with CS before the scanner update whilst all other volunteers where scanned with CS after the scanner update. Therefore it is clear the scanner update has influenced either the scaling of the reconstruction method. 6 Conclusion Just like Pagé Et. Al. showed, the amplitude of the flow curves drops for all modalities if the acceleration rate increases. However the drop in flow amplitude for CS is limited to ml/s between R6 and R8. Since R10 lies between R6 and R8 the difference of ml/s is a maximum difference. The image quality does decrease if R increases. For SENSE the decrease in image quality means that it is no longer possible to draw the correct ROI. For k-t PCA the image quality remains acceptable but the max flow amplitude decreases significantly. For CS the image quality also remains acceptable and the decrease in max flow amplitude is limited to ml/s up to R = 10. This is roughly 7 % of the maximum flow for R6 whilst the drop in amplitude for SENSE and k-t PCA are respectively 33% and 20%. Therefore CS shows promising results for reliably shortening the scan time of a 4D flow MRI. However there are still other limiting factors such as the gating efficiency. It is advisable to do some further research into 4d flow MRI but also into the gating efficiency. 20

22 7 References Heiddel W, Deaths from cardiovasculair disease increase globally while mortality rates decrease, Healthdata, , retrieved Gottwald L. M, Peper E. S, Zang Q, Pronk V, Coolen B. F, Strijkers G. J, Nederveen A. J, Compressed Sensing accelerated 4D flow MRI using a pseudo spiral Cartesian sampling technique with random undersamling in time, Unknown, 2017 Lustig M, Donoh D, Santos J, Compressed sensing MRI, IEEE SIGNAL PROCESSING MAGAZINE,] MARCH 2008, MRI questions, retrieved Pagé G, Bettoni J, Salsac A.V, Balédent O, Influence of the k-t Principal Component Analysis acceleration factor on the accuracy of flow measurement in 4D PC-MRI, Proc. Intl. Soc. Mag. Reson. Med. 25 (2017) Pedersen H, Kozerke S, Ringgaard S, Nehrke K, and Kim W, k-t PCA: Temporally Constrained k-t BLAST Reconstruction Using Principal Component Analysis, Magnetic Resonance in Medicine, 2009, 62, Pruessmann K. P, Weiger M, Scheidegger M. B, Boesiger P, SENSE: Sensitivity Encoding for Fast MRI, Magnetic Resonance in Medicine, 1999, 42,

Compressed Sensing for Rapid MR Imaging

Compressed 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 information

Clinical Importance. Aortic Stenosis. Aortic Regurgitation. Ultrasound vs. MRI. Carotid Artery Stenosis

Clinical 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 information

High dynamic range magnetic resonance flow imaging in the abdomen

High 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 information

Sparse 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 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 information

Accelerated MRI Techniques: Basics of Parallel Imaging and Compressed Sensing

Accelerated 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 information

Redundancy Encoding for Fast Dynamic MR Imaging using Structured Sparsity

Redundancy 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 information

Exam 8N080 - Introduction MRI

Exam 8N080 - Introduction MRI Exam 8N080 - Introduction MRI Friday January 23 rd 2015, 13.30-16.30h For this exam you may use an ordinary calculator (not a graphical one). In total there are 6 assignments and a total of 65 points can

More information

Classification of Subject Motion for Improved Reconstruction of Dynamic Magnetic Resonance Imaging

Classification 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 information

Role of Parallel Imaging in High Field Functional MRI

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 information

G Practical Magnetic Resonance Imaging II Sackler Institute of Biomedical Sciences New York University School of Medicine. Compressed Sensing

G 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 information

Motion Artifacts and Suppression in MRI At a Glance

Motion 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 information

Combination 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 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 information

Module 5: Dynamic Imaging and Phase Sharing. (true-fisp, TRICKS, CAPR, DISTAL, DISCO, HYPR) Review. Improving Temporal Resolution.

Module 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 information

Compressed Sensing Reconstructions for Dynamic Contrast Enhanced MRI

Compressed 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 information

White 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 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 information

A 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 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

PERFORMANCE ANALYSIS BETWEEN TWO SPARSITY-CONSTRAINED MRI METHODS: HIGHLY CONSTRAINED BACKPROJECTION (HYPR) AND

PERFORMANCE ANALYSIS BETWEEN TWO SPARSITY-CONSTRAINED MRI METHODS: HIGHLY CONSTRAINED BACKPROJECTION (HYPR) AND PERFORMANCE ANALYSIS BETWEEN TWO SPARSITY-CONSTRAINED MRI METHODS: HIGHLY CONSTRAINED BACKPROJECTION (HYPR) AND COMPRESSED SENSING (CS) FOR DYNAMIC IMAGING A Thesis by NIBAL ARZOUNI Submitted to the Office

More information

Constrained Reconstruction of Sparse Cardiac MR DTI Data

Constrained 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 information

Fast Imaging Trajectories: Non-Cartesian Sampling (1)

Fast 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

Institute of Cardiovascular Science, UCL Centre for Cardiovascular Imaging, London, United Kingdom, 2

Institute 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 information

K-Space Trajectories and Spiral Scan

K-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 information

Brilliance CT Big Bore.

Brilliance CT Big Bore. 1 2 2 There are two methods of RCCT acquisition in widespread clinical use: cine axial and helical. In RCCT with cine axial acquisition, repeat CT images are taken each couch position while recording respiration.

More information

Weighted-CS for reconstruction of highly under-sampled dynamic MRI sequences

Weighted-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 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

(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 information

Deep 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) 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 information

Single 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 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 information

Image Acquisition Systems

Image Acquisition Systems Image Acquisition Systems Goals and Terminology Conventional Radiography Axial Tomography Computer Axial Tomography (CAT) Magnetic Resonance Imaging (MRI) PET, SPECT Ultrasound Microscopy Imaging ITCS

More information

An Iterative Approach for Reconstruction of Arbitrary Sparsely Sampled Magnetic Resonance Images

An 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 information

Deviceless respiratory motion correction in PET imaging exploring the potential of novel data driven strategies

Deviceless respiratory motion correction in PET imaging exploring the potential of novel data driven strategies g Deviceless respiratory motion correction in PET imaging exploring the potential of novel data driven strategies Presented by Adam Kesner, Ph.D., DABR Assistant Professor, Division of Radiological Sciences,

More information

Analysis of Functional MRI Timeseries Data Using Signal Processing Techniques

Analysis 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 information

Evaluations 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 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 information

Magnetic Resonance Angiography

Magnetic Resonance Angiography Magnetic Resonance Angiography Course: Advance MRI (BIOE 594) Instructors: Dr Xiaohong Joe Zhou Dr. Shadi Othman By, Nayan Pasad Phase Contrast Angiography By Moran 1982, Bryan et. Al. 1984 and Moran et.

More information

Slide 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. 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 information

Respiratory Motion Estimation using a 3D Diaphragm Model

Respiratory Motion Estimation using a 3D Diaphragm Model Respiratory Motion Estimation using a 3D Diaphragm Model Marco Bögel 1,2, Christian Riess 1,2, Andreas Maier 1, Joachim Hornegger 1, Rebecca Fahrig 2 1 Pattern Recognition Lab, FAU Erlangen-Nürnberg 2

More information

Development of fast imaging techniques in MRI From the principle to the recent development

Development 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 information

MR Advance Techniques. Vascular Imaging. Class III

MR Advance Techniques. Vascular Imaging. Class III MR Advance Techniques Vascular Imaging Class III 1 Vascular Imaging There are several methods that can be used to evaluate the cardiovascular systems with the use of MRI. MRI will aloud to evaluate morphology

More information

6 credits. BMSC-GA Practical Magnetic Resonance Imaging II

6 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 information

Parallel Imaging. Marcin.

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 information

NIH Public Access Author Manuscript Med Phys. Author manuscript; available in PMC 2009 March 13.

NIH 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 information

RADIOMICS: potential role in the clinics and challenges

RADIOMICS: potential role in the clinics and challenges 27 giugno 2018 Dipartimento di Fisica Università degli Studi di Milano RADIOMICS: potential role in the clinics and challenges Dr. Francesca Botta Medical Physicist Istituto Europeo di Oncologia (Milano)

More information

Image reconstruction using compressed sensing for individual and collective coil methods.

Image 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 information

Module 4. K-Space Symmetry. Review. K-Space Review. K-Space Symmetry. Partial or Fractional Echo. Half or Partial Fourier HASTE

Module 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 information

UNIVERSITY OF SOUTHAMPTON

UNIVERSITY OF SOUTHAMPTON UNIVERSITY OF SOUTHAMPTON PHYS2007W1 SEMESTER 2 EXAMINATION 2014-2015 MEDICAL PHYSICS Duration: 120 MINS (2 hours) This paper contains 10 questions. Answer all questions in Section A and only two questions

More information

Enhao 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. 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 information

The Near Future in Cardiac CT Image Reconstruction

The Near Future in Cardiac CT Image Reconstruction SCCT 2010 The Near Future in Cardiac CT Image Reconstruction Marc Kachelrieß Institute of Medical Physics (IMP) Friedrich-Alexander Alexander-University Erlangen-Nürnberg rnberg www.imp.uni-erlangen.de

More information

VD-AUTO-SMASH Imaging

VD-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 information

Chapter 3 Set Redundancy in Magnetic Resonance Brain Images

Chapter 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 information

CS 229 Project report: Extracting vital signs from video

CS 229 Project report: Extracting vital signs from video CS 229 Project report: Extracting vital signs from video D.Deriso, N. Banerjee, A. Fallou Abstract In both developing and developed countries, reducing the cost of medical care is a primary goal of science

More information

Edge-Preserving Denoising for Segmentation in CT-Images

Edge-Preserving Denoising for Segmentation in CT-Images Edge-Preserving Denoising for Segmentation in CT-Images Eva Eibenberger, Anja Borsdorf, Andreas Wimmer, Joachim Hornegger Lehrstuhl für Mustererkennung, Friedrich-Alexander-Universität Erlangen-Nürnberg

More information

Zigzag Sampling for Improved Parallel Imaging

Zigzag 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 information

Dynamic Contrast enhanced MRA

Dynamic Contrast enhanced MRA Dynamic Contrast enhanced MRA Speaker: Yung-Chieh Chang Date : 106.07.22 Department of Radiology, Taichung Veterans General Hospital, Taichung, Taiwan 1 Outline Basic and advanced principles of Diffusion

More information

New Technology Allows Multiple Image Contrasts in a Single Scan

New 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 information

MRI Physics II: Gradients, Imaging

MRI 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 information

Compressed Sensing And Joint Acquisition Techniques In Mri

Compressed Sensing And Joint Acquisition Techniques In Mri Wayne State University Wayne State University Theses 1-1-2013 Compressed Sensing And Joint Acquisition Techniques In Mri Rouhollah Hamtaei Wayne State University, Follow this and additional works at: http://digitalcommons.wayne.edu/oa_theses

More information

29 th NATIONAL RADIO SCIENCE CONFERENCE (NRSC 2012) April 10 12, 2012, Faculty of Engineering/Cairo University, Egypt

29 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 information

Optimal Sampling Geometries for TV-Norm Reconstruction of fmri Data

Optimal 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 information

Multi-Dimensional Flow-Preserving Compressed Sensing (MuFloCoS) for Time-Resolved Velocity-Encoded Phase Contrast MRI

Multi-Dimensional Flow-Preserving Compressed Sensing (MuFloCoS) for Time-Resolved Velocity-Encoded Phase Contrast MRI 400 IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 34, NO. 2, FEBRUARY 2015 Multi-Dimensional Flow-Preserving Compressed Sensing (MuFloCoS) for Time-Resolved Velocity-Encoded Phase Contrast MRI Jana Hutter*,

More information

CHAPTER 2. Morphometry on rodent brains. A.E.H. Scheenstra J. Dijkstra L. van der Weerd

CHAPTER 2. Morphometry on rodent brains. A.E.H. Scheenstra J. Dijkstra L. van der Weerd CHAPTER 2 Morphometry on rodent brains A.E.H. Scheenstra J. Dijkstra L. van der Weerd This chapter was adapted from: Volumetry and other quantitative measurements to assess the rodent brain, In vivo NMR

More information

This Time. fmri Data analysis

This Time. fmri Data analysis This Time Reslice example Spatial Normalization Noise in fmri Methods for estimating and correcting for physiologic noise SPM Example Spatial Normalization: Remind ourselves what a typical functional image

More information

EE123 Digital Signal Processing

EE123 Digital Signal Processing EE123 Digital Signal Processing Lecture 24 Compressed Sensing III M. Lustig, EECS UC Berkeley RADIOS https://inst.eecs.berkeley.edu/~ee123/ sp15/radio.html Interfaces and radios on Wednesday -- please

More information

THE STEP BY STEP INTERACTIVE GUIDE

THE 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 information

Correction of Partial Volume Effects in Arterial Spin Labeling MRI

Correction of Partial Volume Effects in Arterial Spin Labeling MRI Correction of Partial Volume Effects in Arterial Spin Labeling MRI By: Tracy Ssali Supervisors: Dr. Keith St. Lawrence and Udunna Anazodo Medical Biophysics 3970Z Six Week Project April 13 th 2012 Introduction

More information

Compressed sensing and undersampling k-space

Compressed sensing and undersampling k-space Boston University OpenBU Theses & Dissertations http://open.bu.edu Boston University Theses & Dissertations 2015 Compressed sensing and undersampling k-space Rahimpour, Yashar https://hdl.handle.net/2144/16257

More information

MRI image formation 8/3/2016. Disclosure. Outlines. Chen Lin, PhD DABR 3. Indiana University School of Medicine and Indiana University Health

MRI image formation 8/3/2016. Disclosure. Outlines. Chen Lin, PhD DABR 3. Indiana University School of Medicine and Indiana University Health MRI image formation Indiana University School of Medicine and Indiana University Health Disclosure No conflict of interest for this presentation 2 Outlines Data acquisition Spatial (Slice/Slab) selection

More information

Compressed Sensing Algorithm for Real-Time Doppler Ultrasound Image Reconstruction

Compressed 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 information

SPM8 for Basic and Clinical Investigators. Preprocessing. fmri Preprocessing

SPM8 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 information

Dynamic Autocalibrated Parallel Imaging Using Temporal GRAPPA (TGRAPPA)

Dynamic 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 information

Magnetic Resonance Elastography (MRE) of Liver Disease

Magnetic Resonance Elastography (MRE) of Liver Disease Magnetic Resonance Elastography (MRE) of Liver Disease Authored by: Jennifer Dolan Fox, PhD VirtualScopics Inc. jennifer_fox@virtualscopics.com 1-585-249-6231 1. Overview of MRE Imaging MRE is a magnetic

More information

XI Signal-to-Noise (SNR)

XI 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 information

7D Cardiac Flow MRI: Techniques & Automation of Reconstruction

7D Cardiac Flow MRI: Techniques & Automation of Reconstruction Cleveland State University EngagedScholarship@CSU ETD Archive 2012 7D Cardiac Flow MRI: Techniques & Automation of Reconstruction Michael G. Ambrosia Cleveland State University Follow this and additional

More information

Iterative CT Reconstruction Using Curvelet-Based Regularization

Iterative 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 information

Spiral keyhole imaging for MR fingerprinting

Spiral 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 information

A Study of Nonlinear Approaches to Parallel Magnetic Resonance Imaging

A 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 information

Functional MRI in Clinical Research and Practice Preprocessing

Functional 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 information

EPSRC Centre for Doctoral Training in Industrially Focused Mathematical Modelling

EPSRC Centre for Doctoral Training in Industrially Focused Mathematical Modelling EPSRC Centre for Doctoral Training in Industrially Focused Mathematical Modelling More Accurate Optical Measurements of the Cornea Raquel González Fariña Contents 1. Introduction... 2 Background... 2 2.

More information

Improved 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 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 information

CHAPTER 9: Magnetic Susceptibility Effects in High Field MRI

CHAPTER 9: Magnetic Susceptibility Effects in High Field MRI Figure 1. In the brain, the gray matter has substantially more blood vessels and capillaries than white matter. The magnified image on the right displays the rich vasculature in gray matter forming porous,

More information

FOREWORD TO THE SPECIAL ISSUE ON MOTION DETECTION AND COMPENSATION

FOREWORD 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 information

Lucy Phantom MR Grid Evaluation

Lucy Phantom MR Grid Evaluation Lucy Phantom MR Grid Evaluation Anil Sethi, PhD Loyola University Medical Center, Maywood, IL 60153 November 2015 I. Introduction: The MR distortion grid, used as an insert with Lucy 3D QA phantom, is

More information

Collaborative Sparsity and Compressive MRI

Collaborative Sparsity and Compressive MRI Modeling and Computation Seminar February 14, 2013 Table of Contents 1 T2 Estimation 2 Undersampling in MRI 3 Compressed Sensing 4 Model-Based Approach 5 From L1 to L0 6 Spatially Adaptive Sparsity MRI

More information

Blood Particle Trajectories in Phase-Contrast-MRI as Minimal Paths Computed with Anisotropic Fast Marching

Blood Particle Trajectories in Phase-Contrast-MRI as Minimal Paths Computed with Anisotropic Fast Marching Blood Particle Trajectories in Phase-Contrast-MRI as Minimal Paths Computed with Anisotropic Fast Marching Michael Schwenke 1, Anja Hennemuth 1, Bernd Fischer 2, Ola Friman 1 1 Fraunhofer MEVIS, Institute

More information

Denoising the Spectral Information of Non Stationary Image using DWT

Denoising the Spectral Information of Non Stationary Image using DWT Denoising the Spectral Information of Non Stationary Image using DWT Dr.DolaSanjayS 1, P. Geetha Lavanya 2, P.Jagapathi Raju 3, M.Sai Kishore 4, T.N.V.Krishna Priya 5 1 Principal, Ramachandra College of

More information

FMRI Pre-Processing and Model- Based Statistics

FMRI Pre-Processing and Model- Based Statistics FMRI Pre-Processing and Model- Based Statistics Brief intro to FMRI experiments and analysis FMRI pre-stats image processing Simple Single-Subject Statistics Multi-Level FMRI Analysis Advanced FMRI Analysis

More information

ANALYSIS OF PULMONARY FIBROSIS IN MRI, USING AN ELASTIC REGISTRATION TECHNIQUE IN A MODEL OF FIBROSIS: Scleroderma

ANALYSIS OF PULMONARY FIBROSIS IN MRI, USING AN ELASTIC REGISTRATION TECHNIQUE IN A MODEL OF FIBROSIS: Scleroderma ANALYSIS OF PULMONARY FIBROSIS IN MRI, USING AN ELASTIC REGISTRATION TECHNIQUE IN A MODEL OF FIBROSIS: Scleroderma ORAL DEFENSE 8 th of September 2017 Charlotte MARTIN Supervisor: Pr. MP REVEL M2 Bio Medical

More information

Abbie M. Diak, PhD Loyola University Medical Center Dept. of Radiation Oncology

Abbie M. Diak, PhD Loyola University Medical Center Dept. of Radiation Oncology Abbie M. Diak, PhD Loyola University Medical Center Dept. of Radiation Oncology Outline High Spectral and Spatial Resolution MR Imaging (HiSS) What it is How to do it Ways to use it HiSS for Radiation

More information

Respiratory 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 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 information

Satoshi Tamura Department of Electrical, Electronic and Computer Engineering, Gifu University, Gifu City, Japan

Satoshi Tamura Department of Electrical, Electronic and Computer Engineering, Gifu University, Gifu City, Japan A fast, noniterative approach for accelerated high-temporal resolution cine-cmr using dynamically interleaved streak removal in the power-spectral encoded domain with low-pass filtering (DISPEL) and modulo-prime

More information

Breast MRI Accreditation Program Clinical Image Quality Guide

Breast MRI Accreditation Program Clinical Image Quality Guide Breast MRI Accreditation Program Clinical Image Quality Guide Introduction This document provides guidance on breast MRI clinical image quality and describes the criteria used by the ACR Breast MRI Accreditation

More information

M R I Physics Course

M 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 information

UvA-DARE (Digital Academic Repository) Motion compensation for 4D PET/CT Kruis, M.F. Link to publication

UvA-DARE (Digital Academic Repository) Motion compensation for 4D PET/CT Kruis, M.F. Link to publication UvA-DARE (Digital Academic Repository) Motion compensation for 4D PET/CT Kruis, M.F. Link to publication Citation for published version (APA): Kruis, M. F. (2014). Motion compensation for 4D PET/CT General

More information

Separate Magnitude and Phase Regularization via Compressed Sensing

Separate Magnitude and Phase Regularization via Compressed Sensing IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 31, NO. 9, SEPTEMBER 2012 1713 Separate Magnitude and Phase Regularization via Compressed Sensing Feng Zhao*, Douglas C. Noll, Senior Member, IEEE, Jon-Fredrik

More information

Image Registration. Prof. Dr. Lucas Ferrari de Oliveira UFPR Informatics Department

Image 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 information

Sampling, Ordering, Interleaving

Sampling, 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 information

A novel noise removal using homomorphic normalization for multi-echo knee MRI

A novel noise removal using homomorphic normalization for multi-echo knee MRI A novel noise removal using homomorphic normalization for multi-echo knee MRI Xuenan Cui 1a),HakilKim 1b), Seongwook Hong 1c), and Kyu-Sung Kwack 2d) 1 School of Information and Communication Engineering,

More information

Assignment 2. Due Feb 3, 2012

Assignment 2. Due Feb 3, 2012 EE225E/BIOE265 Spring 2012 Principles of MRI Miki Lustig Assignment 2 Due Feb 3, 2012 1. Read Nishimura Ch. 3 2. Non-Uniform Sampling. A student has an assignment to monitor the level of Hetch-Hetchi reservoir

More information

Blind Sparse Motion MRI with Linear Subpixel Interpolation

Blind 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 information

Fast Dynamic MRI for Radiotherapy

Fast Dynamic MRI for Radiotherapy 1 Fast Dynamic MRI for Radiotherapy KE SHENG, PH.D., FAAPM DEPARTMENT OF RADIATION ONCOLOGY UNIVERSITY OF CALIFORNIA, LOS ANGELES Disclosures I receive research grants from Varian Medical Systems I am

More information

k-t BLAST Reconstruction From Non-Cartesian k-t Space Sampling

k-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 information

Lab 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 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 information

High dynamic range fusion of magnetic resonance flow imaging data

High dynamic range fusion of magnetic resonance flow imaging data High dynamic range fusion of magnetic resonance flow imaging data Christopher M. Sandino Stanford University EE 367 Final Project Report sandino@stanford.edu Abstract Time-resolved, volumetric phase-contrast

More information