Burst imaging simulation
Published:
An ultrafast excitation and acquisition technique.
Blog posts
2024
Simulatenous Mulit-slice reconstruction tutorial 5
Published:
The slice-leakage is an artefact where the information from one slice is unintentionally transmitted to another slice during reconstruction. Any information from slice A that appears in slice B at the end of the reconstruction is considered slice leakage. The current method involves using a Monte Carlo simulation to impose unique frequency modulations on each slice. After reconstructions, the slice leakage can be determined by quantifying the frequency modulations on each individual slice. For example, if we add a 4 Hz modulation on slice A and a 6 Hz modulation on slice B, any 6 Hz component found on slice A after the reconstruction could be used to indicate slice leakage.
NUFFT tutorial
Published:
Challenges arise when reconstructing from non-uniform sampling patterns in the Fourier domain. Performing a direct Fourier transform in such cases incurs a quadratic computation cost of O($N^2$), rendering it impractical for standard practices at high resolutions. This requires transforming non-Cartesian samples back into a uniformly spaced Cartesian grid, enabling the utilization of more efficient FFT algorithms. Typical regridding procedures involve convolving the acquired data with a predefined gridding kernel weighted by the pre-calculated density compensation function, followed by a resampling process onto the Cartesian grid. The ideal gridding kernel is a sinc function, as its Fourier transform results in a rectangular function. However, since the optimal convolution function is of infinite extent, implementing it directly is impractical. To address this, the convolution kernel is truncated and windowed. A popular choice for this purpose is the use of Kaiser-Bessel functions. These methodologies are commonly known as “gridding” and is a special case of non-uniform Fourier transforms (NUFFT).
bSSFP simulation tutorial 2
Published:
Here, I present various variations of bSSFP: standard bSSFP, fluctuating equilibrium, and alternating steady states. Each variation manipulates sequence parameters such as TR, flip angles, and RF phase to generate unique profiles optimized for specific tasks.
bSSFP simulation tutorial 1
Published:
In bSSFP, the entire spin evolution can be modeled mathematically through bloch equations, which can be broken down into three components: excitation E, relaxation R, and precession P. These three components are matrix operators representing the sequence events. Immediately after the RF pulse R, the spin undergoes T1 and T2 relaxation E and off-resonance precession P. The bSSFP sequence achieves the steady state by repetitively executing these operations, while ensuring that the gradients are refocused within each TR.
2023
Partial Fourier reconstruction tutorial 3
Published:
Partial Fourier techniques aim to recover missing data using only half of the k-space. This reconstruction leverages the conjugate symmetry property of k-space, although it is sensitive to phase errors.
Partial Fourier reconstruction tutorial 2
Published:
Partial Fourier techniques aim to recover missing data using only half of the k-space. This reconstruction leverages the conjugate symmetry property of k-space, although it is sensitive to phase errors.
Partial Fourier reconstruction tutorial 1
Published:
Partial Fourier techniques aim to recover missing data using only half of the k-space. This reconstruction leverages the conjugate symmetry property of k-space, although it is sensitive to phase errors.
Low Rank reconstruction tutorial
Published:
The low rank based method will recover the missing data by enforcing self-consistency among neighboring k-space points in Cartesian space when minimizing the rank of the structured Hankel matrix. The self-consistency refers to the annihilation relationship being satisfied for all locations in k-space. Directly solving rank problems is computationally challenging and falls under the category of NP-hard problems. To address this, the non-convex rank function is replaced with its convex relaxation, which is replacing the non-convex rank function with its convex approximation, known as the nuclear norm.
Simulatenous Mulit-slice reconstruction tutorial 4
Published:
The slice grappa tries to reduce the total artefact without placing constraints on the intraslice and interslice artifacts. So the artefacts can be arbitrarily large(like vectors pointing in opposite directions can cancel out each other, also called artefact cancellation). The split slice trades interslice artefact with intraslice artefact (lower interslice artefact and higher intraslice artefact) by forming a “correlation matrix” that contrains the interslice leakage to be minimum.
2022
Simulatenous Mulit-slice reconstruction tutorial 3
Published:
What if we combine SMS and inplane undersampling ?
Simulatenous Mulit-slice reconstruction tutorial 2
Published:
Slice GRAPPA is used exclusively for SMS, ensuring that the output dimensions match the input dimensions. The kernel is trained for each slice, allowing it to project the SMS data onto each individual slice. Similarly, the slice GRAPPA kernel is obtained through a least squares fit.
Simulatenous Mulit-slice reconstruction tutorial 1
Published:
Simultaneous Multislice (SMS) is a data acceleration technique that acquires multiple slices at the same time, rather than skipping lines within an acquisition. This method can be visualized as capturing all images within an extended field of view and then undersampling with a slice factor in a simulation. In practice, the signals from all slices are acquired together and superimposed linearly, causing multiple slices to overlap in the final image. The objective of the reconstruction process is to separate these slices with the highest possible quality.
EPI distortion correction tutorial 3
Published:
Basic information about EPI has been shared in the previous post. Here, I share a reconstruction pipeline to correct the distortion using a third party package called FSL.
EPI distortion correction tutorial 2
Published:
Basic information about EPI has been shared in the previous post. Here, I share another straightforward method for correcting EPI field distortions, utilizing nonrigid image registration.
EPI distortion correction tutorial 1
Published:
Echo-planar imaging (EPI) offers rapid imaging capabilities by capturing an entire k-space data set in a single acquisition. However, this technique is prone to geometric distortions, which significantly degrade image quality. These distortions are primarily caused by field inhomogeneities, leading to voxel shifts. The voxel shifts are especially pronounced in the phase-encode direction. In this context, we explore the blip up and blip down approaches. I present a straightforward correction method based on line integrals to address these issues.
EPI PEC-SENSE phase correction tutorial
Published:
Due to the alternating direction of the EPI readout lines, hardware imperfections such as timing delays, eddy currents, and gradient coil heating can cause misalignment of the forward and reverse lines in k-space. This misalignment manifests itself in the images as a Nyquist (N/2) ghost in the phase encode direction and sinusoidal modulation of the object in the frequency encode direction.
EPI referenceless phase correction tutorial
Published:
Due to the alternating direction of the EPI readout lines, hardware imperfections such as timing delays, eddy currents, and gradient coil heating can cause misalignment of the forward and reverse lines in k-space. This misalignment manifests itself in the images as a Nyquist (N/2) ghost in the phase encode direction and sinusoidal modulation of the object in the frequency encode direction.
Joint GRAPPA reconstruction tutorial
Published:
bSSFP suffers from unique banding artifacts, which are commonly corrected by using multiple phase cycles. These banding artifacts are a form of spatial modulation. While it is tricky to exploit this with SENSE (you will need coil maps and bssfp profile maps), this additional redundancy can be easily utilized by GRAPPA. To achieve this, we simply extend the GRAPPA kernel to include the phase cycle dimensions.
CG SENSE reconstruction tutorial
Published:
I really enjoy the CG SENSE framework. It enables everything to be formulated as a forward model in the form of a least squares problem. This also allows for the application of various regularization techniques.
SENSE reconstruction tutorial
Published:
SENSE is a parallel imaging technique used to accelerate data acquisition in MRI.
GRAPPA reconstruction tutorial
Published:
GRAPPA is a widely used parallel imaging technique designed for accelerating data acquisition in MRI. It takes advantage of two key facts: first, MRI data is collected from multiple receivers, and second, k-space can be linearly interpolated due to its shift invariance properties.