Ultra-Wideband Millimeter-Wave Dielectric Characteristics of Freshly Excised Normal and Malignant Human Skin Tissues
Ultra-Wideband Millimeter-Wave Dielectric Characteristics of Freshly Excised Normal and Malignant Human Skin Tissues
6, JUNE 2018
Abstract—Millimeter waves have recently gained atten- nonionizing electromagnetic waves. These methods are based
tion for the evaluation of skin lesions and the detection of on the inherent contrast between the electrical parameters of
skin tumors. Such evaluations heavily rely on the dielectric malignant and normal tissues [3]–[5]. It has been shown that
contrasts existing between normal and malignant skin tis-
sues at millimeter-wave frequencies. However, current stud- cancer changes the water content as well as the biochemistry
ies on the dielectric properties of normal and diseased skin (e.g., metal concentration) of tissues [6]–[8], and hence results
tissues at these frequencies are limited and inconsistent. In in a change in their dielectric properties.
this study, a comprehensive dielectric spectroscopy study Microwave imaging techniques have been well developed
is conducted for the first time to characterize the ultra-
over the past years for the detection of breast cancer, lung
wideband dielectric properties of freshly excised normal
and malignant skin tissues obtained from skin cancer pa- cancer, and brain stroke in the frequency range of 300 MHz–
tients having undergone Mohs micrographic surgeries at 10 GHz [9]–[12]. As the frequency increases to the millimeter-
Hackensack University Medical Center. Measurements are wave regime (30–300 GHz), shorter wavelengths allow for the
conducted using a precision slim-form open-ended coaxial achievement of higher spatial resolutions at the cost of reduced
probe in conjunction with a millimeter-wave vector network penetration depths. This makes millimeter waves less suitable
analyzer over the frequency range of 0.5–50 GHz. A one-
pole Cole–Cole model is fitted to the complex permittivity for imaging tissues placed deep within the body. However, the
dataset of each sample. Statistically considerable contrasts penetration depth of these waves varies from 600 μm to 1.2
are observed between the dielectric properties of malignant mm into the body [13], making them very effective for sensing
and normal skin tissues over the ultra-wideband millimeter- pathological changes in different skin layers or in the outer tis-
wave frequency range considered. sue layers of excised organs. It is those layers (just below the
Index Terms—Contrast study, dielectric measurements, surface of the skin or organs) which contain certain epithelial
millimeter-wave measurements, malignant skin tissue, nor- cells where a majority of primary tumors originate [14], [15].
mal skin tissue, skin cancer. As the frequency further increases into the THz regime, electro-
I. INTRODUCTION magnetic waves barely penetrate the surface of the tissues [16]
and therefore have limited ability to detect early-stage tumors
KIN cancer is the most common and fastest growing of
S all cancer types with approximately 3.5 million new cases
diagnosed each year in the US [1]. As all kinds of skin cancer
which reside in deeper tissue layers.
The success and efficacy of using millimeter waves for detect-
ing skin tumors largely relies on the dielectric contrasts between
can be easily treated and managed if diagnosed at an early stage normal and malignant skin tissues at these frequencies. How-
[2], it is imperative that dermatologists be astute in the early ever, a definitive knowledge of such contrasts does not currently
diagnosis and treatment of suspected malignancies. exist as no consistent studies have been performed on evalu-
Recently, there has been considerable interest in the de- ating the millimeter-wave dielectric properties of normal and
tection and management of different types of cancer using malignant skin tissues. There is only one study that measures
the reflection characteristics of basal cell carcinoma skin tissues
Manuscript received April 10, 2017; revised August 7, 2017; accepted in 15 skin cancer patients at the frequency of 42 GHz using a
August 30, 2017. Date of publication September 6, 2017; date of current customized WR-22 waveguide [17]. Although large differences
version May 18, 2018. This work was supported by the U.S. National
Science Foundation under Grant CAREER 1554402. (Corresponding were observed between the reflection properties of tumorous
author: Negar Tavassolian.) and normal skin, the study was performed qualitatively and in a
A. Mirbeik-Sabzevari is with the Department of Electrical and Com- limited and small scale. In addition, they used a relatively large
puter Engineering, Stevens Institute of Technology.
R. Ashinoff is with the Hackensack University Medical Center and also and flange-posing waveguide probe which does not provide a
with the New York University/Langone Medical Center. firm and consistent contact with the skin tissue. In this work,
N. Tavassolian is with the Department of Electrical and Computer we systematically evaluate the dielectric properties of a large
Engineering, Stevens Institute of Technology, Hoboken, NJ 07030 USA
(e-mail: [email protected]). number of ‘freshly-excised’ normal and malignant skin tissue
Digital Object Identifier 10.1109/TBME.2017.2749371 samples over an ultra-wide frequency range of 0.5–50 GHz.
0018-9294 © 2017 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission.
See http://www.ieee.org/publications standards/publications/rights/index.html for more information.
Authorized licensed use limited to: University of Applied Sciences Karlsruhe c/o KIT Library. Downloaded on December 23,2021 at 10:41:29 UTC from IEEE Xplore. Restrictions apply.
MIRBEIK-SABZEVARI et al.: CHARACTERISTICS OF FRESHLY EXCISED NORMAL AND MALIGNANT OF HUMAN SKIN TISSUES 1321
Authorized licensed use limited to: University of Applied Sciences Karlsruhe c/o KIT Library. Downloaded on December 23,2021 at 10:41:29 UTC from IEEE Xplore. Restrictions apply.
1322 IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 65, NO. 6, JUNE 2018
TABLE I TABLE II
DETAILS OF TISSUE COLLECTION CONDITIONS A FULL LIST OF TISSUES MEASURED EX-VIVO
Authorized licensed use limited to: University of Applied Sciences Karlsruhe c/o KIT Library. Downloaded on December 23,2021 at 10:41:29 UTC from IEEE Xplore. Restrictions apply.
MIRBEIK-SABZEVARI et al.: CHARACTERISTICS OF FRESHLY EXCISED NORMAL AND MALIGNANT OF HUMAN SKIN TISSUES 1323
III. METHODS
A. Data Fitting
A one-pole Cole–Cole model [23] is applied to all dielectric
parameters obtained from measurements. Cole–Cole models are
normally used to represent physics-based, frequency-dependent
dielectric properties of biological tissues over wide frequency
bands [24]. It has been shown that a one-pole Cole–Cole model
is sufficient to fit human tissue characteristics [25]. The follow-
ing representation is accordingly considered for each of the 152
datasets: Fig. 2. Histology distribution of the 146 datasets measured for this
study.
εr σs
εr (ω) = εr (ω) − jεr (ω) = εr o (ω) + +
1 + ( jωτ ) 1−α jωε0
(1) Group 4 contains all datasets obtained from malignant SCC
where ω is the angular frequency, εr (ω) and εr (ω) are the real tissues. Fig. 2 shows a histogram of the histology distributions
and imaginary parts of the frequency-dependent dielectric per- of the four groups. We separate negative (normal) BCC tissues
mittivity, and εr o , εr , τ , σs , and α are the Cole–Cole model from negative (normal) SCC tissues to eliminate any possible
parameters which are determined by minimizing the following sources of variabilities in the dielectric properties of the
cost function: corresponding samples. As presented in Table II, the collection
N εr (ωi )−εr c (ωi ) N εr ωi )−εrc (ωi ) process (i.e., body location and time of excision) is different
i=1 (εr (ωi )) + i=1 (εr (ωi )) for each acquired sample. Therefore, we compare malignant
c= (2) BCC/SCC datasets with normal tissues obtained from either
N
the same sample or an adjacent (bordering) sample. This kind
Here, N = 142 is the number of frequency spots of the mea-
of paired comparison is employed to eliminate any variation
sured data in the frequency range of 0.5–50 GHz (frequency
between the samples in each group.
step of 358 MHz), εr (ωi ) and εr (ωi ) are the experimental val-
Averaged dispersion curves are obtained for each group by
ues measured at the frequency of ωi , and εr c (ωi ) and εrc (ωi )
averaging the permittivity values for all samples in the group
are the values generated by (1). The fitting procedure is ac-
at 142 equally-spaced frequency points, followed by apply-
complished in MATLAB [26] using the Levenberg–Marquardt
ing a one-pole Cole–Cole model. The standard deviations of
algorithm.
the averaged curves at approximately 10, 20, 30, and 40 GHz
The dielectric datasets of the samples that cannot be fitted
are also calculated to specify the extent of the curves in each
with a Cole–Cole model (i.e., the fit error between the model
group.
and measured data across the entire range of frequencies exceeds
1) Comparison Between Dielectric Properties of Normal
a pre-determined threshold of 15%), are excluded from further
and Malignant Skin Tissues: Multivariate analysis of vari-
analysis. The threshold value in our fitting procedure was cho-
ance (MANOVA) [28] is employed to evaluate the dielectric
sen to have scaled Euclidean distances between the measured
properties of the samples. MANOVA is a generalized form of
data and the fitted dielectric curves (averaged over frequency
univariate analysis of variance (ANOVA) and compares the mul-
as specified by (2)) well below those distances considered in
tidimensional (here two-dimensional (2D)) distributions of the
similar studies of biological human tissues [5], [24], [27]. A
responses of the tissue groups rather than one dimension at a
total number of 6 datasets were excluded based on this crite-
time. It therefore provides a more robust analysis when the re-
rion, leaving 146 datasets for statistical analysis. The overall
sponse has several dependent vectors. The 2D responses in our
distribution of the remaining datasets are: 98 normal datasets,
case comprise the real and imaginary parts of dielectric permit-
25 BCC datasets, and 23 SCC datasets. Experimental errors
tivities.
such as imperfect contact of the probe aperture with the tissue
MANOVA as described above is used to compare the dielec-
surface and/or unsatisfactory sizes of tissues under assessment
tric properties of normal and malignant samples at roughly 11
can contribute to the failed fits.
equally-spaced frequency spots across 1–50 GHz. Two sepa-
rate analyses are conducted to compare BCC-normal and SCC-
B. Statistical Analyses normal tissues. As mentioned earlier, the two categories in each
For statistical analyses, four groups of datasets are formed analysis are paired. In other words, for each observation in either
according to the three diagnosis types of the tissues col- of the SCC or BCC cancer groups, there is a normal observation
lected. Group 1 contains all measured datasets obtained from in the corresponding normal groups which is obtained from an
negatively-diagnosed BCC tissues, Group 2 contains all datasets adjacent tissue. This adjacent tissue can be either in the same
obtained from negatively-diagnosed SCC tissues, Group 3 sample, another adjacent sample, or an average of both in cases
contains all datasets obtained from malignant BCC tissues, and where both datasets are available.
Authorized licensed use limited to: University of Applied Sciences Karlsruhe c/o KIT Library. Downloaded on December 23,2021 at 10:41:29 UTC from IEEE Xplore. Restrictions apply.
1324 IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 65, NO. 6, JUNE 2018
Fig. 3. Examples of successful Cole–Cole fits to three representative experimental datasets. (a) Real part, (b) imaginary part of dielectric
permittivity as a function of frequency for a normal (negative BCC) sample; (c) real part, (d) imaginary part of dielectric permittivity as a function of
frequency for a normal (negative SCC) sample; (e) real part, (f) imaginary part of dielectric permittivity as a function of frequency for a BCC cancer
tissue sample; (g) real part, (h) imaginary part of dielectric permittivity as a function of frequency for an SCC cancer tissue sample. o: measured
data, solid line: Cole–Cole fit.
2) Comparison Between Patient-to-Patient and Sample- datasets in each subgroup correspond to the same location on the
to-Sample (Within Patient) Variabilities: Two separate nested body). ANOVA will assess the differences between the means of
analyses of variance (ANOVA) tests are conducted at roughly the subgroups (variability between patients) as well as the spread
11 equally-spaced frequency spots across 1–50 GHz (1 GHz, of the observations around the means (sample-to-sample vari-
5 GHz, 10 GHz, 15 GHz, . . . , 50 GHz) to examine whether ability within patients). In this way, the variations between the
the variability in the dielectric permittivity between patients is samples from different patients are compared with the variations
larger than the variability between samples obtained from each between the samples obtained from the same patient (adjacent
patient. One analysis is conducted on the real parts of the dielec- ones).
tric permittivities and the other on the imaginary parts. Group 1
(negative BCC) is considered for this analysis as it has a larger IV. RESULTS AND DISCUSSION
number of datasets (50 datasets from 23 patients) compared to
A. Data Fitting and Exclusion
the other three groups and therefore more precisely satisfies the
normal distribution assumption required for ANOVA analyses. Figs. 3 and 4 demonstrate the physics-based criterion em-
To apply ANOVA, 23 subgroups are formed with each subgroup ployed to ensure that the dielectric measurements met the
containing the datasets obtained from one patient (therefore all Cole-Cole model consistency requirement. Fig. 3 shows four
Authorized licensed use limited to: University of Applied Sciences Karlsruhe c/o KIT Library. Downloaded on December 23,2021 at 10:41:29 UTC from IEEE Xplore. Restrictions apply.
MIRBEIK-SABZEVARI et al.: CHARACTERISTICS OF FRESHLY EXCISED NORMAL AND MALIGNANT OF HUMAN SKIN TISSUES 1325
Authorized licensed use limited to: University of Applied Sciences Karlsruhe c/o KIT Library. Downloaded on December 23,2021 at 10:41:29 UTC from IEEE Xplore. Restrictions apply.
1326 IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 65, NO. 6, JUNE 2018
The small variability observed is due to the different locations (SD = 2.20), while the average imaginary part for SCC cancer
of body from which the datasets are obtained. As mentioned samples is 12.84 (SD = 1.63), which is 16% lower. Similar
before, each group of normal samples included a variety of dif- trends are observed at other frequencies. The largest difference
ferent body locations. If the exact dielectric characteristics are between the average dielectric properties of normal and SCC
needed for each specific location on the body, separate curves cancer samples is approximately 111% (for the real part near
can be obtained followed by another fitting process to produce 50 GHz).
Cole–Cole parameters accordingly. 3) Comparison Between Patient-to-Patient and Sample-
The results obtained here also indicate that over the entire to-Sample (Within Patient) Variabilities: Two ANOVA tests
frequency range of 0.5–50 GHz, malignant BCC (Group 3) were conducted on the 23 subgroups formed from Group 1
has higher real and imaginary parts of dielectric permittivity to evaluate patient-to-patient and sample-to-sample variabili-
compared to normal skin tissue (Group 1). On average, the real ties. Statistically significant differences ( p < 0.05) were ob-
part of the dielectric permittivity is 6% higher near 10 GHz. This served for both the real and imaginary parts of the dielectric
difference exists at higher frequencies as well, with a maximum permittivities at all the 10 frequency spots considered across the
percentage difference of 15% near 28 GHz. The imaginary part 5–50 GHz range. For instance, patient-to-patient and sample-
of the dielectric permittivity is 9% higher, on average, near 10 to-sample standard deviations (SDs) near 40 GHz for the real
GHz for malignant tissues. This difference becomes larger at parts of the dielectric permittivity were 2.85 and 1.16, re-
higher frequencies, with a maximum percentage difference of spectively. For the imaginary part, the SDs were 3.32 and
12% near 50 GHz. 1.85. This means that as expected, there is a greater vari-
SCC tumor tissues (Group 4) also have a higher real part of ability between the patients’ subgroups as compared to the
dielectric permittivity in comparison with their corresponding samples within each patient subgroups. This is expected as
normal skin tissues (Group 2). However, the imaginary part tissues are excised from different body locations for different
for malignant tissues is less than that of normal tissues over patients. However, as explained in Section III-B, this alteration
the entire range of 0.5–50 GHz. For instance, near 10 GHz, does not affect the analysis conducted between normal and ma-
the average real part of the dielectric permittivity is 34.95 for lignant datasets as paired categories were considered for that
normal samples (Group 2), while this value is 41.12 for SCC analysis.
samples (Group 4). The maximum difference between the real
parts of the SCC and normal samples is 111% near 50 GHz. For C. Discussion
the imaginary part, malignant SCC tissues demonstrate lower The results obtained in our study, as demonstrated in Fig. 5,
values over the entire frequency range of 0.5–50 GHz with a indicate that on average, malignant BCC tissues (Group 3)
maximum difference of 49% near 10 GHz. demonstrate a slightly higher dielectric permittivity (both real
2) Comparison Between Normal and Malignant Tissue
and imaginary parts) compared to normal skin tissues (Group
Dielectric Properties: We found statistically significant dif-
1) over the entire frequency range of 0.5–50 GHz. To analyze
ferences between the dielectric properties of normal and BCC and elaborate on these results, we need to understand the di-
cancer samples (p < 0.05) at all the 11 frequency spots consid- electric behavior of the tissues as well as the main components
ered across the 1–50 GHz frequency range. Near 40 GHz, the comprising the skin tissues.
average real part of the dielectric permittivity for normal tissue In [29], three dispersion mechanisms (α, β, γ ) are introduced
samples (Group 1) is 15.74 with a standard deviation (SD) of to characterize the anomalous electric properties of biomaterials.
2.70. At the same frequency, the average real part of the di- The β-dispersion is linked to the cellular structure of biological
electric permittivity for BCC cancer samples (Group 3) is 17.04 materials with low-frequency (up to several MHz) properties
(SD = 1.87), which is 8.5% higher. Similarly, near 40 GHz, the caused by cell membranes [29]. The α-dispersion is attributed
average imaginary part of the dielectric permittivity for normal to the surface admittance which is lateral to the membrane sur-
tissue samples is 14.89 (SD = 3.88), while the average imagi- face. Alpha dispersion can be from below 1 Hz up to 100 kHz
nary part for BCC cancer samples is 16.24 (SD = 1.54), which [29]. The α- and β-dispersion mechanisms fall below the low-
is 9% higher. We see similar trends at other frequencies. The est frequency in the measured data of our study (∼ 500 MHz).
largest difference between the averaged dielectric properties of Hence, they are not considered in our analyses.
normal and BCC cancer tissues is approximately 15% (for the The γ -dispersion is the dielectric dispersion due to the dipolar
real part near 28 GHz). relaxation of water and the capacitive charge of the cellular
Statistically significant differences were also found between membranes in the tissues [29]. The relaxation frequency of the
the dielectric properties of normal and SCC cancer samples γ -dispersion can be roughly estimated from the frequency at
( p < 0.05) at all the 11 frequency spots considered across the which the maximum loss (the imaginary part of the dielectric
1–50 GHz frequency range. Near 40 GHz, the average real part permittivity) is obtained. It is calculated by [29]:
of the dielectric permittivity for normal tissue samples (Group
2) is 15.31 with a standard deviation of 2.59. The average real 1
fc = (3)
part of the dielectric permittivity for SCC cancer tissues (Group 2π τ
4) is 30.28 (SD = 3.17) at the same frequency, which is 97% For the normal and malignant BCC, this frequency is cal-
higher. Similarly, near 40 GHz, the average imaginary part of culated from Equation (3) and Table III to be approximately
the dielectric permittivity for normal tissue samples is 15.27 20 GHz and 21 GHz, respectively. On the other hand, the re-
Authorized licensed use limited to: University of Applied Sciences Karlsruhe c/o KIT Library. Downloaded on December 23,2021 at 10:41:29 UTC from IEEE Xplore. Restrictions apply.
MIRBEIK-SABZEVARI et al.: CHARACTERISTICS OF FRESHLY EXCISED NORMAL AND MALIGNANT OF HUMAN SKIN TISSUES 1327
laxation frequency of water is 19 ± 1 GHz according to the where εbm is the permittivity of the dry biological structural
parametric model of [29] in the range of tissue temperatures in components (chosen to be 3 to represent an effective permittivity
our study. The relaxation frequency agreements between BCC of protein and lipid molecules [39]), εw is the permittivity of
(normal and malignant) tissues and pure water suggest that the water (taken from the Cole-Cole parametric model at 25 °C
relaxation frequency of the γ -dispersion for the normal and ma- [40]) and is the volume fraction of the water component in
lignant BCC skin tissues may be largely influenced by that of the tissue. The 10% increase in the real and imaginary parts
pure water, and relatively independent of other skin constituents. of the dielectric permittivity obtained from (4) is consistent
To confirm this hypothesis, we need to understand the main with the results of our study (according to Fig. 5, the real and
constituents of skin tissues as well as the molecular alterations imaginary parts of the dielectric permittivity are 8.5% and 9%
in skin cancer tissues compared to normal tissues. Living or- higher in BCC tissues compared to their corresponding normal
ganisms are primarily composed of water, organic molecules, tissues, respectively). At other frequencies also, the differences
and ions [31]. Water constitutes 70–85% of the total body mass. predicted by (4) agree well with the data obtained in our study.
It constitutes 60–72% of the skin mass [32]. It is known that In contrast to the case of BCC tissues, for malignant SCC tis-
water in biological materials exists in two states: bulk (free) and sues we observed that the real parts of the dielectric permittivity
hydrated (bound) [33]. Hydrated water interacts closely with were larger on average than their corresponding normal ones,
other molecules as well as with bulk water. In the skin tissue, while the imaginary parts were lower than their corresponding
the majority of water molecules seem to be bound only to each normal tissues. The relaxation frequencies of the γ -dispersion
other, similar to the structure of bulk liquid water, rather than for the normal and malignant SCC tissues are calculated from
interacting with other macromolecules [34] and are therefore Equation (3) and Table III to be approximately 20 GHz and
considered to be in the free state. We also note that the dielectric 51 GHz, respectively. The relaxation frequency for the malig-
response of bound water is below 1 GHz and is therefore mostly nant SCC group is far greater than that of water. This indicates
out of the frequency range of interest [35]. that water may not be the main contributor to the dielectric
The remainder of the skin cellular mass is composed of properties of malignant SCC tissues and that they may be more
carbon-containing molecules, which can be divided into four dependent on other skin components.
classes: proteins, nucleic acids, carbohydrates (sugars), and In [41], Raman spectroscopy was employed to differenti-
lipids [29]. Proteins are polymers of amino acids. Proteins con- ate between normal skin and SCC tissues in-vitro. The Raman
stitute 10–25% of the cell mass, the largest contribution after spectra in the normal tissues showed the presence of vibra-
water [32]. DNA and RNA are polymers of nucleotides. DNA tional bands in protein-specific band peaks with higher inten-
and RNA represent 0.3–1% and 0.8–6% of the total cell mass, sities compared to the carcinoma spectra. Another shift region
respectively [32]. Glucose (sugar) level also varies from 0.1–1% in the Raman spectra showed a difference of intensity in the
of the total cell mass [32]. samples of squamous cell carcinoma, which were attributed to
Several Raman spectroscopy studies have been performed to nucleic acids (DNAs). Another Raman spectroscopy study [42]
investigate the molecular structures of different types of skin showed similar changes. In [43] also, variations were observed
cancer. In [36], the nature of water content and structural alter- in the band patterns corresponding to glucose. However, none of
ations in skin malignancy (specifically BCC) were investigated these studies quantified the changes observed in the constituents.
using Raman spectroscopy. An increase of ∼ 15% in the free There is also no study in the literature on the determination of
(bulk) water content was found in malignant BCC skin tissues the water content in malignant SCC tissues.
compared to normal skin tissues. In another study [37], other According to [44], changes in the complex permittivity (both
molecular structure alterations (including proteins and lipids) of real and imaginary values) and a shift in the relaxation frequency
BCC skin cancer tissues were researched. It was seen that the of tissues are expected when the protein and glucose concentra-
Raman intensity of the lipid-specific band peaks were increased tions are increased. These all show that biochemical alterations
due to an increase in the lipid concentration. Also, no significant in malignant SCC tissues, including changes in protein, nucleic
changes were observed in the intensities of protein and nucleic acid, and glucose concentrations play the most important role
acids bands. Both increases in water and lipid concentrations in their dielectric properties. Further studies are required to in-
result in higher dielectric permittivies (both real and imaginary vestigate the water content of SCC tissues as well as quantify
parts) [35], [38]. the effects of each contributor on the millimeter-wave dielectric
These changes are in full agreement with our findings. To characteristics of malignant SCC tissues.
demonstrate this, we assume a binary mixture of water and As also mentioned in Section I, there has been no prior study
other skin components with a concentration of 60% for water to in the literature which consistently evaluates the millimeter-
represent the skin tissue [32]. A 15% increase in the free water wave dielectric properties of normal and malignant skin tissues.
content in this mixture (corresponding to the BCC cancer tissue In the only available study, in-vivo measurements of reflection
[36]) will result in approximately 10% increase in both real and coefficients of normal and malignant BCC tissues were per-
imaginary parts of the dielectric permittivity at 40 GHz per the formed on the faces of 15 skin cancer patients at 42 GHz [17].
Maxwell–Wagner equation [39]: Our results agree with the results of [17]. The magnitude of
the reflection coefficients from BCC samples are lower on av-
(2εbm + εw ) − 2φ(εbm − εw ) erage than the normal ones although the differences in our case
εmix = εbm (4)
(2εbm + εw ) + φ(εbm − εw ) are less than the ones presented in [17] (due to a difference in
Authorized licensed use limited to: University of Applied Sciences Karlsruhe c/o KIT Library. Downloaded on December 23,2021 at 10:41:29 UTC from IEEE Xplore. Restrictions apply.
1328 IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 65, NO. 6, JUNE 2018
ACKNOWLEDGMENT
The authors would like to thank the Mohs Surgery and Derma-
tology Surgery Unit at Hackensack University Medical Center
for the assistance of the personnel.
REFERENCES
[1] American Cancer Society. (2017). Cancer Facts and Figures.
[Online]. Available: https://old.cancer.org/acs/groups/content/@editorial/
documents/document/acspc-048738.pdf
[2] American Cancer Society. (2017). Cancer Survival Rates for Melanoma
Fig. 6. (a) Real and (b) imaginary parts of the dielectric permittivities Skin Cancer, by Stage. [Online]. Available: https://www.cancer.org/
corresponding to negative BCC (Group 1), negative SCC (Group 2), and cancer/melanoma-skin-cancer/detection-diagnosis-staging/survival-
wet skin provided by Gabriel et al. [24]. The literature data agree well rates-for-melanoma-skin-cancer-by-stage.html
(error <10%) with the dielectric properties of our normal BCC and SCC [3] A. P. O’Rourke et al., “Dielectric properties of human normal, malignant
tissues over the entire range of 0.5–50 GHz. and cirrhotic liver tissue: In vivo and ex vivo measurements from 0.5 to
20 GHz using a precision open-ended coaxial probe,” Phys. Med. Biol.,
vol. 52, pp. 4707–4719, 2007.
the measurement setup and tissue conditions). The dielectric [4] P. T. Nguyen et al., “Three-dimensional microwave hyperthermia for
properties of normal BCC and SCC tissues derived in our study breast cancer treatment in a realistic environment using particle swarm
were also compared with the literature values obtained from optimization,” IEEE Trans. Biomed. Eng., vol. 64, no. 6, pp. 1335–1344,
Jun. 2017.
in-vivo human experiments. Dielectric properties provided by [5] M. Lazebnik et al., “A large-scale study of the ultrawideband microwave
Gabriel et al. [24] for wet skin agreed well (error <10%) with dielectric properties of normal, benign and malignant breast tissues ob-
the dielectric properties of our normal BCC and SCC tissues up tained from cancer surgeries,” Phys. Med. Biol., vol. 52, pp. 6093–6115,
2007
to the frequency of 50 GHz as shown in Fig. 6. The data for wet [6] V. Suntzeff and C. Carruthers, “The water content in the epidermis of mice
skin at the 0.5–20-GHz range were obtained from experiments undergoing carcinogenesis by methylcholanthrene,” Cancer Res., vol. 6,
[24], and for the 20–50-GHz range were extrapolated from the pp. 574–577, 1946.
[7] K. G. Stern and R. Willheim, The Biochemistry of Malignant Tumors.
measurement data available for frequencies up to 20 GHz in Brooklyn, NY, USA: Reference Press, 1943.
[24] using the Cole–Cole parametric model. [8] L. Mulay et al., “Trace-metal analysis of cancerous and non-cancerous
human tissues,” J. Nat. Cancer Inst., vol. 47, pp. 1–13, 1971.
V. CONCLUSION AND FUTURE DIRECTIONS [9] E. Fear et al., “Enhancing breast tumor detection with near-field imaging,”
IEEE Microw. Mag., vol. 3, no. 1, pp. 48–56, Mar. 2002.
This paper reports on the ultra-wideband (0.5–50 GHz) di- [10] E. Porter et al., “An early clinical study of time-domain microwave radar
for breast health monitoring,” IEEE Trans. Biomed. Eng., vol. 63, no. 3,
electric properties of 101 freshly-excised normal and malignant pp. 530–539, Mar. 2016.
skin tissue samples (resulting in 152 datasets) obtained from skin [11] S. Semenov, “Microwave tomography: Review of the progress towards
cancer patients having undergone Mohs micrographic surgeries clinical applications,” Philos. Trans. A Math. Phys. Eng. Sci., vol. 367,
pp. 3021–3042, 2009.
at Hackensack University Medical Center. A one-pole Cole– [12] A. Zamani et al., “Lung cancer detection using frequency-domain mi-
Cole model is fitted to the complex permittivity datasets of each crowave imaging,” Electron. Lett., vol. 51, pp. 740–741, 2015.
sample, followed by applying a physics-based criterion to ex- [13] S. Alekseev et al., “Millimeter wave dosimetry of human skin,” Bioelec-
tromagnetics, vol. 29, pp. 65–70, 2007.
clude the datasets that cannot be fitted by the model from further [14] J. Rosai and L. Ackerman, Surgical Pathology. New York, NY, USA:
analysis. Mosby, 2004.
A multivariate analysis of variance (MANOVA) is employed [15] A. Uong and L. Zon, “Melanocytes in development and cancer,” J. Cell.
Physiol., vol. 222, pp. 38–41, 2010.
to evaluate the effects of malignancy on the dielectric properties [16] C. S. Joseph et al., “Continuous wave terahertz transmission imaging of
of the samples. Our analyses indicate that statistically significant nonmelanoma skin cancers,” Lasers Surg. Med., vol. 43, pp. 457–462,
differences exist between the dielectric properties of both BCC 2011.
[17] A. Taeb et al., “Millimetre-wave waveguide reflectometers for early de-
and SCC skin cancer tissues and their corresponding normal tection of skin cancer,” IET Microw. Antenna, vol. 7, pp. 1182–1186,
skin tissues and the difference is larger for the case of SCC 2013.
skin cancer. We also provide evidence and discussions into the [18] Hackensack Medical Center. (2017). Dermatologic Surgery Unit.
[Online]. Available: http://www.hackensackumc.org/our-services/surgic
obtained results. We suggest that in addition to a change in al-services/plastic-reconstructive-surgery/mohs-surgery-dermatology-
water content, which is the main cause for the case of BCC skin surgery
Authorized licensed use limited to: University of Applied Sciences Karlsruhe c/o KIT Library. Downloaded on December 23,2021 at 10:41:29 UTC from IEEE Xplore. Restrictions apply.
MIRBEIK-SABZEVARI et al.: CHARACTERISTICS OF FRESHLY EXCISED NORMAL AND MALIGNANT OF HUMAN SKIN TISSUES 1329
[19] Keysight Technologies, Inc. (2017). N1501A Dielectric Probe Kit. [31] B. Alberts et al., Molecular Biology of the Cell. New York, NY, USA:
[Online]. Available: http://www.keysight.com/en/pd-2492144-pn- Garland, 2002.
N1501A/dielectric-probe-kit?cc=US&lc=eng [32] F. A. Duck, Physical Properties of Tissue. Bath, U.K.: Academic, 1992.
[20] A. Mirbeik-Sabzevari and N. Tavassolian, “Characterization and valida- [33] Z. Wang et al., “Hydration of fat-free body mass: Review and critique
tion of the slim-form open-ended coaxial probe for the dielectric char- of a classic body-composition constant,” Amer. J. Clin. Nutr., vol. 69,
acterization of biological tissues at millimeter-wave frequencies,” IEEE pp. 833–841, 1992.
Microw. Compon. Lett., to be published. [34] O. F. Nielsen et al., “Optical devices and diagnostic in materials science,”
[21] S. S. Stuchly et al., “A new aperture admittance model for open-ended Proc. SPIE, vol. 4098, pp. 160–168, 2000.
waveguides,” IEEE Trans. Microw. Theory Technol., vol. 42, pp. 192–198, [35] K. R. Foster and H. P. Schwan, “Dielectric properties of tissues and bi-
1994. ological materials: A critical review.” Crit. Rev. Biomed. Eng., vol. 17,
[22] Stevens Institute of Technology. (2017). The Office of Sponsored Pro- pp. 25–104, 1989.
grams (OSP). [Online]. Available: https://www.stevens.edu/research- [36] M. Gniadecka et al., “Water content and structure in malignant and benign
entrepreneurship/faculty-research/sponsored-programs skin tumours,” J. Mol. Struct., vol. 661, pp. 405–410, 2003.
[23] K. S. Cole and R. H. Cole, “Dispersion and absorption in dielectrics—I. [37] M. Gniadecka et al., “Melanoma diagnosis by Raman spectroscopy and
Alternating current characteristics,” J. Chem. Phys., vol. 9, pp. 341–352, neural networks: Structure alterations in proteins and lipids in intact cancer
1941. tissue,” J. Invest. Dermatol., vol. 122, pp. 443–449, 2004.
[24] S. Gabriel et al., “The dielectric properties of biological tissues: III. Para- [38] S. R. Smith and K. R. Foster, “Dielectric properties of low-water-content
metric models for the dielectric spectrum of tissues,” Phys. Med. Biol, tissues,” Phys. Med. Biol., vol. 30, pp. 965–973, 1985.
vol. 41, pp. 2271–2293, 1996. [39] R. Pethig, Dielectric and Electronic Properties of Biological Materials.
[25] M. Lazebnik et al., “Ultra-wideband temperature-dependent dielectric New York, NY, USA: Wiley, 1979.
properties of animal liver tissue in the microwave frequency range,” Phys. [40] W. J. Ellison, “Permittivity of pure water, at standard atmospheric pressure,
Med. Biol., vol. 51, pp. 1941–1955, 2006. over the frequency range 0–25 THz and the Temperature Range 0–100 °C,”
[26] The MathWorks, Inc. MATLAB. 2017. [Online]. Available: https://www. J. Phys. Chem. Ref. Data, vol. 36, pp. 1–8, 2007.
mathworks.com/products/matlab.html [41] R. de M. F. Pereira et al., “Diagnosis of squamous cell carcinoma of human
[27] A. Peyman et al., “Variation in dielectric properties due to pathological skin by Raman spectroscopy,” Proc. SPIE. vol. 5326, pp. 106–112, 2004.
changes in human liver,” Bioelectromagnetics, vol. 36, pp. 603–612, 2015. [42] H. Lui et al., “Real-time Raman spectroscopy for in vivo skin cancer
[28] B. S. Everitt and G. Dunn, Applied Multivariate Data Analysis. London, diagnosis,” Cancer Res., vol. 72, pp. 2491–2500, 2012.
U.K.: Arnold, 2000. [43] A. C. Lieber et al., “Raman microspectroscopy for skin cancer detection
[29] H. P. Schwan, “Electrical properties of tissue and cell suspensions,” in in vitro,” J. Biomed. Opt., vol. 13, pp. 0240131–0240138, 2008.
Proc. IEEE Conf. Eng. Med. Biol. Soc., 1994, pp. A70–A71. [44] M. Zhadobov et al., “Complex permittivity of representative biological
[30] P. M. Meaney et al., “A coaxial dielectric probe technique for distin- solutions in the 2–67 GHz range,” Bioelectromagnetics, vol. 33, pp. 346–
guishing tooth enamel from dental resin,” Adv. Biomed. Eng. Med. Phys., 355, 2012.
vol. 3, pp. 8–17, 2015.
Authorized licensed use limited to: University of Applied Sciences Karlsruhe c/o KIT Library. Downloaded on December 23,2021 at 10:41:29 UTC from IEEE Xplore. Restrictions apply.