Registration of prone and supine CT colonography scans using correlation optimized warping and canonical correlation analysis

Corresponding Supine and Prone Colon Visualization Using Eigenfunction Analysis and Fold Modeling

We present a method for registration and visualization of corresponding supine and prone virtual colonoscopy scans based on eigenfunction analysis and fold modeling. In virtual colonoscopy, CT scans are acquired with the patient in two positions, and their registration is desirable so that physicians can corroborate findings between scans. Our algorithm performs this registration efficiently through the use of Fiedler vector representation (the second eigenfunction of the Laplace-Beltrami operator). This representation is employed to first perform global registration of the two colon positions. The registration is then locally refined using the haustral folds, which are automatically segmented using the 3D level sets of the Fiedler vector. The use of Fiedler vectors and the segmented folds presents a precise way of visualizing corresponding regions across datasets and visual modalities. We present multiple methods of visualizing the results, including 2D flattened rendering and the corresponding 3D endoluminal views. The precise fold modeling is used to automatically find a suitable cut for the 2D flattening, which provides a less distorted visualization. Our approach is robust, and we demonstrate its efficiency and efficacy by showing matched views on both the 2D flattened colons and in the 3D endoluminal view. We analytically evaluate the results by measuring the distance between features on the registered colons, and we also assess our fold segmentation against 20 manually labeled datasets. We have compared our results analytically to previous methods, and have found our method to achieve superior results. We also prove the hot spots conjecture for modeling cylindrical topology using Fiedler vector representation, which allows our approach to be used for general cylindrical geometry modeling and feature extraction.

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.

Endoluminal surface registration for CT colonography using haustral fold matching

•

Novel haustral fold matching algorithm.

•

Achieves 96.1% mean accuracy over 1743 reference points in 17 CTC datasets.

•

New initialisation to non-rigid intensity-based surface registration method.

•

Full method shows 6.0 mm mean error.

•

Use of initialisation shows significant improvement (p < 0.001).

Computed Tomographic (CT) colonography is a technique used for the detection of bowel cancer or potentially precancerous polyps. The procedure is performed routinely with the patient both prone and supine to differentiate fixed colonic pathology from mobile faecal residue. Matching corresponding locations is difficult and time consuming for radiologists due to colonic deformations that occur during patient repositioning.

We propose a novel method to establish correspondence between the two acquisitions automatically. The problem is first simplified by detecting haustral folds using a graph cut method applied to a curvature-based metric applied to a surface mesh generated from segmentation of the colonic lumen. A virtual camera is used to create a set of images that provide a metric for matching pairs of folds between the prone and supine acquisitions. Image patches are generated at the fold positions using depth map renderings of the endoluminal surface and optimised by performing a virtual camera registration over a restricted set of degrees of freedom. The intensity difference between image pairs, along with additional neighbourhood information to enforce geometric constraints over a 2D parameterisation of the 3D space, are used as unary and pair-wise costs respectively, and included in a Markov Random Field (MRF) model to estimate the maximum a posteriori fold labelling assignment.

The method achieved fold matching accuracy of 96.0% and 96.1% in patient cases with and without local colonic collapse. Moreover, it improved upon an existing surface-based registration algorithm by providing an initialisation. The set of landmark correspondences is used to non-rigidly transform a 2D source image derived from a conformal mapping process on the 3D endoluminal surface mesh. This achieves full surface correspondence between prone and supine views and can be further refined with an intensity based registration showing a statistically significant improvement (p < 0.001), and decreasing mean error from 11.9 mm to 6.0 mm measured at 1743 reference points from 17 CTC datasets.

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.

Predicting polyp location on optical colonoscopy from CT colonography by minimal-energy curve modeling of the colonoscope path

The ability to accurately locate a polyp found on computed tomographic colonography (CTC) at subsequent optical colonoscopy (OC) is an important task in colorectal cancer screening. We present a method to more accurately match polyp locations at CTC and OC. A colonoscope was modeled as a flexible tube with negligible stretch and minimal strain. The path of the colonoscope was estimated using a minimal-energy curve method. The energy function was defined and optimized by a subdivision scheme. The prediction of polyp locations at OC from CTC was converted to an optimization problem. The prediction performance was evaluated on 134 polyps by comparing the predicted with the true polyp locations at OC. The method can accurately predict polyp locations at OC to within ±0.5 colonoscope mark (5 cm) for more than 58% of polyps and to within ±1 colonoscope mark (10 cm) for more than 96% of polyps, significantly improving upon previously published methods. This method can be easily incorporated into routine OC practice and allow the colonoscopist to begin the examination by targeting locations of potential polyps found at CTC.

Matching 3-D prone and supine CT colonography scans using graphs

In this paper, we propose a new registration method for prone and supine computed tomographic colonography scans using graph matching. We formulate 3-D colon registration as a graph matching problem and propose a new graph matching algorithm based on mean field theory. In the proposed algorithm, we solve the matching problem in an iterative way. In each step, we use mean field theory to find the matched pair of nodes with highest probability. During iterative optimization, one-to-one matching constraints are added to the system in a step-by-step approach. Prominent matching pairs found in previous iterations are used to guide subsequent mean field calculations. The proposed method was found to have the best performance with smallest standard deviation compared with two other baseline algorithms called the normalized distance along the colon centerline (NDACC) ( p = 0.17) with manual colon centerline correction and spectral matching ( p < 1e-5). A major advantage of the proposed method is that it is fully automatic and does not require defining a colon centerline for registration. For the latter NDACC method, user interaction is almost always needed for identifying the colon centerlines.

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.

Registration under topological change for CT colonography

Computed tomography (CT) colonography is a minimally invasive screening technique for colorectal polyps, in which X-ray CT images of the distended colon are acquired, usually in the prone and supine positions of a single patient. Registration of segmented colon images from both positions will be useful for computer-assisted polyp detection. We have previously presented algorithms for registration of the prone and supine colons when both are well distended and there is a single connected lumen. However, due to inadequate bowel preparation or peristalsis, there may be collapsed segments in one or both of the colon images resulting in a topological change in the images. Such changes make deformable registration of the colon images difficult, and at present, there are no registration algorithms that can accommodate them. In this paper, we present an algorithm that can perform volume registration of prone/supine colon images in the presence of a topological change. For this purpose, 3-D volume images are embedded as a manifold in a 4-D space, and the manifold is evolved for nonrigid registration. Experiments using data from 24 patients show that the proposed method achieves good registration results in both the shape alignment of topologically different colon images from a single patient and the polyp location estimation between supine and prone colon images.