TY - JOUR
T1 - Computer-aided detection and quantification of cavitary tuberculosis from CT scans
AU - Xu, Ziyue
AU - Bagci, Ulas
AU - Kubler, Andre
AU - Luna, Brian
AU - Jain, Sanjay
AU - Bishai, William R.
AU - Mollura, Daniel J.
N1 - Funding Information:
This research is supported by CIDI, the intramural research program of the National Institute of Allergy and Infectious Diseases (NIAID) and the National Institute of Biomedical Imaging and Bioengineering (NIBIB). W. R. Bishai acknowledges R01 AI 079590 and the support of Howard Hughes Medical Institute (HHMI). S. Jain acknowledges the Grant No. DP2 OD006492-01. A. Kubler is funded by internal grant from Imperial College London, and acknowledges Dr. J. S. Friedland and Dr. P. T. G. Elkington. FIG. 1. Cavitation examples from CT scans. (a) Cavity without surrounding consolidation; (b) cavity surrounded by consolidation; (c) dense pathologies with airway and cavities. Airways and cavities are pointed out by dashed and solid arrows, respectively. FIG. 2. Flowchart of the cavity analysis algorithm. FIG. 3. Example of the airway and cavity detection algorithm. (a) Roughly identified candidate groups, generated by thresholding and connected component analysis, with bounding boxes. (b) and (c) Two views showing locations of different structures within the body region. (d) Detected airway (1) and cavities (2) are rendered. Note that a possible false positive (3) can be seen if lung segmentation is not conducted as a preprocessing step. FIG. 4. The result of grayscale morphological reconstruction (right column) and vesselness computation (middle column) for a CT scan (left column) from a rabbit infected with TB. FIG. 5. Airway extraction result (a) without and (b) with local refinement. FIG. 6. CT scan and a zoomed region for a rabbit without (a) and with (b) the breathing controlled mechanism. Breathing artifact was significantly reduced in (b). FIG. 7. Accuracy of the SVM airway and cavity detection algorithm with different feature sets. FIG. 8. (a)–(d) Examples for airway and cavity extraction. FIG. 9. Quantitative evaluation of cavity segmentation as compared with observer 1. FIG. 10. Quantitative evaluation of cavity segmentation as compared with observer 2. FIG. 11. DSC results for proposed method with regard to manual segmentations. FIG. 12. Longitudinal change in cavity volume and surface area for three example subjects. FIG. 13. Automated distance measurement. (a), (c), and (e) DT overlaid by cavity (bulb structure) and airway (tubular structures pointed by arrows) segmentations from three views; (b), (d), and (f) are magnified display of (a), (c), and (e) showing the nearest airway location to cavity. FIG. 14. Manual distance measurements for two images (a) and (b) with magnified display where the minimum distance is measured. FIG. 15. Examples for airway extraction. Week 3, 4, 5 for Rabbits C and D; week 3, 4, 5, 7 for Rabbit I; and week 4 for Rabbit A are shown in the figure. As illustrated, cavities are within close vicinity for Rabbits C, D, and I; while it is away from segmented airway for Rabbit A due to broken structure caused by imaging artifacts. FIG. 16. Minimum distance to airway measurements for four rabbits along time. As illustrated, the distance is relatively stable for an individual rabbit while different between rabbits. FIG. 17. Histogram of minimum distances between cavity and airway in 36 cases. FIG. 18. Airway segmentation under motion artifact on rabbit images without the breathing controlled mechanism. (a) and (d) Two rabbit images under motion artifacts and blurring with zoomed view shown in (b) and (e). (c) 3D rendering of a successful segmentation showing airway with surrounding consolidation in (a) and (b). (f) 3D rendering of a broken case where the airway structure is not fully extracted as pointed by the yellow arrow due to motion artifacts and blurring, in which case the appearances of airway (pointed by solid arrow) and nonairway (pointed by dashed arrow) are similar. FIG. 19. Airway and cavity segmentation of a human CT scan.
PY - 2013/11
Y1 - 2013/11
N2 - Purpose: To present a computer-aided detection tool for identifying, quantifying, and evaluating tuberculosis (TB) cavities in the infected lungs from computed tomography (CT) scans. Methods: The authors' proposed method is based on a novel shape-based automated detection algorithm on CT scans followed by a fuzzy connectedness (FC) delineation procedure. In order to assess interaction between cavities and airways, the authors first roughly identified air-filled structures (airway, cavities, esophagus, etc.) by thresholding over Hounsfield unit of CT image. Then, airway and cavity structure detection was conducted within the support vector machine classification algorithm. Once airway and cavities were detected automatically, the authors extracted airway tree using a hybrid multiscale approach based on novel affinity relations within the FC framework and segmented cavities using intensity-based FC algorithm. At final step, the authors refined airway structures within the local regions of FC with finer control. Cavity segmentation results were compared to the reference truths provided by expert radiologists and cavity formation was tracked longitudinally from serial CT scans through shape and volume information automatically determined through the authors' proposed system. Morphological evolution of the cavitary TB were analyzed accordingly with this process. Finally, the authors computed the minimum distance between cavity surface and nearby airway structures by using the linear time distance transform algorithm to explore potential role of airways in cavity formation and morphological evolution. Results: The proposed methodology was qualitatively and quantitatively evaluated on pulmonary CT images of rabbits experimentally infected with TB, and multiple markers such as cavity volume, cavity surface area, minimum distance from cavity surface to the nearest bronchial-tree, and longitudinal change of these markers (namely, morphological evolution of cavities) were determined precisely. While accuracy of the authors' cavity detection algorithm was 94.61%, airway detection part of the proposed methodology showed even higher performance by 99.8%. Dice similarity coefficients for cavitary segmentation experiments were found to be approximately 99.0% with respect to the reference truths provided by two expert radiologists (blinded to their evaluations). Moreover, the authors noted that volume derived from the authors' segmentation method was highly correlated with those provided by the expert radiologists (R2 = 0.99757 and R 2 = 0.99496, p < 0.001, with respect to the observer 1 and observer 2) with an interobserver agreement of 98%. The authors quantitatively confirmed that cavity formation was positioned by the nearby bronchial-tree after exploring the respective spatial positions based on the minimum distance measurement. In terms of efficiency, the core algorithms take less than 2 min on a linux machine with 3.47 GHz CPU and 24 GB memory. Conclusion: The authors presented a fully automatic method for cavitary TB detection, quantification, and evaluation. The performance of every step of the algorithm was qualitatively and quantitatively assessed. With the proposed method, airways and cavities were automatically detected and subsequently delineated in high accuracy with heightened efficiency. Furthermore, not only morphological information of cavities were obtained through the authors' proposed framework, but their spatial relation to airways, and longitudinal analysis was also provided to get further insight on cavity formation in tuberculosis disease. To the authors' best of knowledge, this is the first study in computerized analysis of cavitary tuberculosis from CT scans.
AB - Purpose: To present a computer-aided detection tool for identifying, quantifying, and evaluating tuberculosis (TB) cavities in the infected lungs from computed tomography (CT) scans. Methods: The authors' proposed method is based on a novel shape-based automated detection algorithm on CT scans followed by a fuzzy connectedness (FC) delineation procedure. In order to assess interaction between cavities and airways, the authors first roughly identified air-filled structures (airway, cavities, esophagus, etc.) by thresholding over Hounsfield unit of CT image. Then, airway and cavity structure detection was conducted within the support vector machine classification algorithm. Once airway and cavities were detected automatically, the authors extracted airway tree using a hybrid multiscale approach based on novel affinity relations within the FC framework and segmented cavities using intensity-based FC algorithm. At final step, the authors refined airway structures within the local regions of FC with finer control. Cavity segmentation results were compared to the reference truths provided by expert radiologists and cavity formation was tracked longitudinally from serial CT scans through shape and volume information automatically determined through the authors' proposed system. Morphological evolution of the cavitary TB were analyzed accordingly with this process. Finally, the authors computed the minimum distance between cavity surface and nearby airway structures by using the linear time distance transform algorithm to explore potential role of airways in cavity formation and morphological evolution. Results: The proposed methodology was qualitatively and quantitatively evaluated on pulmonary CT images of rabbits experimentally infected with TB, and multiple markers such as cavity volume, cavity surface area, minimum distance from cavity surface to the nearest bronchial-tree, and longitudinal change of these markers (namely, morphological evolution of cavities) were determined precisely. While accuracy of the authors' cavity detection algorithm was 94.61%, airway detection part of the proposed methodology showed even higher performance by 99.8%. Dice similarity coefficients for cavitary segmentation experiments were found to be approximately 99.0% with respect to the reference truths provided by two expert radiologists (blinded to their evaluations). Moreover, the authors noted that volume derived from the authors' segmentation method was highly correlated with those provided by the expert radiologists (R2 = 0.99757 and R 2 = 0.99496, p < 0.001, with respect to the observer 1 and observer 2) with an interobserver agreement of 98%. The authors quantitatively confirmed that cavity formation was positioned by the nearby bronchial-tree after exploring the respective spatial positions based on the minimum distance measurement. In terms of efficiency, the core algorithms take less than 2 min on a linux machine with 3.47 GHz CPU and 24 GB memory. Conclusion: The authors presented a fully automatic method for cavitary TB detection, quantification, and evaluation. The performance of every step of the algorithm was qualitatively and quantitatively assessed. With the proposed method, airways and cavities were automatically detected and subsequently delineated in high accuracy with heightened efficiency. Furthermore, not only morphological information of cavities were obtained through the authors' proposed framework, but their spatial relation to airways, and longitudinal analysis was also provided to get further insight on cavity formation in tuberculosis disease. To the authors' best of knowledge, this is the first study in computerized analysis of cavitary tuberculosis from CT scans.
KW - airway tree
KW - cavitary tuberculosis
KW - computer aided detection
KW - fuzzy connectedness
KW - segmentation
UR - http://www.scopus.com/inward/record.url?scp=84889677853&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84889677853&partnerID=8YFLogxK
U2 - 10.1118/1.4824979
DO - 10.1118/1.4824979
M3 - Article
C2 - 24320475
AN - SCOPUS:84889677853
SN - 0094-2405
VL - 40
JO - Medical physics
JF - Medical physics
IS - 11
M1 - 113701
ER -