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

Multi-domain features for reducing false positives in automated detection of clustered microcalcifications in digital breast tomosynthesis

PURPOSE

In digital breast tomosynthesis (DBT) imaging, a microcalcification (MC) cluster may span across different slices and blurring exists in the out-of-focus slices. We developed a radiomics approach to extract features from focus slice and combine multiple spatial domains to reduce false positives (FPs) in an automated pipeline of detecting MC clusters.

METHODS

We performed a retrospective study on a cohort of 290 Chinese women patients with a total of 580 DBT volumes. We developed an automated MC detection pipeline that consists of two stages: an initial detection to identify a set of MC candidates that may include many FPs, followed by a radiomics-based classification model to identify and reduce the FPs. We extract both two-dimensional (2D) and three-dimensional (3D) radiomics features from multiple spatial domains, including a focus slice, projection image, and tomographic volume. A linear discriminant classifier was used coupled with a sequential forward feature selection procedure. The free-response operating characteristics (FROC) curve and partial area under the FROC curve (pAUC) in the FP rate range of 0 to 2 per DBT volume were used to evaluate the model's performance.

RESULTS

At a sensitivity of 90%, the FP rate was reduced from 1.3 to 0.2 per DBT volume after applying the multi-domain-based classification on the initial detections. The multi-domain yielded a significantly higher pAUC compared to the initial detection (increase of pAUC = 0.2278, P < 0.0001), focus slice (increase of pAUC = 0.0345, P = 0.0152), project image (increase of pAUC = 0.1043, P < 0.0001), and tomographic volume (increase of pAUC = 0.0791, P = 0.0032).

CONCLUSION

The radiomic features extracted from the three domains may provide complementary information and their integration can significantly reduce FPs in automated detection of MCs in DBT volumes on a large Chinese women population.

Simplified computer-aided detection scheme of microcalcification clusters in digital breast tomosynthesis images

A computer-aided detection (CADe) algorithm for clustered microcalcifications (MCs) in reconstructed digital breast tomosynthesis (DBT) images is suggested. The MC-like objects were enhanced by a Hessian-based 3D calcification response function, and a signal-to-noise ratio (SNR) enhanced image was also generated to screen the MC clustering seed objects. A connected component segmentation method was used to detect the cluster seed objects, which were considered as potential clustering centers of MCs. Bounding cubes for the accepted clustering seed candidate were generated and the overlapping cubes were combined and examined. After the MC clustering and false-positive (FP) reduction step, the average number of FPs was estimated to be 0.87 per DBT volume with a sensitivity of 90.5%.

Influence of phase connectivity on the relationship among capillary pressure, fluid saturation, and interfacial area in two-fluid-phase porous medium systems

Multiphase flows in porous medium systems are typically modeled at the macroscale by applying the principles of continuum mechanics to develop models that describe the behavior of averaged quantities, such as fluid pressure and saturation. These models require closure relations to produce solvable forms. One of these required closure relations is an expression relating the capillary pressure to fluid saturation and, in some cases, other topological invariants such as interfacial area and the Euler characteristic (or average Gaussian curvature). The forms that are used in traditional models, which typically consider only the relationship between capillary pressure and saturation, are hysteretic. An unresolved question is whether the inclusion of additional morphological and topological measures can lead to a nonhysteretic closure relation. Relying on the lattice Boltzmann (LB) method, we develop an approach to investigate equilibrium states for a two-fluid-phase porous medium system, which includes disconnected nonwetting phase features. A set of simulations are performed within a random close pack of 1964 spheres to produce a total of 42 908 distinct equilibrium configurations. This information is evaluated using generalized additive models to quantitatively assess the degree to which functional relationships can explain the behavior of the equilibrium data. The variance of various model estimates is computed, and we conclude that, except for the limiting behavior close to a single fluid regime, capillary pressure can be expressed as a deterministic and nonhysteretic function of fluid saturation, interfacial area between the fluid phases, and the Euler characteristic. To our knowledge, this work is unique in the methods employed, the size of the data set, the resolution in space and time, the true equilibrium nature of the data, the parametrizations investigated, and the broad set of functions examined. The conclusion of essentially nonhysteretic behavior provides support for an evolving class of two-fluid-phase flow in porous medium systems models.

Three-Dimensional Computer-Aided Detection of Microcalcification Clusters in Digital Breast Tomosynthesis

