1
|
Plummer JW, Hussain R, Bdaiwi AS, Costa ML, Willmering MM, Parra-Robles J, Cleveland ZI, Walkup LL. Analytical corrections for B 1 -inhomogeneity and signal decay in multi-slice 2D spiral hyperpolarized 129 Xe MRI using keyhole reconstruction. Magn Reson Med 2024. [PMID: 38297511 DOI: 10.1002/mrm.30028] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2023] [Revised: 01/08/2024] [Accepted: 01/09/2024] [Indexed: 02/02/2024]
Abstract
PURPOSE Hyperpolarized xenon MRI suffers from heterogeneous coil bias and magnetization decay that obscure pulmonary abnormalities. Non-physiological signal variability can be mitigated by measuring and mapping the nominal flip angle, and by rescaling the images to correct for signal bias and decay. While flip angle maps can be calculated from sequentially acquired images, scan time and breath-hold duration are doubled. Here, we exploit the low-frequency oversampling of 2D-spiral and keyhole reconstruction to measure flip angle maps from a single acquisition. METHODS Flip angle maps were calculated from two images generated from a single dataset using keyhole reconstructions and a Bloch-equation-based model suitable for hyperpolarized substances. Artifacts resulting from acquisition and reconstruction schemes (e.g., keyhole reconstruction radius, slice-selection profile, spiral-ordering, and oversampling) were assessed using point-spread functions. Simulated flip angle maps generated using keyhole reconstruction were compared against the paired-image approach using RMS error (RMSE). Finally, feasibility was demonstrated for in vivo xenon ventilation imaging. RESULTS Simulations demonstrated accurate flip angle maps and B1 -inhomogeneity correction can be generated with only 1.25-fold central-oversampling and keyhole reconstruction radius = 5% (RMSE = 0.460°). These settings also generated accurate flip angle maps in a healthy control (RSME = 0.337°) and a person with cystic fibrosis (RMSE = 0.404°) in as little as 3.3 s. CONCLUSION Regional lung ventilation images with reduced impact of B1 -inhomogeneity can be acquired rapidly by combining 2D-spiral acquisition, Bloch-equation-based modeling, and keyhole reconstruction. This approach will be especially useful for breath-hold studies where short scan durations are necessary, such as dynamic imaging and applications in children or people with severely compromised respiratory function.
Collapse
Affiliation(s)
- J W Plummer
- Center for Pulmonary Imaging Research, Division of Pulmonary Medicine, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, USA
- Department of Biomedical Engineering, University of Cincinnati, Cincinnati, OH, USA
| | - R Hussain
- Center for Pulmonary Imaging Research, Division of Pulmonary Medicine, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, USA
| | - A S Bdaiwi
- Center for Pulmonary Imaging Research, Division of Pulmonary Medicine, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, USA
- Department of Biomedical Engineering, University of Cincinnati, Cincinnati, OH, USA
| | - M L Costa
- Center for Pulmonary Imaging Research, Division of Pulmonary Medicine, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, USA
- Department of Biomedical Engineering, University of Cincinnati, Cincinnati, OH, USA
| | - M M Willmering
- Center for Pulmonary Imaging Research, Division of Pulmonary Medicine, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, USA
| | - J Parra-Robles
- Center for Pulmonary Imaging Research, Division of Pulmonary Medicine, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, USA
| | - Z I Cleveland
- Center for Pulmonary Imaging Research, Division of Pulmonary Medicine, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, USA
- Department of Biomedical Engineering, University of Cincinnati, Cincinnati, OH, USA
- Imaging Research Center, Department of Radiology, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, USA
- Department of Pediatrics, University of Cincinnati, Cincinnati, OH, USA
| | - L L Walkup
- Center for Pulmonary Imaging Research, Division of Pulmonary Medicine, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, USA
- Department of Biomedical Engineering, University of Cincinnati, Cincinnati, OH, USA
- Imaging Research Center, Department of Radiology, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, USA
- Department of Pediatrics, University of Cincinnati, Cincinnati, OH, USA
| |
Collapse
|
2
|
Niedbalski PJ, Willmering MM, Robertson SH, Freeman MS, Loew W, Giaquinto RO, Ireland C, Pratt RG, Dumoulin CL, Woods JC, Cleveland ZI. Mapping and correcting hyperpolarized magnetization decay with radial keyhole imaging. Magn Reson Med 2019; 82:367-376. [PMID: 30847967 DOI: 10.1002/mrm.27721] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2018] [Revised: 01/17/2019] [Accepted: 02/11/2019] [Indexed: 12/11/2022]
Abstract
PURPOSE Hyperpolarized (HP) media enable biomedical imaging applications that cannot be achieved with conventional MRI contrast agents. Unfortunately, quantifying HP images is challenging, because relaxation and radio-frequency pulsing generate spatially varying signal decay during acquisition. We demonstrate that, by combining center-out k-space sampling with postacquisition keyhole reconstruction, voxel-by-voxel maps of regional HP magnetization decay can be generated with no additional data collection. THEORY AND METHODS Digital phantom, HP 129 Xe phantom, and in vivo 129 Xe human (N = 4 healthy; N = 2 with cystic fibrosis) imaging was performed using radial sampling. Datasets were reconstructed using a postacquisition keyhole approach in which 2 temporally resolved images were created and used to generate maps of regional magnetization decay following a simple analytical model. RESULTS Mean, keyhole-derived decay terms showed excellent agreement with the decay used in simulations (R2 = 0.996) and with global attenuation terms in HP 129 Xe phantom imaging (R2 > 0.97). Mean regional decay from in vivo imaging agreed well with global decay values and displayed spatial heterogeneity that matched expected variations in flip angle and oxygen partial pressure. Moreover, these maps could be used to correct variable signal decay across the image volume. CONCLUSIONS We have demonstrated that center-out trajectories combined with keyhole reconstruction can be used to map regional HP signal decay and to quantitatively correct images. This approach may be used to improve the accuracy of quantitative measures obtained from hyperpolarized media. Although validated with gaseous HP 129 Xe in this work, this technique can be generalized to any hyperpolarized agent.
Collapse
Affiliation(s)
- Peter J Niedbalski
- Center for Pulmonary Imaging Research, Division of Pulmonary Medicine, Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio
| | - Matthew M Willmering
- Center for Pulmonary Imaging Research, Division of Pulmonary Medicine, Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio
| | - Scott H Robertson
- Clinical Imaging Physics Group, Department of Radiology, Duke University Medical Center, Durham, North Carolina
| | - Matthew S Freeman
- Center for Pulmonary Imaging Research, Division of Pulmonary Medicine, Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio
| | - Wolfgang Loew
- Imaging Research Center, Department of Radiology, Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio
| | - Randy O Giaquinto
- Imaging Research Center, Department of Radiology, Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio
| | - Christopher Ireland
- Imaging Research Center, Department of Radiology, Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio.,Department of Biomedical Engineering, University of Cincinnati, Cincinnati, Ohio
| | - Ronald G Pratt
- Imaging Research Center, Department of Radiology, Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio
| | - Charles L Dumoulin
- Imaging Research Center, Department of Radiology, Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio.,Department of Biomedical Engineering, University of Cincinnati, Cincinnati, Ohio.,Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati, Ohio
| | - Jason C Woods
- Center for Pulmonary Imaging Research, Division of Pulmonary Medicine, Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio.,Imaging Research Center, Department of Radiology, Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio.,Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati, Ohio
| | - Zackary I Cleveland
- Center for Pulmonary Imaging Research, Division of Pulmonary Medicine, Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio.,Imaging Research Center, Department of Radiology, Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio.,Department of Biomedical Engineering, University of Cincinnati, Cincinnati, Ohio.,Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati, Ohio
| |
Collapse
|
3
|
Letelier K, Urbina J, Andía M, Tejos C, Irarrazaval P, Prieto C, Uribe S. Accelerating dual cardiac phase images using undersampled radial phase encoding trajectories. Magn Reson Imaging 2016; 34:1017-25. [PMID: 27067473 DOI: 10.1016/j.mri.2016.03.017] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2015] [Revised: 03/25/2016] [Accepted: 03/27/2016] [Indexed: 10/22/2022]
Abstract
A three-dimensional dual-cardiac-phase (3D-DCP) scan has been proposed to acquire two data sets of the whole heart and great vessels during the end-diastolic and end-systolic cardiac phases in a single free-breathing scan. This method has shown accurate assessment of cardiac anatomy and function but is limited by long acquisition times. This work proposes to accelerate the acquisition and reconstruction of 3D-DCP scans by exploiting redundant information of the outer k-space regions of both cardiac phases. This is achieved using a modified radial-phase-encoding trajectory and gridding reconstruction with uniform coil combination. The end-diastolic acquisition trajectory was angularly shifted with respect to the end-systolic phase. Initially, a fully-sampled 3D-DCP scan was acquired to determine the optimal percentage of the outer k-space data that can be combined between cardiac phases. Thereafter, prospectively undersampled data were reconstructed based on this percentage. As gold standard images, the undersampled data were also reconstructed using iterative SENSE. To validate the method, image quality assessments and a cardiac volume analysis were performed. The proposed method was tested in thirteen healthy volunteers (mean age, 30years). Prospectively undersampled data (R=4) reconstructed with 50% combination led high quality images. There were no significant differences in the image quality and in the cardiac volume analysis between our method and iterative SENSE. In addition, the proposed approach reduced the reconstruction time from 40min to 1min. In conclusion, the proposed method obtains 3D-DCP scans with an image quality comparable to those reconstructed with iterative SENSE, and within a clinically acceptable reconstruction time.
Collapse
Affiliation(s)
- Karis Letelier
- Biomedical Imaging Center, Pontificia Universidad Católica de Chile, Santiago, Chile; Electrical Engineering Department, Faculty of Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile
| | - Jesus Urbina
- Biomedical Imaging Center, Pontificia Universidad Católica de Chile, Santiago, Chile
| | - Marcelo Andía
- Biomedical Imaging Center, Pontificia Universidad Católica de Chile, Santiago, Chile; Radiology Department, School of Medicine, Pontificia Universidad Católica de Chile, Santiago, Chile
| | - Cristián Tejos
- Biomedical Imaging Center, Pontificia Universidad Católica de Chile, Santiago, Chile; Electrical Engineering Department, Faculty of Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile
| | - Pablo Irarrazaval
- Biomedical Imaging Center, Pontificia Universidad Católica de Chile, Santiago, Chile; Electrical Engineering Department, Faculty of Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile
| | - Claudia Prieto
- Division of imaging sciences, King's College London, London, UK
| | - Sergio Uribe
- Biomedical Imaging Center, Pontificia Universidad Católica de Chile, Santiago, Chile; Radiology Department, School of Medicine, Pontificia Universidad Católica de Chile, Santiago, Chile.
| |
Collapse
|
4
|
Ruppert K. Biomedical imaging with hyperpolarized noble gases. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2014; 77:116701. [PMID: 25360484 DOI: 10.1088/0034-4885/77/11/116701] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Hyperpolarized noble gases (HNGs), polarized to approximately 50% or higher, have led to major advances in magnetic resonance (MR) imaging of porous structures and air-filled cavities in human subjects, particularly the lung. By boosting the available signal to a level about 100 000 times higher than that at thermal equilibrium, air spaces that would otherwise appear as signal voids in an MR image can be revealed for structural and functional assessments. This review discusses how HNG MR imaging differs from conventional proton MR imaging, how MR pulse sequence design is affected and how the properties of gas imaging can be exploited to obtain hitherto inaccessible information in humans and animals. Current and possible future imaging techniques, and their application in the assessment of normal lung function as well as certain lung diseases, are described.
Collapse
|
5
|
Subashi E, Moding EJ, Cofer GP, MacFall JR, Kirsch DG, Qi Y, Johnson GA. A comparison of radial keyhole strategies for high spatial and temporal resolution 4D contrast-enhanced MRI in small animal tumor models. Med Phys 2013; 40:022304. [PMID: 23387766 DOI: 10.1118/1.4774050] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE Dynamic contrast-enhanced (DCE) MRI has been widely used as a quantitative imaging method for monitoring tumor response to therapy. The simultaneous challenges of increasing temporal and spatial resolution in a setting where the signal from the much smaller voxel is weaker have made this MR technique difficult to implement in small-animal imaging. Existing protocols employed in preclinical DCE-MRI acquire a limited number of slices resulting in potentially lost information in the third dimension. This study describes and compares a family of four-dimensional (3D spatial + time), projection acquisition, radial keyhole-sampling strategies that support high spatial and temporal resolution. METHODS The 4D method is based on a RF-spoiled, steady-state, gradient-recalled sequence with minimal echo time. An interleaved 3D radial trajectory with a quasi-uniform distribution of points in k-space was used for sampling temporally resolved datasets. These volumes were reconstructed with three different k-space filters encompassing a range of possible radial keyhole strategies. The effect of k-space filtering on spatial and temporal resolution was studied in a 5 mM CuSO(4) phantom consisting of a meshgrid with 350-μm spacing and in 12 tumors from three cell lines (HT-29, LoVo, MX-1) and a primary mouse sarcoma model (three tumors∕group). The time-to-peak signal intensity was used to assess the effect of the reconstruction filters on temporal resolution. As a measure of heterogeneity in the third dimension, the authors analyzed the spatial distribution of the rate of transport (K(trans)) of the contrast agent across the endothelium barrier for several different types of tumors. RESULTS Four-dimensional radial keyhole imaging does not degrade the system spatial resolution. Phantom studies indicate there is a maximum 40% decrease in signal-to-noise ratio as compared to a fully sampled dataset. T1 measurements obtained with the interleaved radial technique do not differ significantly from those made with a conventional Cartesian spin-echo sequence. A bin-by-bin comparison of the distribution of the time-to-peak parameter shows that 4D radial keyhole reconstruction does not cause significant temporal blurring when a temporal resolution of 9.9 s is used for the subsamples of the keyhole data. In vivo studies reveal substantial tumor heterogeneity in the third spatial dimension that may be missed with lower resolution imaging protocols. CONCLUSIONS Volumetric keyhole imaging with projection acquisition provides a means to increase spatiotemporal resolution and coverage over that provided by existing 2D Cartesian protocols. Furthermore, there is no difference in temporal resolution between the higher spatial resolution keyhole reconstruction and the undersampled projection data. The technique allows one to measure complex heterogeneity of kinetic parameters with isotropic, microscopic spatial resolution.
Collapse
Affiliation(s)
- Ergys Subashi
- Center for In Vivo Microscopy, Department of Radiology, Duke University Medical Center, Durham, NC 27710, USA
| | | | | | | | | | | | | |
Collapse
|
6
|
Mistry NN, Pollaro J, Song J, De Lin M, Johnson GA. Pulmonary perfusion imaging in the rodent lung using dynamic contrast-enhanced MRI. Magn Reson Med 2008; 59:289-97. [PMID: 18228577 DOI: 10.1002/mrm.21353] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
With the development of various models of pulmonary disease, there is tremendous interest in quantitative regional assessment of pulmonary function. While ventilation imaging has been addressed to a certain extent, perfusion imaging for small animals has not kept pace. In humans and large animals perfusion can be assessed using dynamic contrast-enhanced (DCE) MRI with a single bolus injection of a gadolinium (Gd)-based contrast agent. But the method developed for the clinic cannot be translated directly to image the rodent due to the combined requirements of higher spatial and temporal resolution. This work describes a novel image acquisition technique staggered over multiple, repeatable bolus injections of contrast agent using an automated microinjector, synchronized with image acquisition to achieve dynamic first-pass contrast enhancement in the rat lung. This allows dynamic first-pass imaging that can be used to quantify pulmonary perfusion. Further improvements are made in the spatial and temporal resolution by combining the multiple injection acquisition method with Interleaved Radial Imaging and "Sliding window-keyhole" reconstruction (IRIS). The results demonstrate a simultaneous increase in spatial resolution (<200 mum) and temporal resolution (<200 ms) over previous methods, with a limited loss in signal-to-noise-ratio.
Collapse
Affiliation(s)
- Nilesh N Mistry
- Department of Biomedical Engineering, Duke University, Durham, NC 27710, USA
| | | | | | | | | |
Collapse
|
7
|
KUMASHIRO M, MURASE K, ODA K, FUKUSHIGE M, ITO O, NAGAYAMA M, WATANABE Y. Assessment of Time-Resolved, Dynamic, Contrast-Enhanced MRDSA Using Radial Sliding-Window Reconstruction. Magn Reson Med Sci 2008; 7:1-12. [DOI: 10.2463/mrms.7.1] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
|
8
|
Winkelmann S, Schaeffter T, Koehler T, Eggers H, Doessel O. An optimal radial profile order based on the Golden Ratio for time-resolved MRI. IEEE TRANSACTIONS ON MEDICAL IMAGING 2007; 26:68-76. [PMID: 17243585 DOI: 10.1109/tmi.2006.885337] [Citation(s) in RCA: 502] [Impact Index Per Article: 29.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
In dynamic magnetic resonance imaging (MRI) studies, the motion kinetics or the contrast variability are often hard to predict, hampering an appropriate choice of the image update rate or the temporal resolution. A constant azimuthal profile spacing (111.246 degrees), based on the Golden Ratio, is investigated as optimal for image reconstruction from an arbitrary number of profiles in radial MRI. The profile order is evaluated and compared with a uniform profile distribution in terms of signal-to-noise ratio (SNR) and artifact level. The favorable characteristics of such a profile order are exemplified in two applications on healthy volunteers. First, an advanced sliding window reconstruction scheme is applied to dynamic cardiac imaging, with a reconstruction window that can be flexibly adjusted according to the extent of cardiac motion that is acceptable. Second, a contrast-enhancing k-space filter is presented that permits reconstructing an arbitrary number of images at arbitrary time points from one raw data set. The filter was utilized to depict the T1-relaxation in the brain after a single inversion prepulse. While a uniform profile distribution with a constant angle increment is optimal for a fixed and predetermined number of profiles, a profile distribution based on the Golden Ratio proved to be an appropriate solution for an arbitrary number of profiles.
Collapse
|
9
|
de Rochefort L, Maître X, Fodil R, Vial L, Louis B, Isabey D, Croce C, Darrasse L, Apiou G, Caillibotte G, Bittoun J, Durand E. Phase-contrast velocimetry with hyperpolarized3He for in vitro and in vivo characterization of airflow. Magn Reson Med 2006; 55:1318-25. [PMID: 16700024 DOI: 10.1002/mrm.20899] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
This paper describes a technique that combines radial MRI and phase contrast (PC) to map the velocities of hyperpolarized gases ((3)He) in respiratory airways. The method was evaluated on well known geometries (straight and U-shaped pipes) before it was applied in vivo. Dynamic 2D maps of the three velocity components were obtained from a 10-mm slice with an in-plane spatial resolution of 1.6 mm within 1 s. Integration of the in vitro through-plane velocity over the slice matched the input flow within a relative precision of 6.4%. As expected for the given Reynolds number, a parabolic velocity profile was obtained in the straight pipe. In the U-shaped pipe the three velocity components were measured and compared to a fluid-dynamics simulation so the precision was evaluated as fine as 0.025 m s(-1). The technique also demonstrated its ability to visualize vortices and localize characteristic points, such as the maximum velocity and vortex-center positions. Finally, in vivo feasibility was demonstrated in the human trachea during inhalation.
Collapse
Affiliation(s)
- Ludovic de Rochefort
- U2R2M, Unité de Recherche en Résonance Magnétique Médicale, UMR 8081, CNRS-Université Paris-Sud, Le Kremlin-Bicêtre, France.
| | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
10
|
Song HK, Dougherty L. Dynamic MRI with projection reconstruction and KWIC processing for simultaneous high spatial and temporal resolution. Magn Reson Med 2005; 52:815-24. [PMID: 15389936 DOI: 10.1002/mrm.20237] [Citation(s) in RCA: 100] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
A method for dynamic imaging in MRI is presented that enables the acquisition of a series of images with both high temporal and high spatial resolution. The technique, which is based on the projection reconstruction (PR) imaging scheme, utilizes distinct data acquisition and reconstruction strategies to achieve this simultaneous capability. First, during acquisition, data are collected in multiple undersampled passes, with the view angles interleaved in such a way that those of subsequent passes bisect the views of earlier ones. During reconstruction, these views are weighted according to a previously described k-space weighted image contrast (KWIC) technique that enables the manipulation of image contrast by selective filtering. Unlike conventional undersampled PR methods, the proposed dynamic KWIC technique does not suffer from low image SNR or image degradation due to streaking artifacts. The effectiveness of dynamic KWIC is demonstrated in both simulations and in vivo, high-resolution, contrast-enhanced imaging of breast lesions.
Collapse
Affiliation(s)
- Hee Kwon Song
- Department of Radiology, University of Pennsylvania Medical Center, Philadelphia, Pennsylvania 19104, USA.
| | | |
Collapse
|