RSL/MRSL Journal Club 8.27.2010 VIPR/Radial MRI Kitty Moran

1. RSL/MRSL Journal Club 8.27.2010 VIPR/Radial MRI Kitty Moran

2. Outline • VIPR • 2. Aliasing Artifact in Radial MRI

3. Outline 1. 2. • Simulations of Undersampling Artifact

4. Radial/Projection MRI

5. History of Radial MRI Projection Acquisition was the first MRI k-space trajectory (Lauterbur, 1973) Bernstein,et al. Handbook of MRI Pulse Sequences2004

6. History of Radial MRI • In initial MRI scanners rectilinear sampling was utilized because of B0 inhomogeneity and gradient non-linearity • With rectilinear k-space trajectories, the effect of these errors is a shift in the location and amplitude of the point spread function, resulting in geometric and intensity distortion • With PA the effect is not only a shift in location and amplitude but the shape of the point spread function is effected resulting in a blurring artifact Bernstein,et al. Handbook of MRI Pulse Sequences2004

7. Undersampling Cartesian K-space • Aliasing results in wrap around artifacts • Overlap in image space in the direction in which the undersampling occurs • Replications occur at frequency of 1/Δkr where Δkr is the sampling frequency in k-space • To avoid aliasing you want to meet the criteria 1/Δkr = FOV

8. Sampling in Polar Coordinates • Data sampled in the radial (r) and angular (θ) direction • In MRI the image is the inverse Fourier transform of the acquired signal data: • Translating to polar coordinates

9. Sampling in Polar Coordinates

10. Scheffler, Hennig, Reduced Circular Field-of-View Imaging, 2005 • Scheffler and Hennig considered 5 different sampling patterns in their investigation of aliasing artifacts in polar sampling • Sampling in either the kr and/or kθ direction • In the first two examples we consider continuous sampling in one direction and discrete sampling in the other

11. Discrete sampling in kr, continuous sampling in kθ Discrete concentric rings in k-space We assume S(kr,kθ) is circularly symmetric then S(kr,kθ)  S(kr) then

12. Discrete sampling in kr, continuous sampling in kθ Sum over zero order Bessel functions

13. Sum over zero order Bessel functions • Concentric ring lobes surrounding the main peak result in radial aliasing artifacts • Artifacts occur when object size exceeds the circular field of view diameter of 1/Δkr • Object smears circularly • The effective field of view (unaliased region) in polar sampling is a circle of radius 0.5/Δkr

14. Object size 0.375/Δkr • 3*1/Δkr shown • Object size 0.554/Δkr Lauzon, Rutt, MRM 36:940-949 (1996)

15. Discrete sampling in kθ, continuous sampling in kr • Fourier-Bessel integral becomes independent of kr • Corresponding integrals diverge • In the case of continuous radial sampling, the extent of sampling in the radial direction must be limited to a finite value otherwise total measuring time and signal power become infinite

16. Discrete sampling in kθ, continuous sampling in kr

17. Discrete sampling in kθ, continuous sampling in kr Where anθ = 2πrcos(nθΔkθ – θ)

18. PSF consists of radial stripes in the sampling direction that produce streak artifacts • Intensity plot along the angular coordinate θ for a fixed distance r is similar to the Dirichlet kernel for Cartesian sampling

19. For angular undersampling the radial increment Δkr is much smaller than the maximal angular increment Δkθ • Similar to case of continuous radial sampling leads to a star-like PSF • However, small circular flat region at the center • Angular undersampling allows for artifact free imaging as long as the radius of the object does not extend past the radius of the reduced field of view

20. Similar to the first type of continuous angular sampling • PSF shows equidistant ring lobes surrounding the main center peak • Amplitude of successive rings decreases as a function of the radial distance kr

21. “optimal sampling”, increased angular sampling • Circular flat region of the PSF linearly expands to the first ring lobe • Nθ = 2πNr • So intensity profile in the radial direction is composed of two parts • Within the interval (-rΔkr,rΔkr) the profile is approximately given by the Bessel function and therefore is similar to continuous angular sampling • Beyond the first ring lobe, wild oscillations occur

23. ¼ projections ½ projections Radial Imaging * Slide courtesy of Frank Korosec kx ky

24. VIPR • Vastly Undersampled Isotropic Projection Reconstruction • Contrast magnetic resonance angiography • Capture first-pass arterial enhancement • Separate arterial from venous enhancement • Time-resolved methods are desirable to eliminate dependence on specific bolus timing • Trade-off in MRA between spatial resolution and coverage can be greatly reduced by acquiring data with an undersampled 3D projection trajectory in which all three dimensions are symmetrically undersampled Barger, et al.

25. VIPR • Fully sampled 3D PR trajectories are not frequently used • Np = πNr2 • Aliased energy from undersampling the 3D projection trajectory resembles noise more than coherent streaks • Oversampling at the center of k-space and interleaving the 3D projections so that the spatial frequency directions are coarsely sampled every few seconds, allows for the possibility of utilizing a ‘sliding window’ reconstruction • By using a temporal aperture that widens with increasing spatial frequencies, temporal resolution with fast, flexible frame rates can be obtained without any loss in resolution or volume coverage

26. VIPR • 3D spherical coordinate system • Readout direction defined by angle θ from the kz-axis and angle Φ from the ky-axis • Resolution in all three dimensions determined by the maximum k-space radius value (kmax) • Diameter of the full field of view is determined • by the radial sample spacing (Δkr)

