Teniae coli-based circumferential localization system for CT colonography: feasibility study

Progress in Fully Automated Abdominal CT Interpretation

OBJECTIVE

Automated analysis of abdominal CT has advanced markedly over just the last few years. Fully automated assessment of organs, lymph nodes, adipose tissue, muscle, bowel, spine, and tumors are some examples where tremendous progress has been made. Computer-aided detection of lesions has also improved dramatically.

CONCLUSION

This article reviews the progress and provides insights into what is in store in the near future for automated analysis for abdominal CT, ultimately leading to fully automated interpretation.

Feasibility of using the marginal blood vessels as reference landmarks for CT colonography

OBJECTIVE

The purpose of this study was to show the spatial relationship of the colonic marginal blood vessels and the teniae coli on CT colonography (CTC) and the use of the marginal blood vessels for supine-prone registration of polyps and for determination of proper connectivity of collapsed colonic segments.

MATERIALS AND METHODS

We manually labeled the marginal blood vessels on 15 CTC examinations. Colon segmentation, centerline extraction, teniae detection, and teniae identification were automatically performed. For assessment of their spatial relationships, the distances from the marginal blood vessels to the three teniae coli and to the colon were measured. Student t tests (paired, two-tailed) were performed to evaluate the differences among these distances. To evaluate the reliability of the marginal vessels as reference points for polyp correlation, we analyzed 20 polyps from 20 additional patients who underwent supine and prone CTC. The average difference of the circumferential polyp position on the supine and prone scans was computed. Student t tests (paired, two-tailed) were performed to evaluate the supine-prone differences of the distance. We performed a study on 10 CTC studies from 10 patients with collapsed colonic segments by manually tracing the marginal blood vessels near the collapsed regions to resolve the ambiguity of the colon path.

RESULTS

The average distances (± SD) from the marginal blood vessels to the tenia mesocolica, tenia omentalis, and tenia libera were 20.1 ± 3.1 mm (95% CI, 18.5-21.6 mm), 39.5 ± 4.8 mm (37.1-42.0 mm), and 36.9 ± 4.2 mm (34.8-39.1 mm), respectively. Pairwise comparison showed that these distances to the tenia libera and tenia omentalis were significantly different from the distance to the tenia mesocolica (p < 0.001). The average distance from the marginal blood vessels to the colon wall was 15.3 ± 2.0 mm (14.2-16.3 mm). For polyp localization, the average difference of the circumferential polyp position on the supine and prone scans was 9.6 ± 9.4 mm (5.5-13.7 mm) (p = 0.15) and expressed as a percentage of the colon circumference was 3.1% ± 2.0% (2.3-4.0%) (p = 0.83). We were able to trace the marginal blood vessels for 10 collapsed colonic segments and determine the paths of the colon in these regions.

CONCLUSION

The marginal blood vessels run parallel to the colon in proximity to the tenia mesocolica and enable accurate supine-prone registration of polyps and localization of the colon path in areas of collapse. Thus, the marginal blood vessels may be used as reference landmarks complementary to the colon centerline and teniae coli.

A Novel Colon Wall Flattening Model for Computed Tomographic Colonography: Method and Validation

Computed tomographic colonography (CTC) has been developed for screening of colon cancer. Flattening the three-dimensional (3D) colon wall into two-dimensional (2D) image is believed to (1) provide supplementary information to the endoscopic views and further (2) facilitate colon registration, taniae coli (TC) detection, and haustral fold segmentation. Though the previously-used conformal mapping-based flattening methods can preserve the angular geometry, they have the limitations in providing accurate information of the 3D inner colon wall due to the lack of undulating topography. In this paper, we present a novel colon-wall flattening method using a strategy of 2.5D approach. Coupling with the conformal flattening model, the presented new approach builds an elevation distance map to depict the neighborhood characteristics of the inner colon wall. We validated the new method via two CTC applications: TC detection and haustral fold segmentation. Experimental results demonstrated the effectiveness of our strategy for CTC studies.

CT colonography: external clinical validation of an algorithm for computer-assisted prone and supine registration

PURPOSE

To perform external validation of a computer-assisted registration algorithm for prone and supine computed tomographic (CT) colonography and to compare the results with those of an existing centerline method.

MATERIALS AND METHODS

