Skip to main content

Spectral discrimination of breast pathologies in situusing spatial frequency domain imaging



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.


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.


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.


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.


Optimal local management of breast cancer has been hindered by an inability to assess tumor-margin status intraoperatively, predominantly because frozen sections are limited by freezing artifacts in adipose tissue [1, 2], and sensitivities reported for touch-preparation cytology have been inconsistent [3]. Breast-conserving therapy (BCT), which includes local tumor excision followed by moderate-dose radiation therapy, is the standard of care for patients with early invasive breast cancer (stage I and II) and for patients with advanced disease whose tumor burden has been successfully reduced with neoadjuvant chemotherapy [4, 5]. It was the treatment of choice for nearly 75% of the approximately 300,000 new breast cancer patients diagnosed in 2011 [6]. Prospective, randomized trials have demonstrated that BCT survival rates are equivalent to those of mastectomy when surgical margins are negative for residual disease [4, 5], but positive margins have been associated with an increased risk of local recurrence [79] and mortality [10]. Consequently, reexcision of a margin is the standard of care if invasive tumor is detected at the surface of the resected specimen, and also occurs at Dartmouth-Hitchcock Medical Center (DHMC) when ductal carcinoma in situ (DCIS) is found within 2 mm of the inked tissue surface. Wider margins are examined for DCIS because it can be multifocal and interspersed within normal tissues [11]. At DHMC in 2010, invasive cancer was detected at the margin, and in situ cancer was discovered at or within 2 mm of the margin in 20% and 23% (respectively) of patients treated with BCT for invasive carcinomas (n = 129), yielding a 43% reexcision rate (36% actually reexcised). If residual tumor could be detected at one or more margins during the initial surgery, the surgeon could take directed, additional tissue before closing.

Neoplastic processes, from early dysplasia to advanced-stage infiltrating tumors, perturb tissue ultrastructure and thereby alter its optical-scattering spectrum [12, 13]. The scattering spectrum can distinguish tissue pathologies when the optical signal is sampled locally [1418] or filtered by using polarization techniques [1921] to minimize the collection of multiply scattered light. Localized methods, such as optical coherence tomography (OCT) [22], Raman spectroscopy [23], and confocal sampling [24, 25], have been applied to surgical-margin assessment, but they are fundamentally limited in depth sampling by scattering attenuation in tissue. In most realizations, the microscopic field of view (FOV) is too small to evaluate surgical specimens wholly, so they have the same margin-sampling limitations as does conventional histopathology. Resected tissues may include lesions up to 5 cm in diameter, surrounded by a targeted layer of grossly normal breast that can be as thick as 1 cm. Wide-field imaging with localized sampling has recently been realized through multiplexed arrays of probes [17, 18, 26, 27]; these approaches still rely on discrete sampling to form images and thereby incompletely assess disease extent, multifocality, and tumor heterogeneity. Raster-scanning techniques [1618, 28] support high-resolution sampling to assess tumor heterogeneity, but have limited field-of-view and speed to allow scanning the complete surgical specimen in surgical settings. Ideally, the complete surgical specimen would be evaluated (in a noncontact manner) without sacrificing sensitivity to tumor-specific features in the scattering spectrum. Planar spectral imaging techniques have only recently been tailored to surgical resection guidance [2931] and sentinel lymph node (SLN) mapping [30, 32, 33], largely because of the explosive development of molecularly specific NIR probes [34]. However, most methods rely on diffuse light transport, which can be insufficient to resolve important morphologic transformations that have dimensions comparable to the optical wavelength [35], because its spatial resolution is limited by light scattering in tissues [36].

SFDI, a planar-imaging modality pioneered by investigators at the University of California at Irvine and commercialized by Modulated Imaging Inc. for biologic imaging at spatial resolutions between coherent and diffuse optical-imaging techniques, was applied here for wide-field, tissue-type discrimination in nearly 50 surgical breast lesions. SFDI quantitatively resolves subsurface tissue absorption and scattering by analyzing the spatial-modulation transfer function (s-MTF) at multiple NIR wavelengths [37]. Planar, structured light patterns improve signal localization and enable selective depth sampling [38]. Recovered optical parameters are surface-weighted, which may have added value for surgical-margin assessment, where the goal is to detect malignant transformations in the outer millimeters of resected tissue. In this contribution, supervised learning and feature-selection algorithms were implemented to automate spectral discrimination of pathologies in intact, surgical tissues examined with SFDI and to optimize future development of spectroscopic tools for margin assessment.


Spatial frequency domain imaging (SFDI)

