1
|
Chen Z, Contijoch F, Schluchter A, Grady L, Schaap M, Stayman W, Pack J, McVeigh E. Precise measurement of coronary stenosis diameter with CCTA using CT number calibration. Med Phys 2019; 46:5514-5527. [PMID: 31603567 PMCID: PMC7700731 DOI: 10.1002/mp.13862] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2019] [Revised: 09/26/2019] [Accepted: 10/03/2019] [Indexed: 12/14/2022] Open
Abstract
PURPOSE Coronary x-ray computed tomography angiography (CCTA) continues to develop as a noninvasive method for the assessment of coronary vessel geometry and the identification of physiologically significant lesions. The uncertainty of quantitative lesion diameter measurement due to limited spatial resolution and vessel motion reduces the accuracy of CCTA diagnoses. In this paper, we introduce a new technique called computed tomography (CT)-number-Calibrated Diameter to improve the accuracy of the vessel and stenosis diameter measurements with CCTA. METHODS A calibration phantom containing cylindrical holes (diameters spanning from 0.8 mm through 4.0 mm) capturing the range of diameters found in human coronary vessels was three-dimensional printed. We also printed a human stenosis phantom with 17 tubular channels having the geometry of lesions derived from patient data. We acquired CT scans of the two phantoms with seven different imaging protocols. Calibration curves relating vessel intraluminal maximum voxel value (maximum CT number of a voxel, described in Hounsfield Units, HU) to true diameter, and full-width-at-half maximum (FWHM) to true diameter were constructed for each CCTA protocol. In addition, we acquired scans with a small constant motion (15 mm/s) and used a motion correction reconstruction (Snapshot Freeze) algorithm to correct motion artifacts. We applied our technique to measure the lesion diameter in the 17 lesions in the stenosis phantom and compared the performance of CT-number-Calibrated Diameter to the ground truth diameter and a FWHM estimate. RESULTS In all cases, vessel intraluminal maximum voxel value vs diameter was found to have a simple functional form based on the two-dimensional point spread function yielding a constant maximum voxel value region above a cutoff diameter, and a decreasing maximum voxel value vs decreasing diameter below a cutoff diameter. After normalization, focal spot size and reconstruction kernel were the principal determinants of cutoff diameter and the rate of maximum voxel value reduction vs decreasing diameter. The small constant motion had a significant effect on the CT number calibration; however, the motion-correction algorithm returned the maximum voxel value vs diameter curve to that of stationary vessels. The CT number Calibration technique showed better performance than FWHM estimation of diameter, yielding a high accuracy in the tested range (0.8 mm through 2.5 mm). We found a strong linear correlation between the smallest diameter in each of 17 lesions measured by CT-number-Calibrated Diameter (DC ) and ground truth diameter (Dgt ), (DC = 0.951 × Dgt + 0.023 mm, r = 0.998 with a slope very close to 1.0 and intercept very close to 0 mm. CONCLUSIONS Computed tomography-number-Calibrated Diameter is an effective method to enhance the accuracy of the estimate of small vessel diameters and degree of coronary stenosis in CCTA.
Collapse
Affiliation(s)
- Zhennong Chen
- Department of Bioengineering, UC San Diego School of Engineering, La Jolla CA 92037-0412
| | - Francisco Contijoch
- Department of Bioengineering, UC San Diego School of Engineering, La Jolla CA 92037-0412
- Department of Radiology, UC San Diego School of Medicine, La Jolla CA 92123
| | - Andrew Schluchter
- Department of Bioengineering, UC San Diego School of Engineering, La Jolla CA 92037-0412
| | - Leo Grady
- HeartFlow, Inc, Redwood City, CA 94063
| | | | - Web Stayman
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore MD 21205
| | - Jed Pack
- GE Global Research, Niskayuna, NY
| | - Elliot McVeigh
- Department of Bioengineering, UC San Diego School of Engineering, La Jolla CA 92037-0412
- Department of Radiology, UC San Diego School of Medicine, La Jolla CA 92123
- Department of Cardiology, UC San Diego School of Medicine, La Jolla CA 92123
| |
Collapse
|
2
|
Tran TT, Pham VT, Lin C, Yang HW, Wang YH, Shyu KK, Tseng WYI, Su MYM, Lin LY, Lo MT. Empirical Mode Decomposition and Monogenic Signal-Based Approach for Quantification of Myocardial Infarction From MR Images. IEEE J Biomed Health Inform 2018; 23:731-743. [PMID: 29994104 DOI: 10.1109/jbhi.2018.2821675] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Quantification of myocardial infarction on late Gadolinium enhancement cardiovascular magnetic resonance (LGE-CMR) images into heterogeneous infarct periphery (or gray zone) and infarct core plays an important role in cardiac diagnosis, especially in identifying patients at high risk of cardiovascular mortality. However, quantification task is challenging due to noise corrupted in cardiac MR images, the contrast variation, and limited resolution of images. In this study, we propose a novel approach for automatic myocardial infarction quantification, termed DEMPOT, which consists of three key parts: Decomposition of image into intrinsic modes, monogenic phase performing on combined dominant modes, and multilevel Otsu thresholding on the phase. In particular, inspired by the Hilbert-Huang transform, we perform the multidimensional ensemble empirical mode decomposition and 2-D generalization of the Hilbert transform known as the Riesz transform on the MR image to obtain the monogenic phase that is robust to noise and contrast variation. Then, a two-stage algorithm using multilevel Otsu thresholding is accomplished on the monogenic phase to automatically quantify the myocardium into healthy, gray zone, and infarct core regions. Experiments on LGE-CMR images with myocardial infarction from 82 patients show the superior performance of the proposed approach in terms of reproducibility, robustness, and effectiveness.
Collapse
|
3
|
Tsadok Y, Friedman Z, Haluska BA, Hoffmann R, Adam D. Myocardial strain assessment by cine cardiac magnetic resonance imaging using non-rigid registration. Magn Reson Imaging 2015; 34:381-90. [PMID: 26723847 DOI: 10.1016/j.mri.2015.12.035] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2014] [Accepted: 12/20/2015] [Indexed: 11/25/2022]
Abstract
AIMS To evaluate a novel post-processing method for assessment of longitudinal mid-myocardial strain in standard cine cardiac magnetic resonance (CMR) imaging sequences. METHODS AND RESULTS Cine CMR imaging and tagged cardiac magnetic resonance imaging (TMRI) were performed in 15 patients with acute myocardial infarction (AMI) and 15 healthy volunteers served as control group. A second group of 37 post-AMI patients underwent both cine CMR and late gadolinium enhancement (LGE) CMR exams. Speckle tracking echocardiography (STE) was performed in 36 of these patients. Cine CMR, TMRI and STE were analyzed to obtain longitudinal strain. LGE-CMR datasets were analyzed to evaluate scar extent. Comparison of peak systolic strain (PSS) measured from CMR and TMRI yielded a strong correlation (r=0.86, p<0.001). PSS measured from CMR and STE correlated well (r=0.75, p<0.001). A cutoff longitudinal PSS value of -13.14% differentiated non-infarction from any infarcted myocardium, with a sensitivity of 93% and a specificity of 89% (area under curve (AUC) 0.95). PSS value of -9.39% differentiated non-transmural from transmural infarcted myocardium, with a sensitivity of 75% and a specificity of 67% (AUC 0.78). CONCLUSION The present study showed a novel off-line post-processing method for segmental longitudinal strain analysis in mid-myocardium layer based on cine CMR data. The method was found to be highly correlated with strain measurements obtained by TMRI and STE. This tool allows accurate discrimination between different transmurality states of myocardial infarction.
Collapse
Affiliation(s)
- Yossi Tsadok
- Faculty of Biomedical Engineering, Technion-IIT Technion City, Haifa, Israel.
| | - Zvi Friedman
- General Electric Healthcare, Ultrasound, Tirat HaCarmel, Israel
| | - Brian A Haluska
- Cardiovascular Imaging Research Centre, University of Queensland, Brisbane, Australia
| | - Rainer Hoffmann
- Medical Clinic I, University Hospital RWTH Aachen, Aachen, Germany
| | - Dan Adam
- Faculty of Biomedical Engineering, Technion-IIT Technion City, Haifa, Israel
| |
Collapse
|
4
|
Nayak KS, Nielsen JF, Bernstein MA, Markl M, D Gatehouse P, M Botnar R, Saloner D, Lorenz C, Wen H, S Hu B, Epstein FH, N Oshinski J, Raman SV. Cardiovascular magnetic resonance phase contrast imaging. J Cardiovasc Magn Reson 2015; 17:71. [PMID: 26254979 PMCID: PMC4529988 DOI: 10.1186/s12968-015-0172-7] [Citation(s) in RCA: 154] [Impact Index Per Article: 17.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2015] [Accepted: 07/16/2015] [Indexed: 11/10/2022] Open
Abstract
Cardiovascular magnetic resonance (CMR) phase contrast imaging has undergone a wide range of changes with the development and availability of improved calibration procedures, visualization tools, and analysis methods. This article provides a comprehensive review of the current state-of-the-art in CMR phase contrast imaging methodology, clinical applications including summaries of past clinical performance, and emerging research and clinical applications that utilize today's latest technology.
Collapse
Affiliation(s)
- Krishna S Nayak
- Ming Hsieh Department of Electrical Engineering, University of Southern California, 3740 McClintock Ave, EEB 406, Los Angeles, California, 90089-2564, USA.
| | - Jon-Fredrik Nielsen
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, USA.
| | | | - Michael Markl
- Department of Radiology, Northwestern University, Chicago, IL, USA.
| | - Peter D Gatehouse
- Cardiovascular Biomedical Research Unit, Royal Brompton Hospital, London, UK.
| | - Rene M Botnar
- Cardiovascular Imaging, Imaging Sciences Division, Kings's College London, London, UK.
| | - David Saloner
- Department of Radiology and Biomedical Imaging, University of California San Francisco, San Francisco, CA, USA.
| | - Christine Lorenz
- Center for Applied Medical Imaging, Siemens Corporation, Baltimore, MD, USA.
| | - Han Wen
- Imaging Physics Laboratory, National Heart Lung and Blood Institute, Bethesda, MD, USA.
| | - Bob S Hu
- Palo Alto Medical Foundation, Palo Alto, CA, USA.
| | - Frederick H Epstein
- Departments of Radiology and Biomedical Engineering, University of Virginia, Charlottesville, VA, USA.
| | - John N Oshinski
- Departments of Radiology and Biomedical Engineering, Emory University School of Medicine, Atlanta, GA, USA.
| | - Subha V Raman
- Division of Cardiovascular Medicine, The Ohio State University, Columbus, OH, USA.
| |
Collapse
|
5
|
Riordan AJ, Bennink E, Dankbaar JW, Viergever MA, Velthuis BK, Smit EJ, de Jong HWAM. Comparison of partial volume effects in arterial and venous contrast curves in CT brain perfusion imaging. PLoS One 2014; 9:e97586. [PMID: 24858308 PMCID: PMC4032231 DOI: 10.1371/journal.pone.0097586] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2013] [Accepted: 04/21/2014] [Indexed: 11/29/2022] Open
Abstract
Purpose In brain CT perfusion (CTP), the arterial contrast bolus is scaled to have the same area under the curve (AUC) as the venous outflow to correct for partial volume effects (PVE). This scaling is based on the assumption that large veins are unaffected by PVE. Measurement of the internal carotid artery (ICA), usually unaffected by PVE due to its large diameter, may avoid the need for partial volume correction. The aims of this work are to examine i) the assumptions behind PVE correction and ii) the potential of selecting the ICA obviating correction for PVE. Methods The AUC of the ICA and sagittal sinus were measured in CTP datasets from 52 patients. The AUCs were determined by i) using commercial CTP software based on a Gaussian curve-fitting to the time attenuation curve, and ii) by simple integration of the time attenuation curve over a time interval. In addition, frames acquired up to 3 minutes after first bolus passage were used to examine the ratio of arterial and venous enhancement. The impact of selecting the ICA without PVE correction was illustrated by reporting cerebral blood volume (CBV) measurements. Results In 49 of 52 patients, the AUC of the ICA was significantly larger than that of the sagittal sinus (p = 0.017). Measured after the first pass bolus, contrast enhancement remained 50% higher in the ICA just after the first pass bolus, and 30% higher 3 minutes later. CBV measurements were significantly lowered when the ICA was used without PVE correction. Conclusions Contradicting the assumptions underlying PVE correction, contrast in the ICA was significantly higher than in the sagittal sinus, even 3 minutes after the first pass of the contrast bolus. PVE correction might lead to overestimation of CBV if the CBV is calculated using the AUC of the time attenuation curves.
Collapse
Affiliation(s)
- Alan J. Riordan
- Department of Radiology, University Medical Centre, Utrecht, The Netherlands
- * E-mail:
| | - Edwin Bennink
- Department of Radiology, University Medical Centre, Utrecht, The Netherlands
| | - Jan Willem Dankbaar
- Department of Radiology, University Medical Centre, Utrecht, The Netherlands
| | - Max A. Viergever
- Department of Radiology, University Medical Centre, Utrecht, The Netherlands
| | | | - Ewoud J. Smit
- Department of Radiology, University Medical Centre, Utrecht, The Netherlands
| | | |
Collapse
|
6
|
Luo T, Wischgoll T, Kwon Koo B, Huo Y, Kassab GS. IVUS validation of patient coronary artery lumen area obtained from CT images. PLoS One 2014; 9:e86949. [PMID: 24489811 PMCID: PMC3906085 DOI: 10.1371/journal.pone.0086949] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2013] [Accepted: 12/16/2013] [Indexed: 01/02/2023] Open
Abstract
AIMS Accurate computed tomography (CT)-based reconstruction of coronary morphometry (diameters, length, bifurcation angles) is important for construction of patient-specific models to aid diagnosis and therapy. The objective of this study is to validate the accuracy of patient coronary artery lumen area obtained from CT images based on intravascular ultrasound (IVUS). METHODS AND RESULTS Morphometric data of 5 patient CT scans with 11 arteries from IVUS were reconstructed including the lumen cross sectional area (CSA), diameter and length. The volumetric data from CT images were analyzed at sub-pixel accuracy to obtain accurate vessel center lines and CSA. A new center line extraction approach was used where an initial estimated skeleton in discrete value was obtained using a traditional thinning algorithm. The CSA was determined directly without any circular shape assumptions to provide accurate reconstruction of stenosis. The root-mean-square error (RMSE) for CSA and diameter were 16.2% and 9.5% respectively. CONCLUSIONS The image segmentation and CSA extraction algorithm for reconstruction of coronary arteries proved to be accurate for determination of vessel lumen area. This approach provides fundamental morphometric data for patient-specific models to diagnose and treat coronary artery disease.
Collapse
Affiliation(s)
- Tong Luo
- Department of Biomedical Engineering, Indiana University Purdue University Indianapolis, Indianapolis, Indiana, United States of America
| | - Thomas Wischgoll
- Department of Computer Science and Engineering, Wright State University, Fairborn, Ohio, United States of America
| | - Bon Kwon Koo
- Department of Internal Medicine, Seoul National University College of Medicine, Seoul, South Korea
| | - Yunlong Huo
- Department of Biomedical Engineering, Indiana University Purdue University Indianapolis, Indianapolis, Indiana, United States of America
| | - Ghassan S. Kassab
- Department of Biomedical Engineering, Indiana University Purdue University Indianapolis, Indianapolis, Indiana, United States of America
- Department of Surgery, Indiana University Purdue University Indianapolis, Indianapolis, Indiana, United States of America
- Department of Cellular and Integrative Physiology, Indiana University Purdue University Indianapolis, Indianapolis, Indiana, United States of America
| |
Collapse
|
7
|
Riordan AJ, Bennink E, Viergever MA, Velthuis BK, Dankbaar JW, de Jong HWAM. CT brain perfusion protocol to eliminate the need for selecting a venous output function. AJNR Am J Neuroradiol 2013; 34:1353-8. [PMID: 23370476 DOI: 10.3174/ajnr.a3397] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
BACKGROUND AND PURPOSE In CTP, an arterial input function is used for cerebral blood volume measurement. AIFs are often influenced by partial volume effects resulting in overestimated CBV. A venous output function is manually selected to correct for partial volume. This can introduce variability. Our goal was to develop a CTP protocol that enables AIF selection unaffected by partial volume. MATERIALS AND METHODS First, the effects of partial volume on artery sizes/types including the MCA were estimated by using a CTP phantom with 9 protocols (section thicknesses of 1, 1.8, and 5 mm and image resolutions of 0.5, 1, and 1.5 mm). Next, these protocols were applied to clinical CTP studies from 6 patients. The influence of the partial volume effect was measured by comparison of the time-attenuation curves from different artery locations with reference veins. RESULTS AIFs from MCAs were unaffected by partial volume effects when using high image resolution (1 mm) and medium section thickness (1.8 mm). For the clinical data, a total of 104 arteries and 60 veins was selected. The data confirmed that high image resolution and thin section thickness enable selection of MCAs for AIFs free of partial volume influences. In addition, we found that large veins were not insusceptible to partial volume effects relative to large arteries, questioning the use of veins for partial volume correction. CONCLUSIONS A CTP protocol with 1.8-mm section thickness and 1-mm image resolution allows AIF selection unaffected by partial volume effects in MCAs.
Collapse
Affiliation(s)
- A J Riordan
- Department of Radiology, University Medical Center Utrecht, Utrecht, the Netherlands.
| | | | | | | | | | | |
Collapse
|
8
|
New automated Markov-Gibbs random field based framework for myocardial wall viability quantification on agent enhanced cardiac magnetic resonance images. Int J Cardiovasc Imaging 2011; 28:1683-98. [PMID: 22160668 DOI: 10.1007/s10554-011-9991-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/04/2011] [Accepted: 11/29/2011] [Indexed: 10/14/2022]
Abstract
A novel automated framework for detecting and quantifying viability from agent enhanced cardiac magnetic resonance images is proposed. The framework identifies the pathological tissues based on a joint Markov-Gibbs random field (MGRF) model that accounts for the 1st-order visual appearance of the myocardial wall (in terms of the pixel-wise intensities) and the 2nd-order spatial interactions between pixels. The pathological tissue is quantified based on two metrics: the percentage area in each segment with respect to the total area of the segment, and the trans-wall extent of the pathological tissue. This transmural extent is estimated using point-to-point correspondences based on a Laplace partial differential equation. Transmural extent was validated using a simulated phantom. We tested the proposed framework on 14 datasets (168 images) and validated against manual expert delineation of the pathological tissue by two observers. Mean Dice similarity coefficients (DSC) of 0.90 and 0.88 were obtained for the observers, approaching the ideal value, 1. The Bland-Altman statistic of infarct volumes estimated by manual versus the MGRF estimation revealed little bias difference, and most values fell within the 95% confidence interval, suggesting very good agreement. Using the DSC measure we documented statistically significant superior segmentation performance for our MGRF method versus established intensity-based methods (greater DSC, and smaller standard deviation). Our Laplace method showed good operating characteristics across the full range of extent of transmural infarct, outperforming conventional methods. Phantom validation and experiments on patient data confirmed the robustness and accuracy of the proposed framework.
Collapse
|
9
|
King DM, Fagan AJ, Moran CM, Browne JE. Comparative imaging study in ultrasound, MRI, CT, and DSA using a multimodality renal artery phantom. Med Phys 2011; 38:565-73. [DOI: 10.1118/1.3533674] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
10
|
Gratama van Andel HAF, Venema HW, Bol K, Marquering HA, Majoie CB, den Heeten GJ, Grimbergen CA, Streekstra GJ. Model-based measurements of the diameter of the internal carotid artery in CT angiography images. Med Phys 2010; 37:5711-27. [DOI: 10.1118/1.3491808] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
11
|
Kline TL, Zamir M, Ritman EL. Accuracy of microvascular measurements obtained from micro-CT images. Ann Biomed Eng 2010; 38:2851-64. [PMID: 20458628 DOI: 10.1007/s10439-010-0058-7] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2009] [Accepted: 04/26/2010] [Indexed: 10/19/2022]
Abstract
Early changes in branching geometry of microvasculature and its associated impact on the perfusion distribution in diseases, especially those in which different branching generations are affected differently, require the ability to analyze intact vascular trees over a wide range of scales. Micro-CT offers an excellent framework to analyze the microvascular branching geometry. Such an analysis requires methods to be developed that can accurately characterize branching properties, such as branch diameter, length, branching angle, and branch interconnectivity of the microvasculature. The purpose of this article is to report the results of a study of two human intramyocardial coronary vascular tree casts in which the accuracy of micro-CT vascular imaging and its analysis are tested against measurements made through an optical microscope (used as the "gold-standard"). Methods related to image segmentation of the vascular lumen, vessel tree centerline extraction, individual branch segment measurement, and compensating for the non-ideal modulation transfer function of micro-CT scanners are presented. The extracted centerline accurately characterized the hierarchical structure of the vascular tree casts in terms of "parent-branch" relationships which allowed each interbranch segments' dimensions to be compared to the optical measurement method. The comparison results show a close to ideal 1:1 relationship for both length and diameter measurements made by the two methods. Combining the results from both specimens, the standard deviation of the difference between measurement methods was 19 microm for the measurement of interbranch segment diameters (ranging from 12 to 769 microm), and 172 microm for the measurement of interbranch segment lengths (ranging from 14 to 3252 microm). These results suggest that our micro-CT image analysis method can be used to characterize a vascular tree's hierarchical structure, and accurately measure interbranch segment lengths and diameters.
Collapse
Affiliation(s)
- Timothy L Kline
- Department of Physiology and Biomedical Engineering, Mayo Clinic College of Medicine, Rochester, MN 55905, USA
| | | | | |
Collapse
|
12
|
Lin M, Marshall CT, Qi Y, Johnston SM, Badea CT, Piantadosi CA, Johnson GA. Quantitative blood flow measurements in the small animal cardiopulmonary system using digital subtraction angiography. Med Phys 2010; 36:5347-58. [PMID: 19994543 DOI: 10.1118/1.3231823] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
PURPOSE The use of preclinical rodent models of disease continues to grow because these models help elucidate pathogenic mechanisms and provide robust test beds for drug development. Among the major anatomic and physiologic indicators of disease progression and genetic or drug modification of responses are measurements of blood vessel caliber and flow. Moreover, cardiopulmonary blood flow is a critical indicator of gas exchange. Current methods of measuring cardiopulmonary blood flow suffer from some or all of the following limitations--they produce relative values, are limited to global measurements, do not provide vasculature visualization, are not able to measure acute changes, are invasive, or require euthanasia. METHODS In this study, high-spatial and high-temporal resolution x-ray digital subtraction angiography (DSA) was used to obtain vasculature visualization, quantitative blood flow in absolute metrics (ml/min instead of arbitrary units or velocity), and relative blood volume dynamics from discrete regions of interest on a pixel-by-pixel basis (100 x 100 microm2). RESULTS A series of calibrations linked the DSA flow measurements to standard physiological measurement using thermodilution and Fick's method for cardiac output (CO), which in eight anesthetized Fischer-344 rats was found to be 37.0 +/- 5.1 ml/min. Phantom experiments were conducted to calibrate the radiographic density to vessel thickness, allowing a link of DSA cardiac output measurements to cardiopulmonary blood flow measurements in discrete regions of interest. The scaling factor linking relative DSA cardiac output measurements to the Fick's absolute measurements was found to be 18.90 x CODSA = COFick. CONCLUSIONS This calibrated DSA approach allows repeated simultaneous visualization of vasculature and measurement of blood flow dynamics on a regional level in the living rat.
Collapse
Affiliation(s)
- MingDe Lin
- Department of Radiology, Center for In Vivo Microscopy, Duke University Medical Center, Box 3302, Durham, North Carolina 27710, USA
| | | | | | | | | | | | | |
Collapse
|
13
|
Clinical evaluation of angiographic multiple-view 3D reconstruction. Int J Comput Assist Radiol Surg 2009; 4:497-508. [PMID: 20033533 DOI: 10.1007/s11548-009-0361-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2009] [Accepted: 05/11/2009] [Indexed: 10/20/2022]
Abstract
PURPOSE Endovascular treatment of stroke, a leading cause of death in the United States, involves guidance of devices to the intervention site often through tortuous vessels. Typically, these interventions are performed under two- dimensional (2D) fluoroscopy. To facilitate these procedures, we developed and previously presented a multiple-view self-calibration method involving two steps: (1) calibration of the imaging geometry, and (2) reconstruction of the 3D vessel centerline. Only those 2D angiograms obtained during the procedure are used for reconstruction. In this manuscript, we evaluate this technique on a large set (117 cases) of clinical data obtained over a 12-month period. METHODS We evaluated the technique using (1) the RMS difference between the calculated 3D centerlines and the average centerline (before and after application of our method), (2) the difference between the projected 3D centerlines and the 2D indicated centerlines, (3) the translations and rotations calculated by our technique, and (4) intra- and inter-user variations. RESULTS Our approach (1) reduces the RMS 3D differences by a factor of 10, (2) increases the number of projected 3D centerline points lying within 1 mm of the indicated 2D centerline points by over a factor of 2 (from 28 to 71%), (3) provides an assessment of the variations in the gantry geometry as provided by the imaging system, and (4) is insensitive to user variations in indication (<1 mm differences in 3D are seen). CONCLUSIONS These results indicate that this technique will provide more reliable vessel centerlines in the clinical setting without requiring additional acquisitions or increasing dose to the patient.
Collapse
|
14
|
Singh V, Mukherjee L, Dinu PM, Xu J, Hoffmann KR. Limited view CT reconstruction and segmentation via constrained metric labeling. COMPUTER VISION AND IMAGE UNDERSTANDING : CVIU 2008; 112:67-80. [PMID: 19802346 PMCID: PMC2707032 DOI: 10.1016/j.cviu.2008.06.005] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
This paper proposes a new discrete optimization framework for tomographic reconstruction and segmentation of CT volumes when only a few projection views are available. The problem has important clinical applications in coronary angiographic imaging. We first show that the limited view reconstruction and segmentation problem can be formulated as a "constrained" version of the metric labeling problem. This lays the groundwork for a linear programming framework that brings metric labeling classification and classical algebraic tomographic reconstruction (ART) together in a unified model. If the imaged volume is known to be comprised of a finite set of attenuation coefficients (a realistic assumption), given a regular limited view reconstruction, we view it as a task of voxels reassignment subject to maximally maintaining consistency with the input reconstruction and the objective of ART simultaneously. The approach can reliably reconstruct (or segment) volumes with several multiple contrast objects. We present evaluations using experiments on cone beam computed tomography.
Collapse
Affiliation(s)
- Vikas Singh
- Biostatistics & Medical Informatics and Computer Sciences, UW-Madison
| | | | | | | | | |
Collapse
|
15
|
Schmidt MA, Nayak SL, Belli AM, Britten AJ. Accurate geometric calibration in stepping-table digital subtraction angiography. Br J Radiol 2007; 80:835-40. [PMID: 17875596 DOI: 10.1259/bjr/41664978] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022] Open
Abstract
Accurate measurements of vessel dimensions are desirable in many clinical applications. This work uses the known relative motion between X-ray source and the patient in stepping-table digital subtraction angiography (DSA) to provide an accurate geometric calibration for quantitative measurements. The method results in a calibration factor that converts the size of the object measured in pixels on the image to its size in millimetres. The main sources of error relate to: (i) the assessment of relative displacement of a structure in a series of images; (ii) patient motion throughout data acquisition; and (iii) image distortion. Error was evaluated both with a test object consisting of a large grid of ball bearings (2x2 cm spaced) and, in vivo, in five renal DSA examinations performed with identical catheters of known diameter. The calibration factor was calculated with 0.1% accuracy for the test object and at least 2% accuracy in vivo, even with breath holding and pulsatile motion. This demonstrates that the calculation of the calibration factor can be very accurate, and that the method we propose is capable of the submillimetre accuracy required for clinical studies if used in conjunction with an accurate measurement of the vessel size in pixels. In conclusion, accurate geometric measurements can be performed in stepping-table DSA, without the need for external reference objects.
Collapse
Affiliation(s)
- M A Schmidt
- Department of Medical Physics, St George's Hospital, London SW17 0QT, UK.
| | | | | | | |
Collapse
|
16
|
Streekstra GJ, Strackee SD, Maas M, ter Wee R, Venema HW. Model-based cartilage thickness measurement in the submillimeter range. Med Phys 2007; 34:3562-70. [DOI: 10.1118/1.2766759] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
17
|
Hyde DE, Habets DF, Fox AJ, Gulka I, Kalapos P, Lee DH, Pelz DM, Holdsworth DW. Comparison of maximum intensity projection and digitally reconstructed radiographic projection for carotid artery stenosis measurement. Med Phys 2007; 34:2968-74. [PMID: 17822005 DOI: 10.1118/1.2747305] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
Digital subtraction angiography is being supplanted by three-dimensional imaging techniques in many clinical applications, leading to extensive use of maximum intensity projection (MIP) images to depict volumetric vascular data. The MIP algorithm produces intensity profiles that are different than conventional angiograms, and can also increase the vessel-to-tissue contrast-to-noise ratio. We evaluated the effect of the MIP algorithm in a clinical application where quantitative vessel measurement is important: internal carotid artery stenosis grading. Three-dimensional computed rotational angiography (CRA) was performed on 26 consecutive symptomatic patients to verify an internal carotid artery stenosis originally found using duplex ultrasound. These volumes of data were visualized using two different postprocessing projection techniques: MIP and digitally reconstructed radiographic (DRR) projection. A DRR is a radiographic image simulating a conventional digitally subtracted angiogram, but it is derived computationally from the same CRA dataset as the MIP. By visualizing a single volume with two different projection techniques, the postprocessing effect of the MIP algorithm is isolated. Vessel measurements were made, according to the NASCET guidelines, and percentage stenosis grades were calculated. The paired t-test was used to determine if the measurement difference between the two techniques was statistically significant. The CRA technique provided an isotropic voxel spacing of 0.38 mm. The MIPs and DRRs had a mean signal-difference-to-noise-ratio of 30:1 and 26:1, respectively. Vessel measurements from MIPs were, on average, 0.17 mm larger than those from DRRs (P < 0.0001). The NASCET-type stenosis grades tended to be underestimated on average by 2.4% with the MIP algorithm, although this was not statistically significant (P=0.09). The mean interobserver variability (standard deviation) of both the MIP and DRR images was 0.35 mm. It was concluded that the MIP algorithm slightly increased the apparent dimensions of the arteries, when applied to these intra-arterial CRA images. This subpixel increase was smaller than both the voxel size and interobserver variability, and was therefore not clinically relevant.
Collapse
Affiliation(s)
- Derek E Hyde
- Robarts Research Institute, Imaging Research Laboratories, 100 Perth Drive, PO Box 5015, London, Ontario N6A 5K8, Canada.
| | | | | | | | | | | | | | | |
Collapse
|
18
|
Kuhls AT, Patel V, Ionita C, Noël PB, Walczak AM, Rangwala HS, Hoffmann KR, Rudin S. New microangiography system development providing improved small vessel imaging, increased contrast to noise ratios, and multi-view 3D reconstructions. PROCEEDINGS OF SPIE--THE INTERNATIONAL SOCIETY FOR OPTICAL ENGINEERING 2006; 6142. [PMID: 21311731 DOI: 10.1117/12.653654] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Abstract
A new microangiographic system (MA) integrated into a c-arm gantry has been developed allowing precise placement of a MA at the exact same angle as the standard x-ray image intensifier (II) with unchanged source and object position. The MA can also be arbitrarily moved about the object and easily moved into the field of view (FOV) in front of the lower resolution II when higher resolution angiographic sequences are needed. The benefits of this new system are illustrated in a neurovascular study, where a rabbit is injected with contrast media for varying oblique angles. Digital subtraction angiographic (DSA) images were obtained and compared using both the MA and II detectors for the same projection view. Vessels imaged with the MA appear sharper with smaller vessels visualized. Visualization of ~100 μm vessels was possible with the MA whereas not with the II. Further, the MA could better resolve vessel overlap. Contrast to noise ratios (CNR) were calculated for vessels of varying sizes for the MA versus the II and were found to be similar for large vessels, approximately double for medium vessels, and infinitely better for the smallest vessels. In addition, a 3D reconstruction of selected vessel segments was performed, using multiple (three) projections at oblique angles, for each detector. This new MA/II integrated system should lead to improved diagnosis and image guidance of neurovascular interventions by enabling initial guidance with the low resolution large FOV II combined with use of the high resolution MA during critical parts of diagnostic and interventional procedures.
Collapse
Affiliation(s)
- Andrew T Kuhls
- University at Buffalo (State University of New York), Toshiba Stroke Research Center, 3435 Main St., Buffalo, NY USA 14214
| | | | | | | | | | | | | | | |
Collapse
|
19
|
Janssen JP, Koning G, de Koning PJH, Bosch JG, Tuinenburg JC, Reiber JHC. A new approach to contour detection in x-ray arteriograms: the wavecontour. Invest Radiol 2005; 40:514-20. [PMID: 16024989 DOI: 10.1097/01.rli.0000170811.71023.6e] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
OBJECTIVES We sought to develop a novel approach (the Wavecontour) for the detection of contours in vascular x-ray images, designed to eliminate any systematic underestimation or overestimation for vessel sizes in the range of 0.5 to 15 mm and further minimize the influence of the user-defined start points and end points. MATERIALS AND METHODS This method is based on the Wavefront Propagation principle in a 2-stage approach. Two validation experiments were performed: a Plexiglas phantom study (tube sizes ranging from 0.51 to 9.9 mm) and an in vivo patient study (114 patients with various degrees of stenosis). RESULTS The phantom study demonstrated an accuracy of 0.007 mm and a precision of 0.072 mm. The patient study showed a high similarity between the detected and the expert-drawn contours: 93% for a threshold of 1.0 pixel and 81% for a threshold of 0.5 pixels. Furthermore, the contours are robust in complex lesions and are almost independent in the middle part of the segment from the user-defined start point and end point. A variation of only 0.6 pixels exists in the middle 60% of the contours. CONCLUSIONS Our new Wavecontour approach performs very well on phantom images as well as on clinical data over the whole range of 0.5 to 15 mm and results in more robust QCA/QVA analyses.
Collapse
Affiliation(s)
- Johannes P Janssen
- From the Department of Radiology, Leiden University Medical Center, PO Box 9600, 2300 RC, Leiden, The Netherlands.
| | | | | | | | | | | |
Collapse
|
20
|
Divani AA, Tholany CR, Siddiqui AU, AlKawi A, Hussain MS, Kirmani JF, Qureshi AI. Comparison of 1- and 2-Marker Techniques for Calculating System Magnification Factor for Angiographic Measurement of Intracranial Vessels. J Neuroimaging 2005. [DOI: 10.1111/j.1552-6569.2005.tb00336.x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022] Open
|
21
|
Amado LC, Gerber BL, Gupta SN, Rettmann DW, Szarf G, Schock R, Nasir K, Kraitchman DL, Lima JAC. Accurate and objective infarct sizing by contrast-enhanced magnetic resonance imaging in a canine myocardial infarction model. J Am Coll Cardiol 2005; 44:2383-9. [PMID: 15607402 DOI: 10.1016/j.jacc.2004.09.020] [Citation(s) in RCA: 364] [Impact Index Per Article: 19.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/16/2004] [Accepted: 09/03/2004] [Indexed: 01/12/2023]
Abstract
OBJECTIVES To identify an accurate and reproducible method to define myocardial infarct (MI) size, we conducted a study in a closed-chest canine model of acute myocardial infarction, in which MI size was measured using different thresholding techniques and by imaging at different delay times after contrast administration. BACKGROUND The MI size by contrast-enhanced magnetic resonance imaging (CE-MRI) is directly related to long-term prognosis. However, previous measurements were done using nonuniform methods and tended to overestimate nonviable areas. METHODS Thirteen animals underwent 90 min of coronary artery occlusion, followed by reperfusion. The CE-MRI data were acquired within 24 h after reperfusion and compared with triphenyltetrazolium chloride pathology. In the first nine animals, images were obtained approximately 15 min after gadolinium diethylene triamine penta-acetic acid (Gd-DTPA) using an inversion-recovery gradient-echo pulse sequence. To identify the most accurate method, MI size by CE-MRI was measured visually and by semi-automatic thresholding techniques, using different criteria. In four additional animals, images were acquired every 6 min until 30 min after Gd-DTPA. RESULTS Postmortem MI size was 13.5 +/- 2.6% of left ventricular volume. Semi-automatic techniques, using full-width at half-maximum (FWHM) criterion, correlated best with postmortem data (r(2) = 0.94, p < 0.001; results confirmed by Bland-Altman plots). Using FWHM, there was no difference in MI size between different delay times after contrast (15.2 +/- 2.9% to 14.5 +/- 4.2% at 6 and 30 min, respectively; p = NS). CONCLUSIONS When an objective technique is used to define MI size by CE-MRI, accurate infarct size measurements can be obtained from images obtained up to 30 min after contrast administration.
Collapse
Affiliation(s)
- Luciano C Amado
- Department of Medicine, Division of Cardiology, Johns Hopkins University School of Medicine, 600 N. Wolfe Street, Baltimore, MD 2187, USA
| | | | | | | | | | | | | | | | | |
Collapse
|
22
|
Hoffmann KR, Dmochowski J, Nazareth DP, Miskolczi L, Nemes B, Gopal A, Wang Z, Rudin S, Bednarek DR. Vessel size measurements in angiograms: manual measurements. Med Phys 2003; 30:681-8. [PMID: 12722820 DOI: 10.1118/1.1562491] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
Vessel size measurement is perhaps the most often performed quantitative analysis in diagnostic and interventional angiography. Although automated vessel sizing techniques are generally considered to have good accuracy and precision, we have observed that clinicians rarely use these techniques in standard clinical practice, choosing to indicate the edges of vessels and catheters to determine sizes and calibrate magnifications, i.e., manual measurements. Thus, we undertook an investigation of the accuracy and precision of vessel sizes calculated from manually indicated edges of vessels. Manual measurements were performed by three neuroradiologists and three physicists. Vessel sizes ranged from 0.1-3.0 mm in simulation studies and 0.3-6.4 mm in phantom studies. Simulation resolution functions had full-widths-at-half-maximum (FWHM) ranging from 0.0 to 0.5 mm. Phantom studies were performed with 4.5 in., 6 in., 9 in., and 12 in. image intensifier modes, magnification factor = 1, with and without zooming. The accuracy and reproducibility of the measurements ranged from 0.1 to 0.2 mm, depending on vessel size, resolution, and pixel size, and zoom. These results indicate that manual measurements may have accuracies comparable to automated techniques for vessels with sizes greater than 1 mm, but that automated techniques which take into account the resolution function should be used for vessels with sizes smaller than 1 mm.
Collapse
Affiliation(s)
- Kenneth R Hoffmann
- Department of Neurosurgery, Toshiba Stroke Research Center, University at Buffalo, Buffalo, New York 14214-3025, USA.
| | | | | | | | | | | | | | | | | |
Collapse
|