All contributing centers had institutional review board approval; participants provided informed consent. A validation sample of CT colonographic examinations of 51 patients with 68 polyps (6-55 mm) was selected from a publicly available, HIPAA compliant, anonymized archive. No patients were excluded because of poor preparation or inadequate distension. Corresponding prone and supine polyp coordinates were recorded, and endoluminal surfaces were registered automatically by using a computer algorithm. Two observers independently scored three-dimensional endoluminal polyp registration success. Results were compared with those obtained by using the normalized distance along the colonic centerline (NDACC) method. Pairwise Wilcoxon signed rank tests were used to compare gross registration error and McNemar tests were used to compare polyp conspicuity.

RESULTS

Registration was possible in all 51 patients, and 136 paired polyp coordinates were generated (68 polyps) to test the algorithm. Overall mean three-dimensional polyp registration error (mean ± standard deviation, 19.9 mm ± 20.4) was significantly less than that for the NDACC method (mean, 27.4 mm ± 15.1; P = .001). Accuracy was unaffected by colonic segment (P = .76) or luminal collapse (P = .066). During endoluminal review by two observers (272 matching tasks, 68 polyps, prone to supine and supine to prone coordinates), 223 (82%) polyp matches were visible (120° field of view) compared with just 129 (47%) when the NDACC method was used (P < .001). By using multiplanar visualization, 48 (70%) polyps were visible after scrolling ± 15 mm in any multiplanar axis compared with 16 (24%) for NDACC (P < .001).

CONCLUSION

Computer-assisted registration is more accurate than the NDACC method for mapping the endoluminal surface and matching the location of polyps in corresponding prone and supine CT colonographic acquisitions.

Systematic approach: an evidence management strategy for better decision-making

Evidence-based medicine aims to apply best evidences to medical decision making. Although evidence are established using the scientific method, decision making has received less priority. Decision making is a mental process requiring a systematic analysis of evidences for a specific use, leading to a final choice of action. In many scientific disciplines involving decision and control, this process has been significantly improved by making it as a system (a set of interacting entities forming an integrated whole). Hypothesizing that most medical decisions can be described as a system, we present a schematic systematic loop, based on four traditional medical steps (nosology, semeiology, pathophysiology, and therapy), and on the four transitions between these steps. Steps are evaluated by reproducibility, transitions are evaluated by predictability. This leads to formulate eight basic questions for testing the reliability of any loop. We applied this approach to a specific study (EPHESUS) to show its interest in the control of the medical knowledge provided by a study and in indicating where evidence is missing. A systematic approach helps evidence-based medicine in structuring evidences for better medical decisions and in determining which research needs priority.

Haustral fold segmentation with curvature-guided level set evolution

Human colon has complex structures mostly because of the haustral folds. The folds are thin flat protrusions on the colon wall, which complicate the shape analysis for computer-aided detection (CAD) of colonic polyps. Fold segmentation may help reduce the structural complexity, and the folds can serve as an anatomic reference for computed tomographic colonography (CTC). Therefore, in this study, based on a model of the haustral fold boundaries, we developed a level-set approach to automatically segment the fold surfaces. To evaluate the developed fold segmentation algorithm, we first established the ground truth of haustral fold boundaries by experts' drawing on 15 patient CTC datasets without severe under/over colon distention from two medical centers. The segmentation algorithm successfully detected 92.7% of the folds in the ground truth. In addition to the sensitivity measure, we further developed a merit of segmented-area ratio (SAR), i.e., the ratio between the area of the intersection and union of the expert-drawn folds and the area of the automatically segmented folds, to measure the segmentation accuracy. The segmentation algorithm reached an average value of SAR = 86.2%, showing a good match with the ground truth on the fold surfaces. We believe the automatically segmented fold surfaces have the potential to benefit many postprocedures in CTC, such as CAD, taenia coli extraction, supine-prone registration, etc.

Automated teniae coli detection and identification on computed tomographic colonography

PURPOSE

Computed tomographic colonography (CTC) is a minimally invasive technique for colonic polyps and cancer screening. Teniae coli are three bands of longitudinal smooth muscle on the colon surface. Teniae coli are important anatomically meaningful landmarks on human colon. In this paper, the authors propose an automatic teniae coli detection method for CT colonography.

METHODS

