1
|
Rothenberger SM, Zhang J, Markl M, Craig BA, Vlachos PP, Rayz VL. 4D flow MRI velocity uncertainty quantification. Magn Reson Med 2025; 93:397-410. [PMID: 39270010 DOI: 10.1002/mrm.30287] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2024] [Revised: 07/27/2024] [Accepted: 08/20/2024] [Indexed: 09/15/2024]
Abstract
PURPOSE An automatic method is presented for estimating 4D flow MRI velocity measurement uncertainty in each voxel. The velocity distance (VD) metric, a statistical distance between the measured velocity and local error distribution, is introduced as a novel measure of 4D flow MRI velocity measurement quality. METHODS The method uses mass conservation to assess the local velocity error variance and the standardized difference of means (SDM) velocity to estimate the velocity error correlations. VD is evaluated as the Mahalanobis distance between the local velocity measurement and the local error distribution. The uncertainty model is validated synthetically and tested in vitro under different flow resolutions and noise levels. The VD's application is demonstrated on two in vivo thoracic vasculature 4D flow datasets. RESULTS Synthetic results show the proposed uncertainty quantification method is sensitive to aliased regions across various velocity-to-noise ratios and assesses velocity error correlations in four- and six-point acquisitions with correlation errors at or under 3.2%. In vitro results demonstrate the method's sensitivity to spatial resolution, venc settings, partial volume effects, and phase wrapping error sources. Applying VD to assess in vivo 4D flow MRI in the aorta demonstrates the expected increase in measured velocity quality with contrast administration and systolic flow. CONCLUSION The proposed 4D flow MRI uncertainty quantification method assesses velocity measurement error owing to sources including noise, intravoxel phase dispersion, and velocity aliasing. This method enables rigorous comparison of 4D flow MRI datasets obtained in longitudinal studies, across patient populations, and with different MRI systems.
Collapse
Affiliation(s)
- Sean M Rothenberger
- Weldon School of Biomedical Engineering, Purdue University, West Lafayette, Indiana, USA
| | - Jiacheng Zhang
- School of Mechanical Engineering, Purdue University, West Lafayette, Indiana, USA
| | - Michael Markl
- Department of Radiology at the Feinberg School of Medicine, Northwestern University, Chicago, Illinois, USA
| | - Bruce A Craig
- Department of Statistics, Purdue University, West Lafayette, Indiana, USA
| | - Pavlos P Vlachos
- Weldon School of Biomedical Engineering, Purdue University, West Lafayette, Indiana, USA
- School of Mechanical Engineering, Purdue University, West Lafayette, Indiana, USA
| | - Vitaliy L Rayz
- Weldon School of Biomedical Engineering, Purdue University, West Lafayette, Indiana, USA
- School of Mechanical Engineering, Purdue University, West Lafayette, Indiana, USA
| |
Collapse
|
2
|
Degenhardt K, Schmidt S, Aigner CS, Kratzer FJ, Seiter DP, Mueller M, Kolbitsch C, Nagel AM, Wieben O, Schaeffter T, Schulz-Menger J, Schmitter S. Toward accurate and fast velocity quantification with 3D ultrashort TE phase-contrast imaging. Magn Reson Med 2024; 91:1994-2009. [PMID: 38174601 DOI: 10.1002/mrm.29978] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Revised: 11/28/2023] [Accepted: 11/28/2023] [Indexed: 01/05/2024]
Abstract
PURPOSE Traditional phase-contrast MRI is affected by displacement artifacts caused by non-synchronized spatial- and velocity-encoding time points. The resulting inaccurate velocity maps can affect the accuracy of derived hemodynamic parameters. This study proposes and characterizes a 3D radial phase-contrast UTE (PC-UTE) sequence to reduce displacement artifacts. Furthermore, it investigates the displacement of a standard Cartesian flow sequence by utilizing a displacement-free synchronized-single-point-imaging MR sequence (SYNC-SPI) that requires clinically prohibitively long acquisition times. METHODS 3D flow data was acquired at 3T at three different constant flow rates and varying spatial resolutions in a stenotic aorta phantom using the proposed PC-UTE, a Cartesian flow sequence, and a SYNC-SPI sequence as reference. Expected displacement artifacts were calculated from gradient timing waveforms and compared to displacement values measured in the in vitro flow experiments. RESULTS The PC-UTE sequence reduces displacement and intravoxel dephasing, leading to decreased geometric distortions and signal cancellations in magnitude images, and more spatially accurate velocity quantification compared to the Cartesian flow acquisitions; errors increase with velocity and higher spatial resolution. CONCLUSION PC-UTE MRI can measure velocity vector fields with greater accuracy than Cartesian acquisitions (although pulsatile fields were not studied) and shorter scan times than SYNC-SPI. As such, this approach is superior to traditional Cartesian 3D and 4D flow MRI when spatial misrepresentations cannot be tolerated, for example, when computational fluid dynamics simulations are compared to or combined with in vitro or in vivo measurements, or regional parameters such as wall shear stress are of interest.
Collapse
Affiliation(s)
- Katja Degenhardt
- Physikalisch-Technische Bundesanstalt (PTB), Braunschweig and Berlin, Berlin, Germany
- Department of Radiology, Charité - Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany
| | - Simon Schmidt
- Center for Magnetic Resonance Research, University of Minnesota, Minneapolis, Minnesota, USA
- Division of Medical Physics in Radiology, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Christoph S Aigner
- Physikalisch-Technische Bundesanstalt (PTB), Braunschweig and Berlin, Berlin, Germany
| | - Fabian J Kratzer
- Division of Medical Physics in Radiology, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Daniel P Seiter
- Department of Medical Physics, University of Wisconsin, Madison, Wisconsin, USA
| | - Max Mueller
- Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Christoph Kolbitsch
- Physikalisch-Technische Bundesanstalt (PTB), Braunschweig and Berlin, Berlin, Germany
| | - Armin M Nagel
- Division of Medical Physics in Radiology, German Cancer Research Center (DKFZ), Heidelberg, Germany
- Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
| | - Oliver Wieben
- Department of Medical Physics, University of Wisconsin, Madison, Wisconsin, USA
- Department of Radiology, University of Wisconsin Madison, Madison, Wisconsin, USA
| | - Tobias Schaeffter
- Physikalisch-Technische Bundesanstalt (PTB), Braunschweig and Berlin, Berlin, Germany
- School of Imaging Science and Biomedical Engineering, King's College London, London, United Kingdom
- Department of Medical Engineering, Technical University of Berlin, Berlin, Germany
| | - Jeanette Schulz-Menger
- Charité-Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin and Humboldt-Universität zu Berlin, Berlin, Germany
- Working Group on Cardiovascular Magnetic Resonance, Experimental and Clinical Research Center, a joint cooperation between the Charité Medical Faculty and the Max-Delbrueck Center for Molecular Medicine, Berlin, Germany
- DZHK (German Center for Cardiovascular Research), Partner Site Berlin, Berlin, Germany
- Department of Cardiology and Nephrology, HELIOS Hospital Berlin-Buch, Berlin, Germany
| | - Sebastian Schmitter
- Physikalisch-Technische Bundesanstalt (PTB), Braunschweig and Berlin, Berlin, Germany
- Center for Magnetic Resonance Research, University of Minnesota, Minneapolis, Minnesota, USA
- Division of Medical Physics in Radiology, German Cancer Research Center (DKFZ), Heidelberg, Germany
| |
Collapse
|
3
|
Roberts GS, Hoffman CA, Rivera-Rivera LA, Berman SE, Eisenmenger LB, Wieben O. Automated hemodynamic assessment for cranial 4D flow MRI. Magn Reson Imaging 2023; 97:46-55. [PMID: 36581214 PMCID: PMC9892280 DOI: 10.1016/j.mri.2022.12.016] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2022] [Accepted: 12/23/2022] [Indexed: 12/27/2022]
Abstract
Cranial 4D flow MRI post-processing typically involves manual user interaction which is time-consuming and associated with poor repeatability. The primary goal of this study is to develop a robust quantitative velocity tool (QVT) that utilizes threshold-based segmentation techniques to improve segmentation quality over prior approaches based on centerline processing schemes (CPS) that utilize k-means clustering segmentation. This tool also includes an interactive 3D display designed for simplified vessel selection and automated hemodynamic visualization and quantification. The performances of QVT and CPS were compared in vitro in a flow phantom and in vivo in 10 healthy participants. Vessel segmentations were compared with ground-truth computed tomography in vitro (29 locations) and manual segmentation in vivo (13 locations) using linear regression. Additionally, QVT and CPS MRI flow rates were compared to perivascular ultrasound flow in vitro using linear regression. To assess internal consistency of flow measures in vivo, conservation of flow was assessed at vessel junctions using linear regression and consistency of flow along vessel segments was analyzed by fitting a Gaussian distribution to a histogram of normalized flow values. Post-processing times were compared between the QVT and CPS using paired t-tests. Vessel areas segmented in vitro (CPS: slope = 0.71, r = 0.95 and QVT: slope = 1.03, r = 0.95) and in vivo (CPS: slope = 0.61, r = 0.96 and QVT: slope = 0.93, r = 0.96) were strongly correlated with ground-truth area measurements. However, CPS (using k-means segmentation) consistently underestimated vessel areas. Strong correlations were observed between QVT and ultrasound flow (slope = 0.98, r = 0.96) as well as flow at junctions (slope = 1.05, r = 0.98). Mean and standard deviation of flow along vessel segments was 9.33e-16 ± 3.05%. Additionally, the QVT demonstrated excellent interobserver agreement and significantly reduced post-processing by nearly 10 min (p < 0.001). By completely automating post-processing and providing an easy-to-use 3D visualization interface for interactive vessel selection and hemodynamic quantification, the QVT offers an efficient, robust, and repeatable means to analyze cranial 4D flow MRI. This software is freely available at: https://github.com/uwmri/QVT.
Collapse
Affiliation(s)
- Grant S Roberts
- Department of Medical Physics, University of Wisconsin School of Medicine and Public Health, 1111 Highland Avenue #1005, Madison, WI 53705, USA.
| | - Carson A Hoffman
- Department of Medical Physics, University of Wisconsin School of Medicine and Public Health, 1111 Highland Avenue #1005, Madison, WI 53705, USA
| | - Leonardo A Rivera-Rivera
- Department of Medical Physics, University of Wisconsin School of Medicine and Public Health, 1111 Highland Avenue #1005, Madison, WI 53705, USA; Wisconsin Alzheimer's Disease Research Center, University of Wisconsin School of Medicine and Public Health, 600 Highland Avenue, J5/1 Mezzanine, Madison, WI 53792, USA.
| | - Sara E Berman
- Wisconsin Alzheimer's Disease Research Center, University of Wisconsin School of Medicine and Public Health, 600 Highland Avenue, J5/1 Mezzanine, Madison, WI 53792, USA.
| | - Laura B Eisenmenger
- Department of Radiology, University of Wisconsin School of Medicine and Public Health, 600 Highland Avenue, E3/366 Clinical Science Center, Madison, WI 53792, USA.
| | - Oliver Wieben
- Department of Medical Physics, University of Wisconsin School of Medicine and Public Health, 1111 Highland Avenue #1005, Madison, WI 53705, USA; Department of Radiology, University of Wisconsin School of Medicine and Public Health, 600 Highland Avenue, E3/366 Clinical Science Center, Madison, WI 53792, USA.
| |
Collapse
|