We propose computer-aided detection (CADe) algorithm for microcalcification (MC) clusters in reconstructed digital breast tomosynthesis (DBT) images. The algorithm consists of prescreening, MC detection, clustering, and false-positive (FP) reduction steps. The DBT images containing the MC-like objects were enhanced by a multiscale Hessian-based three-dimensional (3D) objectness response function and a connected-component segmentation method was applied to extract the cluster seed objects as potential clustering centers of MCs. Secondly, a signal-to-noise ratio (SNR) enhanced image was also generated to detect the individual MC candidates and prescreen the MC-like objects. Each cluster seed candidate was prescreened by counting neighboring individual MC candidates nearby the cluster seed object according to several microcalcification clustering criteria. As a second step, we introduced bounding boxes for the accepted seed candidate, clustered all the overlapping cubes, and examined. After the FP reduction step, the average number of FPs per case was estimated to be 2.47 per DBT volume with a sensitivity of 83.3%.

Topology preserving thinning of cell complexes

A topology preserving skeleton is a synthetic representation of an object that retains its topology and many of its significant morphological properties. The process of obtaining the skeleton, referred to as skeletonization or thinning, is a very active research area. It plays a central role in reducing the amount of information to be processed during image analysis and visualization, computer-aided diagnosis, or by pattern recognition algorithms. This paper introduces a novel topology preserving thinning algorithm, which removes simple cells-a generalization of simple points-of a given cell complex. The test for simple cells is based on acyclicity tables automatically produced in advance with homology computations. Using acyclicity tables render the implementation of thinning algorithms straightforward. Moreover, the fact that tables are automatically filled for all possible configurations allows to rigorously prove the generality of the algorithm and to obtain fool-proof implementations. The novel approach enables, for the first time, according to our knowledge, to thin a general unstructured simplicial complex. Acyclicity tables for cubical and simplicial complexes and an open source implementation of the thinning algorithm are provided as an additional material to allow their immediate use in the vast number of applications arising in medical imaging and beyond.

Computer-aided detection of clustered microcalcifications in multiscale bilateral filtering regularized reconstructed digital breast tomosynthesis volume

PURPOSE

Develop a computer-aided detection (CADe) system for clustered microcalcifications in digital breast tomosynthesis (DBT) volume enhanced with multiscale bilateral filtering (MSBF) regularization.

METHODS

With Institutional Review Board approval and written informed consent, two-view DBT of 154 breasts, of which 116 had biopsy-proven microcalcification (MC) clusters and 38 were free of MCs, was imaged with a General Electric GEN2 prototype DBT system. The DBT volumes were reconstructed with MSBF-regularized simultaneous algebraic reconstruction technique (SART) that was designed to enhance MCs and reduce background noise while preserving the quality of other tissue structures. The contrast-to-noise ratio (CNR) of MCs was further improved with enhancement-modulated calcification response (EMCR) preprocessing, which combined multiscale Hessian response to enhance MCs by shape and bandpass filtering to remove the low-frequency structured background. MC candidates were then located in the EMCR volume using iterative thresholding and segmented by adaptive region growing. Two sets of potential MC objects, cluster centroid objects and MC seed objects, were generated and the CNR of each object was calculated. The number of candidates in each set was controlled based on the breast volume. Dynamic clustering around the centroid objects grouped the MC candidates to form clusters. Adaptive criteria were designed to reduce false positive (FP) clusters based on the size, CNR values and the number of MCs in the cluster, cluster shape, and cluster based maximum intensity projection. Free-response receiver operating characteristic (FROC) and jackknife alternative FROC (JAFROC) analyses were used to assess the performance and compare with that of a previous study.

RESULTS

Unpaired two-tailed t-test showed a significant increase (p < 0.0001) in the ratio of CNRs for MCs with and without MSBF regularization compared to similar ratios for FPs. For view-based detection, a sensitivity of 85% was achieved at an FP rate of 2.16 per DBT volume. For case-based detection, a sensitivity of 85% was achieved at an FP rate of 0.85 per DBT volume. JAFROC analysis showed a significant improvement in the performance of the current CADe system compared to that of our previous system (p = 0.003).

CONCLUSIONS

MBSF regularized SART reconstruction enhances MCs. The enhancement in the signals, in combination with properly designed adaptive threshold criteria, effective MC feature analysis, and false positive reduction techniques, leads to a significant improvement in the detection of clustered MCs in DBT.

Anatomical annotation on vascular structure in volume rendered images

The precise annotation of vascular structure is desired in computer-assisted systems to help surgeons identify each vessel branch. This paper proposes a method that annotates vessels on volume rendered images by rendering their names on them using a two-pass rendering process. In the first rendering pass, vessel surface models are generated using such properties as centerlines, radii, and running directions. Then the vessel names are drawn on the vessel surfaces. Finally, the vessel name images and the corresponding depth buffer are generated by a virtual camera at the viewpoint. In the second rendering pass, volume rendered images are generated by a ray casting volume rendering algorithm that considers the depth buffer generated in the first rendering pass. After the two-pass rendering is finished, an annotated image is generated by blending the volume rendered image with the surface rendered image. To confirm the effectiveness of our proposed method, we performed a computer-assisted system for the automated annotation of abdominal arteries. The experimental results show that vessel names can be drawn on the corresponding vessel surface in the volume rendered images at a computing cost that is nearly the same as that by volume rendering only. The proposed method has enormous potential to be adopted to annotate the vessels in the 3D medical images in clinical applications, such as image-guided surgery.

Computer-aided detection of clustered microcalcifications in digital breast tomosynthesis: a 3D approach

PURPOSE

To design a computer-aided detection (CADe) system for clustered microcalcifications in reconstructed digital breast tomosynthesis (DBT) volumes and to perform a preliminary evaluation of the CADe system.

METHODS

IRB approval and informed consent were obtained in this study. A data set of two-view DBT of 72 breasts containing microcalcification clusters was collected from 72 subjects who were scheduled to undergo breast biopsy. Based on tissue sampling results, 17 cases had breast cancer and 55 were benign. A separate data set of two-view DBT of 38 breasts free of clustered microcalcifications from 38 subjects was collected to independently estimate the number of false-positives (FPs) generated by the CADe system. A radiologist experienced in breast imaging marked the biopsied cluster of microcalcifications with a 3D bounding box using all available clinical and imaging information. A CADe system was designed to detect microcalcification clusters in the reconstructed volume. The system consisted of prescreening, clustering, and false-positive reduction stages. In the prescreening stage, the conspicuity of microcalcification-like objects was increased by an enhancement-modulated 3D calcification response function. An iterative thresholding and 3D object growing method was used to detect cluster seed objects, which were used as potential centers of microcalcification clusters. In the cluster detection stage, microcalcification candidates were identified using a second iterative thresholding procedure, which was applied to the signal-to-noise ratio (SNR) enhanced image voxels with a positive calcification response. Starting with each cluster seed object as the initial cluster center, a dynamic clustering algorithm formed a cluster candidate by including microcalcification candidates within a 3D neighborhood of the cluster seed object that satisfied the clustering criteria. The number, size, and SNR of the microcalcifications in a cluster candidate and the cluster shape were used to reduce the number of FPs.

RESULTS

The prescreening stage detected a cluster seed object in 94% of the biopsied microcalcification clusters at a threshold of 100 cluster seed objects per DBT volume. After clustering, the detection sensitivity was 90% at 15 marks per DBT volume. After FP reduction, at 85% sensitivity, the average number of FPs estimated using the data set containing microcalcification clusters was 3.8 per DBT volume, and that estimated using the data set free of microcalcification clusters was 3.4. The detection performance for malignant microcalcification clusters was superior to that for benign clusters.

CONCLUSIONS

Our study indicates the feasibility of the 3D approach to the detection of clustered microcalcifications in DBT and that the newly designed enhancement-modulated 3D calcification response function is promising for prescreening. Further work is needed to assess the generalizability of our approach and to improve its performance.

Semi-automatic surface and volume mesh generation for subject-specific biomedical geometries

An overview of surface and volume mesh generation techniques for creating valid meshes to carry out biomedical flows is provided. The methods presented are designed for robust numerical modelling of biofluid flow through subject-specific geometries. The applications of interest are haemodynamics in blood vessels and air flow in upper human respiratory tract. The methods described are designed to minimize distortion to a given domain boundary. They are also designed to generate a triangular surface mesh first and then volume mesh (tetrahedrons) with high quality surface and volume elements. For blood flow applications, a simple procedure to generate a boundary layer mesh is also described. The methods described here are semi-automatic in nature because of the fact that the geometries are complex, and automation of the procedures may be possible if high quality scans are used.