The original CTC slices are first segmented and reconstructed to a 3D colon surface. Then, the 3D colon surface is unfolded using a reversible projection technique. After that the unfolded colon is projected to a 2D height map. The teniae coli are detected using the height map and then reversely projected back to the 3D colon. Since teniae are located at the junctions where the haustral folds meet, the authors apply 2D Gabor filter banks to extract features of haustral folds. The maximum response of the filter banks is then selected as the feature image. The fold centers are then identified based on local maxima and thresholding on the feature image. Connecting the fold centers yields a path of the folds. Teniae coli are extracted as lines running between the fold paths. The authors used the spatial relationship between ileocecal valve (ICV) and teniae mesocolica (TM) to identify the TM, then the teniae omentalis (TO) and the teniae libera (TL) can be identified subsequently.

RESULTS

The authors tested the proposed method on 47 cases of 37 patients, 10 of the patients with both supine and prone CT scans. The proposed method yielded performance with an average normalized root mean square error (RMSE) ( ± standard deviation [95% confidence interval]) of 4.87% ( ± 2.93%, [4.05% 5.69%]).

CONCLUSIONS

The proposed fully-automated teniae coli detection and identification method is accurate and promising for future clinical applications.

Registration of the endoluminal surfaces of the colon derived from prone and supine CT colonography

PURPOSE

Computed tomographic (CT) colonography is a relatively new technique for detecting bowel cancer or potentially precancerous polyps. CT scanning is combined with three-dimensional (3D) image reconstruction to produce a virtual endoluminal representation similar to optical colonoscopy. Because retained fluid and stool can mimic pathology, CT data are acquired with the bowel cleansed and insufflated with gas and patient in both prone and supine positions. Radiologists then match visually endoluminal locations between the two acquisitions in order to determine whether apparent pathology is real or not. This process is hindered by the fact that the colon, essentially a long tube, can undergo considerable deformation between acquisitions. The authors present a novel approach to automatically establish spatial correspondence between prone and supine endoluminal colonic surfaces after surface parameterization, even in the case of local colon collapse.

METHODS

The complexity of the registration task was reduced from a 3D to a 2D problem by mapping the surfaces extracted from prone and supine CT colonography onto a cylindrical parameterization. A nonrigid cylindrical registration was then performed to align the full colonic surfaces. The curvature information from the original 3D surfaces was used to determine correspondence. The method can also be applied to cases with regions of local colonic collapse by ignoring the collapsed regions during the registration.

RESULTS

Using a development set, suitable parameters were found to constrain the cylindrical registration method. Then, the same registration parameters were applied to a different set of 13 validation cases, consisting of 8 fully distended cases and 5 cases exhibiting multiple colonic collapses. All polyps present were well aligned, with a mean (+/- std. dev.) registration error of 5.7 (+/- 3.4) mm. An additional set of 1175 reference points on haustral folds spread over the full endoluminal colon surfaces resulted in an error of 7.7 (+/- 7.4) mm. Here, 82% of folds were aligned correctly after registration with a further 15% misregistered by just onefold.

CONCLUSIONS

The proposed method reduces the 3D registration task to a cylindrical registration representing the endoluminal surface of the colon. Our algorithm uses surface curvature information as a similarity measure to drive registration to compensate for the large colorectal deformations that occur between prone and supine data acquisitions. The method has the potential to both enhance polyp detection and decrease the radiologist's interpretation time.

Supine and prone colon registration using quasi-conformal mapping

In virtual colonoscopy, CT scans are typically acquired with the patient in both supine (facing up) and prone (facing down) positions. The registration of these two scans is desirable so that the user can clarify situations or confirm polyp findings at a location in one scan with the same location in the other, thereby improving polyp detection rates and reducing false positives. However, this supine-prone registration is challenging because of the substantial distortions in the colon shape due to the patient's change in position. We present an efficient algorithm and framework for performing this registration through the use of conformal geometry to guarantee that the registration is a diffeomorphism (a one-to-one and onto mapping). The taeniae coli and colon flexures are automatically extracted for each supine and prone surface, employing the colon geometry. The two colon surfaces are then divided into several segments using the flexures, and each segment is cut along a taenia coli and conformally flattened to the rectangular domain using holomorphic differentials. The mean curvature is color encoded as texture images, from which feature points are automatically detected using graph cut segmentation, mathematic morphological operations, and principal component analysis. Corresponding feature points are found between supine and prone and are used to adjust the conformal flattening to be quasi-conformal, such that the features become aligned. We present multiple methods of visualizing our results, including 2D flattened rendering, corresponding 3D endoluminal views, and rendering of distortion measurements. We demonstrate the efficiency and efficacy of our registration method by illustrating matched views on both the 2D flattened colon images and in the 3D volume rendered colon endoluminal view. We analytically evaluate the correctness of the results by measuring the distance between features on the registered colons.

