Next Article in Journal
An Automatic Detection and Classification System of Five Stages for Hypertensive Retinopathy Using Semantic and Instance Segmentation in DenseNet Architecture
Next Article in Special Issue
Monitoring Bone Density Using Microwave Tomography of Human Legs: A Numerical Feasibility Study
Previous Article in Journal
Prediction of In Vivo Laser-Induced Thermal Damage with Hyperspectral Imaging Using Deep Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Machine Learning to Predict Dielectric Properties of In Vivo Biological Tissue

by
Branislav Gerazov
1,
Daphne Anne Caligari Conti
2,
Laura Farina
3,
Lourdes Farrugia
2,
Charles V. Sammut
2,
Pierre Schembri Wismayer
4 and
Raquel C. Conceição
5,*
1
Faculty of Electrical Engineering and Information Technologies, Ss Cyril and Methodius University in Skopje, 1000 Skopje, North Macedonia
2
Department of Physics, University of Malta, MSD 2080 Msida, Malta
3
Translational Medical Device Lab, National University of Ireland Galway, H91 TK33 Galway, Ireland
4
Department of Anatomy, University of Malta, MSD 2080 Msida, Malta
5
Instituto de Biofísica e Engenharia Biomédica, Faculdade de Ciências da Universidade de Lisboa, 1749-016 Lisboa, Portugal
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(20), 6935; https://doi.org/10.3390/s21206935
Submission received: 24 September 2021 / Revised: 13 October 2021 / Accepted: 16 October 2021 / Published: 19 October 2021

Abstract

:
In this paper we revisited a database with measurements of the dielectric properties of rat muscles. Measurements were performed both in vivo and ex vivo; the latter were performed in tissues with varying levels of hydration. Dielectric property measurements were performed with an open-ended coaxial probe between the frequencies of 500 MHz and 50 GHz at a room temperature of 25 °C. In vivo dielectric properties are more valuable for creating realistic electromagnetic models of biological tissue, but these are more difficult to measure and scarcer in the literature. In this paper, we used machine learning models to predict the in vivo dielectric properties of rat muscle from ex vivo dielectric property measurements for varying levels of hydration. We observed promising results that suggest that our model can make a fair estimation of in vivo properties from ex vivo properties.

Graphical Abstract

1. Introduction

A significant body of fundamental research has been conducted by researchers to measure the dielectric properties of biological tissue (e.g., [1,2,3]). Such studies are of extreme importance to the correct modelling of the electromagnetic behaviour of biological systems, which can be used in a number of applications. Traditionally, the dielectric properties of biological tissues have been studied to address requirements defined by dosimetry studies and often this determined the methodology used. Recently, there has been a great interest in accurate knowledge of the dielectric properties by the microwave medical community [4,5,6], as these facilitate the correct design and development of medical electromagnetic imaging devices [7,8].
Many historical studies have reported dielectric data obtained using different protocols for measuring the dielectric properties of biological tissues. Many studies reported data measured under specified different conditions—such as temperature, frequency or time from excision—(or have failed to report the conditions of experiments) and have also reported different techniques for measuring dielectric properties such as coaxial probe, transmission line, free space, resonant cavity, parallel plate and inductance measurement method [9].
An international networking project, COST Action (MyWAVE, CA17115), has begun with the objective of reviewing current protocols for measuring biological dielectric and thermal properties and converging toward a standard protocol that will allow for data collected across sites to be comparable, eventually merging into a single database. One of the groups involved in this international project has proposed the use of the “Minimum Information Model for Dielectric Measurements of Biological Tissues (MINDER)”, which provides “a framework for recording and storing dielectric data and metadata for measurements conducted on biological tissues”. By ensuring that common metadata is collected across sites, then it is ensured that “the generation of data is repeatable, interoperable, and re-usable” [10,11].
One recurrent discussion within this research community is how to correlate ex vivo with in vivo dielectric properties. Ex vivo properties are easier to measure, whereas in vivo properties are preferred to create realistic electromagnetic models of biological tissue. Additionally, several works in the literature dealt with measurements of tissue dielectric properties in ex vivo conditions; nevertheless, it is important to investigate whether the irreversible variations, occurring in biological tissues after their excision from the host, affect the dielectric properties. Ex vivo measurements on specimens from surgical resections are subject to confounding factors and depend on the frequency range under investigation. At microwave frequencies, the dielectric response is often referred to as the γ-dispersion, which is mainly due to the interaction with water content [12]. Therefore, the variability of hydration in the sample under test significantly influences the resulting measured dielectric properties. It has been shown that in vivo properties can be estimated by ex vivo measurements given that the variation in hydration between the two samples does not exceed 10% [13,14]. Many studies on dielectric properties emphasise the importance of tissue hydration, but very few studies quantify these variations [13,15]. One of the first studies that quantified water content goes back to 1992 by Campbell and Land [16], reporting no sufficient data for comparison. Moreover, they used a resonant technique that requires extensive sample preparation which can compromise the accuracy of the water content evaluated in the sample under test.
Dielectric properties are easier to measure ex vivo, but the in vivo properties represent a more realistic scenario with more accurate measurements of the electromagnetic properties of the biological tissue, as required for the design of microwave medical devices and applications. Thus, it is key to understanding the correlation between the ex vivo and in vivo properties in order to use the former to predict the latter. In this paper, we show preliminary results that demonstrate which information from the recorded metadata can help predict and infer the in vivo dielectric properties from available ex vivo dielectric properties.