27. VIPR

28. VIPR • Stenotic carotid phantom • Separate random noise from undersampled artifact • equal number of excitations and scan time 1,500 projections with 16 signal averages 6,000 projections with 4 signal averages 1,5000 projections with 16 signal averages 12,000 projections with 2 signal averages 24,000 projections with no signal averages

29. VIPR • Plot of relative noise due to increased aliased energy from undersampling • Fewer than 15% of number of projections required for fully sampled k-space leads to significant degradation of image quality • Two carotid phantoms separated by 10 cm • Two fully sampled datasets were subtracted to determine noise not due to undersampling artifact • Noise ratio is noise measured in undersampled datasets relative to noise from subtracted fully sampled datasets

30. VIPR – Time Resolved • Coronal MIPs of time resolved images volumes • Times are time after contrast injection

31. VIPR – Time Resolved • Effect of increasing the width of the top of the temporal filter • Maximum width of the temporal aperture varied from 4 sec to 40 sec • CNR calculated between artery and background

32. VIPR – Isotropic Resolution • Arterial image volume shown 10 seconds after injection • Isotropic resolution allows for reformat in any plane

33. VIPR • Advantages: • No penalty in terms of coverage and spatial resolution • Resolution limits are constrained by noise and artifact • Isotropic resolution provides advantageous in terms of reformatted images • Oversampling center of k-space so more robust to motion • Improved coverage • Limitations: • Artifacts from undersampling decreases CNR

35. HYPR HighlY constrained back Projection Data are acquired in time frames comprised of interleaved and equally spaced k-space projections All of the acquired data is combined and used to form a composite image (either through filtered back projection or gridding and inverse fourier transform) Each projection in a time frame is inverse fourier transformed and back projected in image space Each back projected time point is weighted by the respective projection of the composite image The contributions from all projections in a time frame are summed

36. HYPR

37. HYPR

38. HYPR

39. HYPR

40. In both the polar and Cartesian cases, we have replicates at multiples of 1/Δkr • For Cartesian sampling, the replicates are exact duplicates of the main lobe • For Polar sampling the replications differ from the main lobe in shape and amplitude • In Cartesian sampling, exactly M oscillation periods occur between successive replicates • In polar sampling approximately M oscillation periods occur between successive replicates

41. Explain what artifacts look like in an image

42. Differences between polar and Cartesian undersampling Polar sampling generally produces streaking and radial aliasing arifacts These artifacts appear either inside or outside the circular field of view depending on the size of the object Even for optimal polar sampling, streaks and rings can be observed outside the displayed window In Cartesian sampling, the Dirichlet-shaped PSF of Cartesian sampling just generates an infinite amount of identical copies of the object separated by the dimensions of the rectangular field of view

43. In polar sampling the discretization occurs in both the radial and azimuthal directions We can focus on the radial effects of polar sampling because we effectively have a circularly symmetric ring sampling distribution The Jacobian is the given area that a discrete sample subtends In Cartesian sampling, each sample subtends the same constant area (ΔxΔy), however, in equidistant polar sampling, the sample size changes as a function of radius To normalize the area for each sampling in polar sampling, each sample is weighted by its areal extent which is equivalent to correcting for the sampling density The Fourier-Bessel (Hankel) transform is for circularly symmetric objects

44. Main advantages of Projection Acquisition: Shorter echo times Reduced sensitivity to motion Improved temporal resolution for certain applications Bernstein, pg 898

45. Trajectory: 3D spherical coordinate system Readout direction defined by angle θ from the kz-axis and angle Φ from the ky-axis Resolution in all three dimensions determined by the maximum k-space radius value (kmax) Diameter of the full field of view is determined by the radial sample spacing (Δkr)

46. PA motion artifacts consist of streaks propagating perpendicular to the direction of object motion and are displaced from the image of the moving object itself

47. Aliasing in Projection Acquisitions: Since the first lobe ring occurs at 1/Δkr and the object smears circularly, the inner edge of the smeared replicate occurs at (1/Δkr – rmax) where rmax is the maximal radial extent of the object Thus if rmax is greater than 0.5/Δkr the the replication will overlap the original object The effective field of view (unaliased region) in polar sampling is a circle of radius 0.5/Δkr Object size 0.375/Δkr, 3*1/Δkr shown Object size 0.554/Δkr

48. In both the polar and Cartesian cases, we have replicates at multiples of 1/Δkr For Cartesian sampling, the replicates are exact duplicates of the main lobe For Polar sampling the replications differ from the main lobe in shape and amplitude Oscillations occur due to the finite sampling extent and are thus often referred to as truncation artifact In Cartesian sampling, exactly M oscillation periods occur between successive replicates In polar sampling approximately M oscillation periods occur between successive replicates In the Cartesian case the complex exponentials are purely periodic while the Bessel function’s non-purely periodic nature leads to destructive interference on summation

49. In polar sampling, then the radial sampling requirement is the same as in Cartesian sampling, sampling must be above the Nyquist rate. In other words if the radial spacing in k-space is Δkr, the object must be limited to within a circle of diameter 1/Δkr in the image domain to avoid polar aliasing

50. However, if we assume S(kr,kθ) is circularly symmetric then S(kr,kθ)  S(kr) so