Improving the accuracy of CTC interpretation: computer-aided detection

Computer-aided polyp detection aims to improve the accuracy of the colonography interpretation. The computer searches the colonic wall to look for polyplike protrusions and presents a list of suspicious areas to a physician for further analysis. Computer-aided polyp detection has developed rapidly in the past decade in the laboratory setting and has sensitivities comparable with those of experts. Computer-aided polyp detection tends to help inexperienced readers more than experienced ones and may also lead to small reductions in specificity. In its currently proposed use as an adjunct to standard image interpretation, computer-aided polyp detection serves as a spellchecker rather than an efficiency enhancer.

Normalized distance along the colon centerline: a method for correlating polyp location on CT colonography and optical colonoscopy

OBJECTIVE

The ability to accurately locate a polyp found on CT colonography (CTC) at subsequent optical colonoscopy (OC) is an important part of the successful implementation of CTC for colorectal cancer screening. The purpose of this study was to determine whether a polyp's normalized distance along the colon centerline derived from CTC data can accurately predict its location on OC.

MATERIALS AND METHODS

The polyp population consisted of 152 polyps in 121 patients. CTC polyp findings were verified by same-day segmentally-unblinded OC. Each polyp's normalized distance along the colon centerline was computed by dividing its distance from the anorectal junction measured along the colon centerline by the length of the colon at CTC. The predicted polyp location at OC was computed by multiplying the normalized distance along the colon centerline by the colon length at OC (i.e., the distance to the cecum as determined at full colonoscope insertion). The differences between the true and predicted polyp locations at OC were compared using paired Student's t tests, linear regression, prediction interval assessment, and Bland-Altman analyses.

RESULTS

The differences between the true and predicted polyp locations at OC using the supine and prone CTC-normalized distances along the colon centerline were 2.2 +/- 10.5 cm (mean +/- SD; n = 136) and 1.5 +/- 10.5 cm (n = 135), respectively. The predicted location was within 10 cm of its true location for 71.3% (97/136) to 74.8% (101/135) of polyps and within 20 cm of its true location for 93.3% (126/135) to 93.4% (127/136) of polyps.

CONCLUSION

By computing the normalized distance along the colon centerline of a polyp found at CTC, the location of a polyp at OC can be predicted to within 10 cm (i.e., 1 colonoscope mark) for the majority of polyps.

On concise 3-D simple point characterizations: a marching cubes paradigm

The centerlines of tubular structures are useful for medical image visualization and computer-aided diagnosis applications. They can be effectively extracted by using a thinning algorithm that erodes an object layer by layer until only a skeleton is left. An object point is "simple" and can be safely deleted only if the resultant image is topologically equivalent to the original. Numerous characterizations of 3-D simple points based on digital topology already exist. However, little work has been done in the context of marching cubes (MC). This paper reviews several concise 3-D simple point characterizations in a MC paradigm. By using the Euler characteristic and a few newly observed properties in the context of connectivity-consistent MC, we present concise and more self-explanatory proofs. We also present an efficient method for computing the Euler characteristic locally for MC surfaces. Performance evaluations on different implementations are conducted on synthetic data and multidetector computed tomography examination of virtual colonoscopy and angiography.

Computer-aided detection for CT colonography: update 2007

Computed tomographic colonography (CTC) is an emerging technique for polyp detection in the colon. However, lesion detection can be challenging due to insufficient patient preparation, chosen CT technique or reader imperfection. The primary goal of computer-aided detection (CAD) for CTC is locating possible polyps, and presenting the reader with these polyp candidates. Other goals are sensitivity improvement and reduction of reading time and inter-observer variability. The multistep CAD procedure typically consists of segmentation of the colonic wall (e.g. region growing); selection of intermediate polyp candidates (curvature analysis, sphere fitting, normal analysis, slope density function ...); classification of final candidates for detection and listing suspicious polyps (location, size and volume). Remaining task for the radiologist is the validation or rejection of the polyp candidates. State-of-the-art CAD systems should require minimal or even no user interaction for the extraction of the colonic wall, offer a computation time less than 10-20 min and high sensitivity and specificity for different polyp sizes and shapes, with a low number of false positives. These systems have the potential to increase radiologist's performance and to decrease inter-reader variability. Besides CAD key techniques we also discuss new developments in CAD and describe recent applications facilitating CTC.