1
|
De A, Das N, Saha PK, Comellas A, Hoffman E, Basu S, Chakraborti T. MSO-GP: 3-D segmentation of large and complex conjoined tree structures. Methods 2024; 229:9-16. [PMID: 38838947 DOI: 10.1016/j.ymeth.2024.05.016] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2023] [Revised: 05/03/2024] [Accepted: 05/29/2024] [Indexed: 06/07/2024] Open
Abstract
Robust segmentation of large and complex conjoined tree structures in 3-D is a major challenge in computer vision. This is particularly true in computational biology, where we often encounter large data structures in size, but few in number, which poses a hard problem for learning algorithms. We show that merging multiscale opening with geodesic path propagation, can shed new light on this classic machine vision challenge, while circumventing the learning issue by developing an unsupervised visual geometry approach (digital topology/morphometry). The novelty of the proposed MSO-GP method comes from the geodesic path propagation being guided by a skeletonization of the conjoined structure that helps to achieve robust segmentation results in a particularly challenging task in this area, that of artery-vein separation from non-contrast pulmonary computed tomography angiograms. This is an important first step in measuring vascular geometry to then diagnose pulmonary diseases and to develop image-based phenotypes. We first present proof-of-concept results on synthetic data, and then verify the performance on pig lung and human lung data with less segmentation time and user intervention needs than those of the competing methods.
Collapse
Affiliation(s)
- Arijit De
- Department of Electronics & Telecommunication Engineering, Jadavpur University, Kolkata, India.
| | - Nirmal Das
- Department of Computer Science and Engineering (AIML), Institute of Engineering and Management, Kolkata, India; Department of Computer Science and Engineering, Jadavpur University, Kolkata, India.
| | - Punam K Saha
- Department of Electrical and Computer Engineering & Department of Radiology, University of Iowa, Iowa City, IA 52242, USA.
| | | | - Eric Hoffman
- Department of Internal Medicine, University of Iowa, Iowa City, USA.
| | - Subhadip Basu
- Department of Computer Science and Engineering (AIML), Institute of Engineering and Management, Kolkata, India.
| | - Tapabrata Chakraborti
- Health Sciences Programme, The Alan Turing Institute, London, UK; Department of Medical Physics and Biomedical Engineering, University College London, London, UK.
| |
Collapse
|
2
|
Shape-appearance constrained segmentation and separation of vein and artery in pulsatile tinnitus patients based on MR angiography and flow MRI. Magn Reson Imaging 2019; 61:187-195. [DOI: 10.1016/j.mri.2019.05.026] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2019] [Revised: 05/03/2019] [Accepted: 05/19/2019] [Indexed: 11/21/2022]
|
3
|
Sultana S, Blatt JE, Gilles B, Rashid T, Audette MA. MRI-Based Medial Axis Extraction and Boundary Segmentation of Cranial Nerves Through Discrete Deformable 3D Contour and Surface Models. IEEE TRANSACTIONS ON MEDICAL IMAGING 2017; 36:1711-1721. [PMID: 28422682 DOI: 10.1109/tmi.2017.2693182] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
This paper presents a segmentation technique to identify the medial axis and the boundary of cranial nerves. We utilize a 3-D deformable one-simplex discrete contour model to extract the medial axis of each cranial nerve. This contour model represents a collection of two-connected vertices linked by edges, where vertex position is determined by a Newtonian expression for vertex kinematics featuring internal and external forces, the latter of which include attractive forces toward the nerve medial axis. We exploit multiscale vesselness filtering and minimal path techniques in the medial axis extraction method, which also computes a radius estimate along the path. Once we have the medial axis and the radius function of a nerve, we identify the nerve surface using a two-simplex deformable model, which expands radially and can accommodate any nerve shape. As a result, the method proposed here combines the benefits of explicit contour and surface models, while also achieving a cornerstone for future work that will emphasize shape statistics, static collision with other critical structures, and tree-shape analysis.
Collapse
|
4
|
Payer C, Pienn M, Bálint Z, Shekhovtsov A, Talakic E, Nagy E, Olschewski A, Olschewski H, Urschler M. Automated integer programming based separation of arteries and veins from thoracic CT images. Med Image Anal 2016; 34:109-122. [PMID: 27189777 DOI: 10.1016/j.media.2016.05.002] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2016] [Revised: 04/07/2016] [Accepted: 05/03/2016] [Indexed: 10/24/2022]
Abstract
Automated computer-aided analysis of lung vessels has shown to yield promising results for non-invasive diagnosis of lung diseases. To detect vascular changes which affect pulmonary arteries and veins differently, both compartments need to be identified. We present a novel, fully automatic method that separates arteries and veins in thoracic computed tomography images, by combining local as well as global properties of pulmonary vessels. We split the problem into two parts: the extraction of multiple distinct vessel subtrees, and their subsequent labeling into arteries and veins. Subtree extraction is performed with an integer program (IP), based on local vessel geometry. As naively solving this IP is time-consuming, we show how to drastically reduce computational effort by reformulating it as a Markov Random Field. Afterwards, each subtree is labeled as either arterial or venous by a second IP, using two anatomical properties of pulmonary vessels: the uniform distribution of arteries and veins, and the parallel configuration and close proximity of arteries and bronchi. We evaluate algorithm performance by comparing the results with 25 voxel-based manual reference segmentations. On this dataset, we show good performance of the subtree extraction, consisting of very few non-vascular structures (median value: 0.9%) and merged subtrees (median value: 0.6%). The resulting separation of arteries and veins achieves a median voxel-based overlap of 96.3% with the manual reference segmentations, outperforming a state-of-the-art interactive method. In conclusion, our novel approach provides an opportunity to become an integral part of computer aided pulmonary diagnosis, where artery/vein separation is important.
Collapse
Affiliation(s)
- Christian Payer
- Institute for Computer Graphics and Vision, Graz University of Technology, Austria; Ludwig Boltzmann Institute for Lung Vascular Research, Graz, Austria
| | - Michael Pienn
- Ludwig Boltzmann Institute for Lung Vascular Research, Graz, Austria
| | - Zoltán Bálint
- Ludwig Boltzmann Institute for Lung Vascular Research, Graz, Austria
| | | | - Emina Talakic
- Division of General Radiology, Department of Radiology, Medical University of Graz, Austria
| | - Eszter Nagy
- Division of Pediatric Radiology, Department of Radiology, Medical University of Graz, Austria
| | - Andrea Olschewski
- Ludwig Boltzmann Institute for Lung Vascular Research, Graz, Austria; Experimental Anesthesiology, Department of Anesthesia and Intensive Care Medicine, Medical University of Graz, Austria
| | - Horst Olschewski
- Ludwig Boltzmann Institute for Lung Vascular Research, Graz, Austria; Division of Pulmonology, Department of Internal Medicine, Medical University of Graz, Austria
| | - Martin Urschler
- Institute for Computer Graphics and Vision, Graz University of Technology, Austria; Ludwig Boltzmann Institute for Clinical Forensic Imaging, Graz, Austria; BioTechMed Graz, Austria.
| |
Collapse
|
5
|
Ilicak E, Cetin S, Bulut E, Oguz KK, Saritas EU, Unal G, Çukur T. Targeted vessel reconstruction in non-contrast-enhanced steady-state free precession angiography. NMR IN BIOMEDICINE 2016; 29:532-544. [PMID: 26854004 DOI: 10.1002/nbm.3497] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/15/2015] [Revised: 01/07/2016] [Accepted: 01/08/2016] [Indexed: 06/05/2023]
Abstract
Image quality in non-contrast-enhanced (NCE) angiograms is often limited by scan time constraints. An effective solution is to undersample angiographic acquisitions and to recover vessel images with penalized reconstructions. However, conventional methods leverage penalty terms with uniform spatial weighting, which typically yield insufficient suppression of aliasing interference and suboptimal blood/background contrast. Here we propose a two-stage strategy where a tractographic segmentation is employed to auto-extract vasculature maps from undersampled data. These maps are then used to incur spatially adaptive sparsity penalties on vascular and background regions. In vivo steady-state free precession angiograms were acquired in the hand, lower leg and foot. Compared with regular non-adaptive compressed sensing (CS) reconstructions (CSlow ), the proposed strategy improves blood/background contrast by 71.3 ± 28.9% in the hand (mean ± s.d. across acceleration factors 1-8), 30.6 ± 11.3% in the lower leg and 28.1 ± 7.0% in the foot (signed-rank test, P < 0.05 at each acceleration). The proposed targeted reconstruction can relax trade-offs between image contrast, resolution and scan efficiency without compromising vessel depiction.
Collapse
Affiliation(s)
- Efe Ilicak
- Department of Electrical and Electronics Engineering and the National Magnetic Resonance Research Center, Bilkent University, Ankara, Turkey
| | - Suheyla Cetin
- Faculty of Engineering and Natural Sciences, Sabanci University, Istanbul, Turkey
| | - Elif Bulut
- Hacettepe University Faculty of Medicine, Ankara, Turkey
| | | | - Emine Ulku Saritas
- Department of Electrical and Electronics Engineering and the National Magnetic Resonance Research Center, Bilkent University, Ankara, Turkey
| | - Gozde Unal
- Department of Computer Engineering, Istanbul Technical University, Istanbul, Turkey
| | - Tolga Çukur
- Department of Electrical and Electronics Engineering and the National Magnetic Resonance Research Center, Bilkent University, Ankara, Turkey
| |
Collapse
|
6
|
Data-Dependent Higher-Order Clique Selection for Artery–Vein Segmentation by Energy Minimization. Int J Comput Vis 2015. [DOI: 10.1007/s11263-015-0856-3] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
7
|
Rudyanto RD, Kerkstra S, van Rikxoort EM, Fetita C, Brillet PY, Lefevre C, Xue W, Zhu X, Liang J, Öksüz I, Ünay D, Kadipaşaoğlu K, Estépar RSJ, Ross JC, Washko GR, Prieto JC, Hoyos MH, Orkisz M, Meine H, Hüllebrand M, Stöcker C, Mir FL, Naranjo V, Villanueva E, Staring M, Xiao C, Stoel BC, Fabijanska A, Smistad E, Elster AC, Lindseth F, Foruzan AH, Kiros R, Popuri K, Cobzas D, Jimenez-Carretero D, Santos A, Ledesma-Carbayo MJ, Helmberger M, Urschler M, Pienn M, Bosboom DGH, Campo A, Prokop M, de Jong PA, Ortiz-de-Solorzano C, Muñoz-Barrutia A, van Ginneken B. Comparing algorithms for automated vessel segmentation in computed tomography scans of the lung: the VESSEL12 study. Med Image Anal 2014; 18:1217-32. [PMID: 25113321 DOI: 10.1016/j.media.2014.07.003] [Citation(s) in RCA: 81] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2013] [Revised: 03/01/2014] [Accepted: 07/01/2014] [Indexed: 10/25/2022]
Abstract
The VESSEL12 (VESsel SEgmentation in the Lung) challenge objectively compares the performance of different algorithms to identify vessels in thoracic computed tomography (CT) scans. Vessel segmentation is fundamental in computer aided processing of data generated by 3D imaging modalities. As manual vessel segmentation is prohibitively time consuming, any real world application requires some form of automation. Several approaches exist for automated vessel segmentation, but judging their relative merits is difficult due to a lack of standardized evaluation. We present an annotated reference dataset containing 20 CT scans and propose nine categories to perform a comprehensive evaluation of vessel segmentation algorithms from both academia and industry. Twenty algorithms participated in the VESSEL12 challenge, held at International Symposium on Biomedical Imaging (ISBI) 2012. All results have been published at the VESSEL12 website http://vessel12.grand-challenge.org. The challenge remains ongoing and open to new participants. Our three contributions are: (1) an annotated reference dataset available online for evaluation of new algorithms; (2) a quantitative scoring system for objective comparison of algorithms; and (3) performance analysis of the strengths and weaknesses of the various vessel segmentation methods in the presence of various lung diseases.
Collapse
Affiliation(s)
- Rina D Rudyanto
- Center for Applied Medical Research, University of Navarra, Spain.
| | - Sjoerd Kerkstra
- Diagnostic Image Analysis Group, Radboud University Nijmegen Medical Centre, The Netherlands
| | - Eva M van Rikxoort
- Diagnostic Image Analysis Group, Radboud University Nijmegen Medical Centre, The Netherlands
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | - Marius Staring
- Division of Image Processing (LKEB), Leiden University Medical Center, The Netherlands
| | | | - Berend C Stoel
- Division of Image Processing (LKEB), Leiden University Medical Center, The Netherlands
| | - Anna Fabijanska
- Institute of Applied Computer Science, Lodz University of Technology, Poland
| | - Erik Smistad
- Norwegian University of Science and Technology, Norway
| | - Anne C Elster
- Norwegian University of Science and Technology, Norway
| | | | | | | | | | | | | | - Andres Santos
- Universidad Politécnica de Madrid, Spain; CIBER-BBN, Spain
| | | | - Michael Helmberger
- Graz University of Technology, Institute for Computer Vision and Graphics, Austria
| | - Martin Urschler
- Ludwig Boltzmann Institute for Clinical Forensic Imaging, Graz, Austria
| | - Michael Pienn
- Ludwig Boltzmann Institute for Lung Vascular Research, Graz, Austria
| | - Dennis G H Bosboom
- Diagnostic Image Analysis Group, Radboud University Nijmegen Medical Centre, The Netherlands
| | - Arantza Campo
- Pulmonary Department, Clínica Universidad de Navarra, University of Navarra, Spain
| | - Mathias Prokop
- Diagnostic Image Analysis Group, Radboud University Nijmegen Medical Centre, The Netherlands
| | - Pim A de Jong
- Department of Radiology, University Medical Center, Utrecht, The Netherlands
| | | | | | - Bram van Ginneken
- Diagnostic Image Analysis Group, Radboud University Nijmegen Medical Centre, The Netherlands
| |
Collapse
|
8
|
Tian Y, Chen Q, Wang W, Peng Y, Wang Q, Duan F, Wu Z, Zhou M. A vessel active contour model for vascular segmentation. BIOMED RESEARCH INTERNATIONAL 2014; 2014:106490. [PMID: 25101262 PMCID: PMC4101240 DOI: 10.1155/2014/106490] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/04/2014] [Accepted: 06/12/2014] [Indexed: 11/30/2022]
Abstract
This paper proposes a vessel active contour model based on local intensity weighting and a vessel vector field. Firstly, the energy function we define is evaluated along the evolving curve instead of all image points, and the function value at each point on the curve is based on the interior and exterior weighted means in a local neighborhood of the point, which is good for dealing with the intensity inhomogeneity. Secondly, a vascular vector field derived from a vesselness measure is employed to guide the contour to evolve along the vessel central skeleton into thin and weak vessels. Thirdly, an automatic initialization method that makes the model converge rapidly is developed, and it avoids repeated trails in conventional local region active contour models. Finally, a speed-up strategy is implemented by labeling the steadily evolved points, and it avoids the repeated computation of these points in the subsequent iterations. Experiments using synthetic and real vessel images validate the proposed model. Comparisons with the localized active contour model, local binary fitting model, and vascular active contour model show that the proposed model is more accurate, efficient, and suitable for extraction of the vessel tree from different medical images.
Collapse
Affiliation(s)
- Yun Tian
- College of Information Science & Technology, Beijing Normal University, Beijing 100875, China
| | - Qingli Chen
- Business School, Henan Normal University, Xinxiang 453007, China
| | - Wei Wang
- Department of Obstetrics and Gynecology, Navy General Hospital, Beijing 100048, China
| | - Yu Peng
- School of Design, Communication & Information Technology, University of Newcastle, Callaghan, NSW 2308, Australia
| | - Qingjun Wang
- Department of Radiology, Navy General Hospital, Beijing 100048, China
| | - Fuqing Duan
- College of Information Science & Technology, Beijing Normal University, Beijing 100875, China
| | - Zhongke Wu
- College of Information Science & Technology, Beijing Normal University, Beijing 100875, China
| | - Mingquan Zhou
- College of Information Science & Technology, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
9
|
Rivest-Hénault D, Cheriet M. 3-D curvilinear structure detection filter via structure-ball analysis. IEEE TRANSACTIONS ON IMAGE PROCESSING : A PUBLICATION OF THE IEEE SIGNAL PROCESSING SOCIETY 2013; 22:2849-2863. [PMID: 23335669 DOI: 10.1109/tip.2013.2240005] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
Curvilinear structure detection filters are crucial building blocks in many medical image processing applications, where they are used to detect important structures, such as blood vessels, airways, and other similar fibrous tissues. Unfortunately, most of these filters are plagued by an implicit single structure direction assumption, which results in a loss of signal around bifurcations. This peculiarity limits the performance of all subsequent processes, such as understanding angiography acquisitions, computing an accurate segmentation or tractography, or automatically classifying image voxels. This paper presents a new 3-D curvilinear structure detection filter based on the analysis of the structure ball, a geometric construction representing second order differences sampled in many directions. The structure ball is defined formally, and its computation on a discreet image is discussed. A contrast invariant diffusion index easing voxel analysis and visualization is also introduced, and different structure ball shape descriptors are proposed. A new curvilinear structure detection filter is defined based on the shape descriptors that best characterize curvilinear structures. The new filter produces a vesselness measure that is robust to the presence of X- and Y-junctions along the structure by going beyond the single direction assumption. At the same time, it stays conceptually simple and deterministic, and allows for an intuitive representation of the structure's principal directions. Sample results are provided for synthetic images and for two medical imaging modalities.
Collapse
|
10
|
Filtering and segmentation of 3D angiographic data: Advances based on mathematical morphology. Med Image Anal 2013; 17:147-64. [DOI: 10.1016/j.media.2012.08.004] [Citation(s) in RCA: 53] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2011] [Revised: 07/25/2012] [Accepted: 08/20/2012] [Indexed: 11/23/2022]
|
11
|
Gao Z, Grout RW, Holtze C, Hoffman EA, Saha PK. A new paradigm of interactive artery/vein separation in noncontrast pulmonary CT imaging using multiscale topomorphologic opening. IEEE Trans Biomed Eng 2012; 59:3016-27. [PMID: 22899571 PMCID: PMC4294416 DOI: 10.1109/tbme.2012.2212894] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Distinguishing pulmonary arterial and venous (A/V) trees via in vivo imaging is a critical first step in the quantification of vascular geometry for the purpose of diagnosing several pulmonary diseases and to develop new image-based phenotypes. A multiscale topomorphologic opening (MSTMO) algorithm has recently been developed in our laboratory for separating A/V trees via noncontrast pulmonary human CT imaging. The method starts with two sets of seeds-one for each of A/V trees and combines fuzzy distance transform and fuzzy connectivity in conjunction with several morphological operations leading to locally adaptive iterative multiscale opening of two mutually conjoined structures. In this paper, we introduce the methods for handling "local update" and "separators" into our previous theoretical formulation and incorporate the algorithm into an effective graphical user interface (GUI). Results of a comprehensive evaluative study assessing both accuracy and reproducibility of the method under the new setup are presented and also, the effectiveness of the GUI-based system toward improving A/V separation results is examined. Accuracy of the method has been evaluated using mathematical phantoms, CT images of contrast-separated pulmonary A/V casting of a pig's lung and noncontrast pulmonary human CT imaging. The method has achieved 99% true A/V labeling in the cast phantom and, almost, 92-94% true labeling in human lung data. Reproducibility of the method has been evaluated using multiuser A/V separation in human CT data along with contrast-enhanced CT images of a pig's lung at different positive end-expiratory pressures (PEEPs). The method has achieved, almost, 92-98% agreements in multiuser A/V labeling with ICC for A/V measures being over 0.96-0.99. Effectiveness of the GUI-based method has been evaluated on human data in terms of improvements of accuracy of A/V separation results and results have shown 8-22% improvements in true A/V labeling. Both qualitative and quantitative results found are very promising.
Collapse
Affiliation(s)
- Zhiyun Gao
- Department of Electrical and Computer Engineering, the University of Iowa, Iowa City, IA 52242
| | - Randall W. Grout
- Department of Radiology, The University of Iowa, Iowa City, IA, 52242
| | - Colin Holtze
- Department of Radiology, The University of Iowa, Iowa City, IA, 52242
| | - Eric A. Hoffman
- Department of Radiology and the Department of Biomedical Engineering, The University of Iowa, Iowa City, IA, 52242
| | - Punam K. Saha
- Department of Electrical and Computer Engineering and the Department of Radiology, The University of Iowa, Iowa City, IA 52242, Phone: 319-335-5959, Fax: 319-335-6028
| |
Collapse
|
12
|
Soleimanifard S, Schär M, Hays AG, Weiss RG, Stuber M, Prince JL. VESSEL CENTERLINE TRACKING AND BOUNDARY SEGMENTATION IN CORONARY MRA WITH MINIMAL MANUAL INTERACTION. PROCEEDINGS. IEEE INTERNATIONAL SYMPOSIUM ON BIOMEDICAL IMAGING 2012:1417-1420. [PMID: 23032185 PMCID: PMC3461337 DOI: 10.1109/isbi.2012.6235834] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
Magnetic resonance angiography (MRA) provides a noninvasive means to detect the presence, location and severity of atherosclerosis throughout the vascular system. In such studies, and especially those in the coronary arteries, the vessel luminal area is typically measured at multiple cross-sectional locations along the course of the artery. The advent of fast volumetric imaging techniques covering proximal to mid segments of coronary arteries necessitates automatic analysis tools requiring minimal manual interactions to robustly measure cross-sectional area along the three-dimensional track of the arteries in under-sampled and non-isotropic datasets. In this work, we present a modular approach based on level set methods to track the vessel centerline, segment the vessel boundaries, and measure transversal area using two user-selected endpoints in each coronary of interest. Arterial area and vessel length are measured using our method and compared to the standard Soap-Bubble reformatting and analysis tool in in-vivo non-contrast enhanced coronary MRA images.
Collapse
Affiliation(s)
- Sahar Soleimanifard
- Department of Electrical and Computer Engineering, Johns Hopkins University, Baltimore, MD, USA
| | | | | | | | | | | |
Collapse
|
13
|
Wang Y, Hu D, Liu Y, Li M. Cerebral artery-vein separation using 0.1-Hz oscillation in dual-wavelength optical imaging. IEEE TRANSACTIONS ON MEDICAL IMAGING 2011; 30:2030-2043. [PMID: 21693415 DOI: 10.1109/tmi.2011.2160191] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
We present a novel artery-vein separation method using 0.1-Hz oscillation at two wavelengths with optical imaging of intrinsic signals (OIS). The 0.1-Hz oscillation at a green light wavelength of 546 nm exhibits greater amplitude in arteries than in veins and is primarily caused by vasomotion, whereas the 0.1-Hz oscillation at a red light wavelength of 630 nm exhibits greater amplitude in veins than in arteries and is primarily caused by changes of deoxyhemoglobin concentration. This spectral feature enables cortical arteries and veins to be segmented independently. The arteries can be segmented on the 0.1-Hz amplitude image at 546 nm using matched filters of a modified dual Gaussian model combining with a single Gaussian model. The veins are a combination of vessels segmented on both amplitude images at the two wavelengths using multiscale matched filters of single Gaussian model. Our method can separate most of the thin arteries and veins from each other, especially the thin arteries with low contrast in raw gray images. In vivo OIS experiments demonstrate the separation ability of the 0.1-Hz based segmentation method in cerebral cortex of eight rats. Two validation studies were undertaken to evaluate the performance of the method by quantifying the arterial and venous length based on a reference standard. The results indicate that our 0.1-Hz method is very effective in separating both large and thin arteries and veins regardless of vessel crossover or overlapping to great extent in comparison with previous methods.
Collapse
Affiliation(s)
- Yucheng Wang
- National University of Defense Technology, Changsha 410073, China.
| | | | | | | |
Collapse
|
14
|
Angiographic Image Analysis. ACTA ACUST UNITED AC 2011. [DOI: 10.1007/978-1-4419-9779-1_6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/17/2023]
|
15
|
Post-processing central k-space subtraction for high-resolution arterial peripheral MR angiography. Magn Reson Imaging 2011; 29:835-43. [PMID: 21546189 DOI: 10.1016/j.mri.2011.03.002] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2010] [Revised: 03/04/2011] [Accepted: 03/07/2011] [Indexed: 11/21/2022]
Abstract
Peripheral MR angiography requires high resolution and arterial contrast. Neither can be obtained simultaneously due to the short arterial phase of the contrast agent. To improve temporal resolution, keyhole imaging was developed, which combines high resolution and arterial k-spaces at the time of image acquisition. Here, a related approach is introduced for image post-processing in the Fourier domain. It is demonstrated that simple substitution of the central k-space with low-resolution data leads to severe distortion. Hence, a dedicated calculation scheme is necessary for composite k-space post-processing. A solution is presented for high-resolution arterial peripheral MR angiography that uses subtraction of venous intensities from the central high-resolution k-space. The calculations in the Fourier domain do not require interpolations between the different resolutions. High-resolution steady-state MR angiography, which exhibits contrast-enhanced arteries and veins at an isotropic resolution of 0.65 mm, and standard resolution arterial first-pass MR angiography were combined to obtain images with the resolution of the steady-state images and arterial contrast. Numerical simulations on software phantoms are presented. The operation of the method is demonstrated in five patients.
Collapse
|
16
|
Shang Y, Deklerck R, Nyssen E, Markova A, de Mey J, Yang X, Sun K. Vascular active contour for vessel tree segmentation. IEEE Trans Biomed Eng 2010; 58:1023-32. [PMID: 21138795 DOI: 10.1109/tbme.2010.2097596] [Citation(s) in RCA: 89] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
In this paper, a novel active contour model is proposed for vessel tree segmentation. First, we introduce a region competition-based active contour model exploiting the gaussian mixture model, which mainly segments thick vessels. Second, we define a vascular vector field to evolve the active contour along its center line into the thin and weak vessels. The vector field is derived from the eigenanalysis of the Hessian matrix of the image intensity in a multiscale framework. Finally, a dual curvature strategy, which uses a vesselness measure-dependent function selecting between a minimal principal curvature and a mean curvature criterion, is added to smoothen the surface of the vessel without changing its shape. The developed model is used to extract the liver and lung vessel tree as well as the coronary artery from high-resolution volumetric computed tomography images. Comparisons are made with several classical active contour models and manual extraction. The experiments show that our model is more accurate and robust than these classical models and is, therefore, more suited for automatic vessel tree extraction.
Collapse
Affiliation(s)
- Yanfeng Shang
- Department of Electronics and Informatics, Vrije Universiteit Brussel, IBBT, Brussels 1050, Belgium.
| | | | | | | | | | | | | |
Collapse
|
17
|
Robust CTA lumen segmentation of the atherosclerotic carotid artery bifurcation in a large patient population. Med Image Anal 2010; 14:759-69. [DOI: 10.1016/j.media.2010.05.001] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2009] [Revised: 05/04/2010] [Accepted: 05/04/2010] [Indexed: 11/21/2022]
|
18
|
Chiribiri A, Morton G, Nagel E. Gadofosveset injection for magnetic resonance angiography. ACTA ACUST UNITED AC 2010. [DOI: 10.2217/iim.10.39] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
19
|
Mendrik A, Vonken EJ, van Ginneken B, Smit E, Waaijer A, Bertolini G, Viergever MA, Prokop M. Automatic segmentation of intracranial arteries and veins in four-dimensional cerebral CT perfusion scans. Med Phys 2010; 37:2956-66. [DOI: 10.1118/1.3397813] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
20
|
Saha PK, Gao Z, Alford SK, Sonka M, Hoffman EA. Topomorphologic separation of fused isointensity objects via multiscale opening: separating arteries and veins in 3-D pulmonary CT. IEEE TRANSACTIONS ON MEDICAL IMAGING 2010; 29:840-51. [PMID: 20199919 PMCID: PMC4517589 DOI: 10.1109/tmi.2009.2038224] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
A novel multiscale topomorphologic approach for opening of two isointensity objects fused at different locations and scales is presented and applied to separating arterial and venous trees in 3-D pulmonary multidetector X-ray computed tomography (CT) images. Initialized with seeds, the two isointensity objects (arteries and veins) grow iteratively while maintaining their spatial exclusiveness and eventually form two mutually disjoint objects at convergence. The method is intended to solve the following two fundamental challenges: how to find local size of morphological operators and how to trace continuity of locally separated regions. These challenges are met by combining fuzzy distance transform (FDT), a morphologic feature with a topologic fuzzy connectivity, and a new morphological reconstruction step to iteratively open finer and finer details starting at large scales and progressing toward smaller scales. The method employs efficient user intervention at locations where local morphological separability assumption does not hold due to imaging ambiguities or any other reason. The approach has been validated on mathematically generated tubular objects and applied to clinical pulmonary noncontrast CT data for separating arteries and veins. The tradeoff between accuracy and the required user intervention for the method has been quantitatively examined by comparing with manual outlining. The experimental study, based on a blind seed selection strategy, has demonstrated that above 95% accuracy may be achieved using 25-40 seeds for each of arteries and veins. Our method is very promising for semiautomated separation of arteries and veins in pulmonary CT images even when there is no object-specific intensity variation at conjoining locations.
Collapse
Affiliation(s)
- Punam K. Saha
- Department of Electrical and Computer Engineering and the Department of Radiology, The University of Iowa, Iowa City, IA 52242 USA
| | - Zhiyun Gao
- Department of Electrical and Computer Engineering, The University of Iowa, Iowa City, IA 52242 USA
| | | | - Milan Sonka
- Department of Electrical and Computer Engineering, the Department of Ophthalmology and Visual Sciences, and the Department of Radiation Oncology, The University of Iowa, Iowa City, IA 52242 USA
| | - Eric A. Hoffman
- Department of Radiology and Biomedical Engineering, The University of Iowa, Iowa City, IA 52242 USA
| |
Collapse
|
21
|
Anzidei M, Napoli A, Marincola BC, Nofroni I, Geiger D, Zaccagna F, Catalano C, Passariello R. Gadofosveset-enhanced MR angiography of carotid arteries: does steady-state imaging improve accuracy of first-pass imaging? Comparison with selective digital subtraction angiography. Radiology 2009; 251:457-66. [PMID: 19401574 DOI: 10.1148/radiol.2512081197] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
PURPOSE To evaluate the diagnostic accuracy of gadofosveset-enhanced magnetic resonance (MR) angiography in the assessment of carotid artery stenosis, with digital subtraction angiography (DSA) as the reference standard, and to determine the value of reading first-pass, steady-state, and "combined" (first-pass plus steady-state) MR angiograms. MATERIALS AND METHODS This study was approved by the local ethics committee, and all subjects gave written informed consent. MR angiography and DSA were performed in 84 patients (56 men, 28 women; age range, 61-76 years) with carotid artery stenosis at Doppler ultrasonography. Three readers reviewed the first-pass, steady-state, and combined MR data sets, and one independent observer evaluated the DSA images to assess stenosis degree, plaque morphology and ulceration, stenosis length, and tandem lesions. Interobserver agreement regarding MR angiographic findings was analyzed by using intraclass correlation and Cohen kappa coefficients. Sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) were calculated by using the McNemar test to determine possible significant differences (P < .05). RESULTS Interobserver agreement regarding all MR angiogram readings was substantial. For grading stenosis, sensitivity, specificity, PPV, and NPV were, respectively, 90%, 92%, 91%, and 91% for first-pass imaging; 95% each for steady-state imaging; and 96%, 99%, 99%, and 97% for combined imaging. For evaluation of plaque morphology, respective values were 84%, 86%, 88%, and 82% for first-pass imaging; 98%, 97%, 98%, and 97% for steady-state imaging; and 98%, 100%, 100%, and 97% for combined imaging. Differences between the first-pass, steady-state, and combined image readings for assessment of stenosis degree and plaque morphology were significant (P < .001). CONCLUSION Gadofosveset-enhanced MR angiography is a promising technique for imaging carotid artery stenosis. Steady-state image reading is superior to first-pass image reading, but the combined reading protocol is more accurate.
Collapse
Affiliation(s)
- Michele Anzidei
- Department of Radiological Sciences, University of Rome La Sapienza, Viale Regina Elena 324, 00161 Rome, Italy.
| | | | | | | | | | | | | | | |
Collapse
|
22
|
Qian X, Brennan MP, Dione DP, Dobrucki WL, Jackowski MP, Breuer CK, Sinusas AJ, Papademetris X. A non-parametric vessel detection method for complex vascular structures. Med Image Anal 2009; 13:49-61. [PMID: 18678521 PMCID: PMC2614119 DOI: 10.1016/j.media.2008.05.005] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2007] [Revised: 05/26/2008] [Accepted: 05/30/2008] [Indexed: 10/21/2022]
Abstract
Modern medical imaging techniques enable the acquisition of in vivo high resolution images of the vascular system. Most common methods for the detection of vessels in these images, such as multiscale Hessian-based operators and matched filters, rely on the assumption that at each voxel there is a single cylinder. Such an assumption is clearly violated at the multitude of branching points that are easily observed in all, but the most focused vascular image studies. In this paper, we propose a novel method for detecting vessels in medical images that relaxes this single cylinder assumption. We directly exploit local neighborhood intensities and extract characteristics of the local intensity profile (in a spherical polar coordinate system) which we term as the polar neighborhood intensity profile. We present a new method to capture the common properties shared by polar neighborhood intensity profiles for all the types of vascular points belonging to the vascular system. The new method enables us to detect vessels even near complex extreme points, including branching points. Our method demonstrates improved performance over standard methods on both 2D synthetic images and 3D animal and clinical vascular images, particularly close to vessel branching regions.
Collapse
Affiliation(s)
- Xiaoning Qian
- Department of Diagnostic Radiology, Yale University, New Haven, CT, USA
| | | | | | | | | | | | - Albert J. Sinusas
- Department of Diagnostic Radiology, Yale University, New Haven, CT, USA
- Department of Medicine, Yale University, New Haven, CT, USA
| | - Xenophon Papademetris
- Department of Diagnostic Radiology, Yale University, New Haven, CT, USA
- Department of Biomedical Engineering, Yale University, New Haven, CT, USA
| |
Collapse
|
23
|
Santini F, Patil S, Meckel S, Scheffler K, Wetzel SG. Double-reference cross-correlation algorithm for separation of the arteries and veins from 3D MRA time series. J Magn Reson Imaging 2008; 28:646-54. [DOI: 10.1002/jmri.21499] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022] Open
|
24
|
Gooya A, Liao H, Matsumiya K, Masamune K, Masutani Y, Dohi T. A variational method for geometric regularization of vascular segmentation in medical images. IEEE TRANSACTIONS ON IMAGE PROCESSING : A PUBLICATION OF THE IEEE SIGNAL PROCESSING SOCIETY 2008; 17:1295-1312. [PMID: 18632340 DOI: 10.1109/tip.2008.925378] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
In this paper, a level-set-based geometric regularization method is proposed which has the ability to estimate the local orientation of the evolving front and utilize it as shape induced information for anisotropic propagation. We show that preserving anisotropic fronts can improve elongations of the extracted structures, while minimizing the risk of leakage. To that end, for an evolving front using its shape-offset level-set representation, a novel energy functional is defined. It is shown that constrained optimization of this functional results in an anisotropic expansion flow which is usefull for vessel segmentation. We have validated our method using synthetic data sets, 2-D retinal angiogram images and magnetic resonance angiography volumetric data sets. A comparison has been made with two state-of-the-art vessel segmentation methods. Quantitative results, as well as qualitative comparisons of segmentations, indicate that our regularization method is a promising tool to improve the efficiency of both techniques.
Collapse
Affiliation(s)
- Ali Gooya
- Graduate School of Engineering, The University of Tokyo, Tokyo, Japan.
| | | | | | | | | | | |
Collapse
|
25
|
Abstract
In this chapter, the basic principles of magnetic resonance imaging (MRI) and magnetic resonance spectroscopy (MRS) (Sects. 2.2, 2.3, and 2.4), the technical components of the MRI scanner (Sect. 2.5), and the basics of contrast agents and the application thereof (Sect. 2.6) are described. Furthermore, flow phenomena and MR angiography (Sect. 2.7) as well as diffusion and tensor imaging (Sect. 2.7) are elucidated.
Collapse
|
26
|
Meaney JF, Boyle G, O'Keeffe S. Contrast-enhanced magnetic resonance angiography: Current status, theoretical limitations and future potential. Radiography (Lond) 2007. [DOI: 10.1016/j.radi.2007.09.002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
27
|
Choe TE, Medioni G, Cohen I, Walsh AC, Sadda SR. 2-D registration and 3-D shape inference of the retinal fundus from fluorescein images. Med Image Anal 2007; 12:174-90. [PMID: 18060827 DOI: 10.1016/j.media.2007.10.002] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2007] [Revised: 09/13/2007] [Accepted: 10/01/2007] [Indexed: 11/15/2022]
Abstract
This study presents methods to 2-D registration of retinal image sequences and 3-D shape inference from fluorescein images. The Y-feature is a robust geometric entity that is largely invariant across modalities as well as across the temporal grey level variations induced by the propagation of the dye in the vessels. We first present a Y-feature extraction method that finds a set of Y-feature candidates using local image gradient information. A gradient-based approach is then used to align an articulated model of the Y-feature to the candidates more accurately while optimizing a cost function. Using mutual information, fitted Y-features are subsequently matched across images, including colors and fluorescein angiographic frames, for registration. To reconstruct the retinal fundus in 3-D, the extracted Y-features are used to estimate the epipolar geometry with a plane-and-parallax approach. The proposed solution provides a robust estimation of the fundamental matrix suitable for plane-like surfaces, such as the retinal fundus. The mutual information criterion is used to accurately estimate the dense disparity map. Our experimental results validate the proposed method on a set of difficult fluorescein image pairs.
Collapse
Affiliation(s)
- Tae Eun Choe
- Institute for Robotics and Intelligent Systems, University of Southern California, 3737 Watt way, Los Angeles, CA 90248, USA.
| | | | | | | | | |
Collapse
|
28
|
Li H, Yezzi A. Vessels as 4-D curves: global minimal 4-D paths to extract 3-D tubular surfaces and centerlines. IEEE TRANSACTIONS ON MEDICAL IMAGING 2007; 26:1213-23. [PMID: 17896594 DOI: 10.1109/tmi.2007.903696] [Citation(s) in RCA: 63] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
In this paper, we propose an innovative approach to the segmentation of tubular structures. This approach combines all of the benefits of minimal path techniques such as global minimizers, fast computation, and powerful incorporation of user input, while also having the capability to represent and detect vessel surfaces directly which so far has been a feature restricted to active contour and surface techniques. The key is to represent the trajectory of a tubular structure not as a 3-D curve but to go up a dimension and represent the entire structure as a 4-D curve. Then we are able to fully exploit minimal path techniques to obtain global minimizing trajectories between two user supplied endpoints in order to reconstruct tubular structures from noisy or low contrast 3-D data without the sensitivity to local minima inherent in most active surface techniques. In contrast to standard purely spatial 3-D minimal path techniques, however, we are able to represent a full tubular surface rather than just a curve which runs through its interior. Our representation also yields a natural notion of a tube's "central curve." We demonstrate and validate the utility of this approach on magnetic resonance (MR) angiography and computed tomography (CT) images of coronary arteries.
Collapse
Affiliation(s)
- Hua Li
- School of Electrical and Computer Engineering, Georgia Institute of Technology, Atlanta, GA 30332, USA.
| | | |
Collapse
|
29
|
Kharboutly Z, Treutenaere J, Claude I, Legallais C. Arterio-Venous Fistula: Two Cases Realistic Numerical Blood Flow Simulations. ACTA ACUST UNITED AC 2007; 2007:2980-3. [DOI: 10.1109/iembs.2007.4352955] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
30
|
Tizon X, Lin Q, Hansen T, Borgefors G, Johansson L, Ahlström H, Frimmel H. Identification of the main arterial branches by whole-body contrast-enhanced MRA in elderly subjects using limited user interaction and fast marching. J Magn Reson Imaging 2007; 25:806-14. [PMID: 17348000 DOI: 10.1002/jmri.20848] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022] Open
Abstract
PURPOSE To extract a graph model corresponding to a predefined set of arterial branches from whole-body contrast-enhanced magnetic resonance angiography (CE-MRA) data sets in elderly asymptomatic subjects, a high-incidence group. MATERIALS AND METHODS Maximum intensity projections (MIPs) were used as an interface to place landmarks in the three-dimensional (3D) data sets. These landmarks were linked together using fast marching to form a graph model of the arterial tree. Only vessels of interest were identified. RESULTS We tested our method on 10 subjects. We were able to build a graph model of the main arterial branches that performed well in the presence of vascular pathologies, such as stenosis and aneurysm. The results were rated by an experienced radiologist, with an overall success rate of 80%. CONCLUSION We were able to extract chosen arterial branches in 3D whole-body CE-MRA images with a moderate amount of interaction using a single MIP projection.
Collapse
Affiliation(s)
- Xavier Tizon
- Centre for Image Analysis, SLU/Uppsala University, and Department of Oncology, Radiology and Clinical Immunology, Uppsala University Hospital, Sweden
| | | | | | | | | | | | | |
Collapse
|
31
|
Manniesing R, Viergever MA, Niessen WJ. Vessel axis tracking using topology constrained surface evolution. IEEE TRANSACTIONS ON MEDICAL IMAGING 2007; 26:309-16. [PMID: 17354637 DOI: 10.1109/tmi.2006.891503] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/14/2023]
Abstract
An approach to 3-D vessel axis tracking based on surface evolution is presented. The main idea is to guide the evolution of the surface by analyzing its skeleton topology during evolution, and imposing shape constraints on the topology. For example, the intermediate topology can be processed such that it represents a single vessel segment, a bifurcation, or a more complex vascular topology. The evolving surface is then reinitialized with the newly found topology. Reinitialization is a crucial step since it creates probing behavior of the evolving front, encourages the segmentation process to extract the vascular structure of interest and reduces the risk on leaking of the curve into the background. The method was evaluated in two computed tomography angiography applications: 1) extracting the internal carotid arteries including the region in which they traverse through the skull base, which is challenging due to the proximity of bone structures and overlap in intensity values; 2) extracting the carotid bifurcations including many cases in which they are severely stenosed and contain calcifications. The vessel axis was found in 90% (18/20 internal carotids in ten patients) and 70% (14/20 carotid bifurcations in a different set of ten patients) of the cases.
Collapse
|
32
|
Brieva J, Galvez M, Toumoulin C. Coronary extraction and stenosis quantification in X-ray angiographic imaging. CONFERENCE PROCEEDINGS : ... ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL CONFERENCE 2007; 2004:1714-7. [PMID: 17272035 DOI: 10.1109/iembs.2004.1403515] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
This work describes a method to quantify stenosis in X-ray coronary angiography. Vascular edge extraction is first performed based on a deformable spline algorithm. It makes use of directional S-Gabor filters to build an external energy field that is then used in a snake optimisation scheme. A string matching technique is then applied to match the contour points and obtain a trace between the matched points. This trace allows then computing the vessel diameter and deriving quantitative stenosis measurements. Experimental results are presented on simulated data and real images.
Collapse
Affiliation(s)
- J Brieva
- Departamento de Electrónica, Campus Ciudad de México, México City, Mexico.
| | | | | |
Collapse
|
33
|
van Velsen EFS, Niessen WJ, de Weert TT, de Monyé C, van der Lugt A, Meijering E, Stokking R. Evaluation of an improved technique for lumen path definition and lumen segmentation of atherosclerotic vessels in CT angiography. Eur Radiol 2006; 17:1738-45. [PMID: 17077978 DOI: 10.1007/s00330-006-0469-x] [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] [Received: 01/06/2006] [Revised: 08/09/2006] [Accepted: 08/28/2006] [Indexed: 10/24/2022]
Abstract
Vessel image analysis is crucial when considering therapeutical options for (cardio-) vascular diseases. Our method, VAMPIRE (Vascular Analysis using Multiscale Paths Inferred from Ridges and Edges), involves two parts: a user defines a start- and endpoint upon which a lumen path is automatically defined, and which is used for initialization; the automatic segmentation of the vessel lumen on computed tomographic angiography (CTA) images. Both parts are based on the detection of vessel-like structures by analyzing intensity, edge, and ridge information. A multi-observer evaluation study was performed to compare VAMPIRE with a conventional method on the CTA data of 15 patients with carotid artery stenosis. In addition to the start- and endpoint, the two radiologists required on average 2.5 (SD: 1.9) additional points to define a lumen path when using the conventional method, and 0.1 (SD: 0.3) when using VAMPIRE. The segmentation results were quantitatively evaluated using Similarity Indices, which were slightly lower between VAMPIRE and the two radiologists (respectively 0.90 and 0.88) compared with the Similarity Index between the radiologists (0.92). The evaluation shows that the improved definition of a lumen path requires minimal user interaction, and that using this path as initialization leads to good automatic lumen segmentation results.
Collapse
Affiliation(s)
- Evert F S van Velsen
- Department of Radiology, Biomedical Imaging Group Rotterdam, Erasmus MC, University Medical Center Rotterdam, The Netherlands
| | | | | | | | | | | | | |
Collapse
|
34
|
Zhang HF, Maslov K, Li ML, Stoica G, Wang LV. In vivo volumetric imaging of subcutaneous microvasculature by photoacoustic microscopy. OPTICS EXPRESS 2006; 14:9317-23. [PMID: 19529315 DOI: 10.1364/oe.14.009317] [Citation(s) in RCA: 57] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/18/2023]
Abstract
Photoacoustic microscopy was developed to achieve volumetric imaging of the anatomy and functions of the subcutaneous microvasculature in both small animals and humans in vivo with high spatial resolution and high signal-to-background ratio. By following the skin contour in raster scanning, the ultrasonic transducer maintains focusing in the region of interest. Furthermore, off-focus lateral resolution is improved by using a synthetic-aperture focusing technique based on the virtual point detector concept. Structural images are acquired in both rats and humans, whereas functional images representing hemoglobin oxygen saturation are acquired in rats. After multiscale vesselness filtering, arterioles and venules in the image are separated based on the imaged oxygen saturation levels. Detailed structural information, such as vessel depth and spatial orientation, are revealed by volume rendering.
Collapse
|
35
|
Jiang M, Ji Q, McEwen BF. Model-based automated extraction of microtubules from electron tomography volume. ACTA ACUST UNITED AC 2006; 10:608-17. [PMID: 16871731 DOI: 10.1109/titb.2006.872042] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
We propose a model-based automated approach to extracting microtubules from noisy electron tomography volume. Our approach consists of volume enhancement, microtubule localization, and boundary segmentation to exploit the unique geometric and photometric properties of microtubules. The enhancement starts with an anisotropic invariant wavelet transform to enhance the microtubules globally, followed by a three-dimensional (3-D) tube-enhancing filter based on Weingarten matrix to further accentuate the tubular structures locally. The enhancement ends with a modified coherence-enhancing diffusion to complete the interruptions along the microtubules. The microtubules are then localized with a centerline extraction algorithm adapted for tubular objects. To perform segmentation, we novelly modify and extend active shape model method. We first use 3-D local surface enhancement to characterize the microtubule boundary and improve shape searching by relating the boundary strength with the weight matrix of the searching error. We then integrate the active shape model with Kalman filtering to utilize the longitudinal smoothness along the microtubules. The segmentation improved in this way is robust against missing boundaries and outliers that are often present in the tomography volume. Experimental results demonstrate that our automated method produces results close to those by manual process and uses only a fraction of the time of the latter.
Collapse
Affiliation(s)
- Ming Jiang
- Department of Electrical, Computer, and Systems Engineering, Rensselaer Polytechnic Institute, Troy, NY 12180, USA.
| | | | | |
Collapse
|
36
|
Jiang M, Ji Q, McEwen BF. Automated extraction of fine features of kinetochore microtubules and plus-ends from electron tomography volume. IEEE TRANSACTIONS ON IMAGE PROCESSING : A PUBLICATION OF THE IEEE SIGNAL PROCESSING SOCIETY 2006; 15:2035-48. [PMID: 16830922 DOI: 10.1109/tip.2006.877054] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
Kinetochore microtubules (KMTs) and the associated plus-ends have been areas of intense investigation in both cell biology and molecular medicine. Though electron tomography opens up new possibilities in understanding their function by imaging their high-resolution structures, the interpretation of the acquired data remains an obstacle because of the complex and cluttered cellular environment. As a result, practical segmentation of the electron tomography data has been dominated by manual operation, which is time consuming and subjective. In this paper, we propose a model-based automated approach to extracting KMTs and the associated plus-ends with a coarse-to-fine scale scheme consisting of volume preprocessing, microtubule segmentation and plus-end tracing. In volume preprocessing, we first apply an anisotropic invariant wavelet transform and a tube-enhancing filter to enhance the microtubules at coarse level for localization. This is followed with a surface-enhancing filter to accentuate the fine microtubule boundary features. The microtubule body is then segmented using a modified active shape model method. Starting from the segmented microtubule body, the plus-ends are extracted with a probabilistic tracing method improved with rectangular window based feature detection and the integration of multiple cues. Experimental results demonstrate that our automated method produces results comparable to manual segmentation but using only a fraction of the manual segmentation time.
Collapse
Affiliation(s)
- Ming Jiang
- Department of Electrical, Computer and System Engineering, Rensselaer Polytechnic Institute, Troy, NY 12180, USA
| | | | | |
Collapse
|
37
|
Manniesing R, Velthuis BK, van Leeuwen MS, van der Schaaf IC, van Laar PJ, Niessen WJ. Level set based cerebral vasculature segmentation and diameter quantification in CT angiography. Med Image Anal 2006; 10:200-14. [PMID: 16263325 DOI: 10.1016/j.media.2005.09.001] [Citation(s) in RCA: 84] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2004] [Revised: 03/09/2005] [Accepted: 09/16/2005] [Indexed: 12/28/2022]
Abstract
A level set based method is presented for cerebral vascular tree segmentation from computed tomography angiography (CTA) data. The method starts with bone masking by registering a contrast enhanced scan with a low-dose mask scan in which the bone has been segmented. Then an estimate of the background and vessel intensity distributions is made based on the intensity histogram which is used to steer the level set to capture the vessel boundaries. The relevant parameters of the level set evolution are optimized using a training set. The method is validated by a diameter quantification study which is carried out on phantom data, representing ground truth, and 10 patient data sets. The results are compared to manually obtained measurements by two expert observers. In the phantom study, the method achieves similar accuracy as the observers, but is unbiased whereas the observers are biased, i.e., the results are 0.00+/-0.23 vs. -0.32+/-0.23 mm. Also, the method's reproducibility is slightly better than the inter-and intra-observer variability. In the patient study, the method is in agreement with the observers and also, the method's reproducibility -0.04+/-0.17 mm is similar to the inter-observer variability 0.06+/-0.17 mm. Since the method achieves comparable accuracy and reproducibility as the observers, and since the method achieves better performance than the observers with respect to ground truth, we conclude that the level set based vessel segmentation is a promising method for automated and accurate CTA diameter quantification.
Collapse
Affiliation(s)
- R Manniesing
- Department of Radiology, Image Sciences Institute, University Medical Center Utrecht, Heidelberglaan 100, Room E01.335, 3584 CX Utrecht, The Netherlands.
| | | | | | | | | | | |
Collapse
|
38
|
Vascular Centerline Extraction in 3D MR Angiograms for Phase Contrast MRI Blood Flow Measurement. Int J Comput Assist Radiol Surg 2006. [DOI: 10.1007/s11548-006-0005-0] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
39
|
Hood MN, Ho VB. Contrast Agents: Innovations and Potential Applications for Body MR Angiography. Magn Reson Imaging Clin N Am 2005; 13:189-203, vii. [PMID: 15760765 DOI: 10.1016/j.mric.2004.12.011] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
Affiliation(s)
- Maureen N Hood
- Department of Radiology and Radiological Sciences, Uniformed Services University of the Health Sciences, 4301 Jones Bridge Road, Bethesda, MD 20814-4799, USA.
| | | |
Collapse
|
40
|
Bjørnerud A, Johansson L. The utility of superparamagnetic contrast agents in MRI: theoretical consideration and applications in the cardiovascular system. NMR IN BIOMEDICINE 2004; 17:465-477. [PMID: 15526351 DOI: 10.1002/nbm.904] [Citation(s) in RCA: 89] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
This review will discuss the in vivo physical chemical relaxation properties of superparamagnetic iron oxide particles. Various parameters such as size, magnetization, compartmentalization and water exchange effects and how these alter the behavior of the iron oxide particles in an in vitro vs an in vivo situation with special reference to the cardiovascular system will be exemplified. Furthermore, applications using iron oxide particles for vascular, perfusion and viability imaging as well as assessment of the inflammatory status of a given tissue will be discussed.
Collapse
Affiliation(s)
- Atle Bjørnerud
- Department of Radiology, Rikshospitalet University Hospital, N-0027 Oslo, Norway.
| | | |
Collapse
|
41
|
Chen J, Amini AA. Quantifying 3-D vascular structures in MRA images using hybrid PDE and geometric deformable models. IEEE TRANSACTIONS ON MEDICAL IMAGING 2004; 23:1251-1262. [PMID: 15493693 DOI: 10.1109/tmi.2004.834612] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
The aim of this paper is to present a hybrid approach to accurate quantification of vascular structures from magnetic resonance angiography (MRA) images using level set methods and deformable geometric models constructed with 3-D Delaunay triangulation. Multiple scale filtering based on the analysis of local intensity structure using the Hessian matrix is used to effectively enhance vessel structures with various diameters. The level set method is then applied to automatically segment vessels enhanced by the filtering with a speed function derived from enhanced MRA images. Since the goal of this paper is to obtain highly accurate vessel borders, suitable for use in fluid flow simulations, in a subsequent step, the vessel surface determined by the level set method is triangulated using 3-D Delaunay triangulation and the resulting surface is used as a parametric deformable model. Energy minimization is then performed within a variational setting with a first-order internal energy; the external energy is derived from 3-D image gradients. Using the proposed method, vessels are accurately segmented from MRA data.
Collapse
Affiliation(s)
- Jian Chen
- Cardiovascular Image Analysis Laboratory, Washington University School of Medicine, Box 8086, 660 S. Euclid Ave., St. Louis, MO 63110, USA
| | | |
Collapse
|