The shape of breast cancer.

PURPOSE: Little is known about the three-dimensional shape of breast cancer. Implicit to approaches that localize the center of the tumor for breast-conserving surgery (BCS) of non-palpable cancers is the assumption that breast cancers are spherical about a central point, which may not be accurate. METHODS: Pre-operative supine breast MRI images were obtained of 83 breast cancer patients undergoing partial mastectomy using supine MRI-guided resection techniques. Three-dimensional (3D) tumor models were derived after radiologists outlined tumor edges on successive MRI slices. Ideal resection volumes were determined by adding 1 cm in every dimension to the actual tumor volume. Geometrically defined parameters were used to define tumor shapes and associations between clinical variables and shapes were examined. RESULTS: Seventy-five patients had invasive cancer. Breast cancers were categorized into four tumor shapes: 34% of tumors were discoidal, 29% segmental, 19% spherical, and 18% irregular. If hypothetical spherical excisions were performed, non-spherical cases would excise 143% more tissue than the ideal resection volume. When the 3D shape of each tumor was provided to the surgeon during MR-guided BCS, the percentage of tissue overexcised in non-spherical cases was significantly less (143% vs. 66%, p < 0.001). CONCLUSIONS: Information obtained from a supine MRI can be used to generate 3D tumor models and rapidly classify breast tumor shapes. The vast majority of invasive cancers and DCIS are not spherical. Knowledge of tumor shape may allow surgeons to excise breast cancer more precisely.

A Randomized Prospective Trial of Supine MRI-Guided Versus Wire-Localized Lumpectomy for Breast Cancer.

BACKGROUND: Wire-localized excision of non-palpable breast cancer is imprecise, resulting in positive margins 15-35% of the time. METHODS: Women with a confirmed diagnosis of non-palpable invasive breast cancer (IBC) or ductal carcinoma in situ (DCIS) were randomized to a new technique using preoperative supine magnetic resonance imaging (MRI) with intraoperative optical scanning and tracking (MRI group) or wire-localized (WL group) partial mastectomy. The main outcome measure was the positive margin rate. RESULTS: In this study, 138 patients were randomly assigned. Sixty-six percent had IBC and DCIS, 22% had IBC, and 12% had DCIS. There were no differences in patient or tumor characteristics between the groups. The proportion of patients with positive margins in the MRI-guided surgery group was half that observed in the WL group (12 vs. 23%; p = 0.08). The specimen volumes in the MRI and WL groups did not differ significantly (74 ± 33.9 mL vs. 69.8 ± 25.1 mL; p = 0.45). The pathologic tumor diameters were underestimated by 2 cm or more in 4% of the cases by MRI and in 9% of the cases by mammography. Positive margins were observed in 68% and 58% of the cases underestimated by 2 cm or more using MRI and mammography, respectively, and in 15% and 14% of the cases not underestimated using MRI and mammography, respectively. CONCLUSIONS: A novel system using supine MRI images co-registered with intraoperative optical scanning and tracking enabled tumors to be resected with a trend toward a lower positive margin rate compared with wire-localized partial mastectomy. Margin positivity was more likely when imaging underestimated pathologic tumor size.

Cherenkov imaging for linac beam shape analysis as a remote electronic quality assessment verification tool.

PURPOSE: A remote imaging system tracking Cherenkov emission was analyzed to verify that the linear accelerator (linac) beam shape could be quantitatively measured at the irradiation surface for Quality Audit (QA). METHODS: The Cherenkov camera recorded 2D dose images delivered on a solid acrylonitrile butadiene styrene (ABS) plastic phantom surface for a range of square beam sizes, and 6 MV photons. Imaging was done at source to surface distance (SSD) of 100 cm and compared to GaF film images and linac light fields of the same beam sizes, ranging over 5 × 5 cm2 up to 20 × 20 cm2 . Line profiles of each field were compared in both X and Y jaw directions. Each measurement was repeated on two different Clinac2100 machines. An interreader comparison of the beam width interpretation was completed using procedures commonly employed for beam to light field coincidence verification. Cherenkov measurements are also done for beams of complex treatment plan and isocenter QA. RESULTS: The Cherenkov image widths matched with the measured GaF images and light field images, with accuracy in the range of ±1 mm standard deviation. The differences between the measurements were minor and within tolerance of geometrical requirement of standard linac QA procedures conducted by human setup verification, which had a similar error range. The measurement made by the remote imaging system allowed for beam shape extraction of radiation fields at the SSD location of the beam. CONCLUSIONS: The proposed Cherenkov image acquisition system provides a valid way to remotely confirm radiation field sizes and provides similar information to that obtained from the linac light field or GaF film estimates of the beam size. The major benefit of this approach is that with a fixed installation of the camera, testing could be done completely under software control with automated image analysis, potentially simplifying conventional QA procedures with appropriate calibration of boundary definitions, and the natural extension to capturing dynamic treatment beamlets at SSD could have future value, such as verification of beam plans with complex beam shapes, like IMRT or "star-shot" QA for the isocenter.

