Reduction of false positives by internal features for polyp detection in CT-based virtual colonoscopy
An illustration of two different kernel functions—Uniform kernel (left) and Gaussian kernel (right).
(a)An illustration of a colonic polyp that is fully inside the lumen and has a typical elliptical shape. (b) In a CT image, the geometrical feature of a polyp (indicated by the white rectangle) becomes blurred. (c) The zoomed display of the polyp that shows some kinds of “elliptical” shape for both the detected portion inside the lumen and the inner portion inside the wall. (d) In a CT image, a polyp only shows part of its whole border inside the lumen. The solid curve indicates the outer border in the lumen; the dashed curve indicates the inner border behind the surface that cannot be detected by the shape analysis on the surface directly.
(a) 2D illustration of three rays emitted from a given voxel (indicated by the cross) in a candidate. (b) The CT density profile along , i.e., the middle arrow toward the bottom right in (a). (c) The step-like profile after the Harr wavelet transformation and filtering. (d) The coding sequence of the step-like profile. The above is the initial coding sequence. After merging some small parts, we get the final coding as shown in the bottom. (e) Some special variance pattern can be detected using the coding sequence and the exact position is identified by the first- and second-order derivatives of the original profile. The detected point is indicated by the arrow. Meanwhile, we define a search range (indicated by the dotted vertical line) to reduce the possibility of finding a wrong edge. (f) The detected three new border points (the crosses) along the three rays, respectively. Among them, two are the outer border points and one is the inner border point.
An illustration of Harr wavelet transformation in the edge detector.
(a-1) A polyp (6 mm in diameter) in a CT image. The solid curve indicates all the outer border points that are grouped into the initial candidate. (a-2) The ellipsoid (dashed line) constructed using only the outer border points in the left image. This ellipsoid may only cover a portion of the polyp. (a-3) The ellipsoid generated using both the outer and the inner border points. All the inner border points are detected using the ray-driven edge finder. (b-1) Another polyp candidate (10 mm in diameter) with an irregular surface. Its outer border points are grouped into two different candidates. (b-2) Two different ellipsoids are generated. Either of them covers the whole polyp volume. (b-3) A new ellipsoid generated is generated by merging the two candidates. (c) Another example for a candidate of 7 mm in diameter and its corresponding ellipsoid. (d) Another candidate of 15 mm in diameter that covers only a small section of the whole border of the polyp, and its ellipsoid generated from the small solid curve using the ray-driven edge finder. (e) An example of two polyp-like colon folds that cannot be eliminated by the shape analysis and become polyp candidates, resulting in two ellipsoids. However, due to the lack of inner border information, the upper ellipsoid is spindly. Although the lower ellipsoid shows a similar shape as a true polyp, it has much less inner border information than a true polyp has and, therefore, is likely to be eliminated by the texture information derived from the ellipsoid, as discussed in the next section.
(a) An illustration and definition of the enlarged and shrunk borders for an eROI. (b) The three borders in a CT image. The solid line represents the original border of the eROI; and the two dashed lines represent the shrunk border and the enlarged border of the eROI, respectively.
The 2D/3D plots of the triple-element vectors of voxels from polyps and nonpolyps (or FPs). All voxels are selected by the computer randomly. Pictures (a)–(c) show the 2D plots of the triple-element vectors. Picture (d) shows the 3D plot. The triple-element vectors of polyp voxels show a convergence toward the top right corner (by the circles).
An illustration of the mapping procedure of octsphere parametrization. Picture (a) shows the parametrization of the ellipsoid surface or border. The whole surface is divided into many triangular patches. In picture (b), a CT density profile is drawn, for each patch, along the gradient ray emitted from its central point. The search range is defined by the intersection points of the ray with the enlarged and shrunk borders. If there exists a border in the given search range, this patch is marked [or colored by the red in picture (c)]. Picture (d) shows the mapped image. The gray patches indicate all marked patches, and the black (darker) patches are the unmarked patches. Picture (e) shows a patch pair.
Illustrations of a transformation function (left) and a two-level classifier for the final polyp identification (right).
The FROC curve of the presented CAD scheme for the 153 patient datasets. (a) The FROC curves of the presented CAD scheme for different polyp sizes. (b) The FROC curves for the CAD scheme with different groups of the internal features from all datasets.
The ROC curves resulted from polyp detection using different curvatures. Notation—“Uniform [x] mm” indicates the result of the global curvature with 2 uniform kernel function and a half-curve length of x mm; and “Gaussian-[x]” is the result of the global curvature with Gaussian kernel function and an alpha value of “x” mm. Notation “Uniform” shows the result of a combination of the local and global curvatures. (a) The ROC curves of local and global curvature. (b) The ROC curves from polyp results using the combination of local and global curvature.
Experimental results from the 153 patient datasets after the surface shape-based filtering operation.
Experimental results of the presented CAD scheme after using the internal features.
Experimental results of the presented CAD scheme using different feature groups.
Article metrics loading...
Full text loading...