
==== Front
Neuroradiology
Neuroradiology
Neuroradiology
0028-3940
1432-1920
Springer Berlin Heidelberg Berlin/Heidelberg

38236423
3282
10.1007/s00234-024-03282-6
Advanced Neuroimaging
High-angular resolution diffusion imaging generation using 3d u-net
Suzuki Yuichi 1
Ueyama Tsuyoshi 1
Sakata Kentarou 1
Kasahara Akihiro 1
Iwanaga Hideyuki 1
http://orcid.org/0000-0002-0324-6562
Yasaka Koichiro koyasaka@gmail.com

2
Abe Osamu 12
1 grid.412708.8 0000 0004 1764 7572 Radiology Center, The University of Tokyo Hospital, Tokyo, Japan
2 grid.412708.8 0000 0004 1764 7572 Department of Radiology, The University of Tokyo Hospital, Tokyo, Japan
18 1 2024
18 1 2024
2024
66 3 371387
25 10 2023
28 12 2023
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Purpose

To investigate the effects on tractography of artificial intelligence-based prediction of motion-probing gradients (MPGs) in diffusion-weighted imaging (DWI).

Methods

The 251 participants in this study were patients with brain tumors or epileptic seizures who underwent MRI to depict tractography. DWI was performed with 64 MPG directions and b = 0 s/mm2 images. The dataset was divided into a training set of 191 (mean age 45.7 [± 19.1] years), a validation set of 30 (mean age 41.6 [± 19.1] years), and a test set of 30 (mean age 49.6 [± 18.3] years) patients. Supervised training of a convolutional neural network was performed using b = 0 images and the first 32 axes of MPG images as the input data and the second 32 axes as the reference data. The trained model was applied to the test data, and tractography was performed using (a) input data only; (b) input plus prediction data; and (c) b = 0 images and the 64 MPG data (as a reference).

Results

In Q-ball imaging tractography, the average dice similarity coefficient (DSC) of the input plus prediction data was 0.715 (± 0.064), which was significantly higher than that of the input data alone (0.697 [± 0.070]) (p < 0.05). In generalized q-sampling imaging tractography, the average DSC of the input plus prediction data was 0.769 (± 0.091), which was also significantly higher than that of the input data alone (0.738 [± 0.118]) (p < 0.01).

Conclusion

Diffusion tractography is improved by adding predicted MPG images generated by an artificial intelligence model.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00234-024-03282-6.

Keywords

Diffusion-weighted imaging
High angular resolution diffusion imaging
Artificial intelligence
Orientation distribution function
Tractography
http://dx.doi.org/10.13039/501100001691 Japan Society for the Promotion of Science 20K08016 Suzuki Yuichi The University of TokyoOpen Access funding provided by The University of Tokyo.

issue-copyright-statement© Springer-Verlag GmbH Germany, part of Springer Nature 2024
==== Body
pmcIntroduction

Diffusion-weighted imaging (DWI) is a unique type of MRI able to non-invasively obtain in vivo information. It is known for its excellent ability to detect acute cerebral infarction, especially when combined with apparent diffusion coefficients (ADCs), which can be acquired with DWI. Thus, it has become an indispensable sequence in clinical practice. Diffusion tensor imaging (DTI) was devised as an application of DWI that uses fractional anisotropy (FA) to obtain typical quantitative values [1]. It is used in various clinical and research settings. Diffusion tensor tractography (DTT) depicts bundles of white matter in the brain by tracing eigenvectors. It is used in various clinical situations, such as preoperative planning for brain tumors, postoperative functional evaluations, brain development evaluations, and tracking the course of degenerative diseases [2–6]. However, DTI and DTT approximate diffusion information within voxels to ellipsoids, making the expression of cross-fibers difficult. High-angular resolution diffusion imaging (HARDI) addresses this by representing cross-fiber regions [7]. Based on this approach, other methods have been developed to increase the amount of angular information using motion-probing gradients (MPG). These include Q-ball imaging (QBI) [8], constrained spherical deconvolution (CSD) [9], and generalized q-sampling imaging (GQI) [10]. Using these HARDI-based techniques has enabled the identification of complex white matter pathways at the junction of the auditory radiation and inferior longitudinal fasciculus [11]. QBI tractography has been found more reliable than DTT in identifying the motor language cortex [12] and the brain origins of stroke. CSD tractography has been shown to have better visualization ability than DTI and its clinical usefulness is high [13].

In recent years, the application of machine learning and deep learning technologies to MRI has progressed. A variety of uses of artificial intelligence (AI) in MRI have been reported, including noise removal [14], brain tumor differentiation [15], contrast agent reduction [16], and cerebral aneurysm detection [17]. Some of these have begun to be used in clinical settings. Of course, the application of AI to DWI is also popular. Specifically, there have been reports on noise and artifact removal [18, 19], super-resolution [20], diffusion profile evaluation [21], tractography [22], and white matter extraction [23]. Using AI in this way, it has become possible not only to improve image quality but also to support diagnoses. However, in recent years, DWI analysis and imaging conditions have become more complex. There has been an increase in analytical methods that require a single b-value but many MPGs [8–10] or multiple b-values, such as diffusion kurtosis imaging [24], neurite orientation dispersion and density imaging [25]. Therefore, imaging time using DWI tends to be longer than that for DTI. It is possible to shorten imaging time using multistage simultaneous excitation technology such as multiband echo planar imaging (EPI) [26–28] and simultaneous multi-slice (SMS) [29–31] but there is a trade-off, with poorer image quality due to insufficient imaging data [32]. Therefore, the reduction of imaging time in DWI using these techniques is limited.

Until now, AI has primarily been used in this field to improve aspects of the acquired images, such as image quality, or to increase the precision with which diseases can be differentiated, as mentioned above. In this study, we took a completely different approach, with the aim of having AI generate predicted DW images. We investigated the extent to which it is possible to shorten scan time and ensure image quality (tractography accuracy) by generating images that would usually be acquired with MRI using AI. We evaluate HARDI data, including AI-generated images using several quantitative parameters, and examine the utility of this approach.

Materials and Methods

Participant data

Participants in this study were 251 cases that had already been imaged to depict QBI tractography as an in-hospital clinical examination before brain tumor resection, gamma knife surgery, or epileptic seizure surgery as retrospective data. These inspections were conducted from February 2016 to March 2021.

The analysis of these images analysis was performed using a protocol approved by the Ethics Committee of the University of Tokyo. All procedures were performed in accordance with the tenets of the 2013 revision of the Declaration of Helsinki. The informed consent requirement was waived by the Ethics Committee due to the retrospective nature of the study.

Imaging parameters and preprocessing

All MR images were acquired using a 3T clinical scanner (MAGNETOM Skyra, Siemens Healthcare, Erlangen, Germany; 45 mT/m max. gradient strength, 200 mT·m−1·s−1max. slew rate) and a commercial 20-channel matrix head coil. The DWI parameters for single-shot EPI were: repetition time/echo time = 8900/89 ms, in-plane acceleration factor (generalized auto-calibrating partially parallel acquisition factor) = 2, slice thickness/gap = 2.5/0 mm, number of slices = 60, field-of-view = 240 × 240 mm2, matrix = 96 × 96, voxel size = 2.5 × 2.5 × 2.5 mm3, b = 3000 s/mm2, multiple MPG = 64 directions and a single b = 0 s/mm2 image using phase-encoding directions along the anteroposterior axis, and scan time = 625 s. We also obtained b0 s/mm2 images with reversed phase-encoding directions along the posteroanterior axis for the distortion correction process. In preprocessing, we performed noise removal [33–35] and ringing artifact removal [36, 37] using MRtrix3 [38], distortion and motion correction using FSL top-up [39] and eddy [40], and then B1 field correction [41] using MRtrix3. And 3D T1-weighted image for anatomical image parameters were: repetition time/echo time = 1900/3.16 ms, slice thickness = 1.25 mm, number of slices = 128, field-of-view = 240 × 240 mm2, matrix = 192 × 192, voxel size = 1.25 × 1.25 × 1.25 mm3, Inversion time; 962ms, and scan time = 206 s.