A Patient-Specific 3D-Printed Form Accurately Transfers Supine MRI-Derived Tumor Localization Information to Guide Breast-Conserving Surgery.

BACKGROUND: Wire-localized excision of nonpalpable breast cancer is imprecise, resulting in positive margins 25-30% of the time. METHODS: Patients underwent preoperative supine magnetic resonance imaging (MRI). A radiologist outlined the tumor edges on consecutive images, creating a three-dimensional (3D) view of its location. Using 3D printing, a bra-like plastic form (the Breast Cancer Locator [BCL]) was fabricated, with features that allowed a surgeon to (1) mark the edges of the tumor on the breast surface; (2) inject blue dye into the breast 1 cm from the tumor edges; and (3) place a wire in the tumor at the time of surgery. RESULTS: Nineteen patients with palpable cancers underwent partial mastectomy after placement of surgical cues using patient-specific BCLs. The cues were in place in <5 min and no adverse events occurred. The BCL accurately localized 18/19 cancers. In the 18 accurately localized cases, all 68 blue-dye injections were outside of the tumor edges. Median distance from the blue-dye center to the pathologic tumor edge was 1.4 cm, while distance from the blue dye to the tumor edge was <5 mm in 4% of injections, 0.5-2.0 cm in 72% of injections, and >2 cm in 24% of injections. Median distance from the tumor center to the BCL-localized wire and to the clip placed at the time of diagnosis was similar (0.49 vs. 0.73 cm) on specimen mammograms. CONCLUSIONS: Information on breast cancer location and shape derived from a supine MRI can be transferred safely and accurately to patients in the operating room using a 3D-printed form.

Effects of breast density and compression on normal breast tissue hemodynamics through breast tomosynthesis guided near-infrared spectral tomography.

Optically derived tissue properties across a range of breast densities and the effects of breast compression on estimates of hemoglobin, oxygen metabolism, and water and lipid concentrations were obtained from a coregistered imaging system that integrates near-infrared spectral tomography (NIRST) with digital breast tomosynthesis (DBT). Image data were analyzed from 27 women who underwent four IRB approved NIRST/DBT exams that included fully and mildly compressed breast acquisitions in two projections­craniocaudal (CC) and mediolateral-oblique (MLO)­and generated four data sets per patient (full and moderate compression in CC and MLO views). Breast density was correlated with HbT (r=0.64, p=0.001), water (r=0.62, p=0.003), and lipid concentrations (r=?0.74, p<0.001), but not oxygen saturation. CC and MLO views were correlated for individual subjects and demonstrated no statistically significant differences in grouped analysis. Comparison of compressed and uncompressed imaging demonstrated a significant decrease in oxygen saturation under compression (58% versus 50%, p=0.04). Mammographic breast density categorization was correlated with measured optically derived properties.

Calibration and optimization of 3D digital breast tomosynthesis guided near infrared spectral tomography.

Calibration of a three-dimensional multimodal digital breast tomosynthesis (DBT) x-ray and non-fiber based near infrared spectral tomography (NIRST) system is challenging but essential for clinical studies. Phantom imaging results yielded linear contrast recovery of total hemoglobin (HbT) concentration for cylindrical inclusions of 15 mm, 10 mm and 7 mm with a 3.5% decrease in the HbT estimate for each 1 cm increase in inclusion depth. A clinical exam of a patient's breast containing both benign and malignant lesions was successfully imaged, with greater HbT was found in the malignancy relative to the benign abnormality and fibroglandular regions (11 µM vs. 9.5 µM). Tools developed improved imaging system characterization and optimization of signal quality, which will ultimately improve patient selection and subsequent clinical trial results.

Molecular dyes used for surgical specimen margin orientation allow for intraoperative optical assessment during breast conserving surgery.

A variety of optical techniques utilizing near-infrared (NIR) light are being proposed for intraoperative breast tumor margin assessment. However, immediately following a lumpectomy excision, the margins are inked, which preserves the orientation of the specimen but prevents optical interrogation of the tissue margins. Here, a workflow is proposed that allows for both NIR optical assessment following full specimen marking using molecular dyes which have negligible absorption and scattering in the NIR. The effect of standard surgical inks in contrast to molecular dyes for an NIR signal is shown. Further, the proposed workflow is demonstrated with full specimen intraoperative imaging on all margins directly after the lumpectomy has been excised and completely marked. This work is an important step in the path to clinical feasibility of intraoperative breast tumor margin assessment using NIR optical methods without having to compromise on the current clinical practice of inking resected specimens for margin orientation.

Digital Breast Tomosynthesis guided Near Infrared Spectroscopy: Volumetric estimates of fibroglandular fraction and breast density from tomosynthesis reconstructions.

A multimodality system combining a clinical prototype digital breast tomosynthesis with its imaging geometry modified to facilitate near-infrared spectroscopic imaging has been developed. The accuracy of parameters recovered from near-infrared spectroscopy is dependent on fibroglandular tissue content. Hence, in this study, volumetric estimates of fibroglandular tissue from tomosynthesis reconstructions were determined. A kernel-based fuzzy c-means algorithm was implemented to segment tomosynthesis reconstructed slices in order to estimate fibroglandular content and to provide anatomic priors for near-infrared spectroscopy. This algorithm was used to determine volumetric breast density (VBD), defined as the ratio of fibroglandular tissue volume to the total breast volume, expressed as percentage, from 62 tomosynthesis reconstructions of 34 study participants. For a subset of study participants who subsequently underwent mammography, VBD from mammography matched for subject, breast laterality and mammographic view was quantified using commercial software and statistically analyzed to determine if it differed from tomosynthesis. Summary statistics of the VBD from all study participants were compared with prior independent studies. The fibroglandular volume from tomosynthesis and mammography were not statistically different (p=0.211, paired t-test). After accounting for the compressed breast thickness, which were different between tomosynthesis and mammography, the VBD from tomosynthesis was correlated with (r =0.809, p<0.001), did not statistically differ from (p>0.99, paired t-test), and was linearly related to, the VBD from mammography. Summary statistics of the VBD from tomosynthesis were not statistically different from prior studies using high-resolution dedicated breast computed tomography. The observation of correlation and linear association in VBD between mammography and tomosynthesis suggests that breast density associated risk measures determined for mammography are translatable to tomosynthesis. Accounting for compressed breast thickness is important when it differs between the two modalities. The fibroglandular volume from tomosynthesis reconstructions is similar to mammography indicating suitability for use during near-infrared spectroscopy.

Sub-diffusive scattering parameter maps recovered using wide-field high-frequency structured light imaging.

This study investigates the hypothesis that structured light reflectance imaging with high spatial frequency patterns [Formula: see text] can be used to quantitatively map the anisotropic scattering phase function distribution [Formula: see text] in turbid media. Monte Carlo simulations were used in part to establish a semi-empirical model of demodulated reflectance ([Formula: see text]) in terms of dimensionless scattering [Formula: see text] and [Formula: see text], a metric of the first two moments of the [Formula: see text] distribution. Experiments completed in tissue-simulating phantoms showed that simultaneous analysis of [Formula: see text] spectra sampled at multiple [Formula: see text] in the frequency range [0.05-0.5] [Formula: see text] allowed accurate estimation of both [Formula: see text] in the relevant tissue range [0.4-1.8] [Formula: see text], and [Formula: see text] in the range [1.4-1.75]. Pilot measurements of a healthy volunteer exhibited [Formula: see text]-based contrast between scar tissue and surrounding normal skin, which was not as apparent in wide field diffuse imaging. These results represent the first wide-field maps to quantify sub-diffuse scattering parameters, which are sensitive to sub-microscopic tissue structures and composition, and therefore, offer potential for fast diagnostic imaging of ultrastructure on a size scale that is relevant to surgical applications.

Structured light scatteroscopy.

A new imaging approach, structured light scatteroscopy (SLS), is demonstrated, which offers rapid wide-field imaging of microscopic morphological variations in bulk tissue surfaces. Elastic scattering of light offers exquisite sensitivity to ultrastructural changes at multiple size scales ranging from nanometers to millimeters, but in bulk tissues the confounding effects of molecular absorption and strong multiple scattering of light often lead to a dramatic reduction in scatter contrast and specificity. It is demonstrated that the SLS using structured high spatial frequency illumination and detection to probe the tissue achieves direct, absorption-independent, high-resolution maps of the scattering response. The scattering response is observed to be dependent on both the wavelength and spatial frequency of choice, indicating a potential for multiscale probing of ultrastructural changes in superficial tissue layers. This methodology can be easily applied in most wide-field imaging systems.

Anthropomorphic breast phantoms with physiological water, lipid, and hemoglobin content for near-infrared spectral tomography.

Breast mimicking tissue optical phantoms with sufficient structural integrity to be deployed as stand-alone imaging targets are developed and successfully constructed with biologically relevant concentrations of water, lipid, and blood. The results show excellent material homogeneity and reproducibility with inter- and intraphantom variability of 3.5 and 3.8%, respectively, for water and lipid concentrations ranging from 15 to 85%. The phantoms were long-lasting and exhibited water and lipid fractions that were consistent to within 5% of their original content when measured 2 weeks after creation. A breast-shaped three-compartment model of adipose, fibroglandular, and malignant tissues was created with water content ranging from 30% for the adipose simulant to 80% for the tumor. Mean measured water content ranged from 30% in simulated adipose to 73% in simulated tumor with the higher water localized to the tumor-like material. This novel heterogeneous phantom design is composed of physiologically relevant concentrations of the major optical absorbers in the breast in the near-infrared wavelengths that should significantly improve imaging system characterization and optimization because the materials have stand-alone structural integrity and can be readily molded into the sizes and shapes of tissues commensurate with clinical breast imaging.

Spectral discrimination of breast pathologies in situ using spatial frequency domain imaging.

INTRODUCTION: Nationally, 25% to 50% of patients undergoing lumpectomy for local management of breast cancer require a secondary excision because of the persistence of residual tumor. Intraoperative assessment of specimen margins by frozen-section analysis is not widely adopted in breast-conserving surgery. Here, a new approach to wide-field optical imaging of breast pathology in situ was tested to determine whether the system could accurately discriminate cancer from benign tissues before routine pathological processing. METHODS: Spatial frequency domain imaging (SFDI) was used to quantify near-infrared (NIR) optical parameters at the surface of 47 lumpectomy tissue specimens. Spatial frequency and wavelength-dependent reflectance spectra were parameterized with matched simulations of light transport. Spectral images were co-registered to histopathology in adjacent, stained sections of the tissue, cut in the geometry imaged in situ. A supervised classifier and feature-selection algorithm were implemented to automate discrimination of breast pathologies and to rank the contribution of each parameter to a diagnosis. RESULTS: Spectral parameters distinguished all pathology subtypes with 82% accuracy and benign (fibrocystic disease, fibroadenoma) from malignant (DCIS, invasive cancer, and partially treated invasive cancer after neoadjuvant chemotherapy) pathologies with 88% accuracy, high specificity (93%), and reasonable sensitivity (79%). Although spectral absorption and scattering features were essential components of the discriminant classifier, scattering exhibited lower variance and contributed most to tissue-type separation. The scattering slope was sensitive to stromal and epithelial distributions measured with quantitative immunohistochemistry. CONCLUSIONS: SFDI is a new quantitative imaging technique that renders a specific tissue-type diagnosis. Its combination of planar sampling and frequency-dependent depth sensing is clinically pragmatic and appropriate for breast surgical-margin assessment. This study is the first to apply SFDI to pathology discrimination in surgical breast tissues. It represents an important step toward imaging surgical specimens immediately ex vivo to reduce the high rate of secondary excisions associated with breast lumpectomy procedures.

Direct identification of breast cancer pathologies using blind separation of label-free localized reflectance measurements.

Breast tumors are blindly identified using Principal (PCA) and Independent Component Analysis (ICA) of localized reflectance measurements. No assumption of a particular theoretical model for the reflectance needs to be made, while the resulting features are proven to have discriminative power of breast pathologies. Normal, benign and malignant breast tissue types in lumpectomy specimens were imaged ex vivo and a surgeon-guided calibration of the system is proposed to overcome the limitations of the blind analysis. A simple, fast and linear classifier has been proposed where no training information is required for the diagnosis. A set of 29 breast tissue specimens have been diagnosed with a sensitivity of 96% and specificity of 95% when discriminating benign from malignant pathologies. The proposed hybrid combination PCA-ICA enhanced diagnostic discrimination, providing tumor probability maps, and intermediate PCA parameters reflected tissue optical properties.

System analysis of spatial frequency domain imaging for quantitative mapping of surgically resected breast tissues.

The feasibility of spatial frequency domain imaging (SFDI) for breast surgical margin assessment was evaluated in tissue-simulating phantoms and in fully intact lumpectomy specimens at the time of surgery. Phantom data was evaluated according to contrast-detail resolution, quantitative accuracy and model-data goodness of fit, where optical parameters were estimated by minimizing the residual sum of squares between the measured modulation amplitude and its solutions, modeled according to diffusion and scaled-Monte Carlo simulations. In contrast-detail phantoms, a 1.25-mm-diameter surface inclusion was detectable for scattering contrast >28%; a fraction of this scattering contrast (7%) was detectable for a 10 mm surface inclusion and at least 33% scattering contrast was detected up to 1.5 mm below the phantom surface, a probing depth relevant to breast surgical margin assessment. Recovered hemoglobin concentrations were insensitive to changes in scattering, except for overestimation at visible wavelengths for total hemoglobin concentrations <15 µM. The scattering amplitude increased linearly with scattering concentration, but the scattering slope depended on both the particle size and number density. Goodness of fit was comparable for the diffusion and scaled-Monte Carlo models of transport in spatially modulated, near-infrared reflectance acquired from 47 lumpectomy tissues, but recovered absorption parameters varied more linearly with expected hemoglobin concentration in liquid phantoms for the scaled-Monte Carlo forward model. SFDI could potentially reduce the high secondary excision rate associated with breast conserving surgery; its clinical translation further requires reduced image reconstruction time and smart inking strategies.

Scanning in situ spectroscopy platform for imaging surgical breast tissue specimens.

A non-contact localized spectroscopic imaging platform has been developed and optimized to scan 1 x 1 cm² square regions of surgically resected breast tissue specimens with ~150-micron resolution. A color corrected, image-space telecentric scanning design maintained a consistent sampling geometry and uniform spot size across the entire imaging field. Theoretical modeling in ZEMAX allowed estimation of the spot size, which is equal at both the center and extreme positions of the field with ~5% variation across the designed waveband, indicating excellent color correction. The spot sizes at the center and an extreme field position were also measured experimentally using the standard knife-edge technique and were found to be within ~8% of the theoretical predictions. Highly localized sampling offered inherent insensitivity to variations in background absorption allowing direct imaging of local scattering parameters, which was validated using a matrix of varying concentrations of Intralipid and blood in phantoms. Four representative, pathologically distinct lumpectomy tissue specimens were imaged, capturing natural variations in tissue scattering response within a given pathology. Variations as high as 60% were observed in the average reflectance and relative scattering power images, which must be taken into account for robust classification performance. Despite this variation, the preliminary data indicates discernible scatter power contrast between the benign vs malignant groups, but reliable discrimination of pathologies within these groups would require investigation into additional contrast mechanisms.

Characterization of materials for optimal near-infrared and x-ray imaging of the breast.

Development of a detector case for complete co-registration of images in a non-fiber-based combined near-infrared spectral tomography and digital breast tomosynthesis, required analysis to find materials that could support a breast under full mammographic compression without affecting the x-ray images or the quality of the near infrared measurements. Several possible solutions were considered, and many types of plastics were tested in the development of the detector case. Light channeling within the detector case changed the data obtained in resin and agarose phantoms, lowering recovered absorption values. Additional developments focusing on blocking stray light were successful and permitted a normal subject imaging exam.

A digital x-ray tomosynthesis coupled near infrared spectral tomography system for dual-modality breast imaging.

A Near Infrared Spectral Tomography (NIRST) system has been developed and integrated into a commercial Digital Breast Tomosynthesis (DBT) scanner to allow structural and functional imaging of breast in vivo. The NIRST instrument uses an 8-wavelength continuous wave (CW) laser-based scanning source assembly and a 75-element silicon photodiode solid-state detector panel to produce dense spectral and spatial projection data from which spectrally constrained 3D tomographic images of tissue chromophores are produced. Integration of the optical imaging system into the DBT scanner allows direct co-registration of the optical and DBT images, while also facilitating the synergistic use of x-ray contrast as anatomical priors in optical image reconstruction. Currently, the total scan time for a combined NIRST-DBT exam is ~50s with data collection from 8 wavelengths in the optical scan requiring ~42s to complete. The system was tested in breast simulating phantoms constructed using intralipid and blood in an agarose matrix with a 3 cm x 2 cm cylindrical inclusion at 1 cm depth from the surface. Diffuse image reconstruction of total hemoglobin (HbT) concentration resulted in accurate recovery of the lateral size and position of the inclusion to within 6% and 8%, respectively. Use of DBT structural priors in the NIRST reconstruction process improved the quantitative accuracy of the HbT recovery, and led to linear changes in imaged versus actual contrast, underscoring the advantages of dual-modality optical imaging approaches. The quantitative accuracy of the system can be further improved with independent measurements of scattering properties through integration of frequency or time domain data.

Scatter spectroscopic imaging distinguishes between breast pathologies in tissues relevant to surgical margin assessment.

PURPOSE: A new approach to spectroscopic imaging was developed to detect and discriminate microscopic pathologies in resected breast tissues; diagnostic performance of the prototype system was tested in 27 tissues procured during breast conservative surgery. EXPERIMENTAL DESIGN: A custom-built, scanning in situ spectroscopy platform sampled broadband reflectance from a 150-µm-diameter spot over a 1 × 1 cm(2) field using a dark field geometry and telecentric lens; the system was designed to balance sensitivity to cellular morphology and imaging the inherent diversity within tissue subtypes. Nearly 300,000 broadband spectra were parameterized using light scattering models and spatially dependent spectral signatures were interpreted using a cooccurrence matrix representation of image texture. RESULTS: Local scattering changes distinguished benign from malignant pathologies with 94% accuracy, 93% sensitivity, 95% specificity, and 93% positive and 95% negative predictive values using a threshold-based classifier. Texture and shape features were important to optimally discriminate benign from malignant tissues, including pixel-to-pixel correlation, contrast and homogeneity, and the shape features of fractal dimension and Euler number. Analysis of the region-based diagnostic performance showed that spectroscopic image features from 1 × 1 mm(2) areas were diagnostically discriminant and enabled quantification of within-class tissue heterogeneities. CONCLUSIONS: Localized scatter-imaging signatures detected by the scanning spectroscopy platform readily distinguished benign from malignant pathologies in surgical tissues and showed new spectral-spatial signatures of clinical breast pathologies.

Dark-field scanning in situ spectroscopy platform for broadband imaging of resected tissue.

A dark-field geometry spectral imaging system is presented to raster scan thick tissue samples in situ in 1.5 cm square sections, recovering full spectra from each 100 µm diameter pixel. This spot size provides adequate resolution for wide field scanning, while also facilitating scatter imaging without requiring sophisticated light-tissue transport modeling. The system is demonstrated showing accurate estimation of localized scatter parameters and the potential to recover absorption-based contrast from broadband reflectance data measured from 480 nm up to 750 nm in tissue phantoms. Results obtained from xenograft pancreas tumors show the ability to quantitatively image changes in localized scatter response in this fast-imaging geometry. The polychromatic raster scan design allows the rapid scanning necessary for use in surgical/clinical applications where timely decisions are required about tissue pathology.

Semi-automated segmentation and classification of digital breast tomosynthesis reconstructed images.

Digital breast tomosynthesis (DBT) is a limited-angle tomographic x-ray imaging technique that reduces the effect of tissue superposition observed in planar mammography. An integrated imaging platform that combines DBT with near infrared spectroscopy (NIRS) to provide co-registered anatomical and functional imaging is under development. Incorporation of anatomic priors can benefit NIRS reconstruction. In this work, we provide a segmentation and classification method to extract potential lesions, as well as adipose, fibroglandular, muscle and skin tissue in reconstructed DBT images that serve as anatomic priors during NIRS reconstruction. The method may also be adaptable for estimating tumor volume, breast glandular content, and for extracting lesion features for potential application to computer aided detection and diagnosis.