2. Materials

This work investigated the ex vivo dielectric properties of biological tissue samples (i.e., rat muscle) and to studied how the variation in tissue hydration may help predict in vivo dielectric properties.
Spectroscopy measurements were conducted using an Agilent Slim Form 85,070 probe connected to a R&S ZVA-50 Vector Network Analyzer (VNA) in a frequency range between 500 MHz and 50 GHz in a linear scale. The VNA measured the complex reflection coefficient, S11, which was converted to complex permittivity, ε*, using the Agilent software 85070E (Agilent, 2012).
Each tissue sample underwent a dehydration process using a loss-on drying method. The dielectric properties of each tissue sample were measured in vivo in situ, when freshly excised and after dehydration, and then correlated with the percentage loss of water content. Various techniques were studied to determine the water content of the samples: loss-on-drying (LOD) method, centrifugation and freezing, with LOD being the most reliable method for biological tissues and, thus, used for the purpose of this study. In the meantime, intermediate measurements were made, obtaining a set of dielectric data and their corresponding water content. Levels of dehydration between 10% and 70% were achieved. The dielectric measurement data of biological tissues in vivo and ex vivo used in this study are reported in [13] which outline the experimental protocol adopted for the measurement collection. We would like to emphasise that the in vivo and ex vivo samples from each animal were the same.

3. Data Processing and Discussion of Results

We modelled the ex vivo dielectric data in several machine learning models. The database comprised a large number of observations (each measurement of complex permittivity) which were characterised by a number of input features, for example, the frequencies at which the measurements were recorded and the level of hydration for a tissue. We knew the ground truth of the in vivo measurement of complex permittivity for each ex vivo observation. We used a nested cross-validation methodology to train and test our model. In the outer cross-validation loop, we took each rat as a test set and conducted an inner cross-validation loop with the rest of the rats. In the inner loop, we performed a randomised search of the model’s hyperparameters using a leave-one-out group (rat) cross-validation. The best performing hyperparameters were then used to train a model in the whole inner loop set, which was then evaluated on the test set in the outer loop. Finally, the score for each model was calculated as an average score across the outer loop iterations. We present the regression calculation of the in vivo complex permittivity of the validation set and compare these with the ground truth measured results.
We completed the following five tasks:
  • Task 1—load the measurement data set in Python;
  • Task 2—explore and visualise the data;
  • Task 3—structure the data for regression model training and evaluation;
  • Task 4—choose a set of machine learning regression models to evaluate;
  • Task 5—train and evaluate the models with the data in a nested cross validation loop tuning their hyperparameters.
The data set comprised measurements of the relative permittivity ε r of rat muscle tissue given by:
ε r = ε ε 0 ,
ε r = ε r i ε r ,  
where ε* is the actual permittivity of the tissue, and ε0 is the permittivity of vacuum. ε r and ε r are the real and imaginary parts of the relative permittivity, respectively. Relative permittivity was measured in the range of 101 frequency data points from 500 MHz to 50 GHz, with a step of 495 MHz (101 frequency measurements). Both in vivo and ex vivo properties were measured.
For ex vivo properties, permittivity was tracked across several stages of dehydration which were thermally induced and determined by the percentage of weight loss. Four rats were used for multiple measurements, with the ex vivo dehydration levels varying between the rats. For each rat, two thigh muscles were measured in a varying number of randomly chosen locations to obtain a good estimate of the average permittivity of the tissue. In total, 3535 in vivo and 33,229 ex vivo data points at different frequencies, for a total of eight muscle samples, at every dehydration state (in ex vivo), were available in the data set. In the data set we used for this analysis, we discarded some of the data due to the noted calibration problems with the measurement device. The final spread of the data used across the different rats is given in Table 1. The total number of measured data points used was 25,149, taken from 249 measurements (each at 101 frequencies).

