PurposeWe investigate terahertz (THz) polarimetry imaging of seven human breast cancer surgical specimens. The goal is to enhance image contrast between adjacent tissue types of cancer, healthy collagen, and fat in excised breast tumors. Based on the biological perception of random growth of cancer and invasion of surrounding healthy tissues in the breast, we hypothesize that cancerous cells interact with the THz electric field in a different manner compared with healthy cells. This difference can be best captured using multiple polarizations instead of single polarization.ApproachTime domain pulsed signals are experimentally collected from each pixel of the specimen in horizontal–horizontal, vertical–horizontal, vertical–vertical, and horizontal–vertical polarizations. The time domain pulses are transformed to the frequency domain to obtain the power spectra and 16 Mueller matrix images. The whole-slide pathology imaging was used to interpret and label all images.ResultsThe results of the cross and co-polarization power spectrum images demonstrated a strong dependency on the tissue orientation with respect to the emitted and detected electric fields. At the 130-deg rotation angle of the scanned samples, the detector showed the strongest reflected signal in cross-polarization. Furthermore, the Mueller matrix images consistently demonstrated patterns in fresh and block tissues confirming the differentiation between tissue types in breast tumor specimens.ConclusionsTHz polarimetry imaging shows a potential for improving image contrast in excised tumor tissues compared with single polarization imaging. Cross-polarization signals demonstrated smaller amplitudes compared with co-polarized signals. However, averaging the signal during measurements has tremendously improved the image. Furthermore, in post-processing, averaging the frequency domain images and the Mueller matrix elements with respect to frequency has led to better image contrast. Some patterns in the Mueller matrix images were difficult to interpret leading to the necessity of more investigation of the Mueller matrix and its physiological interpretation of breast tumor tissues.
Purpose: We investigate the enhancement in terahertz (THz) images of freshly excised breast tumors upon treatment with an optical clearance agent. The hyperspectral imaging and spectral classifications are used to quantitatively demonstrate the image enhancement. Glycerol solution with 60% concentration is applied to excised breast tumor specimens for various time durations to investigate the effectiveness on image enhancement.
Approach: THz reflection spectroscopy is utilized to obtain the absorption coefficient and the index of refraction of untreated and glycerol-treated tissues at each frequency up to 3 THz. Two classifiers, spectral angular mapping (SAM) based on several kernels and Euclidean minimum distance (EMD) are implemented to evaluate the effectiveness of the treatment. The testing raw data is obtained from five breast cancer specimens: two untreated specimens and three specimens treated with glycerol solution for 20, 40, or 60 min. All tumors used in the testing data have healthy tissues adjacent to cancerous ones consistent with the challenge faced in lumpectomy surgeries.
Results: The glycerol-treated tissues showed a decrease in the absorption coefficients compared with untreated tissues, especially as the period of treatment increased. Although the sensitivity metric of the classifier presented higher values in the untreated tissues compared with the treated ones, the specificity and accuracy metrics demonstrated higher values for the treated tissues compared with the untreated ones.
Conclusions: The biocompatible glycerol solution is a potential optical clearance agent in THz imaging while keeping the histopathology imaging intact. The SAM technique provided a good classification of cancerous tissues despite the small amount of cancer in the training data (only 7%). The SAM exponential kernel and EMD presented classification accuracy of ∼80 % to 85% compared with linear and polynomial kernels that provided accuracy ranging from 70% to 80%. Overall, glycerol treatment provides a potential improvement in cancer classification in freshly excised breast tumors.
Purpose: The objective of this study is to quantitatively evaluate terahertz (THz) imaging for differentiating cancerous from non-cancerous tissues in mammary tumors developed in response to injection of N-ethyl-N-nitrosourea (ENU) in Sprague Dawley rats.
Approach: While previous studies have investigated the biology of mammary tumors of this model, the current work is the first study to employ an imaging modality to visualize these tumors. A pulsed THz imaging system is utilized to experimentally collect the time-domain reflection signals from each pixel of the rat’s excised tumor. A statistical segmentation algorithm based on the expectation-maximization (EM) classification method is implemented to quantitatively assess the obtained THz images. The model classification of cancer is reported in terms of the receiver operating characteristic (ROC) curves and the areas under the curves.
Results: The obtained low-power microscopic images of 17 ENU-rat tumor sections exhibited the presence of healthy connective tissue adjacent to cancerous tissue. The results also demonstrated that high reflection THz signals were received from cancerous compared with non-cancerous tissues. Decent tumor classification was achieved using the EM method with values ranging from 83% to 96% in fresh tissues and 89% to 96% in formalin-fixed paraffin-embedded tissues.
Conclusions: The proposed ENU breast tumor model of Sprague Dawley rats showed a potential to obtain cancerous tissues, such as human breast tumors, adjacent to healthy tissues. The implemented EM classification algorithm quantitatively demonstrated the ability of THz imaging in differentiating cancerous from non-cancerous tissues.
Terahertz imaging and spectroscopy characterization of freshly excised breast cancer tumors are presented in the range 0.15 to 3.5 THz. Cancerous breast tissues were obtained from partial or full removal of malignant tumors while healthy breast tissues were obtained from breast reduction surgeries. The reflection spectroscopy to obtain the refractive index and absorption coefficient is performed on experimental data at each pixel of the tissue, forming tomographic images. The transmission spectroscopy of the refractive index and absorption coefficient are retrieved from experimental data at few tissue points. The average refractive index and absorption coefficients for cancer, fat, and collagen tissue regions are compared between transmission and reflection modes. The reflection mode offers the advantage of retrieving the electrical properties across a significantly greater number of points without the need for sectioning or altering the freshly excised tissue as in the transmission mode. The terahertz spectral power images and the tomographic images demonstrated good qualitative comparison with pathology.
This paper investigates terahertz (THz) imaging and classification of freshly excised murine xenograft breast cancer tumors. These tumors are grown via injection of E0771 breast adenocarcinoma cells into the flank of mice maintained on high-fat diet. Within 1 h of excision, the tumor and adjacent tissues are imaged using a pulsed THz system in the reflection mode. The THz images are classified using a statistical Bayesian mixture model with unsupervised and supervised approaches. Correlation with digitized pathology images is conducted using classification images assigned by a modal class decision rule. The corresponding receiver operating characteristic curves are obtained based on the classification results. A total of 13 tumor samples obtained from 9 tumors are investigated. The results show good correlation of THz images with pathology results in all samples of cancer and fat tissues. For tumor samples of cancer, fat, and muscle tissues, THz images show reasonable correlation with pathology where the primary challenge lies in the overlapping dielectric properties of cancer and muscle tissues. The use of a supervised regression approach shows improvement in the classification images although not consistently in all tissue regions. Advancing THz imaging of breast tumors from mice and the development of accurate statistical models will ultimately progress the technique for the assessment of human breast tumor margins.
Photoconductive antennas (PCAs) have been extensively utilized for the generation and detection of both pulsed broadband and single frequency continuous wave terahertz (THz) band radiation. These devices form the basis of many THz imaging and spectroscopy systems, which have demonstrated promising applications in various industries and research fields. The development of THz PCA technology through the last 30 years is reviewed. The key modalities of improving device performance are identified, and literature is reviewed to summarize the progress made in these areas. The goal of this review is to provide a collection of all relevant literature to bring researchers up to date on the current state and remaining challenges of THz PCA technology.
This work focuses on pulsed terahertz imaging for the application of surgical margin assessment of breast cancer. Various phantom tissue types and orientations are tested here to refine imaging methodology that can detect breast cancer up to 0.5-1.0 mm from the edge of the sample. The depth of the cancer within the sample is estimated using time of flight analysis of the reflected peaks in the pulsed time domain signal. Breast tissue phantoms have been designed to resemble fresh infiltrating ductal carcinoma, fibroglandular tissue, and fatty tissue of the breast to accomplish this work.
The goal of this work was to develop phantoms that match the refractive indices and absorption coefficients between 0.15 and 2.0 THz of the freshly excised tissues commonly found in breast tumors. Since a breast cancer tumor can contain fibrous and fatty tissues alongside the cancerous tissues, a phantom had to be developed for each. In order to match the desired properties of the tissues, oil in water emulsions were solidified using the proven phantom component TX151. The properties of each potential phantom were verified through THz time-domain spectroscopy on a TPS Spectra 3000. Using this method, phantoms for fibrous and cancerous tissue were successfully developed while a commercially available material was found which matched the optical properties of fatty tissue.
This paper evaluates image enhancement and visualization techniques for pulsed terahertz (THz) images of tissue samples. Specifically, our research objective is to effectively differentiate between heterogeneous regions of breast tissues that contain tumors diagnosed as triple negative infiltrating ductal carcinoma (IDC). Tissue slices and blocks of varying thicknesses were prepared and scanned using our lab’s THz pulsed imaging system. One of the challenges we have encountered in visualizing the obtained images and differentiating between healthy and cancerous regions of the tissues is that most THz images have a low level of details and narrow contrast, making it difficult to accurately identify and visualize the margins around the IDC. To overcome this problem, we have applied and evaluated a number of image processing techniques to the scanned 3D THz images. In particular, we employed various spatial filtering and intensity transformation techniques to emphasize the small details in the images and adjust the image contrast. For each of these methods, we investigated how varying filter sizes and parameters affect the amount of enhancement applied to the images. Our experimentation shows that several image processing techniques are effective in producing THz images of breast tissue samples that contain distinguishable details, making further segmentation of the different image regions promising.
This work seeks to obtain the properties of paraffin-embedded breast cancer tumor tissues using transmission imaging and spectroscopy. Formalin-fixed and paraffin-embedded breast tumors are first sectioned into slices of 20 μm and 30 μm and placed between two tsurupica slides. The slides are then scanned in a pulsed terahertz system using transmission imaging. The tissue regions in adjacent pathology section are compared to the transmission imaging scan in order to define a region of points over which to average the electrical properties results from the scan.
This work aims to enact a quick and reasonable estimation of breast cancer margin thickness using time of flight analysis of embedded breast cancer tissue. A pulsed terahertz system is used to obtain reflection imaging scans from breast cancer tumors that are formalin-fixed and embedded in paraffin blocks. Time of flight analysis is then used to compare the reflection patterns seen within the block to pathology sections and paraffin-embedded sections that are taken throughout the depth of the tumor in order to estimate the three-dimensional boundaries of the tumor.
In this work a new plasmonic thin-film based terahertz photoconductive antenna is proposed. The computational method utilized to design the antenna is outlined, as well as the steps and preliminary results for the fabrication and characterization of the device. The model predicted over two orders of magnitude increase in the peak photocurrent as compared to a conventional device design, while slightly reducing the width of the induced current pulse. This indicates that the proposed design will be effective as a high efficiency terahertz emitter. In addition to the computational modeling, preliminary results demonstrating the proposed fabrication processes and experimental characterization are presented. It is demonstrated that when using a pyroelectric detector to quantify the output terahertz power it is important to first quantify the power of the IR photons generated by thermal relaxation in the device.
The presented reconstruction algorithm is based on merging the fast multipole method (FMM) for the forward solver with the rapidly convergent descent method for the cost function minimization algorithm (by Fletcher and Powell 1964 and 1970). Parametric results are presented showing the potential of the proposed computational algorithm.
KEYWORDS: Scattering, Dielectrics, Polarimetry, 3D acquisition, Monte Carlo methods, Land mines, Dielectric polarization, Target detection, Free space, Electrical engineering
This work presents the Mueller matrix elements for scattering from dielectric targets buried beneath 2-D random rough surfaces (3-D scattering problem). The fully polarimetric scattering matrix S is computed for hundreds of computer generated random rough surface realizations and hence the Mueller matrix elements are obtained. The numerical results show that if one relies only on the co- or cross-polarized intensities (i.e. the magnitude of the four elements of the polarimetric scattering matrix S); it is very difficult to sense the buried objects. However investigating all the sixteen Mueller matrix elements greatly helps in detecting these objects.
In this work, radar images of dielectric targets buried beneath 2-D rough ground surfaces are simulated. These images are based on the amplitude and phase of all polarization of incident and scattered waves. Both bistatic and backscatter intensity-based images are considered in this work to compare the image quality in each case. All simulations are obtained using the integral equation based technique; the Steepest Descent Fast Multiple method (SDFMM). The numerical results show that images based on the Covariance matrix are very unique when the surrounding illuminated areas have almost the same rough surface profile. However, these images deteriorate as the surface profile greatly changes in the surrounding ground areas.
In realistic landmine fields, the anti-personnel plastic mine is often buried nearby a clutter-object under the ground. The presence of a second object buried near the mine under a two-dimensional (2-D) rough ground can easily obscure the target and/or cause a false alarm. The separation distance between the AP mine and clutter-object plays a significant role on the probability of true or false alarm in this situation. A rigorous electromagnetic model has been developed to analyze the scattering mechanism between the target and the clutter-object, between the target and the rough ground, between clutter-object and the rough ground and the multiple scattering between different spots on the rough ground itself. The new rigorous model is based on the classical electromagnetic equivalence theorem leading to producing six new integral equations. Using the Method of Moment (MoM), the new integral equations are transformed into a linear system of equations to be solved for the unknown electric and magnetic currents on the surface of three scatterers; 1) rough ground, 2)target and 3)clutter-object. The MoM impedance matrix completely represents every interaction between these three scatterers. The superior Steepest Descent Fast Multipole Method (SDFMM) is used to tremendously speed up the computations of the unknown MoM surface currents.
KEYWORDS: Finite-difference time-domain method, Signal processing, Target detection, Land mines, General packet radio service, Dielectrics, Reflection, Surface roughness, Receivers, Monte Carlo methods
In ground penetrating radar (GPR) antipersonnel mine sensing, in which the target is small, shallow and often of low dielectric contrast, detection is challenging. One of the difficulties is that it is hard to distinguish the target signal from the omnipresent random rough ground reflection clutter. In this work, a Monte Carlo computational simulation using 2-dimensional (2-D) transverse magnetic (TM) finite difference time domain (FDTD) with multiple rough surfaces is implemented to investigate single TNT target buries in dispersive soil. Based on the effects of the random rough surface on an impulse GPR signal and the knowledge of wave propagation differences in different media - air, soil, and TNT - a special background average process using physics based signal processing (PBSP) is performed to remove the ground clutter signal. This procedure first involves shifting and scaling multiple time signals from target-free random rough ground to establish the nominal (average) ground reflection pulse shape. Next, this nominal pulse shape is correlated in time with each trial signal, then shifted and scaled to match the ground surface clutter of that trial signal. Subtracting this shifted scaled clutter signal from the trial signal ideally leaves the target signal (with some additional multiple scattering between the target and ground surface). The PBSP algorithm reapplied in cases for which surface scattering occurs at multiple points. The statistical results of PBSP surface clutter removal indicate that the detection performance degrades with increasing surface roughness and decreasing burial depth. Hypothesis testing on the processed results proved to be successful in a detection and estimation point of view. This paper presents the detection performances in terms of Receiver Operating Characteristics (ROC) for various ground surface roughness and target burial depth cases. Also demonstrated is the performance improvements expected from multiple views: indicating that a multi-bistatic configuration appears to be superior to multistatic transmitter/receiver geometry with minimum combinations.
Scattering of electromagnetic waves from multilayered random rough surfaces is crucial for subsurface sensing applications. A multiple interaction method of moments (MoM) model is used in this work to analyze scattering from two-dimensional multilayered random rough ground (3-D scattering problem) especially when the underground layer is deeply buried under the air/ground interface. The presented model removes a barrier and enables the application of the Steepest Descent Fast Multipole Method (SDFMM) to certain 3-D non-quasi-planar structures. The conventional SDFMM has been used to analyze electromagnetic wave scattering from quasi-planar structures where the scatterer's height is a fraction of a free-space wavelength. The presented model is based on multiple interactions mechanism between the air/ground interface and the buried underground layer. The basic idea of the proposed multiple scattering model is to decompose the non-quasi-planar multilayered ground into two quasi-planar scatterers where the conventional SDFMM can be applied separately to each one. The interactions between the sub-quasi-planar scatterers are calculated using the electromagnetic vector potentials near-field expressions. This model is tested and validated with the MoM on a variety of geometries. The results show that the strongest signature of the buried scatterer is mainly due to the first multiple interaction mechanism (ground-object-ground) while the contributions from repeating this mechanism become insignificant even for lossless and/or slightly lossy underground.
KEYWORDS: Mining, Interfaces, Dielectrics, Near field, Scattering, Gaussian beams, Monte Carlo methods, Optical spheres, Electromagnetism, General packet radio service
The Steepest Descent Fats Multilevel Multiple Method (SDFMM) is used to analyze the distorting effects of random rough ground surfaces on scattered electromagnetic waves from buried TNT mines. The SDFMM method is an integral equation- based fast algorithm that is well suited for 2D penetrable rough surfaces in the frequency domain, and it is used to calculate the unknown surface currents on both the rough ground and the buried target as well. In this study all interactions between the rough interface and the buried target are taken into account. The scattered near field E- patterns of an incident Gaussian beam are calculated at different locations above the mean plane of the dielectric rough interface. The receiver locations are chosen to simulate GPR measurement protocols. The dimensions and burial depth of the TNT mine are smaller than the free space wavelength with material slightly different from the surrounding soil. The average and the standard deviation of the scattered fields for just the target are calculated and results showed that the presence of the rough interface tremendously distorts the target signal even for the small roughness parameters. Moreover, results showed the degradation of signal as the TNT mine is located away from incident beam. This knowledge can significantly contribute to inventing better sensing systems for less false alarm detection strategies.
The finite difference time domain technique, FDTD, is used to calculate the scattered field din the n era zone from 1D random rough surfaces. Different statistics for the random surface will be assumed in this work. First, the random rough surface will be characterized by one-scale roughness with Gaussian distribution for the heights and Gaussian auto-correlation function. In the second part, the surface will be assumed to have two-scale roughness with the same Gaussian statistics as before. The statistics of the scattered fields are calculated in this work using Monte Carlo simulations. Numerical results comparing scattered fields from one-scale roughness and two-scale roughness are shown. The results obtained indicate that the distortion in the scattered signals is primarily due to the small-scale roughness while the two-scale roughness causes more time delay In the scattered signals. Different rough surface parameters will be used to quantify their effect on the statistics of scattered signals.
The like and cross polarized single and double scattered fields are derived using a full wave approach. This approach is based on the complete expansion of the electromagnetic fields, the imposition of exact boundary conditions and the conversion of Maxwell's equations into generalized telegraphists equations for the scattered wave amplitudes. Thus, the zero order iterative solutions for the generalized telegraphists equations yield the primary electromagnetic (source) fields impressed upon the rough surface. The first and second order iterative solutions to the generalized telegropherts' equations yield the single and double scattered fields. This can be clearly demonstrated by taking the geometric optics limit of the full wave solutions. To obtain the corresponding like and cross polarized scatter cross sections, as in the case of scattering from one dimensional rough surfaces, it is necessary to account for contributions from the quasi parallel double scatter paths as well as the quasi antiparallel double scatter paths. However, for scattering from two dimensional rough surfaces, these paths are not restricted to the plane of incidence. The full wave solutions for the double scattered fields are expressed as six dimensional integrals, that account for the complete wave spectra of scattered fields and the coordinate variables at a pair of points on the rough surface. These expressions are used to obtain the multidimensional integrals for the like and cross polarized cross sections. To make these solutions tractable for computational purposes, a high frequency approximation of the full wave double scatter cross sections are expressed as four dimensional integrals involving scatter wave vector variables. These results can be evaluated in significantly less time than standard numerical solutions of the integral equations. Moreover, the physical interpretation of the results shed light on the impact, of different statistical parameters of the random rough surfaces, upon the backscatter enhancement.
Access to the requested content is limited to institutions that have purchased or subscribe to SPIE eBooks.
You are receiving this notice because your organization may not have SPIE eBooks access.*
*Shibboleth/Open Athens users─please
sign in
to access your institution's subscriptions.
To obtain this item, you may purchase the complete book in print or electronic format on
SPIE.org.
INSTITUTIONAL Select your institution to access the SPIE Digital Library.
PERSONAL Sign in with your SPIE account to access your personal subscriptions or to use specific features such as save to my library, sign up for alerts, save searches, etc.