A compact SFDI system (purchased from Modulated Imaging Inc., Irvine, CA) illuminated breast surgical-resection specimens with a harmonically modulated, planar source at four NIR wavelengths (658, 730, 850, and 970 nm) [37]. Structured light patterns were projected onto each tissue surface at 30 spatial frequencies uniformly distributed between 0 mm-1 and 0.33 mm-1 by using high-power light-emitting diodes (LEDs), a projection system, and a digital micro-mirror device. The projector and camera subsystems were described in previous publications [37, 39] and were fully integrated in a portable platform mounted on a z-axis post. In total, 360 images were acquired per tissue (30 spatial frequencies × 3 phase offsets × 4 wavelengths) in approximately 10 minutes. Data were simultaneously acquired over the full field in a noncontact geometry, in which the acquisition field of view (FOV) was determined by magnification of the illumination and collection optics, here optimized to image a 5.5-inch × 7.5-inch area. A 12-bit CCD-based camera, co-registered with the projector, collected diffusely reflected light at a 696 × 520 pixel resolution with 2 × 2 pixel binning (full pixel resolution of the camera is 1,392 × 1,040). Fundamentally, however, the detected spatial resolution is limited by the physics of light transport in tissue, and contrast-detail analysis was previously performed to evaluate the minimized size of detectable scattering contrast in tissue-simulating phantoms [40]. A three-phase demodulation scheme was then used to extract the frequency-dependent modulation amplitude and to remove noise and ambient light from the imaged field [41]. The spatial frequency of the illumination pattern integrates optical parameterization with effective sampling depth, in which higher spatial frequencies attenuate more rapidly and thereby sample superficial tissue volumes. The effective sampling depth is also determined by the tissue optical properties. Theoretically, probing depths of 1 to 8 mm are expected in tissues; experimentally, we demonstrated detection of scattering inclusions up to depths of 1.5 mm in tissue-simulating phantoms by using the higher spatial frequencies used in this study [40]. At each imaging session, the modulation amplitude of a Siloxane titanium dioxide (TiO2) reflectance standard with known optical properties was imaged to calibrate spatial nonuniformity in the illumination and imaging systems. A 1-mm uniform height offset was erroneously, but systematically, introduced between measures of the tissue and reference standard. The offset did not affect tissue discrimination and resulted in negligible scaling of the modulation amplitude because reflectance decays with distance, h, according to an inverse square law (the distance between the CCD lens and tissue was ~120 mm h 2 2 / h 1 2 = 119 m m 2 / 120 m m 2 1 [42].

Optical parameters were estimated by minimizing the residual sum of squares between the measured modulation amplitude and its forward simulation, according to a scaled Monte Carlo model of light transport [37]. Monte Carlo solutions were stored in a look-up table (LUT) to reduce the computational burden associated with iterative estimation of spectral parameters [39]. Subsequently, spectral parameters were fit according to absorption by the endogenous tissue chromophores oxygenated-hemoglobin, deoxygenated-hemoglobin, and water, and the wavelength-dependence of light scattering, and produced estimates of scattering amplitude (A), scattering slope (b), total hemoglobin concentration (HbT), percentage hemoglobin oxygenation (%O 2 ), and percentage water (%H 2 O). Data-acquisition and analysis methods are discussed more comprehensively in a companion article [39].

Imaging protocol for the surgical specimens

The ability of SFDI to distinguish histopathology, here used as the diagnostic gold standard, was evaluated at the cut surfaces of 47 surgical breast tissues. In this HIPAA-compliant, prospective study, approved by the Institutional Review Board for the protection of human subjects at Dartmouth, written informed consent was not required for participants, although an Information Sheet detailing the study was provided with an opt-out provision. All of the study patients were female. A breakdown of the pathologic characteristics of the imaged cancers is provided in Table 1. In addition to these 27 biopsy-proven cancers, 11 fibroadenomas and nine normal (including fibrocystic disease) breast cases were imaged. At Dartmouth, the breast-resection specimens are inked with six pre-assigned colors to preserve orientation within the breast, so that, if it is found by pathologic evaluation that the margins are involved by tumor, the patient is called back to undergo a second resection of those specific margins rather than the entire surgical cavity.

Table 1 Summary of clinical and pathology data for study participants with a cancer diagnosis

The fresh, inked specimens, from patients who did not decline this research use of their tissues, were transported immediately to the Department of Pathology. By standard protocol, all of the specimens were “bread loafed” in their entirety into consecutive slices up to 0.5 cm in thickness, and, for this study, one face of one slice representing the lesion and its relation to the nearest margin was imaged. No surgical margin was imaged en face, given the required orientation inking by the surgeon. The exact area of tissue imaged in one slice was outlined with four inked pins pushed through the tissue block. The tissue was returned to Pathology after imaging, within 30 minutes of resection, for standard histologic processing. The imaged slices were formalin-fixed with the pins in situ before being removed for usual processing. The tiny inked pinholes survive processing and are evident in the histologic sections of that imaged tissue, facilitating optimal registration with the histologic findings. The effects of tissue shrinkage after processing and sectioning distortion were minimized by ensuring optimal formalin concentration and fixation time in accordance with College of American Pathologist (CAP) best practices. Specimen imaging did not affect procedure time in the operating room, processing time in the Department of Pathology, or the content and time to verification of the final pathology report.

The findings in histologic sections representing the imaged areas were included in the final pathology report. Figure 1 illustrates co-registration between photographs of the cut face of one slice in the specimen, with the lesions outlined in diagnosis-defined colors (row 1), corresponding histopathology (row 2), and spectral images (rows 3 to 7). An experienced breast pathologist (WAW) outlined tissue areas containing target breast lesions on digital photographs of each specimen, guided by co-registered hematoxylin and eosin (H&E)-stained tissue sections from the imaged surface. Affine transformations were manually performed in Photoshop to convert the outlined tissue areas into masks for region-based image analysis in MATLAB (The MathWorks Inc., Natick, MA, USA).

Figure 1

Representative spectral parameter maps for tissue subtypes. Spectral parameter maps corresponding to the pathology subtypes: normal (including fibrocystic disease) (red outline), fibroadenoma (blue outline), DCIS, invasive cancer and partially treated invasive cancer after neoadjuvant chemotherapy (all black outline), and fat (yellow outline or label). Row 1 is a tissue photograph of the cut face of one slice of the specimen with the lesion; row 2 is the corresponding histology; row 3 is the scattering-amplitude maps; row 4 is the scattering slope maps; row 5 is the hemoglobin concentration maps; row 6 is the percentage oxygenated hemoglobin maps; and row 7 is the percentage water maps.

In total, 48 specimens from 47 patients were consecutively imaged in this study; one was excluded from analysis because surgical inks bled into the primary tissue field. Surgical inks, sometimes observed at the tissue edge, were avoided during region of interest (ROI) selection and masked according to the R2 coefficient of determination (pixels excluded if R2 < 0.95). Spatial frequency- and wavelength-dependent reflectance measurements were evaluated on a pixel-by-pixel basis in 59 regions of interest (>265,000 spectral pixels). Discrimination was assessed between benign and malignant pathologies and between the benign pathology subtypes (normal and fibroadenoma), and the malignant pathology subtypes (DCIS, invasive cancer, and partially treated invasive cancer after neoadjuvant chemotherapy), as listed in Table 2. In this analysis, fibrocystic disease was grouped with normal breast pathologies. Tumors treated with neoadjuvant chemotherapy before surgical resection (n = 5) were imaged to observe how chemotherapy alters tissue optical properties, the better to inform diffuse tomographic monitoring of treatment response to therapy [43].

Table 2 Summary of tissue subtypes imaged

Histopathology and immunohistochemical correlates

Spectral maps were associated with morphologic and immunohistochemical markers identified in adjacent, stained sections of the tissue, cut in the exact geometry imaged in situ. Possible correlates between immunohistochemical measures and spectral parameters were evaluated according to the Pearson correlation coefficient. CD31 (platelet endothelial cell adhesion molecule) was used to stain immunohistochemically the preexisting endothelial cells as a general indicator of normal vasculature, even though malignant vessels have been reported to retain this antigen [44]. CD105 (Endoglin), thought to antagonize the inhibitory effects of transforming growth factor-beta (TGF-β) on cell proliferation and capillary formation [44], was used to stain immunohistochemically a transmembrane glycoprotein expressed predominantly on the surface of highly proliferating endothelial cells, as an indicator of tumor angiogenesis. Cytokeratin CK5D3 was used as a marker for breast epithelium (benign and malignant). All tissue was fixed in 10% buffered formalin (Biochemical Science Inc, Swedesboro, NJ, USA), dehydrated through graded alcohols, and paraffin embedded. Tissue sections (4 μm) were coated with adhesive (Sta-on; Surgipath Medical Industries, Richmond, IL, USA), mounted on glass slides, and stained with H&E for initial review. Tissue sections were air-dried for at least 30 minutes and then loaded onto a Leica Bond Max automated immunostainer. Here, the sections were baked (30 minutes at 60°C), dewaxed for 30 minutes at 72°C, rinsed with alcohol, and washed in Bond wash buffer. Antigen retrieval was accomplished with Bond epitope Retrieval 2, Ar9640 (pH 8.9 to 9.1) for 20 minutes at 100°C; followed by cooling for 12 minutes and a rinse in Bond wash buffer. All dewaxed, antigen-retrieval, and detection reagents were supplied by Leica Microsystems (Bannockburn, IL, USA). The primary antibodies were incubated for 15 minutes: CK5D3 at 1:100 (Biogenix, Fremont, CA, USA), CD31 at 1:50 (Dako, Carpinteria, CA, USA), and CD105 at 1:60 (Vector, Burlingame, CA, USA). Diaminobenzidine was applied for visualization with a hematoxylin counterstain.

Quantitation of immunomarkers

Using Media Cybernetics™ Image Pro (Media Cybernetics, Bethesda, MD, USA) and automated stage control bundled software (Surveyor© Automated Specimen Scanning stage control bundled software (Objective Imaging Ltd., Cambridge, UK), whole H&E-stained slides and CD105 and CD31-immunostained slides, representing the imaged areas of lesional breast tissue, were digitally scanned at high resolution and montaged. Mean vessel density (MVD) was quantified according to the areas of CD31-positive (preexisting vasculature) or CD105-positive (tumor-induced) vessels, and segmented in pseudo-color, as a percentage of the total slide area. Mean vessel area (MVA) was determined from the combined areas of CD31-positive or CD105-positive blood vessels, segmented in pseudo-color, and measured in square micrometers. The epithelial component of each imaged area was defined as the percentage of the total slide area showing positive immunostaining for CK5D3 and segmented in pseudo-color. On the same slide, the morphologically distinctive white spaces representing lipocytes (emptied of lipid as a result of fixation and tissue processing) were also segmented in pseudo-color as a percentage of the total slide. The remaining percentage of the slide represented connective tissue stroma. The diagnostic ROIs were selected by WAW, a board-certified surgical pathologist with 15 years of expertise in breast pathology and image analysis. The image thresholding and processing were performed by MCS, a technologist with 6 years of expertise in the Image Pro software.

Automated discrimination of breast pathologies

A nearest-neighbor learning algorithm was implemented to automate spectral discrimination of benign and malignant pathologies, and all pathology subtypes listed in Table 1; performance was evaluated by using a threefold cross-validation [45]. All data were randomly divided into three nonoverlapping sets, with an equal number of reflectance spectra per diagnostic class per set. Two of these sets were used for training, and the third set was used for validation. Training pixels were associated with a known diagnosis, according to the pathologist’s demarcation of lesions; the diagnosis of each validation point was also identified by the pathologist, but remained blind to the classifier to evaluate its performance. Classification error was taken to be the percentage of misclassified pixels in the validation set, where a misclassification means that the diagnosis assigned by the classifier did not match the diagnosis provided by the pathologist. Performance metrics were quantified 3 times for all possible permutations of the training and validation sets, and reported values were the average of these three executions. The nearest-neighbor algorithm assigned each unclassified parameter set to the majority diagnosis of its k nearest neighbors found in the training space by an efficient k-dimensional tree-search algorithm [46]. A whitening transformation was applied to all spectral parameters before diagnostic discrimination to prevent amplitude weighting. Additionally, outliers were removed from the training set according to their interquartile fractions (comprising 5% of the total data set). Query points were never marked as outliers, because this information was not known a priori.

Receiver operating characteristic (ROC) analysis was used to optimize classifier sensitivity and bias as a function of nearest-neighbor number, k. Confidence intervals (α = 0.05) for the sensitivity and specificity were computed according to the Yates χ2 interval [47].

Optical parameter ranking

An iterative search algorithm, sequential floating forward selection (SFFS), was implemented to rank the contribution of each spectral parameter to tissue-type discrimination. The Bhattacharyya statistical distance, J ij , was used to measure the separation between two diagnostic classes (i,j) [48, 49] and was generalized to all spectral parameters (n = 5) by:

J = i = 1 5 j = 1 5 P i P j J ij

Here, Pi represents the prior probability of class i determined by its fraction of pixels in the training set. The search algorithm ranked all parameters in order of the contribution to diagnostic-class separability.


Optical biomarkers of breast pathologies

Examples of optical maps co-registered with digital photographs of the surgically resected tissues are shown in Figure 1 for each lesion type evaluated. Boxplots in Figure 2 quantify the extensive distribution of optical and immunohistochemistry markers for all reflectance measures evaluated on a pixel-by-pixel basis (>265,000 measurements); mean values per patient are plotted as green dots that overlie the boxplots. Median values were compared to identify optical signatures specific to each diagnostic class, and interquartile fractions reflect the inherent heterogeneity observed within and between tissues. Median values were significantly different within 95% confidence limits for nonoverlapping notches. Distributions of the scattering amplitude and scattering slope showed less variance than functional absorption parameters, like hemoglobin concentration and oxygenation status, because the morphologic transformations that determine light scattering occur on submicron and even subnanometer length scales [50, 51]. Absorption parameters typically signal broader functions, and their distributions were characterized by greater variance because of surgery-induced artifacts like vascular disruption, eradication of the tissue’s oxygen supply, and eventual tissue dehydration. Efforts were made to image all tissues within the same time window and to avoid sampling of cauterization and hemorrhage, but limited artifacts are plausible. Surgical inks at the edge of the resected specimens may have also contributed to outliers with artificially high hemoglobin concentrations. This variability was also detected immunohistochemically in the subsequently fixed tissues, as shown in Figure 2f. However, SFDI assess the full tissue surface, accounting for inherent biologic heterogeneity, so that on average, functional biomarkers may strengthen tissue-type discrimination.

Figure 2

Diagnostic distributions of spectral parameters and immunohistochemistry correlates. (a-e) Boxplots of recovered spectral parameters for all tissues per pathology subtype, indicated by color: normal (including fibrocystic disease) (green), fibroadenoma (blue), ductal carcinoma in situ (pink), invasive cancer (red), and partially treated invasive cancer after neoadjuvant chemotherapy (orange). (f) Corresponding boxplots of immunohistochemical measures of percentage of stroma, epithelium, and fat, and CD31-positive, CD105-positive mean vascular density (MVD) and area (MVA). Box plot red bars indicate the population median; green dots indicate the mean value per patient ROI; and red crosses indicate outliers.

The distributions in Figure 2 show that observed hallmarks of normal or fibrocystic pathologies included a heightened scattering amplitude and slope, as compared with malignant pathologies. Lower scattering amplitude was observed in fibroadenomas, which contain clearly encapsulated and abundant stromal cellularity, confirming diagnostic trends established in previous studies of localized spectroscopic scattering from breast tissues [16]. Higher hemoglobin and reduced oxygenation values were observed in invasive cancers, as compared with the benign pathologies. DCIS was spectrally distinguished from other pathologies by decreased oxygenation and increased water levels, suggesting heightened metabolism and a proportional propensity for more surrounding fibrous stroma. However, a small DCIS patient population (n = 4) and spectral parameterization limited interpretation of this signature (fits were based on the extinction of pure, unbound water and did not include absorption by adipose tissue). Neoadjuvant chemotherapy increased the scattering slope when comparing its untreated, invasive counterparts, regardless of response to therapy. An increase in scattering slope suggests cell shrinkage, given its increased propensity for Rayleigh-type scattering.

Increased oxygenation levels were also induced in tumors partially responsive to chemotherapy and could be valuable indicators for therapy monitoring. Many aggressive tumors thrive in hypoxic environments [52], so increased oxygenation levels suggest a return to the normal phenotype.

Pairwise correlations between spectral parameters and quantitative immunohistochemistry are summarized in Table 3. Collagen in the extracellular stroma demonstrated a strong, positive correlation with scattering slope because the small fibrils that form collagen fibers further enhance Rayleigh-type scattering [53]. In contrast, a strong negative correlation was observed between tissue epithelial content and the scattering slope. Malignant tissues expressed greater levels of angiogenic vessels and preexisting vascular endothelium. This net growth in vascularity was detected spectroscopically by increased total hemoglobin, but no optical distinction was made between CD105-specific and CD31-specific vasculature. Immunohistochemical measures were sampled in discrete locations about the tissue surface, limiting its complete representation of ROIs.

Table 3 Spectral-immunohistochemical correlates

Pathology discrimination in resected breast tissues

Spatial frequency- and wavelength-dependent reflectance measures were applied to the discriminant classifier on a pixel-by-pixel basis, and separation of all tissue subtypes was achieved with 82% accuracy. Figure 3a shows that diagnostic sensitivity increased with nearest neighbor number at the cost of specificity or bias. ROC curves revealed that 11 and nine nearest neighbors optimized separation of benign and malignant pathologies and all pathology subtypes listed in Table 1, respectively. The confusion matrix in Figure 3b explicitly reflects the true and predicted diagnosis as a percentage of total number of diagnosed pixels. True positives per tissue type are along the diagonal, and misclassifications are off-axis; showing that, for example, 16% of treated invasive tumors were misclassified as normal. Figure 3c through e shows a false-color map of the multiparametric diagnosis co-registered with pathology for a malignant and benign tissue to illustrate application of the classifier to a patient diagnosis. Diagnostic performance is summarized in Table 4, which reports the sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and accuracy observed per tissue type. Discrimination of benign from malignant pathologies was highly specific (93%) and reasonably sensitive (79%), although sensitivity was limited in underrepresented pathologies like DCIS and residual invasive cancers after neoadjuvant chemotherapy (61% to 66%). Discrimination between benign and malignant pathologies was more accurate (88%) because these distribution sizes were more equally represented in the classifier training set. Ultimately, the surgeon can moderately regulate diagnostic sensitivity and specificity by varying the number of nearest neighbors used by the classifier to generate a diagnosis, as shown in Figure 3a. The NPV (89%) between benign and malignant pathologies was higher than the PPV (86%). In practical terms, surgeons are more interested in a high NPV when high confidence exists that no residual cancer has been left behind. A lower PPV, correlating with a false-positive call on any margin, would mean that the surgeon would take a little more tissue at the appropriate margin at the time of the first surgery. Even if that margin does not subsequently turn out to be positive, this can be justified clinically by the improved cosmetic outcome, the accuracy of immediate margin reexcision, the reduced costs of a single surgical procedure, and patient well-being by avoiding a second, separate reexcision procedure. High NPVs were also reported for underrepresented tissue types, but these values appear inflated when the number of true negatives significantly outweighed the number of true positives.

Figure 3

Optimization and performance of the nearest-neighbor classifier for diagnostic discrimination. (a) Nearest-neighbor number optimization by using the receiver operating characteristic curve for discrimination between all pathology subtypes. (b) A confusion matrix showing the true and predicted diagnosis for all spectroscopic measures according to the nearest-neighbor classifier, presented as a percentage of the total number of diagnosed pixels (N). A multiparametric diagnostic map for a malignant (row 1) and benign (row 2) tissue is illustrated (c through e); the patient histology is shown in column (c), a photograph of the tissue imaged by SFDI is shown in column (d), and the patient-specific diagnostic map generated by the classifier is shown in column (e).

Table 4 Summary of diagnostic performance

Diagnostic value of optical parameters

The SFFS algorithm ranked the region-averaged (1) scattering amplitude, (2) percentage water, (3) total hemoglobin concentration, (4) the scattering slope, and (5) percentage oxygen, as most significant to tissue discrimination in that order. The scattering amplitude was most valuable to pathology discrimination, even through qualitatively, the scattering slope appeared to localize better with suggestive fibroglandular lesions. The scattering slope was potentially undervalued by the limitations of pathology co-registration. Percentage oxygenated hemoglobin was least valuable to a diagnosis, mainly because of ex vivo fluctuations in hemoglobin oxygenation status. Bydlon [54] showed that percentage of oxygenated hemoglobin varies nonlinearly within 30 minutes of resection, but other spectral parameters, like the scattering coefficients and total hemoglobin concentration, remained relatively stable after excision.


To achieve the best survival outcomes for patients with breast cancer, the goal for BCT is to completely resect the tumor with negative margins and simultaneously take a minimum of tissue to preserve the shape of the breast [4, 5]. At Dartmouth and elsewhere, inking the specimen with multiple colors to document better its orientation within the breast ensures that if a margin were found to be involved by tumor pathologic examination, a more directed approach to additional tissue removal at a second resection procedure is associated with lower recurrence rates and better cosmesis [55]. This reexcision would be better performed at the time of the first surgery rather than waiting at least 48 hours, as is currently the case. Although a preliminary pathological analysis can be performed while the surgeon waits in the surgical suite, these current methods are far less accurate [13] and not amenable to surgical workflow or operating room efficiency. As a result, the initial surgery is concluded, and patients recover before the final pathological margin assessment. Effective and efficient patient care demands an intraoperative technique that minimizes the need for reexcisions and maximizes the cosmetic outcomes by assisting the surgeon in removing the minimum normal tissue around the tumor necessary to ensure a negative margin within a single procedure.

NIR optical parameters imaged by a new planar imaging modality, SFDI, were combined with a discriminant classifier to detect and differentiate breast pathologies in situ. Spectroscopic scattering was emphasized because absorption biomarkers varied markedly with the surgical procedure itself, and a feature-ranking algorithm identified the scattering amplitude as most significant to tissue-type discrimination. All images were interpreted according to histopathology for direct clinical relevance. Spectral parameters applied to a discriminant classifier distinguished benign from malignant pathologies with 88% accuracy and the pathology subtypes: normal (includes fibrocystic disease), fibroadenoma, DCIS, invasive cancer, and partially treated invasive cancer after neoadjuvant chemotherapy), were identified with 82% accuracy. Discrimination of benign from malignant pathologies was highly specific (93%) and reasonably sensitive (79%). The precision of spectral-image co-registration with pathology may have also limited the accuracy of discrimination; ROIs were chosen conservatively here because this study was the first to explore SFDI optical-parameter contrast in resected breast tissues. Microscopic co-registration may not be necessary for surgical-margin assessment, in which the ultimate clinical goal is to determine if residual cancer is present at the surface of a resected tissue. The primary advantage of SFDI is near-real-time assessment of intact, surgical specimens, so that decisions about additional margin sampling can be made at the time of primary surgery without damaging the sample. Clinical workflow and the sample remain intact, so that margin status may be confirmed postoperatively by conventional histopathology. Although its modest PPV (0.86 when discriminating between benign and malignant pathologies) could result in a wider primary excision margin, information is gained only at the time of surgery, and a secondary excision would still be avoided. More important, SFDI identifies breast pathologies with a high NPV (0.89 when discriminating between benign and malignant pathologies and 0.91 to 0.99 when discriminating between all pathology subtypes), ensuring that tissues left unresected are truly negative for residual disease. Ultimately, the surgeon has control over the tradeoff between sensitivity and specificity according to the number of nearest neighbors used by the classifier, as shown in the ROC curve in Figure 3. Therefore, decisions about the width of excisional margin could be made in a patient-specific manner.

Application of SFDI to breast surgical margin assessment may be realized soon, given its development as a commercial package and the discriminatory power demonstrated here. Its planar illumination scheme rapidly imaged 48 surgical tissues, and only one was rejected from analysis because ink (applied postoperatively) contaminated the primary imaged field. Image interpretation was highly efficient when using an LUT for parameter optimization and the k-dimensional search tree in the discriminant classifier. Further reductions in data-acquisition time are possible by limiting the number of spatial frequencies sampled; Cuccia et al.[37] demonstrated optical parameter recovery in tissue simulating phantoms with just two spatial frequencies. Recent hardware modifications by Modulated Imaging, Inc., enabled acquisition of 15 spatial frequencies at four wavelengths in less than 2 minutes; this has been tested in four fully intact lumpectomy specimens before gross sectioning by pathology [39].

Inking strategies coordinated with the surgeon will be necessary to translate SFDI further to the intraoperative setting. Although inking is necessary for conventional histopathology and to validate new diagnostic techniques, it is not required for intraoperative margin assessment. Ideally, the surgeon would apply the colored inks subsequent to imaging, but before removal of tissue from the imaging platform. The surgeon would be responsible for moving the specimen to and from the imaging platform and additionally could use sutures to safeguard knowledge of its orientation until inks are applied. In that way, feedback on margin status is provided at the time of primary surgery, and the specimen remains intact for conventional histopathology, maintaining the clinical gold standard. In this study, bread-loafed sections of the resected tissue were imaged to enable accurate co-registration with pathology for quantitative assessment of SFDI diagnostic performance, as compared with the clinical gold standard. SFDI images a surface area much larger than the microscopic FOV sampled by a pathologist, so it does not have the sampling limitations encountered in conventional histopathology. An initial demonstration of imaging fully intact surgical margins en face, before inking, is reported in our companion systems-analysis article [40]. Edge artifacts caused by surface profile changes were sometimes observed in intact specimens, but a three-phase amplitude demodulation scheme has been implemented to correct for surface profile changes [42]. Future development of SFDI for surgical-margin assessment should focus on enhancing tumor-specific scattering contrast, improving methods for spectral image co-registration with pathology, and automating image processing for near-real-time diagnostic feedback at the time of surgery. Increasing the density of visible-NIR sampling may also improve depth resolution and the quantification of significant absorption signatures like water.


SFDI combines quantitative spectroscopy with depth sampling appropriate for margin assessment and optimally balances the tradeoffs between wide-field acquisition and signal localization. Most important, it performs diagnostic assessment rapidly and on intact surgical specimens, providing intraoperative feedback to the surgeon. Submillimeter probing volumes limit microscopic resolution, but are here demonstrated to be diagnostically powerful and clinically pragmatic. Planar imaging of NIR optical parameters, in contrast to probing discrete tissue regions by using fiberoptics, comprehensively images known biologic heterogeneity within pathology subtypes. In combination with a discriminant classifier, this wealthy data set readily provides a tissue-type diagnosis. SFDI distinguished benign from malignant pathologies within 47 surgically resected breast tissues with 88% accuracy and discriminated between all pathology subtypes with 82% accuracy. Discrimination of benign from malignant pathologies was highly specific (93%) and reasonably sensitive (79%). Although spectral absorption and scattering features were essential components of the classifier, scattering exhibited lower variance, and the scattering amplitude contributed most to tissue-type separation. The scattering slope was also exquisitely sensitive to stromal and epithelial distributions measured with quantitative immunohistochemistry. Certain clinical challenges, mainly orienting spectroscopic images with pathology and the patient, and integrating data acquisition with immediate processing at the time of surgery, remain to be effected, but here we demonstrated that SFDI can rapidly discriminate between microscopic pathologies at the surface of lumpectomy tissues. This work validates the continued translation of SFDI toward evaluation of surgical specimens immediately ex vivo and at the time of primary surgery to reduce significantly the secondary excision rate currently associated with breast lumpectomy procedures.



Percentage water


Percentage hemoglobin oxygen


Scattering amplitude


Scattering slope


Breast-conserving therapy/surgery


College of American Pathologists


Charge-coupled device




Platelet endothelial cell adhesion molecule


Ductal carcinoma in situ


Dartmouth-Hitchcock Medical Center


Estrogen-receptor protein




Fluorescence in situ hybridization


Field of view


Hematoxylin and eosin


Total concentration of hemoglobin


Health Insurance Portability and Accountability Act




Infiltrating lobular carcinoma


Intermediate grade


Invasive cancer


Treated invasive cancer (residual)


Nearest neighbor number


Light-emitting diodes


Look-up table


Mean vessel area


Mean vessel density








Negative predictive value


Optical coherence tomography




Positive predictive value


Progesterone-receptor protein




Receiver operator characteristic


Region of interest


Spatial frequency domain imaging


Sequential floating-forward selection


Sentinel lymph node


Spatial-modulation transfer function


Transforming growth factor


Titanium dioxide


Wendy A. Wells.


  1. 1.

    Ferreiro J, Gisvold J, Bostwick D: Accuracy of frozen section diagnosis of mammographically detected breast biopsies: results of 1,490 consecutive cases. Am J Surg Pathol. 1995, 19: 1267-1271. 10.1097/00000478-199511000-00006.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Tinnemans J, Wobbes T, Holland R: Mammographic and histopathologic correlation of non-palpable lesions of the breast and reliability of frozen section diagnosis. Surg Gynecol Obstet. 1987, 165: 523-529.

    CAS  PubMed  Google Scholar 

  3. 3.

    Saarela A, Paloneva T, Rissanen T, Kiminiemi H: Determinants of positive histologic margins and residual tumor after lumpectomy for early breast cancer: a prospective study with special reference to touch preparation cytology. J Surg Oncol. 1997, 66: 248-253. 10.1002/(SICI)1096-9098(199712)66:4<248::AID-JSO5>3.0.CO;2-B.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Veronesi U, Cascinelli N, Mariani L, Greco M, Saccozzi R, Luini A, Aguilar M, Marubini E: Twenty-year follow-up of a randomized study comparing breast-conserving surgery with radical mastectomy for early breast cancer. N Engl J Med. 2002, 347: 1227-1232. 10.1056/NEJMoa020989.

    Article  PubMed  Google Scholar 

  5. 5.

    Fisher B, Anderson S, Bryant J, Margolese RG, Deutsch M, Fisher ER, Jeong J, Wolmark N: Twenty-year follow-up of a randomized trial comparing total mastectomy, lumpectomy, and lumpectomy plus irradiation for the treatment of invasive breast cancer. N Engl J Med. 2002, 347: 1233-1241. 10.1056/NEJMoa022152.

    Article  PubMed  Google Scholar 

  6. 6.

    American Cancer Society: Cancer Facts & Figures 2011. 2011, Atlanta: American Cancer Society

    Google Scholar 

  7. 7.

    Scopa CD, Aroukatos P, Tsamandas AC, Aletra C: Evaluation of margin status in lumpectomy specimens and residual breast carcinoma. Breast J. 2006, 12: 150-153. 10.1111/j.1075-122X.2006.00223.x.

    Article  PubMed  Google Scholar 

  8. 8.

    Schnitt SJ, Abner A, Gelman R, Connolly JL, Recht A, Duda RB, Eberlein TJ, Mayzel K, Silver B, Harris JR: The relationship between microscopic margins of resection and the risk of local recurrence in patients with breast cancer treated with breast-conserving surgery and radiation therapy. Cancer. 1994, 74: 1746-1751. 10.1002/1097-0142(19940915)74:6<1746::AID-CNCR2820740617>3.0.CO;2-Y.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Spivack B, Khanna MM, Tafra L, Juillard G, Giuliano AE: Margin status and local recurrence after breast-conserving surgery. Arch Surg Chicago. 1994, 129: 952-956. 10.1001/archsurg.1994.01420330066013.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Clarke M, Collins R, Darby S, Davies C, Elphinstone P, Evans E, Godwin J, Gray R, Hicks C, James S, MacKinnon E, McGale P, McHugh T, Peto R, Taylor C, Wang Y, Early Breast Cancer Trialists' Collaborative Group (EBCTCG): , et al: Effects of radiotherapy and of differences in the extent of surgery for early breast cancer on local recurrence and 15-year survival: an overview of the randomised trials. Lancet. 2005, 366: 2087-2106.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Singletary S: Surgical margins in patients with early-stage breast cancer treated with breast conservation therapy. Am J Surg. 2002, 184: 383-393. 10.1016/S0002-9610(02)01012-7.

    Article  PubMed  Google Scholar 

  12. 12.

    Michels R, Foschum F, Kienle A: Optical properties of fat emulsions. Optics Exp. 2008, 16: 5907-5925. 10.1364/OE.16.005907.

    CAS  Article  Google Scholar 

  13. 13.

    Mourant JR, Boyer J, Hielscher AH, Bigio IJ: Influence of the scattering phase function on light transport measurements in turbid media performed with small source-detector separations. Optics Lett. 1996, 21: 546-548. 10.1364/OL.21.000546.

    CAS  Article  Google Scholar 

  14. 14.

    van Veen RLP, Amelink A, Menke-Pluymers M, van der Pol C, Sterenborg H: Optical biopsy of breast tissue using differential path-length spectroscopy. Physics Med Biol. 2005, 50: 2573-2581. 10.1088/0031-9155/50/11/009.

    Article  Google Scholar 

  15. 15.

    Canpolat M, Mourant JR: Particle size analysis of turbid media with a single optical fiber in contact with the medium to deliver and detect white light. Appl Opt. 2001, 40: 3792-3799. 10.1364/AO.40.003792.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Laughney AM, Krishnaswamy V, Rizzo EJ, Schwab MC, Barth RJ, Pogue BW, Paulsen KD, Wells WA: Scatter spectroscopic imaging distinguishes between breast pathologies in tissues relevant to surgical margin assessment. Clin Cancer Res. 2012, 18: 10-

    Article  Google Scholar 

  17. 17.

    Krishnaswamy V, Laughney AM, Wells WA, Paulsen KD, Pogue BW: Scanning in situ spectroscopy platform for imaging surgical breast tissues. Opt Express. 2013, 21: 2185-2194. 10.1364/OE.21.002185.

    Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Krishnaswamy V, Hoopes PJ, Samkoe KS, O’Hara JA, Hasan T, Pogue BW: Quantitative imaging of scattering changes associated with epithelial proliferation, necrosis, and fibrosis in tumors using microsampling reflectance spectroscopy. J Biomed Optics. 2009, 14: 014004-10.1117/1.3065540.

    Article  Google Scholar 

  19. 19.

    Perelman LT, Backman V, Wallace M, Zonios G, Manoharan R, Nusrat A, Shields S, Seiler M, Lima C, Hamano T, et al: Observation of periodic fine structure in reflectance from biological tissue: a new technique for measuring nuclear size distribution. Phys Rev Lett. 1998, 80: 627-630. 10.1103/PhysRevLett.80.627.

    CAS  Article  Google Scholar 

  20. 20.

    Gurjar RS, Backman V, Perelman LT, Georgakoudi I, Badizadegan K, Itzkan I, Dasari RR, Feld MS: Imaging human epithelial properties with polarized light-scattering spectroscopy. Nat Med. 2001, 7: 1245-1248. 10.1038/nm1101-1245.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Kim YL, Liu Y, Turzhitsky VM, Roy HK, Wali RK, Backman V: Coherent backscattering spectroscopy. Optics Lett. 2004, 29: 1906-1908. 10.1364/OL.29.001906.

    Article  Google Scholar 

  22. 22.

    Nguyen F, Zysk A, Chaney E, Kotynek J, Oliphant U, Bellafiore F, Rowland K, Johnson P, Boppart S: Intraoperative evaluation of breast tumor margins with optical coherence tomography. Cancer Res. 2009, 69: 8790-8796. 10.1158/0008-5472.CAN-08-4340.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Haka AS, Volynskaya Z, Gardecki JA, Nazemi J, Lyons J, Hicks D, Fitzmaurice M, Dasari RR, Crowe JP, Feld MS: In vivo margin assessment during partial mastectomy breast surgery using Raman spectroscopy. Cancer Res. 2006, 66: 3317-3322. 10.1158/0008-5472.CAN-05-2815.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Brown J, Wilke L, Geradts J, Kennedy S, Palmer G, Ramanujam N: Quantitative optical spectroscopy: a robust tool for direct measurement of breast cancer vascular oxygenation and total hemoglobin content in vivo. Cancer Res. 2009, 69: 2919-2926. 10.1158/0008-5472.CAN-08-3370.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Bigio I, Bown S, Briggs G, Kelley C, Lakhani S, Pickard D, Ripley P, Rose I, Saunders C: Diagnosis of breast cancer using elastic-scattering spectroscopy: preliminary clinical results. J Biomed Opt. 2000, 5: 221-228. 10.1117/1.429990.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Brown JQ, Bydlon TM, Kennedy SA, Richards L, Junker MS, Palmer GM, Geradts J, Wilke LG, Ramanujam N: Intraoperative optical breast tissue characterization device for tumor margin assessment. Cancer Res. 2009, 69: 101s-101s.

    Google Scholar 

  27. 27.

    Bydlon T, Kennedy S, Richards L, Brown J, Yu B, Junker M, Gallagher J, Geradts J, Wilke L, Ramanujam N: Performance metrics of an optical spectral imaging system for intra-operative assessment of breast tumor margins. Opt Express. 2010, 18: 8058-8076. 10.1364/OE.18.008058.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Yu C-C, Lau C, O’Donoghue G, Mirkovic J, McGee S, Galindo L, Elackattu A, Stier E, Grillone G, Badizadegan K, et al: Quantitative spectroscopic imaging for non-invasive early cancer detection. Optics Exp. 2008, 16: 16227-16239. 10.1364/OE.16.016227.

    CAS  Article  Google Scholar 

  29. 29.

    Bogaards A, Sterenborg HJCM, Trachtenberg J, Wilson BC, Lilge L: In vivo quantification of fluorescent molecular markers in real time by ratio imaging for diagnostic screening and image-guided surgery. Lasers Surg Med. 2007, 39: 605-613. 10.1002/lsm.20525.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Gioux S, Choi HS, Frangioni JV: Image-guided surgery using invisible near-infrared light: fundamentals of clinical translation. Mol Imaging. 2010, 9: 237-255.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Themelis G, Yoo JS, Soh K-S, Schulz R, Ntziachristos V: Real-time intraoperative fluorescence imaging system using light-absorption correction. J Biomed Optics. 2009, 14: 064012-10.1117/1.3259362.

    Article  Google Scholar 

  32. 32.

    Troyan SL, Kianzad V, Gibbs-Strauss SL, Gioux S, Matsui A, Oketokoun R, Ngo L, Khamene A, Azar F, Frangioni JV: The FLARE((TM)) intraoperative near-infrared fluorescence imaging system: a first-in-human clinical trial in breast cancer sentinel lymph node mapping. Ann Surg Oncol. 2009, 16: 2943-2952. 10.1245/s10434-009-0594-2.

    Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Xu H, Rice BW: In-vivo fluorescence imaging with a multivariate curve resolution spectral unmixing technique. J Biomed Optics. 2009, 14: 064011-10.1117/1.3258838.

    Article  Google Scholar 

  34. 34.

    Ntziachristos V, Bremer C, Weissleder R: Fluorescence imaging with near-infrared light: new technological advances that enable in vivo molecular imaging. Eur Radiol. 2003, 13: 195-208.

    PubMed  Google Scholar 

  35. 35.

    Mourant JR, Freyer JP, Hielscher AH, Eick AA, Shen D, Johnson TM: Mechanisms of light scattering from biological cells relevant to noninvasive optical-tissue diagnostics. Appl Optics. 1998, 37: 3586-3593. 10.1364/AO.37.003586.

    CAS  Article  Google Scholar 

  36. 36.

    Svaasand LO, Spott T, Fishkin JB, Pham T, Tromberg BJ, Berns MW: Reflectance measurements of layered media with diffuse photon-density waves: a potential tool for evaluating deep burns and subcutaneous lesions. Physics Med Biol. 1999, 44: 801-813. 10.1088/0031-9155/44/3/020.

    CAS  Article  Google Scholar 

  37. 37.

    Cuccia DJ, Bevilacqua F, Durkin AJ, Ayers FR, Tromberg BJ: Quantitation and mapping of tissue optical properties using modulated imaging. J Biomed Optics. 2009, 14: 024012-10.1117/1.3088140.

    Article  Google Scholar 

  38. 38.

    Pham TH, Spott T, Svaasand LO, Tromberg BJ: Quantifying the properties of two-layer turbid media with frequency-domain diffuse reflectance. Appl Optics. 2000, 39: 4733-4745. 10.1364/AO.39.004733.

    CAS  Article  Google Scholar 

  39. 39.

    Mayo JG, Laster WR, Andrews CM, Schabel FM: Success and failure in the treatment of solid tumors: 3. “Cure” of metastatic Lewis lung carcinoma with methyl-CCNU (NSC-95442) and surgery-chemotherapy. Cancer Chemother Rep. 1972, 56: 183-195.

    CAS  PubMed  Google Scholar 

  40. 40.

    Laughney AM, Krishnaswamy V, Rice TB, Cuccia DJ, Barth RJ, Tromberg BJ, Paulsen KD, Pogue BW, Wells WA: System analysis of spatial frequency domain imaging for quantitative mapping of surgically resected breast tissues. J Biomed Optics. 2013, 18: 036012-10.1117/1.JBO.18.3.036012.

    Article  Google Scholar 

  41. 41.

    Neil MAA, Juskaitis R, Wilson T: Method of obtaining optical sectioning by using structured light in a conventional microscope. Optics Lett. 1997, 22: 1905-1907. 10.1364/OL.22.001905.

    CAS  Article  Google Scholar 

  42. 42.

    Gioux S, Mazhar A, Cuccia DJ, Durkin AJ, Tromberg BJ, Frangioni JV: Three-dimensional surface profile intensity correction for spatially modulated imaging. J Biomed Optics. 2009, 14: 034045-10.1117/1.3156840.

    Article  Google Scholar 

  43. 43.

    Kukretia S, Cerussi A, Tromberg B, Gratton E: Intrinsic near-infrared spectroscopic markers of breast tumors. Dis Markers. 2008, 25: 281-290. 10.1155/2008/651308.

    Article  Google Scholar 

  44. 44.

    Pakalniskis MG, Wells WA, Schwab MC, Froehlich HM, Jiang S, Li Z, Tosteson TD, Poplack SP, Kaufman PA, Pogue BW, et al: Tumor angiogenesis change estimated by using diffuse optical spectroscopic tomography: demonstrated correlation in women undergoing neoadjuvant chemotherapy for invasive breast cancer?. Radiology. 2011, 259: 365-374. 10.1148/radiol.11100699.

    Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Laughney AM, Krishnaswamy V, Garcia-Allende PB, Conde OM, Wells WA, Paulsen KD, Pogue BW: Automated classification of breast pathology using local measures of broadband reflectance. J Biomed Opt. 2010, 15: 066019-10.1117/1.3516594.

    Article  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Matlab Statistical Toolbox. 2012, Natick MA, USA,: Mathworks

  47. 47.

    Wallis S: Binomial distributions, probability and Wilson’s confidence interval. 2009, University College London,,

    Google Scholar 

  48. 48.

    Garcia-Allende PB, Krishnaswamy V, Hoopes PJ, Samkoe KS, Conde OM, Pogue BW: Automated identification of tumor microscopic morphology based on macroscopically measured scatter signatures. J Biomed Optics. 2009, 14: 034034-10.1117/1.3155512.

    Article  Google Scholar 

  49. 49.

    Gomez-Chova L, Calpe J, Camps-Valls G, Martín JD, Soria E, Vila J, Alonso-Chorda L, Moreno J: Feature selection of hyperspectral data through local correlation and SFFS for crop classification. IEEE Int Geosci Remote Sensing Sympos. 2003, 1: 555-557.

    Google Scholar 

  50. 50.

    Backman V, Gopal V, Kalashnikov M, Badizadegan K, Gurjar R, Wax A, Georgakoudi I, Mueller M, Boone CW, Dasari RR, et al: Measuring cellular structure at submicrometer scale with light scattering spectroscopy. IEEE J Select Topics Quantum Elect. 2001, 7: 887-893. 10.1109/2944.983289.

    CAS  Article  Google Scholar 

  51. 51.

    Subramanian H, Roy HK, Pradhan P, Goldberg MJ, Muldoon J, Brand RE, Sturgis C, Hensing T, Ray D, Bogojevic A, et al: Nanoscale cellular changes in field carcinogenesis detected by partial wave spectroscopy. Cancer Res. 2009, 69: 5357-5363. 10.1158/0008-5472.CAN-08-3895.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Groenhuis RAJ, Ferwerda HA, Tenbosch JJ: Scattering and absorption of turbid materials determined from reflection measurements. 1. Theory. Appl Optics. 1983, 22: 2456-2462. 10.1364/AO.22.002456.

    CAS  Article  Google Scholar 

  53. 53.

    Saidi IS, Jacques SL, Tittel FK: Mie and Rayleigh modeling of visible-light scattering in neonatal skin. Appl Optics. 1995, 34: 7410-7418. 10.1364/AO.34.007410.

    CAS  Article  Google Scholar 

  54. 54.

    Bydlon TM, Barry WT, Kennedy SA, Brown JQ, Gallagher JE, Wilke LG, Geradts J, Ramanujam N: Advancing optical imaging for breast margin assessment: an analysis of excisional time, cautery, and patent blue dye on underlying sources of contrast. PLoS One. 2012, 7: e51418-10.1371/journal.pone.0051418.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Gibson G, Lesnikoski BE, Yoo J, Mott LA, Cady B, Barth RJ: A comparison of ink-directed and traditional whole cavity re-excision for breast lumpectomy specimens with positive margins. Ann Surg Oncol. 2001, 9: 693-704.

    Article  Google Scholar 

Download references


The authors thank Kari Rozenkranz, MD, and Burt Eisenberg, MD, in the Department of Surgery, and Vincent Memoli, MD, Candice Black, DO, Xiaoying Liu, MD, and Laura Tafe, MD, in the Department of Pathology at Dartmouth Hitchcock Medical Center for their help informing patients and procuring and diagnosing the breast surgical specimens. The immunohistochemistry was performed by Rebecca O’Meara, HT (ASCP), in the Department of Pathology Translational Research Shared Resource Laboratory of the Geisel School of Medicine at Dartmouth, the Dartmouth Hitchcock Medical Center, and the Norris Cotton Cancer Center. This work was supported by the NCRR 1R21RR024411-01A1, the NIH PO1 CA80139, and NIH P41EB015890 (Laser Microbeam and Medical Program, LAMMP), and the Department of Defense BC093811.

Author information



Corresponding author

Correspondence to Wendy A Wells.

Additional information

Competing interests

AML, VK, EJR, MCS, RJB, KDP, BWP, and WAW declare no competing interests. BT and DC report patents, owned by the University of California, that are related to the SFDI technology and analysis methods. The University of California has licensed SFDI technologies to two companies, Caliper/Perkin-Elmer, Inc., and Modulated Imaging, Inc. BT and DC are cofounders of Modulated Imaging. This research was completed without financial support of either company.

Authors’ contributions

BT and DC provided the original SFDI technology and analysis methods. AML and VK designed, adapted, and analyzed the performance of the new imaging technology. AML evaluated all images and quantified their diagnostic significance. EJR procured the breast tissue and ensured that the clinical standard-of-care for breast-specimen processing and reporting was maintained at all times. RJB talked to the patients about the study, performed the surgery, and ensured that the specimens were transported from the operating room in a timely manner. AML, WAW, EJR, BWP, and KDP participated in the design of the study and drafted the manuscript. MCS performed the quantitative image analysis of the immunohistochemical correlates. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Laughney, A.M., Krishnaswamy, V., Rizzo, E.J. et al. Spectral discrimination of breast pathologies in situusing spatial frequency domain imaging. Breast Cancer Res 15, R61 (2013).

Download citation


  • Breast-conserving surgery/therapy
  • Near-infrared spectroscopy
  • Spatial frequency domain imaging
  • Diagnostic pathology