Convolutional neural network

Our dataset was divided into training (n = 191), validation (n = 30), and test (n = 30) sets. We employed a 3D U-net learning model (Fig. 1) created using Python 3.9.7 (https://www.python.org/) and tensorflow-gpu 2.8.0 (https://www.tensorflow.org/) on a computer equipped with a Core i9-12900F central processing unit, 128 GB of random access memory, and a GeForce RTX 3090 graphic processing unit. We used the b0 image and the first 32 axes of HARDI data as input data. The remaining 32 axes of HARDI data were used as the teaching data, and the MPG32 axes output was used as the prediction data. Table 1 shows the MPG array using the electrostatic repulsion method [42] in this study.Fig. 1 Pipeline depiction of our artificial intelligence imaging generation method using 3D U-net.  MPG, motion-probing gradient 

Table 1 The MPG array using the electrostatic repulsion method  

MPG is motion-probing gradient and MPG no.0 indicates b0 image

At the time of learning, we calculated the loss function as the mean squared error (MSE) between the teaching and prediction data. We used Adam as the optimizer. The hyperparameters were optimized using the validation data. The hyperparameters for the 3D U-net model were: kernel size, 3; stride, 2; activation function, rectified linear unit. Those for the supervised learning process were: minibatch size, 15; number of epochs, 300. The trained model was applied to the test data to obtain the predicted data.

Image evaluations

First, the structural similarity (SSIM) of the prediction data was calculated for each MPG axis, with the latter half of the reference data as the gold standard. We calculated the SSIMs with the multiSSIM3 function of MATLAB 2020b. And we also calculated the peak signal to noise ratio (PSNR). Next, we calculated the orientation distribution function (ODF) and the similarity of the data for the first 32 MPG axes only (input data). The ODF of the reference data and the ODF of the input plus prediction data were compared. ODF calculations were performed using the Diffusion toolkit (https://trackvis.org/dtk/), and quantitative value comparisons were performed using Jensen-Shannon divergences (JSDs) and angular correlation coefficients (ACCs) [43]. JSDs were used to quantify the similarities between two fiber orientation distributions (FODs) or ODFs. We projected both ODFs onto 181 values distributed equally over a sphere. JSD was defined as:JSD(P,Q)=DKLP,M+DKL(Q,M)2

M(i)=Pi+Q(i)2

where P(i) and Q(i) are the magnitudes of the histological and MRI FODs (or ODFs) along index i (i = 1 …181), and DKL is the Kullback–Leibler divergence:DKL(P,Q)=∑iPilog(P(i)Q(i))

We calculated the correlation of functions over a sphere given the spherical harmonic (SH) expansions of both functions. Given two spherical functions and their SH expansions,U(θ,φ)=∑l=0∞∑m=-llulmYlm(θ,φ)

V(θ,φ)=∑l=0∞∑m=-llvlmYlm(θ,φ)

the ACC of the functions is calculated as:ACC=∑l=1∞∑m=-llulmv∗lm[∑l′=1∞∑m′=-l′l′|ul′m′|2]12∑l′′=1∞∑m′′=-l′′l′′|vl′′m′′|2|12

The calculation areas for the SSIM, JSD, and ACC were masked using the b0 image.

We performed diffusion tensor analysis and generated FA and color FA from the eigenvalues in our software. In addition, regions of interest (ROIs) were set at the cerebral peduncle (ROI1), posterior peduncle of the internal capsule (ROI2), corona radiata (ROI3), and semiovale center (ROI4) in the normal side on the color FA map. The subjects for ROI measurement were 24 patients, excluding 6 patients with bilateral disease, brain stem disease, and multiple metastases. And the average value and standard deviation (SD) for FA and eigenvalue1 (E1) were calculated. The ROIs were manually set with the consent of one radiological technologist (17 years of MRI experience) and one radiologist (Fig. 2).Fig. 2 ROI for measurement FA and E1

Red ellipsoids indicate ROIs in each axial section on the color FA map. ROI1 indicates the cerebral peduncle, ROI2 indicates posterior peduncle of the internal capsule, ROI3 indicates corona radiata and ROI4 indicates semiovale center in the normal side.

Finally, left and right corticospinal tract (CST) and frontal aslant tract (FAT) with Q-ball imaging tractography and generalized q-sampling imaging tractography were visualized, and dice similarity coefficients (DSCs) were calculated for the reference data.

QBI tractography was performed using TrackVis (https://trackvis.org/), and GQI tractography was performed using DSI Studio (https://dsi-studio.labsolver.org/). ROIs of CST were automatically obtained by TractSeg [23] analysis. The seed ROI was the cerebral peduncle, with the primary motor cortex as the target point. The midsagittal section was manually set as the avoidance ROI (Fig. 3).Fig. 3 ROIs for depicting CST

Yellow indicates the cerebral peduncle for seed point, green indicates the primary motor cortex for target point, and blue indicates the midsagittal section for avoidance ROI.

The ROIs for FAT was manually set with the consent of one radiological technologist (17 years of MRI experience) and one radiologist. T1-weighted image was registered to the b0 image using FMRIB's Linear Image Registration Tool; FLIRT (https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FLIRT). The seed point is (supplementary motor area; SMA and pre-SMA), target point is set in the gill cover of inferior frontal gyrus (IFG), the triangular part of IFG, and the orbital part of IFG [44] (Fig. 4).Fig. 4 ROIs for depicting FAT

Green area indicates the SMA and pre-SMA for seed point, blue area indicates the gill cover of inferior frontal gyrus (IFG), the triangular part of IFG, and the orbital part of IFG.

The default conditions in each piece of software were used for settings such as the tracking algorithms. After voxelization of the corticospinal tracts, we calculated the DSC [45] The DSC as:DSC=2∗V(x∩y)Vx+V(y)

where x is a corticospinal tract voxel from the reference data, y is the corresponding voxel from a fiber bundle generated using input or input plus precision data and V is the volume of the relevant voxel.

Statistics

Post hoc analysis with Freidman tests was conducted with a Bonferroni correction (p < 0.0167) applied to the FA and E1.The Wilcoxon signed-rank test (p < 0.05) was applied to the JSD, ACC, and DSC data.

Statistical analyses were performed using EZR version 4.3.1 software (https://www.jichi.ac.jp/saitama-sct/SaitamaHP.files/download.html).

Unless otherwise specified, data were presented as mean ± SD. Effect size [46] was also calculated for each test.

Results

Patient backgrounds

Tables 2 and 3 shows the medical conditions of the participants in each data group. The average ages (± SD) of the patients in the training, validation, and test sets were 45.7 (± 19.1), 41.6 (± 19.1), and 49.6 (± 18.3), respectively. Table 2 Patient conditions in the three groups

Table 3 Patient detail conditions in the test data

For brain diseases, cases with only one patient in the training data were classified as “other.” Similarly, when there was only one case with a given tumor location, it was classified as “other.”

Tumor location and brain disease types with no relevant data or numerical values were marked with a hyphen (-) (Table 2).

Training and validation of the convolutional neural network

The MSE values from the training process are shown in Fig. 5. The values for the training and validation data after 300 epochs were 15.935 and 18.613, respectively.Fig. 5 Loss function results

MSE, mean squared error.

DWI similarity evaluation using the test data

The average SSIM for the test data was 0.964 (± 0.010). The minimum and maximum SSIMs for all data were 0.913 and 0.982. The average PSNR for the test data was 35.41 dB (± 1.59). The minimum and maximum PSNRs for all data were 30.36 and 38.32 (Suppl. data). Image comparisons were made between the reference data and predicted data. For each MPG, the predicted image had similar anisotropic contrast to the reference image (Fig. 6).Fig. 6 Comparison of prediction and reference images from an example case

The case shown is a 27-year-old female with an astrocytoma in the left thalamus.

For each MPG, the prediction image had the same anisotropic contrast as the reference image.

MPG, motion-probing gradient.

There were no obvious outliers from the trend in the average SSIM of each MPG (Fig. 7).Fig. 7 Structural similarities in each motion-probing gradient

The blue dots represent the average structural similarities. The error bars represent the standard deviations for each motion-probing gradient.

MPG, motion-probing gradient; SSIM, structural similarity.

Diffusion profile evaluation of the test data

For both the JSD and the ACC, input plus prediction data was significantly closer to the reference data than input data alone (p < 0.001); that is, the ODF were similar. For the JSD, the input data was 1.25*10–4 and the input plus prediction data was 6.75*10–5. For the ACC, the input data was 9.9954*10–1 and the input plus prediction data was 9.9974*10–1. Figures 8 and 9 show the results in an example case.Fig. 8 JSD and ACC distributions in an example case

Fig. 9 The image results from the JSD and ACC of an example case

Input + predicted data had JSD closer to 0 and ACC closer to 1 than input data only. This showed a high degree of similarity to the diffusion profile of the reference data from the JSD and ACC definition equations. This result was the same for all individual data and overall average.

The case shown is a 75-year-old man with a primary central nervous system lymphoma in the right frontal parietal area. The closer the JSD is to 0, the higher the similarity to the reference data diffusion profile.

Overall, there was greater similarity (profile shifted to the left) with input plus prediction data (orange sticks) than with input data alone (blue sticks).

ACC is a quantitative value indicating that the closer it gets to 1, the higher the similarity. Overall similarity was better (profile shifted to the left) with input plus prediction data (orange sticks) than with input data alone (blue sticks).

ACC, angular correlation coefficient; JSD, Jensen-Shannon divergence.

The case shown is a 71-year-old man with anaplastic astrocytoma in the left temporal lobe. It can be seen that the addition of prediction data improved the quantitative values of the JSD and ACC. However, especially in ACC, even input data alone showed high similarity to reference data, so adding prediction data did not have a large effect.

ACC, angular correlation coefficient; JSD, Jensen-Shannon divergence.

Tensor analysis of the test data

Example images of FA and color FA map are shown in Fig. 10 and Fig. 11. In each figure, the upper row is the images of reference data, the middle row is the images of only input data, and the lower row is the images of input plus prediction data. There was no visually obvious difference.Fig. 10 The image results from the FA of an example case

Fig. 11 The image results from the Color FA of an example case

The case shown is a 37-year-old woman with glioblastoma in the right temporal lobe. The upper row is the images of reference data, the middle row is the images of only input data, and the lower row is the images of input plus prediction data. There was no visually obvious difference.

The case shown is a 37-year-old woman with glioblastoma in the right temporal lobe. The upper row is the images of reference data, the middle row is the images of only input data, and the lower row is the images of input plus prediction data. There was no visually obvious difference.

The average FA and SD in ROI1 was 0.552 (± 0.077) for the reference data, 0.562 (± 0.078) for the only input data, and 0.554 (± 0.080) for the input plus prediction data, respectively. The average FA and SD in ROI2 was 0.521 (± 0.044) for the reference data, 0.512 (± 0.041) for the only input data, and 0.520 (± 0.046) for the input plus prediction data, respectively. The average FA and SD in ROI3 was 0.471 (± 0.044) for the reference data, 0.470 (± 0.052) for the only input data, and 0.471 (± 0.048) for the input plus prediction data, respectively. The average FA and SD in ROI4 was 0.434 (± 0.049) for the reference data, 0.437 (± 0.050) for the only input data, and 0.437 (± 0.051) for the input plus prediction data, respectively. In ROI2, there were significant differences between reference data and only input data (p < 0.001) and between only input data and input plus prediction data (p < 0.01).

The average E1 [µm/ms2] and SD in ROI1 was 0.664 (± 0.11) for the reference data, 0.659 (± 0.11) for the only input data, and 0.655 (± 0.11) for the input plus prediction data, respectively. The average E1 [µm/ms2] and SD in ROI2 was 0.682 (± 0.067) for the reference data, 0.680 (± 0.066) for the only input data, and 0.680 (± 0.071) for the input plus prediction data, respectively. The average E1 [µm/ms2] and SD in ROI3 was 0.684 (± 0.051) for the reference data, 0.683 (± 0.048) for the only input data, and 0.680 (± 0.049) for the input plus prediction data, respectively. The average E1 [µm/ms2] and SD in ROI4 was 0.695 (± 0.057) for the reference data, 0.694 (± 0.057) for the only input data, and 0.695 (± 0.059) for the input plus prediction data, respectively. There were no significant differences among comparison conditions (Suppl. data).

The effect sizes were 0.213 and 0.164 for FA in ROI2 between Reference and Only input, and Only input and Input plus prediction, where there was a significant difference, respectively. Effect sizes in these comparisons were "small (0.1 to 0.3)”. Other effect sizes were less than 0.1 (Suppl data)(Fig. 12).Fig. 12 The results of ROIs measurement

ROI1 indicates the cerebral peduncle, ROI2 indicates posterior peduncle of the internal capsule, ROI3 indicates corona radiata and ROI4 indicates semiovale center in the normal side.

In the (a) FA value comparison, significant differences occurred between only input and reference, only input and input plus prediction. In the (b) E1 value comparison, there was no significant difference.

Tractography evaluation of the test data

For the tractography DSCs, three bundles of CST could not be visualized by QBI tractography using TractVis (right side of subject no. 3, 22, and 30). Therefore, these three were excluded from our comparison, which was made using the remaining 57 cases. The average DSC for the input data alone was 0.697 (± 0.070). The average DSC for the input plus prediction data was 0.715 (± 0.064). The DSC was significantly higher with prediction data than with input data alone (p < 0.05). GQI tractography using DSI Studio also failed to visualize two bundles of CST (sub no. 3, right of 30). Thus, a comparison was made using 58 samples, excluding these two cases. The average DSC for the input data alone was 0.738 (± 0.118). The average DSC for the input plus prediction data was 0.769 (± 0.091). Again, the DSC was significantly higher with prediction data than with input data alone (p < 0.01).

The effect size was 0.263 and 0.288 for QBI and GQI, respectively.

The resulting image examples shown in Fig. 13 and 14 are for the same patient. They illustrate the tendency for the DSC to decrease as a result of increases in mis-tracking and decreases in the extraction area (Figs. 13 and 14). For the case shown in Fig. 13, the DSCs of the left and right corticospinal tract with GQI were 0.693 and 0.740 for the input data, and 0.721 and 0.808 for the input plus prediction data, respectively.Fig. 13 Example tractography result generated using generalized q-sampling imaging

Fig. 14 Example of a tractography result generated using Q-ball imaging

The case shown is a 34-year-old man with a glioma in the left frontal lobe. The values represent the DSCs in the left and right corticospinal tracts. The DSCs of the left and right corticospinal tract were 0.693 and 0.740 for input data, and 0.721 and 0.808 for input plus prediction data, respectively.

DSC, dice similarity coefficient; L, left; R, right.

For the case shown in Fig. 14, the DSCs of the left and right CST with QBI were 0.676 and 0.714 for the input data, and 0.687 and 0.734 for the input plus prediction data, respectively.

The case shown is a 34-year-old man with a glioma in the left frontal lobe. The values represent the DSCs for the left and right corticospinal tract. The DSCs of the left and right corticospinal tract were 0.676 and 0.714 for the input data, and 0.687 and 0.734 for the input plus prediction data, respectively.

DSC, dice similarity coefficient; L, left; R, right.

An increase in the DSC values was observed in 35/57 bundles (61.4%) with QBI and 42/58 bundles (72.4%) with GQI. This indicates that the AI model produced a general improvement in the tractography visualization ability.

As for FAT, three bundles of FAT could not be visualized by QBI tractography using TractVis (both side of subject no. 21, and right side of subject no. 29). Therefore, these three were excluded from our comparison, which was made using the remaining 57 cases. The average DSC for the input data alone was 0.502 (± 0.124). The average DSC for the input plus prediction data was 0.531 (± 0.128). The DSC was significantly higher with prediction data than with input data alone (p < 0.001). GQI tractography using DSI Studio also failed to visualize two bundles of FAT (right side of subject no. 27). Thus, a comparison was made using 59 samples, excluding one case. The average DSC for the input data alone was 0.793 (± 0.101). The average DSC for the input plus prediction data was 0.816 (± 0.093). Again, the DSC was significantly higher with prediction data than with input data alone (p < 0.001).

The effect size was 0.235 and 0.240 for QBI and GQI, respectively.

An increase in the DSC values was observed in 42/57 bundles (73.7%) with QBI and 40/59 bundles (67.8%) with GQI. This indicates that the AI model produced a general improvement in the tractography visualization ability.

Discussion

In this study, we generated unimaged DWI from imaged DWI HARDI data. Research on image generation using AI has evaluated the generation of CT from MR T1-weighted images [47], the generation of methionine positron emission tomography from contrast-enhanced T1-weighted images [48], and the generation of fluid-attenuated inversion recovery from DWI [49]. However, this is the first study to generate predicted images with the same contrast in the same series. We have shown that the second half of a HARDI dataset generated from the first half closely matches the imaged second-half data, with almost no bias in the results due to the MPG axis. Absolute similarity evaluation cannot be performed using SSIM alone, but previous research with 3D U-net on image similarity using other SSIM [50–52] suggests that the images generated by this study have a relatively high degree of similarity. However, the smoothness of the image may have caused the SSIM to show a high value, and further investigation is required. The QBI and GQI used in this study were HARDI-based datasets. HARDI-based datasets provide more information about cross-fibers than DTI, the imaging time is longer for requiring more MPG direction. While it is possible to reduce the imaging time using SMS, it is known that an increase in the SMS factor leads to a deterioration in image quality. Moreover, this type of artifact differs from classical artifacts that appear continuously within or along the slice, making interpretation challenging [31]. Our method can shorten imaging time by generating the second half of the imaged data from the first half. The generation of images in this way can reduce the incidence of artifacts caused by SMS.

DWI data is generally preprocessed for noise removal [33–35], distortion correction, and motion correction [39, 40]. Regarding body movement, the longer the imaging time, the more likely it is that the body will move from its initial position, which will lower the correction accuracy. However, this was an resting-state fMRI study, and it has been reported that 76% of the data in such studies could be collected in the first 10 min of the 20-min imaging time [53]. On this basis, it can be inferred that the longer the scan time, the lower the reliability of the data from the latter half of the scan time. With this in mind, the method presented here could effectively improve the accuracy and reliability of the imaging data through its halving of imaging times.

SSIM was used to evaluate the similarity of the original image, and JSD and ACC were used to evaluate the similarity of the diffusion profile. In this study, the ODF values generated using input plus prediction data were more similar to the ODF of the reference data than the ODF values calculated using input data alone. This was shown by the JSD and ACC values. It has been reported that 60–90% of the brain is chiasm fibers [54]. Furthermore, it has been established that MPGs of about 60 axes are required to visualize tractography using CSD [55] and that visualization performance is improved when using 64-axis data rather than 30-axis data in DTI analyses. It has further been reported that 30 is an insufficient number of MPG axes [11]. Of course, it cannot be said that 64 axes of MPG are sufficient but the effectiveness of adding images generated by this 64-axis MPG AI has been demonstrated.

We performed tensor analysis to calculate quantitative values and we compared FA and E1 using ROI measurements. No difference was visually observed in FA or color FA map (Fig. 10 and 11). Regarding ROI measurements, there was no significant difference between input plus prediction data and reference data. Although a significant difference occurred only in FA in the input data of ROI2, this result can be said to suggest the usefulness of this method.

We evaluated tractography visualization performance using QBI and GQI analysis. In both analysis algorithms, similar to the ODF study, adding the prediction data improved the value obtained using only the input data for reference. From the perspective of tractography depiction, we found that adding AI data improves the similarity of the predicted data to the reference data. A previous study that evaluated the reproducibility of QBI tractography targeting optic radiation found tractography DSCs using scan-rescan data of 0.6–0.75 [56]. The current input plus prediction data result (DSC = 0.715) for CST was within this range, suggesting that the images may be equivalent to scan-rescans.

On the other hand, the FAT result (DSC = 0.531) was out of range (0.6–0.75). This may be due to the effect of the arcuate fasciculus intersecting the FAT. As a result, we considered that FAT's DSC is overall lower than CST's DSC.

Regarding the difference in analysis algorithms, the DSCs produced by GQI analysis were superior to those from QBI. Although the DSCs were not compared, a previous study found no significant difference in the percentage of true fibers in the pyramidal tract but that the QBI calculated a superior proportion to the GQI [57]. And this report found a high proportion of false fiber, which caused the GQI's DSC to be higher than the QBI's DSC. The present authors assume that these results support the current study.

Regarding the improvement of CST tractography visualization ability by AI, a DSC increase of 61.4% for the QBI and 72.4% for the GQI was observed. When evaluated using the DSC classification of a previous study [54], we found that adding prediction data led to improved outcomes when the input data alone had DSC ≥ 0.6 (good or excellent), but no major improvement was seen (DSC of QBI increased by 0.013, GQI increased by 0.017). On the other hand, in one case, there was a decrease in the DSC in GQI when the DSC < 0.6, but the DSC for QBI showed an increase of 0.087 and the DSC for GQI showed an increase of 0.264 (Fig. 15).Fig. 15 DSC improvement through the addition of prediction data in cases with DSC < 0.6 for input data alone. (a) The result in QBI. (b) The result in GQI

The DSC values decreased with the addition of prediction data in a GQI case but increased in all other cases. The improvement in DSCs resulting from the addition of the prediction data was more effective when the input data DSC was lower.

DSC, dice similarity coefficient; CST, corticospinal tract; GQI, generalized q-sampling imaging; L, left; QBI, Q-ball imaging; R, right.

For the improvement of FAT tractography visualization ability by AI, a DSC increase of 73.7% for the QBI and 67.8% for the GQI was observed. Similar to the trend for CST, we found that adding predicted data led to improved results when the DSC of input data alone was 0.6 or higher (good or excellent) [53], but no significant improvement was observed (DSC of QBI increased by 0.010 and GQI increased by 0.024). On the other hand, when DSC < 0.6, the GQI DSC remained almost unchanged (0.307 to 0.300) in only one case. Regarding QBI, DSC worsened in 8 patients (average decrease of 0.07, maximum decrease of 0.172), and DSC improved in 30 patients (average increase of 0.07, maximum increase of 0.206).

This result suggests that cases with lower DSCs and lower accuracy due to disease may benefit from this approach, especially in CST. However, there were cases in which the DSC was < 0.6 even on the healthy side, so further verification is necessary.

Our study had some limitations. The generated images were only verified using one MRI machine and one type of MPG array. Also, the second 50% of the data was generated from the first 50%. Other generation ratios (e.g., generating the last 30% from the first 70%) were not attempted for comparison. Therefore, the relationship between the generation rate and image similarity was not ascertained. For the diffusion profile, only the ODF was evaluated. We had planned to test another algorithm using the FOD generated by CSD analysis but, because the test data was from brain tumor patients, it was necessary to focus only on the white matter as the ROI. The ROI was set using ss3t-CSD, which can extract white matter like even with a single unit of b-value data [58].The patient’s data for this study have many different diagnoses, it is not always the same tracts can be compared. The data set for this study was performed as part of a preoperative examination for neurosurgery, and imaging time was limited. Therefore, compared to previous studies, the resolution in this study is poorer, so the influence of partial volume effects may be stronger. In addition, although the gold standard of this tractography was used as the reference data, it is unclear whether the fiber tracing results are accurate.

Conclusion

We developed an artificial intelligence model based on the 3D U-net architecture capable of predicting the latter 32 axes of HARDI data (MPG64 axes) from the initial 32 axes. The integration of predicted data into the input significantly improved diffusion profiles and tractography. This suggests the potential use of this approach to reduce MRI scan times by incorporating AI-generated predicted data alongside the scanned data.

Supplementary Information

Below is the link to the electronic supplementary material.Supplementary file1 (XLSX 51 KB)

Abbreviations

ACC Angular correlation coefficient

ADC Apparent diffusion coefficient

AI Artificial intelligence

CST Corticospinal tract

DSC Dice similarity coefficient

DTI Diffusion tensor imaging

DTT Diffusion tensor tractography

EPI Echo planar imaging

FAT Frontal aslant tract

FOD Fiber orientation distribution

GQI Generalized q-sampling imaging

HARDI High angular resolution diffusion imaging

JSD Jensen-Shannon Divergence

MPG Motion-probing gradient

ODF Orientation distribution function

QBI Q-ball imaging

ROI Region of interest

SMS Simultaneous multi-slice

SSIM Structural similarity

Funding

Open Access funding provided by The University of Tokyo. This work was supported by KAKENHI (20K08016) in JAPAN.

Data Availability

The authors confirm that the data supporting the findings of this study are available within the article and its supplementary materials.

Declarations

Ethical approval

The participants in this study had already been imaged to depict QBI tractography. The analysis of these images was performed using a protocol approved by the Ethics Committee of 11605-(5) and 2561-(25). All procedures were performed in accordance with the tenets of the 2013 revision of the Declaration of Helsinki.

Informed consent

The informed consent requirement was waived by the Ethics Committee due to the retrospective nature of the study.

Competing Interests

None.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References  

1. Basser PJ Mattiello J LeBihan D MR diffusion tensor spectroscopy and imaging Biophys J 1994 66 1 259 267 10.1016/S0006-3495(94)80775-1 8130344
Basser PJ, Mattiello J, LeBihan D (1994) MR diffusion tensor spectroscopy and imaging. Biophys J 66(1):259–267. 10.1016/S0006-3495(94)80775-18130344 10.1016/S0006-3495(94)80775-1
2. Mori S Crain BJ Chacko VP van Zijl PC Three-dimensional tracking of axonal projections in the brain by magnetic resonance imaging Ann Neurol 1999 45 2 265 269 10.1002/1531-8249(199902)45:2<265::aid-ana21>3.0.co;2-3 9989633
Mori S, Crain BJ, Chacko VP, van Zijl PC (1999) Three-dimensional tracking of axonal projections in the brain by magnetic resonance imaging. Ann Neurol 45(2):265–269. 10.1002/1531-8249(199902)45:2%3c265::aid-ana21%3e3.0.co;2-39989633 10.1002/1531-8249(199902)45:2<265::aid-ana21>3.0.co;2-3
3. Romano A D'Andrea G Minniti G Mastronardi L Ferrante L Fantozzi LM Eur Radiol 2009 19 12 2798 2808 10.1007/s00330-009-1483-6 19533147
Romano A, D’Andrea G, Minniti G, Mastronardi L, Ferrante L, Fantozzi LM et al (2009) Eur Radiol 19(12):2798–2808. 10.1007/s00330-009-1483-619533147 10.1007/s00330-009-1483-6
4. Kovanlikaya I Firat Z Kovanlikaya A Uluğ AM Mutlu Cihangiroglu M John M Assessment of the corticospinal tract alterations before and after resection of brainstem lesions using Diffusion Tensor Imaging (DTI) and tractography at 3T Eur J Radiol 2011 77 3 383 391 10.1016/j.ejrad.2009.08.012 19767164
Kovanlikaya I, Firat Z, Kovanlikaya A, Uluğ AM, Mutlu Cihangiroglu M, John M et al (2011) Assessment of the corticospinal tract alterations before and after resection of brainstem lesions using Diffusion Tensor Imaging (DTI) and tractography at 3T. Eur J Radiol 77(3):383–391. 10.1016/j.ejrad.2009.08.01219767164 10.1016/j.ejrad.2009.08.012
5. Dubois J Hertz-Pannier L Dehaene-Lambertz G Cointepas Y Le Bihan D Assessment of the early organization and maturation of infants' cerebral white matter fiber bundles: a feasibility study using quantitative diffusion tensor imaging and tractography Neuroimage 2006 30 4 1121 32 10.1016/j.neuroimage.2005.11.022 16413790
Dubois J, Hertz-Pannier L, Dehaene-Lambertz G, Cointepas Y, Le Bihan D (2006) Assessment of the early organization and maturation of infants’ cerebral white matter fiber bundles: a feasibility study using quantitative diffusion tensor imaging and tractography. Neuroimage 30(4):1121–32. 10.1016/j.neuroimage.2005.11.02216413790 10.1016/j.neuroimage.2005.11.022
6. Kamagata K Motoi Y Tomiyama H Abe O Ito K Shimoji K Relationship between cognitive impairment and white-matter alteration in Parkinson’s disease with dementia: tract-based spatial statistics and tract-specific analysis Eur Radiol 2013 23 7 1946 1955 10.1007/s00330-013-2775-4 23404139
Kamagata K, Motoi Y, Tomiyama H, Abe O, Ito K, Shimoji K et al (2013) Relationship between cognitive impairment and white-matter alteration in Parkinson’s disease with dementia: tract-based spatial statistics and tract-specific analysis. Eur Radiol 23(7):1946–1955. 10.1007/s00330-013-2775-423404139 10.1007/s00330-013-2775-4
7. Tuch DS, Weisskoff RM, Belliveau JW, Wedeen VJ (1999) High angular resolution diffusion imaging of the human brain. Proceedings of the 7th Annual Meeting of ISMRM. https://cds.ismrm.org/ismrm-1999/PDF2/321.pdf. Accessed 1 Oct 2023
8. Tuch DS Q-ball imaging Magn Reson Med 2004 52 6 1358 1372 10.1002/mrm.20279 15562495
Tuch DS (2004) Q-ball imaging. Magn Reson Med 52(6):1358–1372. 10.1002/mrm.2027915562495 10.1002/mrm.20279
9. Tournier J-D Calamante F Connelly A Robust determination of the fibre orientation distribution in diffusion MRI: non-negativity constrained super-resolved spherical deconvolution Neuroimage 2007 35 4 1459 1472 10.1016/j.neuroimage.2007.02.016 17379540
Tournier J-D, Calamante F, Connelly A (2007) Robust determination of the fibre orientation distribution in diffusion MRI: non-negativity constrained super-resolved spherical deconvolution. Neuroimage 35(4):1459–1472. 10.1016/j.neuroimage.2007.02.01617379540 10.1016/j.neuroimage.2007.02.016
10. Yeh F-C Wedeen VJ Tseng W-Y Generalized q-sampling imaging IEEE Trans Med Imaging 2010 29 9 1626 1635 10.1109/TMI.2010.2045126 20304721
Yeh F-C, Wedeen VJ, Tseng W-Y (2010) Generalized q-sampling imaging. IEEE Trans Med Imaging 29(9):1626–1635. 10.1109/TMI.2010.204512620304721 10.1109/TMI.2010.2045126
11. Berman JI Lanza MR Blaskey L Edgar JC Roberts TPL High angular resolution diffusion imaging probabilistic tractography of the auditory radiation Am J Neuroradiol 2013 34 8 1573 1578 10.3174/ajnr.A3471 23493892
Berman JI, Lanza MR, Blaskey L, Edgar JC, Roberts TPL (2013) High angular resolution diffusion imaging probabilistic tractography of the auditory radiation. Am J Neuroradiol 34(8):1573–1578. 10.3174/ajnr.A347123493892 10.3174/ajnr.A3471
12. Koike T Tanaka S Kin T Suzuki Y Takayanagi S Takami H (2022) Accurate Preoperative Identification of Motor Speech Area as Termination of Arcuate Fasciculus Depicted by Q-Ball Imaging Tractography World Neurosurg 2022 164 e764 e771 10.1016/j.wneu.2022.05.041 35595046
Koike T, Tanaka S, Kin T, Suzuki Y, Takayanagi S, Takami H et al (2022) (2022) Accurate Preoperative Identification of Motor Speech Area as Termination of Arcuate Fasciculus Depicted by Q-Ball Imaging Tractography. World Neurosurg 164:e764–e771. 10.1016/j.wneu.2022.05.04135595046 10.1016/j.wneu.2022.05.041
13. A M Auriat, M R Borich, N J Snow, K P Wadden, L A Boyd (2015) Neuroimage Clin 14:7:771–81. 10.1016/j.nicl.2015.03.007. eCollection 2015.
14. Kidoh M Shinoda K Kitajima M Isogawa K Nambu M Deep Learning Based Noise Reduction for Brain MR Imaging: Tests on Phantoms and Healthy Volunteers Magn Reson Med Sci 2020 19 3 195 206 10.2463/mrms.mp.2019-0018 31484849
Kidoh M, Shinoda K, Kitajima M, Isogawa K, Nambu M et al (2020) Deep Learning Based Noise Reduction for Brain MR Imaging: Tests on Phantoms and Healthy Volunteers. Magn Reson Med Sci 19(3):195–206. 10.2463/mrms.mp.2019-001831484849 10.2463/mrms.mp.2019-0018
15. Yogananda CGB Shah BR Vejdani-Jahromi M Nalawade SS Murugesan GK A novel fully automated MRI-based deep-learning method for classification of IDH mutation status in brain gliomas Neuro Oncol 2020 22 3 402 411 10.1093/neuonc/noz199 31637430
Yogananda CGB, Shah BR, Vejdani-Jahromi M, Nalawade SS, Murugesan GK et al (2020) A novel fully automated MRI-based deep-learning method for classification of IDH mutation status in brain gliomas. Neuro Oncol 22(3):402–411. 10.1093/neuonc/noz19931637430 10.1093/neuonc/noz199
16. Gong E Pauly JM Wintermark M Zaharchuk G Deep learning enables reduced gadolinium dose for contrast-enhanced brain MRI J Magn Reson Imaging 2018 48 2 330 340 10.1002/jmri.25970 29437269
Gong E, Pauly JM, Wintermark M, Zaharchuk G (2018) Deep learning enables reduced gadolinium dose for contrast-enhanced brain MRI. J Magn Reson Imaging 48(2):330–340. 10.1002/jmri.2597029437269 10.1002/jmri.25970
17. Sichtermann T Faron A Sijben R Teichert N Freiherr J Wiesmann M Deep Learning-Based Detection of Intracranial Aneurysms in 3D TOF-MRA AJNR Am J Neuroradiol 2019 40 1 25 32 10.3174/ajnr.A5911 30573461
Sichtermann T, Faron A, Sijben R, Teichert N, Freiherr J, Wiesmann M (2019) Deep Learning-Based Detection of Intracranial Aneurysms in 3D TOF-MRA. AJNR Am J Neuroradiol 40(1):25–32. 10.3174/ajnr.A591130573461 10.3174/ajnr.A5911
18. Tian Q Bilgic B Fan Q Liao C Ngamsombat C DeepDTI: High-fidelity six-direction diffusion tensor imaging using deep learning Neuroimage 2020 219 117017 10.1016/j.neuroimage.2020.117017 32504817
Tian Q, Bilgic B, Fan Q, Liao C, Ngamsombat C et al (2020) DeepDTI: High-fidelity six-direction diffusion tensor imaging using deep learning. Neuroimage 219:117017. 10.1016/j.neuroimage.2020.11701732504817 10.1016/j.neuroimage.2020.117017
19. Tian Q Li Z Fan Q Polimeni JR Bilgic B SDnDTI: Self-supervised deep learning-based denoising for diffusion tensor MRI Neuroimage 2022 253 119033 10.1016/j.neuroimage.2022.119033 35240299
Tian Q, Li Z, Fan Q, Polimeni JR, Bilgic B et al (2022) SDnDTI: Self-supervised deep learning-based denoising for diffusion tensor MRI. Neuroimage 253:119033. 10.1016/j.neuroimage.2022.11903335240299 10.1016/j.neuroimage.2022.119033
20. Qiyuan Tian, Ziyu Li, Qiuyun Fan, Chanon Ngamsombat, Yuxin Hu et al (2021) SRDTI: Deep learning-based super-resolution for diffusion tensor MRI. arXiv:2102.09069. 10.48550/arXiv.2102.09069.
21. Karimi D Vasung L Jaimes C Machado-Rivas F Khan S A machine learning-based method for estimating the number and orientations of major fascicles in diffusion-weighted magnetic resonance imaging Med Image Anal 2021 72 102129 10.1016/j.media.2021.102129 34182203
Karimi D, Vasung L, Jaimes C, Machado-Rivas F, Khan S et al (2021) A machine learning-based method for estimating the number and orientations of major fascicles in diffusion-weighted magnetic resonance imaging. Med Image Anal 72:102129. 10.1016/j.media.2021.10212934182203 10.1016/j.media.2021.102129
22. Ranjeet Ranjan Jha Sudhir K Pathak, Vishwesh Nath, Walter Schneider, B V Rathish Kumar et al Magn Reson Imaging 2022 90 1 16 10.1016/j.mri.2022.03.004 35341904
Ranjeet Ranjan Jha (2022) Sudhir K Pathak, Vishwesh Nath, Walter Schneider, B V Rathish Kumar et al. Magn Reson Imaging 90:1–16. 10.1016/j.mri.2022.03.00435341904 10.1016/j.mri.2022.03.004
23. Richards TJ Anderson KL Anderson JS Fully automated segmentation of the corticospinal tract using the TractSeg algorithm in patients with brain tumors Clin Neurol Neurosurg 2021 210 107001 10.1016/j.clineuro.2021.107001 34749021
Richards TJ, Anderson KL, Anderson JS (2021) Fully automated segmentation of the corticospinal tract using the TractSeg algorithm in patients with brain tumors. Clin Neurol Neurosurg 210:107001. 10.1016/j.clineuro.2021.10700134749021 10.1016/j.clineuro.2021.107001
24. Jensen JH Helpern JA Ramani A Hanzhang Lu Kaczynski K Diffusional kurtosis imaging: the quantification of non-gaussian water diffusion by means of magnetic resonance imaging Magn Reson Med 2005 53 6 1432 1440 10.1002/mrm.20508 15906300
Jensen JH, Helpern JA, Ramani A, Hanzhang Lu, Kaczynski K (2005) Diffusional kurtosis imaging: the quantification of non-gaussian water diffusion by means of magnetic resonance imaging. Magn Reson Med 53(6):1432–1440. 10.1002/mrm.2050815906300 10.1002/mrm.20508
25. Zhang H Schneider T Wheeler-Kingshott CA Alexander DC NODDI: practical in vivo neurite orientation dispersion and density imaging of the human brain Neuroimage 2012 61 4 1000 1016 10.1016/j.neuroimage.2012.03.072 22484410
Zhang H, Schneider T, Wheeler-Kingshott CA, Alexander DC (2012) NODDI: practical in vivo neurite orientation dispersion and density imaging of the human brain. Neuroimage 61(4):1000–1016. 10.1016/j.neuroimage.2012.03.07222484410 10.1016/j.neuroimage.2012.03.072
26. Moeller S Yacoub E Olman CA Auerbach E Strupp J Multiband multislice GE-EPI at 7 tesla, with 16-fold acceleration using partial parallel imaging with application to high spatial and temporal whole-brain fMRI Magn Reson Med 2010 63 5 1144 1153 10.1002/mrm.22361 20432285
Moeller S, Yacoub E, Olman CA, Auerbach E, Strupp J et al (2010) Multiband multislice GE-EPI at 7 tesla, with 16-fold acceleration using partial parallel imaging with application to high spatial and temporal whole-brain fMRI. Magn Reson Med 63(5):1144–1153. 10.1002/mrm.2236120432285 10.1002/mrm.22361
27. Setsompop K Gagoski BA Polimeni JR Witzel T Wedeen VJ Blipped-controlled aliasing in parallel imaging for simultaneous multislice echo planar imaging with reduced g-factor penalty Magn Reson Med 2012 67 5 1210 1224 10.1002/mrm.23097 21858868
Setsompop K, Gagoski BA, Polimeni JR, Witzel T, Wedeen VJ et al (2012) Blipped-controlled aliasing in parallel imaging for simultaneous multislice echo planar imaging with reduced g-factor penalty. Magn Reson Med 67(5):1210–1224. 10.1002/mrm.2309721858868 10.1002/mrm.23097
28. Junqian Xu Moeller S Auerbach EJ Strupp J Smith SM Evaluation of slice accelerations using multiband echo planar imaging at 3 T Neuroimage 2013 83 991 1001 10.1016/j.neuroimage.2013.07.055 23899722
Junqian Xu, Moeller S, Auerbach EJ, Strupp J, Smith SM et al (2013) Evaluation of slice accelerations using multiband echo planar imaging at 3 T. Neuroimage 83:991–1001. 10.1016/j.neuroimage.2013.07.05523899722 10.1016/j.neuroimage.2013.07.055
29. Larkman DJ Hajnal JV Herlihy AH Coutts GA Young IR Ehnholm G Use of multicoil arrays for separation of signal from multiple slices simultaneously excited J Magn Reson Imaging 2001 13 2 313 317 10.1002/1522-2586(200102)13:2<313::aid-jmri1045>3.0.co;2-w 11169840
Larkman DJ, Hajnal JV, Herlihy AH, Coutts GA, Young IR, Ehnholm G (2001) Use of multicoil arrays for separation of signal from multiple slices simultaneously excited. J Magn Reson Imaging 13(2):313–317. 10.1002/1522-2586(200102)13:2%3c313::aid-jmri1045%3e3.0.co;2-w11169840 10.1002/1522-2586(200102)13:2<313::aid-jmri1045>3.0.co;2-w
30. Feinberg DA Setsompop K Ultra-fast MRI of the human brain with simultaneous multi-slice imaging J Magn Reason 2013 229 90 100 10.1016/j.jmr.2013.02.002
Feinberg DA, Setsompop K (2013) Ultra-fast MRI of the human brain with simultaneous multi-slice imaging. J Magn Reason 229:90–100. 10.1016/j.jmr.2013.02.00210.1016/j.jmr.2013.02.002
31. Barth M Breuer F Koopmans PJ Norris DG Poser BA Simultaneous multislice (SMS) imaging techniques Magn Reson Med 2016 75 1 63 81 10.1002/mrm.25897 26308571
Barth M, Breuer F, Koopmans PJ, Norris DG, Poser BA (2016) Simultaneous multislice (SMS) imaging techniques. Magn Reson Med 75(1):63–81. 10.1002/mrm.2589726308571 10.1002/mrm.25897
32. Mitsuda M Suzuki Y Kunimatsu A Kasahara A Watanabe Y Feasibility of Diffusion Tensor Imaging at 15T Using Multi-Band Echo Planar Acquisition Magn Reson Med Sci 2017 16 2 169 175 10.2463/mrms.tn.2015-0159 27599585
Mitsuda M, Suzuki Y, Kunimatsu A, Kasahara A, Watanabe Y et al (2017) Feasibility of Diffusion Tensor Imaging at 15T Using Multi-Band Echo Planar Acquisition. Magn Reson Med Sci 16(2):169–175. 10.2463/mrms.tn.2015-015927599585 10.2463/mrms.tn.2015-0159
33. Veraart J Novikov DS Christiaens D Ades-Aron B Sijbers J Fieremans E Denoising of diffusion MRI using random matrix theory Neuroimage 2016 15 142 394 406 10.1016/j.neuroimage.2016.08.016
Veraart J, Novikov DS, Christiaens D, Ades-Aron B, Sijbers J, Fieremans E (2016) Denoising of diffusion MRI using random matrix theory. Neuroimage 15(142):394–406. 10.1016/j.neuroimage.2016.08.01610.1016/j.neuroimage.2016.08.016
34. Veraart J Fieremans E Novikov DS Diffusion MRI noise mapping using random matrix theory Magn Reson Med 2016 76 5 1582 1593 10.1002/mrm.26059 26599599
Veraart J, Fieremans E, Novikov DS (2016) Diffusion MRI noise mapping using random matrix theory. Magn Reson Med 76(5):1582–1593. 10.1002/mrm.2605926599599 10.1002/mrm.26059
35. Cordero-Grande L Christiaens D Hutter J Price AN Hajnal JV Complex diffusion-weighted image estimation via matrix recovery under general noise models Neuroimage 2019 15 200 391 404 10.1016/j.neuroimage.2019.06.039
Cordero-Grande L, Christiaens D, Hutter J, Price AN, Hajnal JV (2019) Complex diffusion-weighted image estimation via matrix recovery under general noise models. Neuroimage 15(200):391–404. 10.1016/j.neuroimage.2019.06.03910.1016/j.neuroimage.2019.06.039
36. Kellner E Dhital B Kiselev VG Reisert M Gibbs-ringing artifact removal based on local subvoxel-shifts Magn Reson Med 2016 76 5 1574 1581 10.1002/mrm.26054 26745823
Kellner E, Dhital B, Kiselev VG, Reisert M (2016) Gibbs-ringing artifact removal based on local subvoxel-shifts. Magn Reson Med 76(5):1574–1581. 10.1002/mrm.2605426745823 10.1002/mrm.26054
37. Bautista T O'Muircheartaigh J Hajnal JV Tournier, J-D Removal of Gibbs ringing artefacts for 3D acquisitions using subvoxel shifts Proc ISMRM 2021 29 3535
Bautista T, O’Muircheartaigh J, Hajnal JV, Tournier, J-D (2021) Removal of Gibbs ringing artefacts for 3D acquisitions using subvoxel shifts. Proc ISMRM 29:3535
38. Tournier J-D Smith R Raffelt D Tabbara R Dhollander T Pietsch M MRtrix3: A fast, flexible and open software framework for medical image processing and visualization Neuroimage 2019 15 202 116137 10.1016/j.neuroimage.2019.116137
Tournier J-D, Smith R, Raffelt D, Tabbara R, Dhollander T, Pietsch M et al (2019) MRtrix3: A fast, flexible and open software framework for medical image processing and visualization. Neuroimage 15(202):116137. 10.1016/j.neuroimage.2019.11613710.1016/j.neuroimage.2019.116137
39. Andersson JLR Skare S Ashburner J How to correct susceptibility distortions in spin-echo echo-planar images: application to diffusion tensor imaging Neuroimage 2003 20 2 870 888 10.1016/S1053-8119(03)00336-7 14568458
Andersson JLR, Skare S, Ashburner J (2003) How to correct susceptibility distortions in spin-echo echo-planar images: application to diffusion tensor imaging. Neuroimage 20(2):870–888. 10.1016/S1053-8119(03)00336-714568458 10.1016/S1053-8119(03)00336-7
40. Andersson JLR Sotiropoulos SN An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging Neuroimage 2016 15 125 1063 1078 10.1016/j.neuroimage.2015.10.019
Andersson JLR, Sotiropoulos SN (2016) An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging. Neuroimage 15(125):1063–1078. 10.1016/j.neuroimage.2015.10.01910.1016/j.neuroimage.2015.10.019
41. Tustison NJ Avants BB Cook PA Zheng Y Egan A Yushkevich PA N4ITK: improved N3 bias correction IEEE Trans Med Imaging 2010 29 6 1310 1320 10.1109/TMI.2010.2046908 20378467
Tustison NJ, Avants BB, Cook PA, Zheng Y, Egan A, Yushkevich PA et al (2010) N4ITK: improved N3 bias correction. IEEE Trans Med Imaging 29(6):1310–1320. 10.1109/TMI.2010.204690820378467 10.1109/TMI.2010.2046908
42. Jones DK Horsfield MA Simmons A Optimal strategies for measuring diffusion in anisotropic systems by magnetic resonance imaging Magn Reson Med 1999 42 3 515 525 10.1002/(SICI)1522-2594(199909)42:3<515::AID-MRM14>3.0.CO;2-Q 10467296
Jones DK, Horsfield MA, Simmons A (1999) Optimal strategies for measuring diffusion in anisotropic systems by magnetic resonance imaging. Magn Reson Med 42(3):515–52510467296 10.1002/(SICI)1522-2594(199909)42:3<515::AID-MRM14>3.0.CO;2-Q
43. Schilling KG Janve V Gao Y Stepniewska I Landman BA Anderson AW Histological validation of diffusion MRI fiber orientation distributions and dispersion Neuroimage 2018 15 165 200 221 10.1016/j.neuroimage.2017.10.046
Schilling KG, Janve V, Gao Y, Stepniewska I, Landman BA, Anderson AW (2018) Histological validation of diffusion MRI fiber orientation distributions and dispersion. Neuroimage 15(165):200–221. 10.1016/j.neuroimage.2017.10.04610.1016/j.neuroimage.2017.10.046
44 Szmuda T Rogowska M Słoniewski P Abuhaimed A Szmuda M Springer J Frontal aslant tract projections to the inferior frontal gyrus Folia Morphol (Warsz) 17 2017 76 4 574 581 10.5603/FM.a2017.0039
Szmuda T, Rogowska M, Słoniewski P, Abuhaimed A, Szmuda M, Springer J et al (2017) Frontal aslant tract projections to the inferior frontal gyrus. Folia Morphol (Warsz) 17 76(4):574–581. 10.5603/FM.a2017.003910.5603/FM.a2017.0039
45. Dice LR Measures of the Amount of Ecologic Association Between Species Ecology 1945 26 3 297 302 10.2307/1932409
Dice LR (1945) Measures of the Amount of Ecologic Association Between Species. Ecology 26(3):297–302. 10.2307/193240910.2307/1932409
46. Cohen J (1988) Statistical Power Analysis for the Behavioral Sciences. Routledge, London
47. Han X MR-based synthetic CT generation using a deep convolutional neural network method Med Phys 2017 44 4 1408 1419 10.1002/mp.12155 28192624
Han X (2017) MR-based synthetic CT generation using a deep convolutional neural network method. Med Phys 44(4):1408–1419. 10.1002/mp.1215528192624 10.1002/mp.12155
48. Takita H Matsumoto T Tatekawa H Katayama Y Nakajo K Uda T AI-based Virtual Synthesis of Methionine PET from Contrast-enhanced MRI: Development and External Validation Study Radiology 2023 308 2 e223016 10.1148/radiol.223016 37526545
Takita H, Matsumoto T, Tatekawa H, Katayama Y, Nakajo K, Uda T et al (2023) AI-based Virtual Synthesis of Methionine PET from Contrast-enhanced MRI: Development and External Validation Study. Radiology 308(2):e223016. 10.1148/radiol.22301637526545 10.1148/radiol.223016
49. Benzakoun J Deslys M-A Legrand L Hmeydia G Turc G Hassen WB Synthetic FLAIR as a Substitute for FLAIR Sequence in Acute Ischemic Stroke Radiology 2022 303 1 153 159 10.1148/radiol.211394 35014901
Benzakoun J, Deslys M-A, Legrand L, Hmeydia G, Turc G, Hassen WB et al (2022) Synthetic FLAIR as a Substitute for FLAIR Sequence in Acute Ischemic Stroke. Radiology 303(1):153–159. 10.1148/radiol.21139435014901 10.1148/radiol.211394
50. Le J Tian Ye Mendes J Wilson B Ibrahim M DiBella E Deep learning for radial SMS myocardial perfusion reconstruction using the 3D residual booster U-net Magn Reson Imaging 2021 83 178 188 10.1016/j.mri.2021.08.007 34428512
Le J, Tian Ye, Mendes J, Wilson B, Ibrahim M, DiBella E et al (2021) Deep learning for radial SMS myocardial perfusion reconstruction using the 3D residual booster U-net. Magn Reson Imaging 83:178–188. 10.1016/j.mri.2021.08.00734428512 10.1016/j.mri.2021.08.007
51. Chen C Raymond C Speier W Jin X Cloughesy TF Enzmann D Synthesizing MR Image Contrast Enhancement Using 3D High-Resolution ConvNets IEEE Trans Biomed Eng 2023 70 2 401 412 10.1109/TBME.2022.3192309 35853075
Chen C, Raymond C, Speier W, Jin X, Cloughesy TF, Enzmann D et al (2023) Synthesizing MR Image Contrast Enhancement Using 3D High-Resolution ConvNets. IEEE Trans Biomed Eng 70(2):401–412. 10.1109/TBME.2022.319230935853075 10.1109/TBME.2022.3192309
52. Freedman JN Gurney-Champion OJ Nill S Shiarli A-M Bainbridge HE Mandeville HC Rapid 4D-MRI reconstruction using a deep radial convolutional neural network: Dracula Radiother Oncol 2021 159 209 217 10.1016/j.radonc.2021.03.034 33812914
Freedman JN, Gurney-Champion OJ, Nill S, Shiarli A-M, Bainbridge HE, Mandeville HC et al (2021) Rapid 4D-MRI reconstruction using a deep radial convolutional neural network: Dracula. Radiother Oncol 159:209–217. 10.1016/j.radonc.2021.03.03433812914 10.1016/j.radonc.2021.03.034
53. Dosenbach NUF Koller JM Earl EA Miranda-Dominguez O Klein RL Van AN Real-time motion analytics during brain MRI improve data quality and reduce costs Neuroimage 2017 161 80 93 10.1016/j.neuroimage.2017.08.025 28803940
Dosenbach NUF, Koller JM, Earl EA, Miranda-Dominguez O, Klein RL, Van AN et al (2017) Real-time motion analytics during brain MRI improve data quality and reduce costs. Neuroimage 161:80–93. 10.1016/j.neuroimage.2017.08.02528803940 10.1016/j.neuroimage.2017.08.025
54. Jeurissen B Leemans A Tournier J-D Derek K Jones, Jan Sijbers (2013) Investigating the prevalence of complex fiber configurations in white matter tissue with diffusion magnetic resonance imaging Hum Brain Mapp 2013 34 11 2747 2766 10.1002/hbm.22099 22611035
Jeurissen B, Leemans A, Tournier J-D (2013) Derek K Jones, Jan Sijbers (2013) Investigating the prevalence of complex fiber configurations in white matter tissue with diffusion magnetic resonance imaging. Hum Brain Mapp 34(11):2747–2766. 10.1002/hbm.2209922611035 10.1002/hbm.22099
55. Tournier J-D Calamante F Connelly A MRtrix: Diffusion tractography in crossing fiber regions Int J Imaging Syst Technol 2012 22 53 66 10.1002/ima.22005
Tournier J-D, Calamante F, Connelly A (2012) MRtrix: Diffusion tractography in crossing fiber regions. Int J Imaging Syst Technol 22:53–66. 10.1002/ima.2200510.1002/ima.22005
56. Dayan M Kreutzer S Clark CA Tractography of the optic radiation: a repeatability and reproducibility study NMR Biomed 2015 28 4 423 431 10.1002/nbm.3266 25703088
Dayan M, Kreutzer S, Clark CA (2015) Tractography of the optic radiation: a repeatability and reproducibility study. NMR Biomed 28(4):423–431. 10.1002/nbm.326625703088 10.1002/nbm.3266
57. Suo X Guo L Dianxun Fu Ding H Li Y Qin W A Comparative Study of Diffusion Fiber Reconstruction Models for Pyramidal Tract Branches Front Neurosci 2021 9 15 777377 10.3389/fnins.2021.777377
Suo X, Guo L, Dianxun Fu, Ding H, Li Y, Qin W (2021) A Comparative Study of Diffusion Fiber Reconstruction Models for Pyramidal Tract Branches. Front Neurosci 9(15):777377. 10.3389/fnins.2021.77737710.3389/fnins.2021.777377
58. Tournier JD, Mito R, Raffelt D, Connelly A (2019) Improved white matter response function estimation for 3-tissue constrained spherical deconvolution. Proceedings of the 27th Annual Meeting of ISMRM. https://archive.ismrm.org/2019/0555.html. Accessed 1 Oct 2023