3.1. Exploratory Data Analysis

3.1.1. Data Loading, Cleaning and Restructuring

In the first step of our exploration, we loaded the measured data points, originally available across a number of Excel spreadsheets into Pandas (Python Data Analysis Library https://pandas.pydata.org/, accessed on 23 September 2021) dataframes. To do so, we first re-organised the column labels in the spreadsheets. We then loaded the data into Python and eliminated the measurements with known calibration issues. Finally, we restructured all the data into a single master dataframe containing one data point per row, i.e., 25,149 rows.

3.1.2. Number of Measurements and Dehydration Levels

To visualise the spread of the dehydration levels and number of measurements per rat/level, we plotted the scatter plot shown in Figure 1. The dehydration level shown below 0% encodes in vivo measurements. The colour codes for the four rats are red for rat 3, blue for rat 4, green for rat 5 and purple for rat 6. In vivo muscle measurements are represented in coloured squares (there is no muscle designator for the in vivo measurements), whereas ex vivo measurements of the measured muscles 1 and 2 are represented in circles and triangles, respectively. Scatter points were plotted with an added offset to avoid their overlap. We can observe that the number of measurements went from four measurements for each rat–muscle combination of ex vivo without dehydration, up to 11 measurements for rat 3 muscle 1. The average number of measurements per rat–muscle combination was 5.8. We can further observe that the two muscles of rat 3 were measured at three dehydration levels, as were those for rat 5; the muscles for rat 4 were measured at five dehydration levels, and the muscles for rat 6 were measured at eight dehydration levels.

3.1.3. Visualising the Measurements Per Frequencies

To visualise the measurements, we plotted the real and imaginary part of the εr ( ε r and ε r , respectively) for a set frequency as shown in Figure 2. We can observe that the data did capture the general dynamics of permittivity change for increasing dehydration levels. We can also observe that the variability in the measurements also increased with the increase in dehydration. Moreover, we can note that the inter-rat spread of the measurements was rather consistent, although there were some inconsistencies, e.g., measurements at the ~30% dehydration level were only recorded for rats 4 and 6. Finally, we can note that there were no measurements made between the levels of 40% and 50% dehydration for any of the rats. One can also infer the technical challenges behind ensuring specific dehydration levels at higher percentages, as it was not possible to measure at the same dehydration level: ~54% was the highest dehydration possible in rat 3, ~59% and 64% in rat 6, ~62% in rat 4 and ~64% and ~68% in rat 5. In addition, the measurements at higher dehydration levels varied more compared to measurements at lower dehydration levels.

3.1.4. Visualising the Measurements for all Frequencies in 3D

To plot ε r and ε r across all frequencies using a 3D plot, we calculated the mean and standard deviation of both real and imaginary parts, respectively, for each dehydration level, grouping the measurements from different rats and muscles. The resulting averages of ε′ and ε″ are shown in the two surface plots in Figure 3. We can observe that ε’ (the real part) had a more monotonous change with respect to both the frequency and dehydration levels. It decreased steadily for higher frequency and dehydration levels.
The changes in the imaginary part were more complex. Its value started high and dropped to a “valley” for low frequencies (approximately 3 GHz), which was followed by a “plateau” in the midrange of the measured frequencies (approximately 24 GHz) and a final drop for the highest measured frequencies (approximately 40 GHz). This irregularity decreased with the rise in the dehydration level, and the curve flattened for the highest dehydration levels (above 50%). We can also note that for the lowest frequencies, the imaginary part of the permittivity increased with the dehydration level.

3.1.5. Visualising the Measurements for all Frequencies in 2D

We also plotted the average ε r and ε r for all frequencies in a 2D plot, inspired by how spectrograms are represented. The results are shown in Figure 4. To obtain this plot, we first binned all the dehydration levels to bins that were multiples of 10. Since there were no data for the 40% dehydration bin, we then used linear interpolation to obtain the values for it. These plots show the trends we observed in the 3D plots of the data more clearly.

3.2. Regression Problem

3.2.1. Restructuring with Average Nonuniform Samples

To formulate the task of estimating the in vivo permittivity from the measured ex vivo data, we chose to structure the data so that we grouped it by rat and averaged the measurements per dehydration level. Each average measurement is a sample with the target being the in vivo data for the rat measured. Since the targets are estimates of a continuous variable, it can be posed as a machine learning regression problem, in contrast to a classification problem in which the model is trained to estimate a target’s class. In our data re-organisation, we ignored the muscle identifier for the ex vivo measurements, as there was none for the in vivo ones.
Thus, we constructed a six-dimensional feature vector that comprised frequency, dehydration level, and the mean and standard deviation of both ε′ and ε″, extracted from the ex vivo data. The target was, in turn, two-dimensional comprising only the means of the ε′ and ε″ of the in vivo data. This restructuring of the data amounted to a total of 3634 observations.
Even though the data showed that the dynamics were influenced by the dehydration level, the measurement series made for different dehydration levels per rat cannot be directly used in most regression models because they are nonuniformly sampled. Namely, the permittivity was not measured at regular dehydration points across the different rats. Thus, we depended on the regression models to implicitly learn the influence of dehydration of the target values based on the set of training samples provided.

3.2.2. Data Imputation

One way to deal with nonuniformly sampled data is through the use of data imputation, i.e., resampling the data at regular sample points using interpolated values from the original sampling of the data. To assess whether this approach was applicable to our data set, we took the mean data for rat 3 at a single frequency and used various interpolation and extrapolation methods to impute points over the whole range of dehydration levels. In particular, we used a spline interpolator with the following orders: 0 (zero), 1 (linear) and 2 (quadratic) as well as an interpolator imputing the value of the nearest known point (nearest).
The results of this experiment are shown in Figure 5. The plots suggest that the usefulness of interpolation is limited by the access to only a handful of measurements to calculate the means. Unfortunately, due to the experimental nature of this study, we cannot explore imputing further. We do, however, publish these results so that other authors designing a measurement campaign can address this issue.

3.3. Chosen Regression Models

We based our analysis on the regression models implemented in the powerful machine learning package scikit-learn (Machine Learning in Python https://scikit-learn.org/, 13 October 2021). Indeed, scikit-learn provides a large variety of models from which to choose. From them, we chose to use models based on four different categories: linear regression, kernels, trees and tree ensembles, and neural networks. In total 11 models were selected, trained and evaluated.

3.3.1. Linear Regression-Based Models

In our evaluation we incorporated:
  • Linear regression—the ordinary least squares method evaluated for different degrees of the polynomial feature generator;
  • Ridge regression—linear least squares with L2 regularisation of the weight parameters;
  • Lasso regression—linear least squares with L1 regularisation of the weight parameters;
  • Elastic net regression—linear least squares with combined L1 and L2 regularisation of the weight parameters.

3.3.2. Kernel-Based Models

From the available kernel-based models we evaluated:
  • Epsilon-support vector regression (SVR)—an extension of the support vector machines used for classification, the SVR supports various kernels from which we used the radial basis function (RBF) kernel;
  • Kernel ridge regression (KRR)—combines ridge regression with the kernel trick to obtain a model identical in form to SVR but trained with a different loss function. Similarly, we again used the RBF kernel.

3.3.3. Tree and Ensemble of Tree-Based Models

From the tree and ensemble of trees-based regressors we selected the:
  • Decision tree regressor (regtree)—similarly to the decision tree classifier but outputs an estimate of the continuous target value from each leaf; splitting decisions are made in each node based on the best split among all the features;
  • Random forest regressor (randforest)—an ensemble method based on averaging the output of a collection of decision tree regressors, each trained on a random subset of the data set; splitting is done with the best split among a random subset of all the features;
  • Extremely randomised trees regressor (extratree)—an ensemble method based on averaging the output of a collection of decision tree regressors, each trained on a random subset of the data set; splitting is based on the best randomly chosen threshold for a random subset of features;
  • Gradient boosted regression trees (gradboost)—an additive model made from a cascade of weak tree regressors, each trained on the previous one’s output.

3.3.4. Neural Network-Based Models

From the neural network-based models, we used the multilayer perceptron (MLP) regressor. Since scikit-learn allows for building models with multiple hidden layers; we also trained and evaluated deep neural network (DNN)-based regressors.

3.4. Model Training and Evaluation

3.4.1. Feature Augmentation

To increase the number of features available to the linear regression-based classifiers, we added polynomial feature combinations up to the 3rd degree. For example, a 2nd order polynomial feature generator that receives as input a two-dimensional feature vector (a, b), would output the six-dimensional feature vector (1, a, b, a2, ab, b2). This feature augmentation allows linear models to learn a polynomial data manifold.

3.4.2. Data Normalisation

Since some regression models work better when the input features have similar scales, we included a standard scaler in the cross-validation pipeline. The standard scaler removes the mean and scales the features to unit variance as in:
x ^ = x μ σ ,  
where μ is the mean value of feature x, and σ is its standard deviation.
An important caveat is that the data must be scaled using the mean and standard deviation of the training set exclusively. Doing otherwise would leak information from the validation or test sets into the model training.

3.4.3. Hyperparameter Tuning

Random hyperparameter tuning was carried out for each of the selected regression models. For each hyperparameter, we defined a valid range of values and a random sampling function to sample from this range. Most hyperparameters were sampled using an exponential distribution, and some were sampled using a uniform distribution. Finally, some hyperparameter settings were limited to a list of values.

3.4.4. Nested Cross-Validation Loop

Each of the models were trained and evaluated using a “leave one group out” nested cross-validation loop in which the data were grouped according to the rat identifier. In each iteration of the outer loop, one rat was selected to be the test rat. In the inner loop, each of the three remaining rats were taken out to be the validation rat in each iteration. Within each inner loop iteration, the model was trained with various hyperparameters and evaluated on the validation rat, and the hyperparameters that resulted in the best validation score were used to retrain the model for the three rats, i.e., including the validation rat. The newly trained model was then evaluated on the test rat. After the four rats were selected as test rats and the outer loop was completed, the average test score was calculated for the four rats.
The best results were selected based on the R2 coefficient of determination regression score measure defined as:
R 2 = 1 i = 0 N 1 ( y i y ^ i ) 2 i = 0 N 1 ( y i μ ) 2 ,
where N is the number of samples, y is the target value, μ is its average, and y is the target estimate output by the model. The best R2 score is a value close to 1, but the score can be negative. The advantage over using the mean square error (MSE), which is partly in the nominator in the above equation:
MSE = 1 N i = 0 N 1 ( y i y ^ i ) 2 ,  
Is that a model that always predicts the mean target value, disregarding the input features, would have a R2 score of 0.

3.4.5. Regression Models Results

The results obtained using the selected regression models on predicting the real and imaginary part of the in vivo permittivity, ε′ and ε″, based on its ex vivo measurements at various levels of dehydration, are shown in Table 2. The linear regression model was included to have a non-machine learning algorithm to predict in vivo dielectric properties. As expected, due to the pronounced non-linearity of the data, this algorithm performed the worst compared to the other machine learning algorithms. We can observe that most models do well on the set task, with the four best performing models being: extremely randomised trees, random forests, gradient boosted trees and regression trees, all scoring below 1.2 average MSE. From the remaining models, ridge regression is the next in line in terms of performance.

3.5. Most Important Features

To estimate the most important features for estimating the in vivo target values, we took the best performing model—the extremely randomised trees—and repeated the nested cross-validation loop for 1000 iterations, aggregating the feature importance for each of the four models trained within one loop. Because of its stochastic nature, this process resulted with some 4000 trained models and as much feature importance values for each of the features. These are shown in the box plot in Figure 6. We can see that the two most important features are the frequency of the measured value and the mean of ε r .

4. Conclusions

The presented results show that estimating the in vivo permittivity of muscle tissue in rats, based on measurements of its ex vivo permittivity is feasible, and current state-of-the-art regression models provide promising results. The models best suited for the task are ensemble tree-based methods, especially the extremely randomised trees regressor, which yielded the lowest MSE and highest R2 metrics. This study is an important first step towards the development of tools that can predict in vivo dielectric properties from widely available ex vivo dielectric properties.
In future work, we intend to include more metadata in our machine learning model to better estimate the in vivo dielectric properties of biological tissues as well as extending this study to other types of tissues other than rat muscle. Moreover, we will attempt to use a state-of-the-art model that can deal with nonuniformly sampled data, which is the recurrent neural network (RNN). In order to use this data with an RNN-based regressor, the data will have to be restructured in a dehydration series format.

Author Contributions

The experimental data collection study was designed by L.F. (Lourdes Farrugia) and C.V.S., while the data was collected by D.A.C.C. and L.F. (Laura Farina). All animals were handled by P.S.W. The data analysis study in this paper was conceptualised by B.G. and R.C.C. and executed by B.G. The original draft was prepared by B.G. and R.C.C. Review and editing were completed by all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Fundação para a Ciência e a Tecnologia-FCT, FCT/MEC (PIDDAC) under the Strategic Programme UIDB/00645/2020, COST Action TD1301—MiMed with an STSM grant entitled “Characterization of the Dielectric Properties of Biological Tissues and Their Correlation to Tissue Hydration” (COST-STSM-TD1301-100416-071441) and COST Action CA17115 MyWAVE with an STSM grant entitled “Application of Machine Learning to Predict Dielectric Properties of in vivo Biological Tissue” (ECOST-STSM-CA17115-080419-107078).

Institutional Review Board Statement

This paper features secondary analysis of dielectric data previously published in “Characterization of the Dielectric Properties of Biological Tissues and Their Correlation to Tissue Hydration” (https://doi.org/10.1088/2057-1976/aafc1a) [13], the authors of which are also co-authors of this paper. A self-assessment form for this study was submitted for records to the Faculty of Science Research Ethics Committee—University of Malta Form ID. 15062020 9735.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Joines, W.T. Frequency-Dependent Absorption of Electromagnetic Energy in Biological Tissue. IEEE Trans. Biomed. Eng. 1984, BME-31, 17–20. [Google Scholar] [CrossRef] [PubMed]
  2. Campbell, A.M. Measurements and Analysis of the Microwave Dielectric Properties of Tissues. Ph.D. Thesis, University of Glasgow, Glasgow, UK, 1990. [Google Scholar]
  3. Gabriel, S.; Lau, R.W.; Gabriel, C. The Dielectric Properties of Biological Tissues: III. Parametric Models for the Dielectric Spectrum of Tissues. Phys. Med. Biol. 1996, 41, 2271–2293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Gabriel, S.; Lau, R.W.; Gabriel, C. The Dielectric Properties of Biological Tissues: II. Measurements in the Frequency Range 10 Hz to 20 GHz. Phys. Med. Biol. 1996, 41, 2251–2269. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Lazebnik, M.; Popovic, D.; McCartney, L.; Watkins, C.B.; Lindstrom, M.J.; Harter, J.; Sewall, S.; Ogilvie, T.; Magliocco, A.; Breslin, T.M.; et al. A Large-Scale Study of the Ultrawideband Microwave Dielectric Properties of Normal, Benign and Malignant Breast Tissues Obtained from Cancer Surgeries. Phys. Med. Biol. 2007, 52, 6093–6115. [Google Scholar] [CrossRef] [PubMed]
  6. Taflove, A.; Hagness, S.C. Computational Electrodynamics: The Finite-Difference Time-Domain Method, 2nd ed.; Artech House: Boston, MA, USA, 2000. [Google Scholar]
  7. Amin, B.; Elahi, M.A.; Shahzad, A.; Porter, E.; O’Halloran, M. A review of the dielectric properties of the bone for low frequency medical technologies. Biomed. Phys. Eng. Express 2019, 5, 022001. [Google Scholar] [CrossRef]
  8. Gioia, A.L.; Porter, E.; Merunka, I.; Shahzad, A.; Salahuddin, S.; Jones, M.; O’Halloran, M. Open-Ended Coaxial Probe Technique for Dielectric Measurement of Biological Tissues: Challenges and Common Practices. Diagnostics 2018, 8, 40. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Agilent Basics of Measuring the Dielectric Properties of Materials. Available online: www.agilent.com (accessed on 8 October 2021).
  10. Minder—Minimum Information Model for Dielectric Measurements of Biological Tissues. Available online: https://www.bio-minder.com/ (accessed on 9 October 2021).
  11. Porter, E.; Gioia, A.L.; Salahuddin, S.; Decker, S.; Shahzad, A.; Elahi, M.A.; O’Halloran, M.; Beyan, O. Minimum information for dielectric measurements of biological tissues (MINDER): A framework for repeatable and reusable data. Int. J. RF Microw. Comput. Aided Eng. 2017, 28, 1–17. [Google Scholar] [CrossRef] [Green Version]
  12. Schwan, H.P. Electrical properties of tissue and cell suspensions. Adv. Biol. Med. Phys. 1957, 5, 147–209. [Google Scholar] [PubMed]
  13. Pollacco, D.A.; Farina, L.; Sammut, C.V.; Farrugia, L.; Schembri-Wismayer, P. Characterization of the dielectric properties of biological tissues and their correlation to tissue hydration. IEEE TDEI 2018, 25, 2191–2197. [Google Scholar] [CrossRef]
  14. Farrugia, L.; Wismayer, P.S.; Mangion, L.Z.; Sammut, C.V. Accurate in vivo dielectric properties of liver from 500 MHz to 40 GHz and their correlation to ex vivo measurements. Electromagnet. Biol. Med. 2016, 35, 365–373. [Google Scholar] [CrossRef] [PubMed]
  15. Meo, S.D.; Bonello, J.; Farhat, I.; Farrugia, L.; Pasian, M.; Podesta, M.T.; Camilleri Podesta, M.T.; Suleiman, S.; Calleja-Agius, J.; Sammut, C.V. The variability of dielectric permittivity of biological tissues with water content. JEMWA 2021, 31, 1–21. [Google Scholar]
  16. Campbell, A.M.; Land, D.V. Dielectric Properties of Female Human Breast Tissue Measured In Vitro at 3.2 GHz. Phys. Med. Biol. 1992, 37, 193–210. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Spread of dehydration levels and number of measurements per level per rat. The in vivo measurements are shown below the ex vivo measurements. The colour codes for the four rats are red, blue, green and purple for rats labelled 3, 4, 5 and 6, respectively. In vivo muscle measurements are represented in coloured squares (there is no muscle designator for the in vivo measurements), whereas ex vivo measurements of measured muscles 1 and 2 are represented in circles and triangles, respectively. Scatter points were plotted with an added offset to avoid their overlap.
Figure 1. Spread of dehydration levels and number of measurements per level per rat. The in vivo measurements are shown below the ex vivo measurements. The colour codes for the four rats are red, blue, green and purple for rats labelled 3, 4, 5 and 6, respectively. In vivo muscle measurements are represented in coloured squares (there is no muscle designator for the in vivo measurements), whereas ex vivo measurements of measured muscles 1 and 2 are represented in circles and triangles, respectively. Scatter points were plotted with an added offset to avoid their overlap.
Sensors 21 06935 g001
Figure 2. Measurements of the real (top) and imaginary (bottom) part of the permittivity (ε′ and ε″, respectively) for a frequency of 500 MHz. Measurements were performed on 4 rats (labelled as 3, 4, 5 and 6, colours red, blue, green and purple, respectively); in vivo measurements were made on a single muscle (square), and ex vivo measurements were performed on two muscles (1 and 2, circle and triangle shapes, respectively).
Figure 2. Measurements of the real (top) and imaginary (bottom) part of the permittivity (ε′ and ε″, respectively) for a frequency of 500 MHz. Measurements were performed on 4 rats (labelled as 3, 4, 5 and 6, colours red, blue, green and purple, respectively); in vivo measurements were made on a single muscle (square), and ex vivo measurements were performed on two muscles (1 and 2, circle and triangle shapes, respectively).
Sensors 21 06935 g002
Figure 3. Average real (top) and imaginary (bottom) parts of the relative permittivity ( ε r and ε r , respectively) across all frequencies in 3D. Note that the colour bar limits were adjusted for better contrast; the actual ranges are 7–65 for Re( ε r ) and 3–68 for Im( ε r ).
Figure 3. Average real (top) and imaginary (bottom) parts of the relative permittivity ( ε r and ε r , respectively) across all frequencies in 3D. Note that the colour bar limits were adjusted for better contrast; the actual ranges are 7–65 for Re( ε r ) and 3–68 for Im( ε r ).
Sensors 21 06935 g003
Figure 4. Binned and interpolated values of average real (top) and imaginary (bottom) parts of the permittivity ( ε r and ε r , respectively) across all frequencies. Note that the colour bar limits were adjusted for better contrast, the actual ranges are 7–65 for Re( ε r ) and 3–68 for Im( ε r ).
Figure 4. Binned and interpolated values of average real (top) and imaginary (bottom) parts of the permittivity ( ε r and ε r , respectively) across all frequencies. Note that the colour bar limits were adjusted for better contrast, the actual ranges are 7–65 for Re( ε r ) and 3–68 for Im( ε r ).
Sensors 21 06935 g004aSensors 21 06935 g004b
Figure 5. Interpolation permittivity results, average real (top) and imaginary (bottom) parts (ε′ and ε″, respectively), for data imputing for rat 3 at a frequency of 500 MHz for different interpolation methods: nearest neighbour (blue) and spline interpolation with orders of 0 (zero, orange), 1 (linear, green) and 2 (quadratic, red).
Figure 5. Interpolation permittivity results, average real (top) and imaginary (bottom) parts (ε′ and ε″, respectively), for data imputing for rat 3 at a frequency of 500 MHz for different interpolation methods: nearest neighbour (blue) and spline interpolation with orders of 0 (zero, orange), 1 (linear, green) and 2 (quadratic, red).
Sensors 21 06935 g005
Figure 6. Feature importance across 1000 runs of the nested cross-validation loop for the extremely randomised trees regressor.
Figure 6. Feature importance across 1000 runs of the nested cross-validation loop for the extremely randomised trees regressor.
Sensors 21 06935 g006
Table 1. Summary of the number of data points for each rat as measured in vivo, ex vivo (freshly excised) and at different dehydration levels.
Table 1. Summary of the number of data points for each rat as measured in vivo, ex vivo (freshly excised) and at different dehydration levels.
RatIn VivoEx Vivo without DehydrationEx Vivo with DehydrationEx Vivo Dehydration LevelsTotal
rat 3808808414155757
rat 47078085959107474
rat 5808808323264848
rat 6707808555597070
Total3030323218,8873025,149
Table 2. Average metrics obtained with the various regression models: linear regression (linear), ridge regression (ridge), lasso regression (lasso), elastic regression (elastic), support vector machines regressor (SVR), k-nearest neighbour regressor (KNNR), regression tree (regtree), random forest (randforest), extra trees (extratree), gradboost (gradient boost regression) and multilayer perceptron (MLP).
Table 2. Average metrics obtained with the various regression models: linear regression (linear), ridge regression (ridge), lasso regression (lasso), elastic regression (elastic), support vector machines regressor (SVR), k-nearest neighbour regressor (KNNR), regression tree (regtree), random forest (randforest), extra trees (extratree), gradboost (gradient boost regression) and multilayer perceptron (MLP).
ModelNumber of Evaluated HyperparametersMSE
Re(εr)
MSE
Im(εr)
MSE
Average
R2 Score
Re(εr)
R2 Score
Im(εr)
R2 Score
Average
linear10003.033.523.280.9820.6670.825
ridge10001.181.561.370.9930.8550.924
lasso10001.331.591.460.9920.8550.924
elastic10001.821.361.590.9890.8760.933
SVR1001.631.351.490.9910.8770.934
KRR1001.431.461.440.9920.8680.930
regtree10001.250.991.120.9930.9130.953
randforest1001.110.981.050.9940.9130.954
extratree1000.990.940.960.9940.9180.956
gradboost1001.130.931.030.9940.9190.956
MLP1001.631.291.460.9900.8830.937
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gerazov, B.; Caligari Conti, D.A.; Farina, L.; Farrugia, L.; Sammut, C.V.; Schembri Wismayer, P.; Conceição, R.C. Application of Machine Learning to Predict Dielectric Properties of In Vivo Biological Tissue. Sensors 2021, 21, 6935. https://doi.org/10.3390/s21206935

AMA Style

Gerazov B, Caligari Conti DA, Farina L, Farrugia L, Sammut CV, Schembri Wismayer P, Conceição RC. Application of Machine Learning to Predict Dielectric Properties of In Vivo Biological Tissue. Sensors. 2021; 21(20):6935. https://doi.org/10.3390/s21206935

Chicago/Turabian Style

Gerazov, Branislav, Daphne Anne Caligari Conti, Laura Farina, Lourdes Farrugia, Charles V. Sammut, Pierre Schembri Wismayer, and Raquel C. Conceição. 2021. "Application of Machine Learning to Predict Dielectric Properties of In Vivo Biological Tissue" Sensors 21, no. 20: 6935. https://doi.org/10.3390/s21206935

APA Style

Gerazov, B., Caligari Conti, D. A., Farina, L., Farrugia, L., Sammut, C. V., Schembri Wismayer, P., & Conceição, R. C. (2021). Application of Machine Learning to Predict Dielectric Properties of In Vivo Biological Tissue. Sensors, 21(20), 6935. https://doi.org/10.3390/s21206935

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop