A study on the computerized fractal analysis of architectural distortion in screening mammograms

[Radiomics: Big Data Instead of Biopsies in the Future?]

Precision medicine is increasingly pushed forward, also with respect to upcoming new targeted therapies. Individual characterization of diseases on the basis of biomarkers is a prerequisite for this development. So far, biomarkers are characterized clinically, histologically or on a molecular level. The implementation of broad screening methods (“Omics”) and the analysis of big data – in addition to single markers – allow to define biomarker signatures. Next to “Genomics”, “Proteomics”, and “Metabolicis”, “Radiomics” gained increasing interest during the last years. Based on radiologic imaging, multiple radiomic markers are extracted with the help of specific algorithms. These are correlated with clinical, (immuno-) histopathological, or genomic data. Underlying structural differences are based on the imaging metadata and are often not visible and therefore not detectable without specific software. Radiomics are depicted numerically or by graphs. The fact that radiomic information can be extracted from routinely performed imaging adds a specific appeal to this method. Radiomics could potentially replace biopsies and additional investigations. Alternatively, radiomics could complement other biomarkers and thus lead to a more precise, multimodal prediction. Until now, radiomics are primarily used to investigate solid tumors. Some promising studies in head and neck cancer have already been published.

A computer-aided detection of the architectural distortion in digital mammograms using the fractal dimension measurements of BEMD

Achieving a high performance for the detection and characterization of architectural distortion in screening mammograms is important for an efficient breast cancer early detection. Viewing a mammogram image as a rough surface that can be described using the fractal theory is a well-recognized approach. This paper presents a new fractal-based computer-aided detection (CAD) algorithm for characterizing various breast tissues in screening mammograms with a particular focus on distinguishing between architectural distortion and normal breast parenchyma. The proposed approach is based on two underlying assumptions: (i) monitoring the variation pattern of fractal dimension, with the changes of the image resolution, is a useful tool to distinguish textural patterns of breast tissue, (ii) the bidimensional empirical mode decomposition (BEMD) algorithm appropriately generates a multiresolution representation of the mammogram. The proposed CAD has been tested using different validation datasets of mammographic regions of interest (ROIs) extracted from the Digital Database for Screening Mammography (DDSM) database. The validation ROI datasets contain architectural distortion (AD), normal breast tissue, and AD surrounding tissue. The highest classification performance, in terms of area under the receiver operating characteristic curve, of Az = 0.95 was achieved when the proposed approach applied to distinguish 187 architectural distortion depicting regions from 2191 normal breast parenchyma regions. The obtained results validate the underlying hypothesis and demonstrate that effectiveness of capturing the variation of the fractal dimension measurements within an appropriate multiscale representation of the digital mammogram. Results also reveal that this tool has the potential of prescreening other key and common mammographic signs of early breast cancer.

A novel method for quantification of beam's-eye-view tumor tracking performance

PURPOSE

In-treatment imaging using an electronic portal imaging device (EPID) can be used to confirm patient and tumor positioning. Real-time tumor tracking performance using current digital megavolt (MV) imagers is hindered by poor image quality. Novel EPID designs may help to improve quantum noise response, while also preserving the high spatial resolution of the current clinical detector. Recently investigated EPID design improvements include but are not limited to multi-layer imager (MLI) architecture, thick crystalline and amorphous scintillators, and phosphor pixilation and focusing. The goal of the present study was to provide a method of quantitating improvement in tracking performance as well as to reveal the physical underpinnings of detector design that impact tracking quality. The study employs a generalizable ideal observer methodology for the quantification of tumor tracking performance. The analysis is applied to study both the effect of increasing scintillator thickness on a standard, single-layer imager (SLI) design as well as the effect of MLI architecture on tracking performance.

METHODS

The present study uses the ideal observer signal-to-noise ratio (d') as a surrogate for tracking performance. We employ functions which model clinically relevant tasks and generalized frequency-domain imaging metrics to connect image quality with tumor tracking. A detection task for relevant Cartesian shapes (i.e., spheres and cylinders) was used to quantitate trackability of cases employing fiducial markers. Automated lung tumor tracking algorithms often leverage the differences in benign and malignant lung tissue textures. These types of algorithms (e.g., soft-tissue localization - STiL) were simulated by designing a discrimination task, which quantifies the differentiation of tissue textures, measured experimentally and fit as a power-law in trend (with exponent β) using a cohort of MV images of patient lungs. The modeled MTF and NPS were used to investigate the effect of scintillator thickness and MLI architecture on tumor tracking performance.

RESULTS

Quantification of MV images of lung tissue as an inverse power-law with respect to frequency yields exponent values of β = 3.11 and 3.29 for benign and malignant tissues, respectively. Tracking performance with and without fiducials was found to be generally limited by quantum noise, a factor dominated by quantum detective efficiency (QDE). For generic SLI construction, increasing the scintillator thickness (gadolinium oxysulfide - GOS) from a standard 290 μm to 1720 μm reduces noise to about 10%. However, 81% of this reduction is appreciated between 290 and 1000 μm. In comparing MLI and SLI detectors of equivalent individual GOS layer thickness, the improvement in noise is equal to the number of layers in the detector (i.e., 4) with almost no difference in MTF. Further, improvement in tracking performance was slightly less than the square-root of the reduction in noise, approximately 84-90%. In comparing an MLI detector with an SLI with a GOS scintillator of equivalent total thickness, improvement in object detectability is approximately 34-39%.

CONCLUSIONS

We have presented a novel method for quantification of tumor tracking quality and have applied this model to evaluate the performance of SLI and MLI EPID designs. We showed that improved tracking quality is primarily limited by improvements in NPS. When compared to very thick scintillator SLI, employing MLI architecture exhibits the same gains in QDE, but by mitigating the effect of optical Swank noise, results in more dramatic improvements in tracking performance.

Characterization of Architectural Distortion in Mammograms Based on Texture Analysis Using Support Vector Machine Classifier with Clinical Evaluation

Architecture distortion (AD) is an important and early sign of breast cancer, but due to its subtlety, it is often missed on the screening mammograms. The objective of this study is to create a quantitative approach for texture classification of AD based on various texture models, using support vector machine (SVM) classifier. The texture analysis has been done on the region of interest (ROI) selected from the original mammogram. A comprehensive analysis has been done on samples from three databases; out of which, two data sets are from the public domain, and the third data set is for clinical evaluation. The public domain databases are IRMA version of digital database for screening mammogram (DDSM) and Mammographic Image Analysis Society (MIAS). For clinical evaluation, the actual patient's database has been obtained from ACE Healthways, Diagnostic Centre Ludhiana, India. The significant finding of proposed study lies in appropriate selection of the size of ROIs. The experiments have been done on fixed size of ROIs as well as on the ground truth (variable size) ROIs. Best results pertain to an accuracy of 92.94 % obtained in case of DDSM database for fixed-size ROIs. In case of MIAS database, an accuracy of 95.34 % is achieved in AD versus non-AD (normal) cases for ground truth ROIs. Clinically, an accuracy of 88 % was achieved for ACE dataset. The results obtained in the present study are encouraging, as optimal result has been achieved for the proposed study in comparison with other related work in the same area.

Characterizing Architectural Distortion in Mammograms by Linear Saliency

Architectural distortion (AD) is a common cause of false-negatives in mammograms. This lesion usually consists of a central retraction of the connective tissue and a spiculated pattern radiating from it. This pattern is difficult to detect due the complex superposition of breast tissue. This paper presents a novel AD characterization by representing the linear saliency in mammography Regions of Interest (ROI) as a graph composed of nodes corresponding to locations along the ROI boundary and edges with a weight proportional to the line intensity integrals along the path connecting any pair of nodes. A set of eigenvectors from the adjacency matrix is then used to extract discriminant coefficients that represent those nodes with higher salient lines. A dimensionality reduction is further accomplished by selecting the pair of nodes with major contribution for each of the computed eigenvectors. The set of main salient lines is then assembled as a feature vector that inputs a conventional Support Vector Machine (SVM). Experimental results with two benchmark databases, the mini-MIAS and DDSM databases, demonstrate that the proposed linear saliency domain method (LSD) performs well in terms of accuracy. The approach was evaluated with a set of 246 RoI extracted from the DDSM (123 normal tissues and 123 AD) and a set of 38 ROI from the mini-MIAS collections (19 normal tissues and 19 AD) respectively. The classification results showed respectively for both databases an accuracy rate of 89 % and 87 %, a sensitivity rate of 85 % and 95 %, and a specificity rate of 93 % and 84 %. Likewise, the area under curve (A _z ) of the Receiver Operating Characteristic (ROC) curve was 0.93 for both databases.

A new approach for the detection of architectural distortions using textural analysis of surrounding tissue

Aiming at improving the performance of computer-aided detection of architectural distortion (AD) in mammograms, this paper investigates whether textural patterns of AD surrounding tissue (ST) have the potential to detect AD signatures. More specifically, for characterizing the presence of AD; we investigated the application of textural analysis for discriminating between AD surrounding tissue and normal breast parenchyma. We evaluated the underlying hypothesis using a dataset of 2544 regions of interest (ROI) obtained from the Digital Database for Screening Mammography (DDSM). The ROI dataset contained 353 ST regions and 2191 normal parenchyma related regions. The bidimensional empirical mode decomposition (BEMD) algorithm was, first, applied to extract, from each ROI, the 2D intrinsic mode functions (2DIMF) or detail subbands. Then, statistical signatures of IMF layers were computed and used along with the fractal dimension, estimated from the original ROI, for discriminating ST from the normal breast tissue. The statistical analysis of various textural descriptors demonstrated the significant difference between characteristics of AD surrounding tissue and normal breast parenchyma. The highest AD recognition results of Az = 0.869, obtained from the textural analysis of AD surrounding tissue, is very promising and comparable with Az = 0.913 produced from characterizing AD regions.

Radiomics: a new application from established techniques

The increasing use of biomarkers in cancer have led to the concept of personalized medicine for patients. Personalized medicine provides better diagnosis and treatment options available to clinicians. Radiological imaging techniques provide an opportunity to deliver unique data on different types of tissue. However, obtaining useful information from all radiological data is challenging in the era of "big data". Recent advances in computational power and the use of genomics have generated a new area of research termed Radiomics. Radiomics is defined as the high throughput extraction of quantitative imaging features or texture (radiomics) from imaging to decode tissue pathology and creating a high dimensional data set for feature extraction. Radiomic features provide information about the gray-scale patterns, inter-pixel relationships. In addition, shape and spectral properties can be extracted within the same regions of interest on radiological images. Moreover, these features can be further used to develop computational models using advanced machine learning algorithms that may serve as a tool for personalized diagnosis and treatment guidance.

An automated method for detecting architectural distortions on mammograms using direction analysis of linear structures

Architectural distortion is one of the most important findings when evaluating mammograms for breast cancer. Abnormal breast architecture is characterized by the presence of spicules, which are distorted mammary structures that are not accompanied by an increased density or mass. We have been developing an automated method for detecting spiculated architectural distortions by analyzing linear structures extracted by normal curvature. However, some structures that are possibly related to distorted areas are not extracted using this method. The purpose of this study was to develop a new automated method for direction analysis of linear structures to improve detection performance in mammography. The direction of linear structures in each region of interest (ROI) was first determined using a direction filter and a background filter that can define one of eight directions (0°, 22.5°, 45°, 67.5°, 90°, 112.5°, 135°, and 157.5°). The concentration and isotropic indexes were calculated using the determined direction of the linear structures in order to extract the candidate areas. Discriminant analysis was performed to eliminate false positives results. Our database consisted of 168 abnormal images containing 174 distorted areas and 580 normal images. The sensitivity of the new method was 81%. There were 2.6 and 4.2 false positives per image using the new and previous methods, respectively. These findings show that our new method is effective for detecting spiculated architectural distortions.

In the digital era, architectural distortion remains a challenging radiological task

AIM

To compare readers' performance in detecting architectural distortion (AD) compared with other breast cancer types using digital mammography.

MATERIALS AND METHODS

Forty-one experienced breast screen readers (20 US and 21 Australian) were asked to read a single test set of 30 digitally acquired mammographic cases. Twenty cases had abnormal findings (10 with AD, 10 non-AD) and 10 cases were normal. Each reader was asked to locate and rate any abnormalities. Lesion and case-based performance was assessed. For each collection of readers (US; Australian; combined), jackknife free-response receiver operating characteristic (JAFROC), figure of merit (FOM), and inferred receiver operating characteristic (ROC), area under curve (Az) were calculated using JAFROC v.4.1 software. Readers' sensitivity, location sensitivity, JAFROC, FOM, ROC, Az scores were compared between cases groups using Wilcoxon's signed ranked test statistics.

RESULTS

For lesion-based analysis, significantly lower location sensitivity (p=0.001) was shown on AD cases compared with non-AD cases for all reader collections. The case-based analysis demonstrated significantly lower ROC Az values (p=0.02) for the first collection of readers, and lower sensitivity for the second collection of readers (p=0.04) and all-readers collection (p=0.008), for AD compared with non-AD cases.

CONCLUSIONS

The current work demonstrates that AD remains a challenging task for readers, even in the digital era.

An efficient fractal method for detection and diagnosis of breast masses in mammograms

In this paper, we present an efficient fractal method for detection and diagnosis of mass lesion in mammogram which is one of the abnormalities in mammographic images. We used 110 images that were carefully selected by a radiologist, and their abnormalities were also confirmed by biopsy. These images included circumscribed benign, ill-defined, and spiculated malignant masses. Firstly, we discriminated lesions automatically using new fractal dimensions. The results which were examined by different types of breast density showed that the proposed method was able to yield quite satisfactory detection results. Secondly, noting that contours of masses playing the most important role in diagnosis of different mass types, we defined new fractal features based on information extraction from the contours. This information is able to identify the roughness in mass contours and determines the extent of spiculation or smoothness of the masses. In this manner, in classification of the spiculated malignant masses from the circumscribed benign tumors, we achieved highly satisfactory results, i.e., 0.98 measured in terms of area under ROC curve (AUC). In this paper, it is also shown that the roughness in contours is a suitable characteristic feature for diagnosis of ill-defined malignant tumors with AUC equal to 0.94 in their classification. The extracted information was also found to be useful in the classification of early malignancies whereas in the classification of spiculated and ill-defined malignant masses in their early stage from those of benign tumors, we achieved high accuracy of 0.99 and 0.90 for AUC, respectively.

Detection of architectural distortion in prior mammograms via analysis of oriented patterns

We demonstrate methods for the detection of architectural distortion in prior mammograms of interval-cancer cases based on analysis of the orientation of breast tissue patterns in mammograms. We hypothesize that architectural distortion modifies the normal orientation of breast tissue patterns in mammographic images before the formation of masses or tumors. In the initial steps of our methods, the oriented structures in a given mammogram are analyzed using Gabor filters and phase portraits to detect node-like sites of radiating or intersecting tissue patterns. Each detected site is then characterized using the node value, fractal dimension, and a measure of angular dispersion specifically designed to represent spiculating patterns associated with architectural distortion. Our methods were tested with a database of 106 prior mammograms of 56 interval-cancer cases and 52 mammograms of 13 normal cases using the features developed for the characterization of architectural distortion, pattern classification via quadratic discriminant analysis, and validation with the leave-one-patient out procedure. According to the results of free-response receiver operating characteristic analysis, our methods have demonstrated the capability to detect architectural distortion in prior mammograms, taken 15 months (on the average) before clinical diagnosis of breast cancer, with a sensitivity of 80% at about five false positives per patient.

Measures of divergence of oriented patterns for the detection of architectural distortion in prior mammograms

PURPOSE

We propose a method for the detection of architectural distortion in prior mammograms of interval-cancer cases based on the expected orientation of breast tissue patterns in mammograms.

METHODS

The expected orientation of the breast tissue at each pixel was derived by using automatically detected landmarks including the breast boundary, the nipple, and the pectoral muscle (in mediolateral-oblique views). We hypothesize that the presence of architectural distortion changes the normal expected orientation of breast tissue patterns in a mammographic image. The angular deviation of the oriented structures in a given mammogram as compared to the expected orientation was analyzed to detect potential sites of architectural distortion using a measure of divergence of oriented patterns. Each potential site of architectural distortion was then characterized using measures of spicularity and angular dispersion specifically designed to represent spiculating patterns. The novel features for the characterization of spiculating patterns include an index of divergence of spicules computed from the intensity image and Gabor magnitude response using the Gabor angle response; radially weighted difference and angle-weighted difference (AWD) measures of the intensity, Gabor magnitude, and Gabor angle response; and AWD in the entropy of spicules computed from the intensity, Gabor magnitude, and Gabor angle response.

RESULTS

Using the newly proposed features with a database of 106 prior mammograms of 56 interval-cancer cases and 52 mammograms of 13 normal cases, through feature selection and pattern classification with an artificial neural network, an area under the receiver operating characteristic curve of 0.75 was obtained. Free-response receiver operating characteristic analysis indicated a sensitivity of 0.80 at 5.3 false positives (FPs) per patient. Combining the features proposed in the present paper with others described in our previous works led to significant improvement with a sensitivity of 0.80 at 3.7 FPs per patient.

CONCLUSION

The proposed methods can detect architectural distortion in prior mammograms taken 15 months (on the average) before clinical diagnosis of breast cancer, but the FP rate needs to be reduced.

Measures of angular spread and entropy for the detection of architectural distortion in prior mammograms

PURPOSE

Architectural distortion is an important sign of early breast cancer. We present methods for computer-aided detection of architectural distortion in mammograms acquired prior to the diagnosis of breast cancer in the interval between scheduled screening sessions.

METHODS

Potential sites of architectural distortion were detected using node maps obtained through the application of a bank of Gabor filters and linear phase portrait modeling. A total of 4,224 regions of interest (ROIs) were automatically obtained from 106 prior mammograms of 56 interval-cancer cases, including 301 true-positive ROIs, and from 52 mammograms of 13 normal cases. Each ROI was represented by three types of entropy measures of angular histograms composed with the Gabor magnitude response, angle, coherence, orientation strength, and the angular spread of power in the Fourier spectrum, including Shannon's entropy, Tsallis entropy for nonextensive systems, and Rényi entropy for extensive systems.

RESULTS

Using the entropy measures with stepwise logistic regression and the leave-one-patient-out method for feature selection and cross-validation, an artificial neural network resulted in an area under the receiver operating characteristic curve of 0.75. Free-response receiver operating characteristics indicated a sensitivity of 0.80 at 5.2 false positives (FPs) per patient.

CONCLUSION

The proposed methods can detect architectural distortion in prior mammograms taken 15 months (on the average) before clinical diagnosis of breast cancer, with a high sensitivity and a moderate number of FPs per patient. The results are promising and may be improved with additional features to characterize subtle abnormalities and larger databases including prior mammograms.

Three-dimensional in silico breast phantoms for multimodal image simulations

Anatomic simulators have provided researchers with the realistic objects needed to develop and evaluate medical imaging approaches. Today we have new insights into the cellular biology of breast tissues that is driving many new targeted diagnostic and therapeutic approaches, including molecular imaging. We report on our initial efforts to build a scalable framework for the construction of realistic 3-D in silico breast phantoms (ISBP) capable of leveraging existing knowledge and yet adaptable to fully integrate future discoveries. The ISBP frames are developed with scalable anatomical shapes and morphologic features as adapted from a rich literature on this topic. Frames are populated with tissue subtypes essential for imaging and object contrast functions are assigned. These data can be resampled to match the intrinsics scales of various imaging modalities; we explore mammography, sonography and computed tomography. Initial comparisons between simulated and clinical images demonstrate reasonable agreement and provides guidance for future development of a more realistic ISBP. An end-to-end simulation of breast images is described to demonstrate techniques for including stochastic variability and deterministic physical principles on which image formation is based.

Fractal dimension as a measure of altered actin cytoskeleton in MC3T3-E1 cells under simulated microgravity using 3-D/2-D clinostats

Osteoblasts, the bone-forming cells, respond to various mechanical forces, such as stretch and fluid shear force in essentially similar ways. The cytoskeleton, as the load-bearing architecture of the cell, is sensitive to altered inertial forces. Disruption of the cytoskeleton will result in alteration of cellular structure and function. However, it is difficult to quantitatively illustrate cytoskeletal rearrangement because of the complexity of cytoskeletal structure. Usually, the morphological changes in actin organization caused by external stimulus are basically descriptive. In this study, fractal dimensions (D) analysis was used to quantify the morphological changes in the actin cytoskeleton of osteoblast-like cells (MC3T3-E1) under simulated microgravity using 3-D/2-D clinostats. The ImageJ software was used to count the fractal dimension of actin cytoskeleton by box-counting methods. Real-time PCR and immunofluroscent assays were used to further confirm the results obtained by fractal dimension analysis. The results showed significant decreases in D value of actin cytoskeleton, β-actin mRNA expression, and the mean fluorescence intensity of F-actin in osteoblast-like cells after 24 or 48 h of incubation under 3-D/2-D clinorotation condition compared with control. The findings indicate that 3-D/2-D clinorotation affects both actin cytoskeleton architecture and mRNA expression, and fractal may be a promising approach for quantitative analysis of the changes in cytoskeleton in different environments.

Detection of architectural distortion in prior mammograms

We present methods for the detection of sites of architectural distortion in prior mammograms of interval-cancer cases. We hypothesize that screening mammograms obtained prior to the detection of cancer could contain subtle signs of early stages of breast cancer, in particular, architectural distortion. The methods are based upon Gabor filters, phase portrait analysis, a novel method for the analysis of the angular spread of power, fractal analysis, Laws' texture energy measures derived from geometrically transformed regions of interest (ROIs), and Haralick's texture features. With Gabor filters and phase portrait analysis, 4224 ROIs were automatically obtained from 106 prior mammograms of 56 interval-cancer cases, including 301 true-positive ROIs related to architectural distortion, and from 52 mammograms of 13 normal cases. For each ROI, the fractal dimension, the entropy of the angular spread of power, 10 Laws' measures, and Haralick's 14 features were computed. The areas under the receiver operating characteristic curves obtained using the features selected by stepwise logistic regression and the leave-one-ROI-out method are 0.76 with the Bayesian classifier, 0.75 with Fisher linear discriminant analysis, and 0.78 with a single-layer feed-forward neural network. Free-response receiver operating characteristics indicated sensitivities of 0.80 and 0.90 at 5.8 and 8.1 false positives per image, respectively, with the Bayesian classifier and the leave-one-image-out method.

Evaluation of an improved algorithm for producing realistic 3D breast software phantoms: application for mammography

PURPOSE

This work presents an improved algorithm for the generation of 3D breast software phantoms and its evaluation for mammography.

METHODS

The improved methodology has evolved from a previously presented 3D noncompressed breast modeling method used for the creation of breast models of different size, shape, and composition. The breast phantom is composed of breast surface, duct system and terminal ductal lobular units, Cooper's ligaments, lymphatic and blood vessel systems, pectoral muscle, skin, 3D mammographic background texture, and breast abnormalities. The key improvement is the development of a new algorithm for 3D mammographic texture generation. Simulated images of the enhanced 3D breast model without lesions were produced by simulating mammographic image acquisition and were evaluated subjectively and quantitatively. For evaluation purposes, a database with regions of interest taken from simulated and real mammograms was created. Four experienced radiologists participated in a visual subjective evaluation trial, as they judged the quality of the simulated mammograms, using the new algorithm compared to mammograms, obtained with the old modeling approach. In addition, extensive quantitative evaluation included power spectral analysis and calculation of fractal dimension, skewness, and kurtosis of simulated and real mammograms from the database.

RESULTS

The results from the subjective evaluation strongly suggest that the new methodology for mammographic breast texture creates improved breast models compared to the old approach. Calculated parameters on simulated images such as beta exponent deducted from the power law spectral analysis and fractal dimension are similar to those calculated on real mammograms. The results for the kurtosis and skewness are also in good coincidence with those calculated from clinical images. Comparison with similar calculations published in the literature showed good agreement in the majority of cases.

CONCLUSIONS

The improved methodology generated breast models with increased realism compared to the older model as shown in evaluations of simulated images by experienced radiologists. It is anticipated that the realism will be further improved using an advanced image simulator so that simulated images may be used in feasibility studies in mammography.

Computer-aided detection of architectural distortion in prior mammograms of interval cancer

Architectural distortion is an important sign of breast cancer, but because of its subtlety, it is a common cause of false-negative findings on screening mammograms. This paper presents methods for the detection of architectural distortion in mammograms of interval cancer cases taken prior to the detection of breast cancer using Gabor filters, phase portrait analysis, fractal analysis, and texture analysis. The methods were used to detect initial candidates for sites of architectural distortion in prior mammograms of interval cancer and also normal control cases. A total of 4,224 regions of interest (ROIs) were automatically obtained from 106 prior mammograms of 56 interval cancer cases, including 301 ROIs related to architectural distortion, and from 52 prior mammograms of 13 normal cases. For each ROI, the fractal dimension and Haralick's texture features were computed. Feature selection was performed separately using stepwise logistic regression and stepwise regression. The best results achieved, in terms of the area under the receiver operating characteristics curve, with the features selected by stepwise logistic regression are 0.76 with the Bayesian classifier, 0.73 with Fisher linear discriminant analysis, 0.77 with an artificial neural network based on radial basis functions, and 0.77 with a support vector machine. Analysis of the performance of the methods with free-response receiver operating characteristics indicated a sensitivity of 0.80 at 7.6 false positives per image. The methods have good potential in detecting architectural distortion in mammograms of interval cancer cases.

Fractal-feature distance analysis of contrast-detail phantom image and meaning of pseudo fractal dimension and complexity

The purposes of our studies are to examine whether or not fractal-feature distance deduced from virtual volume method can simulate observer performance indices and to investigate the physical meaning of pseudo fractal dimension and complexity. Contrast-detail (C-D) phantom radiographs were obtained at various mAs values (0.5 - 4.0 mAs) and 140 kVp with a computed radiography system, and the reference image was acquired at 13 mAs. For all C-D images, fractal analysis was conducted using the virtual volume method that was devised with a fractional Brownian motion model. The fractal-feature distances between the considered and reference images were calculated using pseudo fractal dimension and complexity. Further, we have performed the C-D analysis in which ten radiologists participated, and compared the fractal-feature distances with the image quality figures (IQF). To clarify the physical meaning of the pseudo fractal dimension and complexity, contrast-to-noise ratio (CNR) and standard deviation (SD) of images noise were calculated for each mAs and compared with the pseudo fractal dimension and complexity, respectively. A strong linear correlation was found between the fractal-feature distance and IQF. The pseudo fractal dimensions became large as CNR increased. Further, a linear correlation was found between the exponential complexity and image noise SD.

Methodology for generating a 3D computerized breast phantom from empirical data

The initial process for creating a flexible three-dimensional computer-generated breast phantom based on empirical data is described. Dedicated breast computed-tomography data were processed to suppress noise and scatter artifacts in the reconstructed image set. An automated algorithm was developed to classify the breast into its primary components. A preliminary phantom defined using subdivision surfaces was generated from the segmented data. To demonstrate potential applications of the phantom, simulated mammographic image data were acquired of the phantom using a simplistic compression model and an analytic projection algorithm directly on the surface model. The simulated image was generated using a model for a polyenergetic cone-beam projection of the compressed phantom. The methods used to create the breast phantom generate resulting images that have a high level of tissue structure detail available and appear similar to actual mammograms. Fractal dimension measurements of simulated images of the phantom are comparatively similar to measurements from images of real human subjects. A realistic and geometrically defined breast phantom that can accurately simulate imaging data may have many applications in breast imaging research.

Computerized characterization of prostate cancer by fractal analysis in MR images

PURPOSE

To explore the potential of computerized characterization of prostate MR images by extracting the fractal features of texture and intensity distributions as indices in the differential diagnosis of prostate cancer.

MATERIALS AND METHODS

MR T2-weighted images (T2WI) of 55 patients with pathologic results detected by ultrasound guided biopsy were collected and then divided in two groups, 27 with prostate cancer (PCa) and 28 with no histological abnormality. Texture fractal dimension (TFD) and histogram fractal dimension (HFD) were calculated to analyze complexity features of regions of Interest (ROIs) selected from the peripheral zone. Two-sample t-tests were performed to evaluate group differences for both parameters. Receiver operating characteristic (ROC) analysis was used to estimate the performance of TFD and HFD for discriminating PCa.

RESULTS

Significant differences were found in both TFD and HFD between the two patient groups. The areas under the ROC curves of TFD and HFD were 0.691 and 0.966, respectively, in distinguishing prostatic carcinoma from normal peripheral zone. As characterized by the fractal indices, cancerous prostatic tissue exhibited smoother texture and lower variation in intensity distribution than normal prostatic tissue.

CONCLUSION

The study suggests that TFD and HFD depict the changes in texture and intensity distribution associated with prostate cancer on T2WI. Both TFD and HFD provide promising quantitative indices for cancer identification. HFD performs better than TFD offering a more robust MR-based indicator in the diagnosis of prostatic carcinoma.

Characterization and classification of tumor lesions using computerized fractal-based texture analysis and support vector machines in digital mammograms

OBJECTIVE

This paper presents a detailed study of fractal-based methods for texture characterization of mammographic mass lesions and architectural distortion. The purpose of this study is to explore the use of fractal and lacunarity analysis for the characterization and classification of both tumor lesions and normal breast parenchyma in mammography.

MATERIALS AND METHODS

We conducted comparative evaluations of five popular fractal dimension estimation methods for the characterization of the texture of mass lesions and architectural distortion. We applied the concept of lacunarity to the description of the spatial distribution of the pixel intensities in mammographic images. These methods were tested with a set of 57 breast masses and 60 normal breast parenchyma (dataset1), and with another set of 19 architectural distortions and 41 normal breast parenchyma (dataset2). Support vector machines (SVM) were used as a pattern classification method for tumor classification.

RESULTS

Experimental results showed that the fractal dimension of region of interest (ROIs) depicting mass lesions and architectural distortion was statistically significantly lower than that of normal breast parenchyma for all five methods. Receiver operating characteristic (ROC) analysis showed that fractional Brownian motion (FBM) method generated the highest area under ROC curve (A ( z ) = 0.839 for dataset1, 0.828 for dataset2, respectively) among five methods for both datasets. Lacunarity analysis showed that the ROIs depicting mass lesions and architectural distortion had higher lacunarities than those of ROIs depicting normal breast parenchyma. The combination of FBM fractal dimension and lacunarity yielded the highest A ( z ) value (0.903 and 0.875, respectively) than those based on single feature alone for both given datasets. The application of the SVM improved the performance of the fractal-based features in differentiating tumor lesions from normal breast parenchyma by generating higher A ( z ) value.

CONCLUSION

FBM texture model is the most appropriate model for characterizing mammographic images due to self-affinity assumption of the method being a better approximation. Lacunarity is an effective counterpart measure of the fractal dimension in texture feature extraction in mammographic images. The classification results obtained in this work suggest that the SVM is an effective method with great potential for classification in mammographic image analysis.

A pilot study of architectural distortion detection in mammograms based on characteristics of line shadows

OBJECTIVE

We present herein a novel algorithm for architectural distortion detection that utilizes the point convergence index with the likelihood of lines (e.g., spiculations) relating to architectural distortion.

MATERIALS AND METHODS

Validation was performed using 25 computed radiography (CR) mammograms, each of which has an architectural distortion with radiating spiculations. The proposed method comprises five steps. First, the lines were extracted on mammograms, such as spiculations of architectural distortion as well as lines in the mammary gland. Second, the likelihood of spiculation for each extracted line was calculated. In the third step, point convergence index weighted by this likelihood was evaluated at each pixel to enhance distortion only. Fourth, local maxima of the index were extracted as candidates for the distortion, then classified based on nine features in the last step.

RESULTS

Point convergence index without the proposed likelihood generated 84.48/image false-positives (FPs) on average. Conversely, the proposed index succeeded in decreasing this number to 12.48/image on average when sensitivity was 100%. After the classification step, number of FPs was reduced to 0.80/image with 80.0% sensitivity.

CONCLUSION

Combination of the likelihood of lines with point convergence index is effective in extracting architectural distortion with radiating spiculations.