1. Introduction
Changes in blood vessel pulsation and cerebrospinal fluid dynamics cause cyclic deformation of the brain [
1,
2,
3,
4,
5,
6]. Such brain motion can be altered by neurological pathologies such as hydrocephalus [
7,
8,
9], Chiari I malformation [
10,
11,
12], idiopathic intracranial hypertension [
13,
14], and age-related diseases in small cerebral vessels [
15,
16,
17]. Furthermore, the glymphatic system [
18], responsible for the macroscopic waste clearance of the central nervous system, is thought to depend on CSF pulsation in the perivascular pathways to remove metabolic waste from the brain [
19,
20]. Quantifying the cardiac- and CSF-induced (pulsatile) brain motion
in vivo may therefore provide valuable insights into neurological pathologies and their effects on the brain’s biomechanical response, making it a useful tool for research, and it has the potential to assist in the diagnosis and management of such diseases.
Various MRI techniques are available that enable the visualization and quantification of pulsatile brain motion. Phase contrast MRI (PC-MRI) [
3,
21,
22,
23] uses bipolar gradients to measure tissue velocity over the cardiac cycle in different directions, which can lead to longer scanning times. To estimate the displacement, the measured velocities need to be integrated over time, which can introduce additional errors. Complementary Spatial Modulation of Magnetization (CSPAMM) [
3,
24,
25] measures tissue displacement by tagging or nulling the signal in the region of interest. This technique has demonstrated an ability to measure brain motion in the cranial–caudal direction [
25]. However, due to this technique’s low spatial resolution, it may pose difficulty in detecting small brain motion in 3D within a clinically feasible scan time. Displacement encoding with stimulated echoes (DENSE) MRI [
5,
6,
26,
27,
28] is a technique that can measure tissue displacement with high spatial and temporal resolution and is sensitive to small displacements (0.01 mm). However, it requires encoding in different directions and suffers from a low SNR, both of which translate to long scan times.
Amplified MRI (aMRI) [
29,
30] is a relatively new method that can visualize pulsatile brain tissue motion over the cardiac cycle by amplifying subtle temporal intensity changes in cardiac gated ‘cine’ MRI data. aMRI offers several advantages over other techniques, including pronounced tissue contrast (when coupled with the balanced steady-state free precession (bSSFP) MRI sequence), a short scan time (it does not require displacement encoding in multiple directions), and a higher spatial and temporal resolution compared to DENSE-MRI.
Additionally, the ability of aMRI to reveal mesoscopic brain deformations makes it a promising tool for assessing various neurological disorders, including Chiari Malformation [
30,
31], aneurysms (aFlow [
32,
33]), concussion [
34], and acute hemorrhagic stroke [
35]. Recently, aMRI was extended to three dimensions, enabling the visualization of brain motion in all three planes (3D aMRI [
36,
37]). Compared to 2D acquisition, 3D aMRI enhances image quality and improves the visualization of pulsatile brain motion by capturing out-of-plane motion (third direction). The high spatial and temporal resolution of the underlying cine acquisition, along with the conspicuous brain tissue and CSF contrast enabled by the underlying bSSFP sequence, allows for the observation of the brain’s biomechanical response in exquisite detail.
However, 3D aMRI still lacks the ability to quantify the sub-voxel motion field in physical units. Our study introduces a novel 3D quantitative aMRI (3D q-aMRI) post-processing algorithm that enables the visualization and quantification of pulsatile sub-voxel brain motion by applying it to 3D bSSFP cine data. The new 3D q-aMRI is an extension of the original 3D aMRI algorithm and is based on solving the optical flow equation over the coefficients of the 3D complex linear steerable pyramid decomposition [
38,
39,
40].
To validate the developed algorithm, we tested 3D q-aMRI on a 3D digital phantom to assess its behavior and limitations. in vivo validation was performed to fine-tune the algorithm hyperparameters and to test the ability of the algorithm to quantify sub-voxel displacement motion observed in 3D aMRI. The results from 3D q-aMRI were qualitatively compared to phase-contrast MRI (PC-MRI), the conventional technique for capturing brain motion. The repeatability and reproducibility of 3D q-aMRI were evaluated in six healthy volunteers using a scan–rescan approach. Finally, preliminary 3D q-aMRI data were acquired on nine participants: four healthy controls and five with confirmed neurodegenerative disorders. Our aim was to test the potential diagnostic value of 3D q-aMRI in identifying disease-induced biomechanical differences.
2. Materials and Methods
2.1. Human Subjects
Experiments were conducted under ethical approval from The New Zealand Health and Disability Ethics Committee and Stanford University. Thirteen healthy adult volunteers were scanned (eight males and five females between 21 and 37 years old), with eight repeated scans acquired on six of these volunteers within the same session. Data were also acquired from nine subjects: four healthy controls (aged 55–75 years; one male and three females) and five individuals with confirmed neurodegenerative disorders. These included two patients with mild cognitive impairment (MCI) due to dementia (a 70-year-old female and a 74-year-old male), two patients with dementia due to Alzheimer’s disease (AD) (aged 62 years and 85 years), and one 76-year-old male with Lewy body dementia (LBD).
2.2. MRI Acquisition
A 3T MRI scanner (SIGNA Premier; GE Healthcare, Milwaukee, WI, USA) was used with an AIR™ 48-channel head coil. Three-dimensional volumetric cardiac-gated cine (bSSFP/FIESTA) MRI datasets were acquired as follows: Sagittal plane, FOV = 24 × 24 cm2, matrix size = 200 × 200, TR/TE/flip-angle = 2.8 ms/1.4 ms/25°, acceleration factor = 8, resolution = 1.2 mm isotropic, peripheral pulse gating with retrospective binning to 20 cardiac phases, 120 slices for whole brain coverage, and a scan time of 2:30 min.
For segmentation purposes, we also collected T1 MPRAGE volumetric data from the six volunteers who were scanned for repeatability and reproducibility analyses. The following parameters were used: Sagittal plane, FOV = 24 × 24 cm2, matrix size = 200 × 200, TR/TE/flip-angle = 2500 ms/2.4 ms/8°, resolution = 1.2 mm isotropic.
For comparison to past studies using the common phase contrast MRI (PC-MRI) approach, we also collected a cine PC-MRI dataset on two volunteers, with the following parameters: FOV = 24 × 24 cm2, velocity encoding = 0.5 cm/s, matrix size = 128 × 128 (zero-filling to 200 × 200), slice thickness = 4 mm, and TR/TE/flip-angle = 27 ms/14 ms/10°. Velocity encoding was performed in the right–left direction in the axial plane, and in the superior–inferior direction for the sagittal and coronal planes.
PET/MR scans were performed using a simultaneous time-of-flight (TOF)-enabled PET/MR scanner (3T SIGNA Premier; GE Healthcare, Milwaukee, WI), with high sensitivity (23.3 cps/kBq). MRI data were collected simultaneously with the PET data from both healthy controls and participants with dementia using the same aMRI parameters as previously mentioned.
2.3. Motion Estimation
3D q-aMRI is a natural extension of the original 3D aMRI algorithm (
Figure 1). As such, the algorithm starts by decomposing the brain volumetric ‘cine’ MRI data using the 3D steerable pyramid decomposition [
36,
37]. The 3D steerable pyramid is a complex linear multiscale and multi-orientation decomposition in which steerable filters are used in a recursive scheme (low-pass filtering and downsampling) [
41,
42]. The scales basis functions are bandpass filters in the frequency domain. They are calculated in polar coordinates as a multiplication of low- and high-pass filters. The low-pass and high-pass filters for each scale are given by the following equations:
where
and
s represents the level scale. The bandpass filter of each level is given by the following:
The angular (orientation) filters are 3D cones oriented along the six vertices of cuboctahedron (
Table 1) and satisfy the following equation in the frequency domain:
where
are the directions of the axes of symmetry of the six basis filters
. The 3D steerable filters are constructed as follows:
and the final image response is calculated as follows:
where
s and
j are the scaling factor and orientation direction respectively;
is the Fourier transform of the image; and
is the inverse Fourier transform. The decomposition outputs a complex number with amplitude (
) and phase (
) at each scale and orientation.
The phases of the decomposition encode information about the sub-voxel motion, similar to the Fourier shift theorem. This information is exploited to achieve both magnification (visualization) and quantification of the voxel displacement field. The phases are separated from the amplitude component and temporarily bandpass filtered to isolate the desired frequency range, and to remove any DC component.
The bandpass phases proceed through two distinct processing pathways: (a) the original 3D aMRI algorithm [
36], which includes phase denoising using an amplitude-weighted Gaussian smoothing filter, multiplication by the amplification factor
, and reconstruction that results in motion magnification, and (b) the new voxel displacement quantification algorithm. The quantification algorithm is based on solving the optical flow equation, which assumes the brightness constancy constraint (that the brightness of a moving point in an image remains constant over time). However, there is a subtle distinction from the standard optical flow approach that we use here. The equation is instead solved over the filtered phases of the coefficients of the 3D complex-valued steerable pyramid decomposition. Here, we are guided by the work from Fleet and Jepson [
38], who showed the general case whereby the temporal evolution of (spatial) contours with constant phase in the image subbands (such as used in a complex wavelet decomposition) correspond to the motion field. Mathematically, this can be written as follows:
Here,
are the phases of the 3D complex steerable pyramid, and
are the voxel displacement fields we would like to estimate. We can expand the right-hand side using the Taylor series around
to obtain the following:
where
and
are higher-order terms in the Taylor expansion. The main assumption here is that the motion is at the sub-voxel level, and as a result, the higher-order terms in the Taylor expansion are negligible. Therefore, the solution can be approximated using only the linear term:
We solve this linear equation using a weighted least squares objective function over the various orientations and scales of the 3D steerable pyramid decomposition [
39] as follows:
where
W is a Gaussian window (which assumes that the motion field is smooth around a given voxel);
and
are the phase and amplitude of the output of the filters in the 3D steerable pyramid decomposition; and
is the phase difference between the frames
and
. Note that the extent of the Gaussian window is calculated according to the standard deviation
value using the following equation:
and that throughout this work, we use a Gaussian window with
, which equates to a support window of
pixels. In addition, we notice that phases with orientations corresponding to edges—which are more reliable for motion estimation—will have a larger weighting term
.
Furthermore, conventional discrete derivative operators (e.g., the intermediate or the central kernels) are sensitive to noise when estimating derivatives in high spatial frequency images, such as in the filter responses of the higher levels of the 3D steerable pyramid. Therefore, to increase the algorithm’s robustness to noise and improve the estimated motion field accuracy, the phases’ spatial derivatives
were estimated using the following approach [
38]:
where
is the image response for the
filter (complex number) and
is the imaginary operator.
still contains components with high spatial frequencies, but it can be demodulated as follows, with the image response expressed as:
where
is the peak tuning frequency of the corresponding subband’s filter. Taking the derivative of
results in the following:
where
is the demodulated filter response, which contains only low spatial frequencies. Finally, we notice that traditional derivative filters are often inappropriate for multi-dimensional problems [
43]; therefore,
was approximated using the following derivative filters:
where the
Derivative filter is applied in the desired direction, and the
PreFilter is applied in the two remaining directions. In comparison to the traditional intensity-based optical flow algorithms, which estimate motion larger than one voxel, this approach of using image local phase information can estimate sub-voxel motion on the order of 0.01 pixel size [
39,
44].
2.4. Digital Phantom Simulation
To evaluate the performance of 3D q-aMRI, a 3D digital phantom was used to mimic the subtle deformation, intensity, and contrast of the lateral ventricles observed in 3D aMRI (
Supplemental Video S1). The phantom consists of a 3D cylinder with an initial height
and an initial radius
and was simulated using MATLAB (MathWorks, Natick, MA, version 2022b). The two ends of the 3D cylinder were deformed according to the equation
, resulting in cyclic compression and tension along its axis (
Figure 2a). We made the following assumptions:
The phantom volume remains constant over time, resulting in expansion and compression along the cylinder’s radius.
The radius follows a parabolic curve, , where is the maximum radius and z is the position along the axis of the cylinder of height . These assumptions result in the following equation for a and :
Here
, which was selected as in our previous 3D aMRI work [
37]. This choice of
ensures that the phantom undergoes sub-voxel motions.
The estimated motion field was evaluated against the true motion field using the following metrics:
- (1)
Pearson’s Linear Correlation
- (2)
- (3)
Ninety-ninth percentile of the error distribution
The following were evaluated: (1) the algorithm’s accuracy as a function of the magnitude of the true motion (
Figure 2b); (2) the algorithm’s robustness to various Gaussian noise levels (
Table 2); and (3) the effects of the algorithm’s hyperparameters (number of pyramid levels in the least squares, temporal filtering, and Gaussian window
parameter) and various temporal frequency motions on the estimated motion field (
Supplemental Figures S1–S4).
2.5. In Vivo Validation
First, it is important to note that we undertook the following standardization procedures:
- (1)
Since the native image size of the acquisitions determines the extent and resolution of the 3D steerable pyramid filters, we ensured that all in vivo datasets were zero-padded to a uniform dimension of 256 × 256 × 256.
- (2)
Accordingly, for the analysis of different image resolutions, we used the corresponding values to match the physical extent of the Gaussian smoothing filter. For example, to compare a 2.4 mm isotropic resolution (128 × 128 × 128) versus an isotropic 1.2 mm isotropic resolution (256 × 256 × 256) dataset (the latter with ), we first zero-padded the 2.4 mm dataset and used to ensure the filter has the same extent as the 1.2 mm dataset in physical units.
- (3)
To generate a sharp frequency response of the temporal filter that will pass the desired frequency band, the number of cardiac phases was extended to 5 periods—which is a valid procedure since the motion of the cine data is periodic. This was performed since the native number of cardiac phases of 20 may generate side lobes in the filter frequency spectrum—which in turn may propagate errors in the voxel displacement estimation due to noise.
To validate the ability of 3D q-aMRI to accurately quantify the observed motion in 3D aMRI, we followed three steps. First, 3D aMRI was used to amplify the signal in the 4D
in vivo data. Then, the voxel displacement field was calculated using 3D q-aMRI, and a second amplified video was created by warping the first volume in the data according to an amplified version of the estimated motion field. Next, normalized temporal standard deviation maps were calculated for both amplified videos (
Figure 3), and these were compared qualitatively.
Lastly, we also tested and optimized the algorithm’s hyperparameters (number of pyramid levels, temporal filtering, and Gaussian window parameter) for the in vivo data. Lastly, the effect of varying image resolution on the estimated voxel displacement field was tested.
2.6. In Vivo Analysis
An
in vivo analysis was conducted to provide insights into the performance and potential applications of 3D q-aMRI in biomedical research and clinical practice. First, we qualitatively compared 3D q-aMRI against ground truth phase contrast (PC)-MRI data. The velocities were integrated over time to extract the displacement field. The motion profile was also evaluated in the through-plane direction through the cerebral aqueduct. Second, the repeatability and reproducibility of 3D q-aMRI were assessed using 4 scans and 4 rescans (for a total of 8 scans) in 6 healthy volunteers. For each volunteer, the following analyses were conducted: initially, the MPRAGE data were registered to the first ‘cine’ bSSFP scan using rigid body registration in MATLAB (MathWorks, Natick, MA, USA, 2022). Second, the registered MPRAGE data were used to segment [
45] the following brain regions: lateral ventricles, the 3rd ventricle, the 4th ventricle, the brainstem, and the cerebellum. Next, the voxel displacement field was extracted for each scan using 3D q-aMRI. Each scan (of the eight) was registered to the first scan using rigid body registration, and the extracted transformation was applied to each of the voxel displacement fields. Lastly, the average motion profile in the superior and inferior directions throughout the cardiac cycle was extracted for each of the segmented brain regions. In addition, statistical analysis was conducted as follows: For each brain region and subject, the distance between the average signal (computed over the eight scans) and each individual measurement was computed using dynamic time warping. Dynamic Time Warping (DTW) is a technique used to measure the similarity between two time series by warping their time axes to achieve the best match between corresponding elements. This alignment minimizes the distance between the sequences. The calculated distances were used to calculate the Intraclass Correlation Coefficient (ICC). Finally, 3D q-aMRI data were acquired on patients with dementia and a healthy age-matched control to assess the potential diagnostic value in identifying disease-induced biomechanical differences. Note that this is an early preliminary test, suggesting that, regardless of the results of these scans, further research will be necessary to determine the technique’s potential as a diagnostic tool.
4. Discussion
This work introduces a novel 3D quantitative aMRI (3D q-aMRI) post-processing algorithm that enables the visualization and quantification of the sub-voxel pulsatile brain motion in physical units. The algorithm is based on solving the optical flow equation over the phase coefficients of the 3D complex steerable pyramid decomposition and as such is a natural extension to the original 3D aMRI algorithm. By extending the original 3D aMRI to 3D q-aMRI and applying it to volumetric in vivo data, this method offers the capability to visualize and quantify the 3D piston-like motion of the brain and change in ventricular shape as the brain deforms through the cardiac cycle.
The phantom simulations suggest that when properly tuned, the 3D q-aMRI is capable of quantifying motion at a scale greater than 0.005 of a pixel size (the theoretical limit) for an SNR as low as 25. This is achieved while maintaining a high correlation (0.95), low error (13%), and a relatively small error distribution tail (99% of the error values fall below 30%). In general, phase-based methods, such as those used in 3D q-aMRI, are relatively robust to noise, as demonstrated by phantom simulations. However, it is important to note that when the data contain high levels of noise, the algorithm is likely to produce significant errors. In addition, in practice, there are likely confounding factors such as field inhomogeneity, partial volume effects, and various sources of ‘noise’, which were not comprehensively represented in the phantom simulation. These factors are likely to limit 3D q-aMRI’s ability to accurately estimate motions as small as 0.005 pixels in size. Yet, as demonstrated in
Figure 8b, the motion profile observed in the cerebral aqueduct is continuous and smooth. This suggests that the algorithm has successfully detected coherent signal changes corresponding to voxel displacement at a scale of 0.01 pixels.
A noteworthy observation from the phantom simulation is that the algorithm performs very accurately, even for very low-contrast data, such as what is seen in the brain ventricles. As long as some level of texture exists in the data and the voxels exhibit subtle signal changes corresponding to motion, the algorithm appears capable of accurately estimating the voxel displacement field-even in the presence of considerable noise (low SNR). This claim comes with the condition that in a small neighborhood, the motion field remains relatively constant-as this factor influences the algorithm’s accuracy. Thankfully, this assumption seems valid for the human brain but will require further evaluation in cases that may present extreme brain motion patterns. It is worth mentioning that some texture almost always exists in the in vivo data. For example, the lateral ventricles, which are supposed to contain only CSF (same T1 and T2 property, and as such the same contract) contain some texture. Finally, the results suggest that 3D q-aMRI can adequately quantify motion across a range of temporal frequencies.
The study’s
in vivo validation demonstrates that 3D q-aMRI successfully quantifies (in the form of a voxel displacement field) the observed signal in 3D aMRI. It is important to distinguish between two types of quantification: (1) quantification of the 3D aMRI unamplified signal and (2) comparison of 3D q-aMRI with true quantitative methods, such as PC-MRI. This validation procedure supports the first type of quantification. The algorithm hyperparameters were tuned empirically according to the phantom simulation observations. It is important to note that brain motion induced by other physiological signals, such as respiration and vasomotion, may require different hyperparameter tuning. However, for cardiac-induced motion, we observed that most of the motion occurs in the first two levels of the 3D steerable pyramid. This suggests that the pulsatile brain motion is composed of local deformations rather than global translation (rigid body motion). In addition, most of the motion is concentrated within the one to four heart rate harmonics band. This suggests that the cardiac input, which is periodic, contains higher harmonics in addition to the main heart rate (fundamental harmonic) and/or that the brain, as a nonlinear system, is excited by higher harmonic modes in response to the cardiac input. In addition, we noticed that for a Gaussian filter with
, the estimated motion field profile is almost unchanged. This observation helps to determine a threshold for how much smoothing is needed to remove noise without oversmoothing the estimated field. This observation is similar to our previous work on 3D aMRI [
36,
37], where we visually observed that for a Gaussian filter with
, the amplified data contain minimal noise, and the amplitude of the brain motion is smaller. Finally, we determined the range of image isotropic resolution (1.2–1.8 mm) for which the estimated motion field remains relatively consistent. Larger motion displacement was observed in low-resolution data. This is likely to be due predominantly to errors in low-resolution images, which are plagued by partial volume effects. Our analysis here emphasizes the importance of standardizing 3D q-aMRI for clinical applications, including various parameters, such as image resolution and dimension, Gaussian filter (
value), and temporal filtering frequency band, which have an impact on the extracted brain motion field (mainly the amplitude). Notice that other true quantitative methods such as DENSE-MRI are also subject to variation in tissue displacement estimation as a function of image resolution [
47].
It is worth mentioning that, like many other MRI techniques, 3D q-aMRI is prone to rigid body motions, which can be categorized into coherent and non-coherent motions. Noncoherent motion results in spatial blurring in the ‘cine’ MRI data, which in turn results in an error in the estimated voxel displacement field. Coherent motion-such as bulk head motion induced by blood pulsation-is detectable by 3D q-aMRI and can be corrected by subtracting the estimated voxel displacement field in the skull (which is assumed to be static) from the estimated displacement field.
There is an important limitation of 3D q-aMRI that must be highlighted. The algorithm is unable to estimate motion accurately beyond 1.5 pixels. This is primarily because small motions are assumed for the linearization of the phase constancy equation (Equation (
8)). When dealing with heavily pulsating vessels like the basilar artery, which can move through several pixels in a cine dataset, these dynamic structures can introduce errors in the estimated motion field. This underscores the importance of thoroughly examining both the raw and amplified data in conjunction with the displacement field maps.
Previous research has illustrated the intricate nature of the brain’s biomechanics, showing the interplay between the CSF and the pulsatile brain motion. The general characteristics of pulsatile brain motion and the estimated voxel displacement field extracted by 3D q-aMRI resemble the displacement field extracted by PC-MRI and other studies [
5,
6]. The predominant displacement of tissue and fluid occurs in the cranial–caudal direction in the sagittal and coronal planes. Additionally, the axial plane shows expansion/contraction motion, and the most significant tissue displacement is concentrated around the midbrain, brainstem, cerebellar tonsil, and hypothalamus regions. Due to the reduced perfusion of brain tissue toward the cranium boundary, minimal displacement is observed in the frontal lobe, parietal lobe, occipital lobe, temporal lobe, and posterior cerebellum.
It is noteworthy that our choice of the 3D q-aMRI hyperparameters resulted in a motion/flow profile through the cerebral aqueduct, which matched the profile previously reported by [
46], which used PC-MRI. While this is a promising finding that supports the current choice of the hyperparameters used here, the 3D q-aMRI flow profile was an order of magnitude (mm/s) smaller than that reported in the literature [
46,
48]. This is an interesting finding. Although 3D q-aMRI does not directly measure CSF flow like PC-MRI, it may provide complementary information about tissue compliance alongside flow measurements. It seems like the cerebral aqueduct pulsatile motion may result from the pressure of the CSF flowing through it- and that regional brain displacement and CSF flow may be tightly coupled. This observation suggests that 3D q-aMRI may be capable of detecting the CSF pulsatile flow profile in various brain compartments, such as the subarachnoid and perivascular space. Again, while the flow profile itself is not quantifiable, it could lead to an important method of indirectly assessing CSF homeostasis, and further testing should be performed to confirm this.
The in vivo validation of our study underscores the potential of 3D q-aMRI to serve as a reliable quantitative tool in clinical studies, as it can generate ’heat maps’ of pulsatile brain motion, which may assist in assessing disease-induced biomechanical differences. Our study observed that the most substantial brain motion occurs in the superior–inferior direction, which is consistent with previous studies. This motion is likely attributed to the pulsatile nature of the CSF flow, which is strongly influenced by the cardiac cycle.
Additionally, our results demonstrate that 3D q-aMRI is both repeatable and reproducible. The intra-correlation coefficient (ICC) analysis suggests significant differences between within-subject and across-subject measurements, indicating the potential to detect individual differences even within this healthy volunteer cohort. However, in some repeat scans, we did observe some variation in the amplitude of the voxel displacement profile (
Supplemental Figure S6). Upon examining these ‘cine’ data, we identified global temporal flickering (
Supplemental Video S9). The source of the flickering is still under investigation, but we suspect it may be related to the high acceleration factor (eight in our case) and the reconstruction scheme. So far, 3D q-aMRI has been applied to product sequences, such as the GE ’cine’ FIESTA. Our future work aims to jointly optimize the acquisition and reconstruction of 3D q-aMRI to further enhance its repeatability. It is also worth noting that heart rate exerts a significant influence on the estimated motion magnitude. We believe that this is a direct physiological consequence of cardiac pulsatility. Our data is cardiac-gated, and motion is assumed to result from blood flow and pulsations transferred through the larger arteries within a cardiac cycle. Hence, any alterations in the amplitude or frequency of these pulsations-such as natural heart rate variations-are likely to impact the timescale available for transferring these pulsations carried by the pulsatile blood flow. Monitoring the effects of heart rate variations during the scan and understanding their impact on data interpretation are topics that will be explored in future studies.
It is also important to recognize that while 3D q-aMRI demonstrated reliable detection of sub-voxel brain motion, other MRI techniques (CSPAMM, PC-MRI, and DENSE-MRI) are also capable of detecting this type of motion. With its high temporal resolution, high spatial resolution, good tissue contrast, and short scan time, 3D q-aMRI is a practical clinical alternative method. For example, the acquisition time for the PC-MRI data in this study was 4.5 min (for only three slices, with a single encoding direction, and with a much lower resolution of 1.9 × 1.9 × 4 mm), and currently, the total acquisition time for a DENSE-MRI scan in a 3T scanner with 2.2 mm isotropic resolution is 16.5 min [
5]. In addition, with a longer scanning time (15 min), 3D q-aMRI can achieve 0.8 mm isotropic spatial resolution without compromising temporal resolution and SNR. For comparison, the acquisition time for DENSE to achieve this resolution will be more than twenty times longer. This advantage might enable the detection of pulsatile motion of small features such as the perivascular spaces.
However, it is important to highlight that 3D q-aMRI estimates the voxel displacement field based on the signal intensity temporal changes in bSSFP ‘cine’ MRI. While the measured signal in DENSE-MRI and PC-MRI corresponds to the genuine spin displacement, in bSSFP, there are other physiological sources that may impact the signal intensity and cause 3D q-aMRI to interpret it as motion (like any other optical flow method)-even if they do not necessarily correlate to true displacement. As such, the method should not be considered as a true displacement quantification method like DENSE-MRI. There is no doubt that the extracted voxel displacement field is highly correlated and resembles physical units with tissue displacement, as seen in our PC-MRI analysis and [
46,
48,
49], but further work is also necessary to characterize the true physiology behind the 3D q-aMRI signal variation.
Preliminary data from our study indicate that 3D q-aMRI has the potential for identifying biomechanical differences induced by disease and could serve as a diagnostic tool for detecting early signs of neurodegenerative disorders. Our preliminary study suggests that abnormal pulsatile brain motion might correlate with changes associated with cognitive decline; however, it is important to note that these findings are based on a small number of healthy volunteers, and further studies with larger and more diverse populations are necessary to validate these results and establish the clinical utility of 3D q-aMRI.