Next Article in Journal
Placenta-Targeted Nanoparticles Loaded with PFKFB3 Overexpression Plasmids Enhance Angiogenesis and Placental Function
Previous Article in Journal
SelfCoLearn: Self-Supervised Collaborative Learning for Accelerating Dynamic MR Imaging
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An End-to-End Deep Learning Approach for Quantitative Microwave Breast Imaging in Real-Time Applications

1
Dipartimento di Scienze Motorie e del Benessere, University of Napoli Parthenope, Via Medina 40, 80133 Napoli, Italy
2
Centro Direzionale, Dipartimento di Ingegneria, University of Napoli Parthenope, 80143 Napoli, Italy
*
Author to whom correspondence should be addressed.
Bioengineering 2022, 9(11), 651; https://doi.org/10.3390/bioengineering9110651
Submission received: 12 September 2022 / Revised: 11 October 2022 / Accepted: 28 October 2022 / Published: 4 November 2022

Abstract

:
(1) Background: In this paper, an artificial neural network approach for effective and real-time quantitative microwave breast imaging is proposed. It proposes some numerical analyses for the optimization of the network architecture and the improvement of recovery performance and processing time in the microwave breast imaging framework, which represents a fundamental preliminary step for future diagnostic applications. (2) Methods: The methodological analysis of the proposed approach is based on two main aspects: firstly, the definition and generation of a proper database adopted for the training of the neural networks and, secondly, the design and analysis of different neural network architectures. (3) Results: The methodology was tested in noisy numerical scenarios with different values of SNR showing good robustness against noise. The results seem very promising in comparison with conventional nonlinear inverse scattering approaches from a qualitative as well as a quantitative point of view. (4) Conclusion: The use of quantitative microwave imaging and neural networks can represent a valid alternative to (or completion of) modern conventional medical imaging techniques since it is cheaper, safer, fast, and quantitative, thus suitable to assist medical decisions.

1. Introduction

1.1. Motivation

Breast cancer is one of the commonest types of cancer which affects women, and its early detection is of vital importance for successful treatments [1,2,3]. Among the commonest clinical imaging and diagnostic modalities employed for this aim it is worth mentioning X-ray mammography, magnetic resonance imaging (MRI), ultrasound scanning and nuclear medicine [4].
X-ray mammography is based on the use of very-high-frequency radiations (10 PHz–10 EHz) which travels along the scanned area. Thus, it exploits ionizing radiations and it represents an uncomfortable exam due to compression of the breast to carry out the diagnosis. Some other limitations are related to the limited dynamic range, low contrast and grainy image, for which it is difficult to visualise very subtle lesions in women who have implants or surgical scars [5]. Another disadvantage is related to its poor spatial resolution and the requirement of a large storage place.
Conversely, due to the use of ionizing radiation to perform the investigation, the use of breast MRI is more recommended for women who are high risk, representing a valuable alternative to X-ray mammography since it does not involve any ionizing radiation exposure. This medical exam produces good spatial resolution images but has some limitations related to the low specificity which results in further tests and biopsies, which drives higher costs [6]. Besides being expensive, this imaging modality is not portable, is slow in the acquisition process and is unsuitable for patients with metallic devices implanted.
Breast ultrasound imaging, which is based on the use of ultrasonic mechanical waves for non-invasive diagnostics, is used as a follow-up test for abnormalities found by mammogram and provides some guidelines for biopsy inspections. This medical examination produces both qualitative and quantitative diagnostic information with a good image quality [7]. The main drawbacks are related to the operator-dependent nature of this exam and to the poor resolution of the image and the low contrast.
Complementary information can be inferred from nuclear medicine, whose functional images are based on the molecular properties of the tissues and on the injected radioactive substance. This imaging modality uses ionizing radiation, as in the mammography case, and it is very expensive. Different from the other imaging modalities, it can investigate the physiological function of the system, but has limited resolution and slow imaging time [8].
The drawbacks and limitations of the aforementioned diagnostic methodologies have motivated the research community to develop new imaging techniques and modalities to realise early, reliable and inexpensive diagnostics. In this framework, the use of microwaves for breast imaging and cancer detection has received significant attention, since they might provide better sensitivity and safety due to their non-ionizing nature [9,10].
Microwave tomography exploits microwave signals to investigate breast tissues and provide quantitative permittivity and conductivity maps (the dielectric properties) of the imaged tissues. This imaging technique offers several advantages over other classical methods since it is not expensive, is non-ionizing, and is comfortable with respect to the treatment [11]. The capability of these systems to perform the imaging is related to the fact that different types of tissues, including both normal and malignant ones, have different electrical properties in the microwave frequency range, both in terms of permittivity and conductivity values [12]. Thus, by virtue of these differences, microwave imaging might represent a good candidate to discriminate between healthy and cancer tissues, and the clinical benefits deriving from its use for imaging purposes are supported by several articles in the scientific literature [9,11,13,14,15].
In particular, microwave systems can provide complementary information on the investigated regions which can be merged with the one deriving from traditional exams to support medical decisions. Nevertheless, there are still some technical challenges and open topics to be solved, which involve both hardware and software [11,15]. Among them, it is worth mentioning the challenges related to the choice of an effective coupling medium to maximise the microwave signal transmitted into the body, the design of proper microwave antennas to allow a high enough number of sensors in a small area to investigate the region of interest, and the reduction of the impact of mutual coupling between close antennas and electromagnetic interference with other radio-frequency devices.
Moreover, another main limitation, especially for quantitative microwave imaging, is related to the well-known nonlinearity and ill-posedness issues and their consequences affecting these kinds of inverse problems, such as the presence of local minima [16], and whose difficulty further increases when dealing with three-dimensional modelling, which is computationally demanding. In particular, the impact of the ill-posedness can be partially mitigated by exploiting some a priori information [17,18,19].
It is worth noting that, since the use of microwave imaging in the medical community is relatively recent compared to other state-of-the-art approaches, the authors intended this work to explore the potentialities that the proposed neural-network-based methodology allows us to reach in terms of imaging performance in the context of microwave breast imaging. For these motivations, in this manuscript we did not report a comparison with other conventional imaging modalities nor some recoveries including the malignant-tissue case, which will be properly considered in future work. The main aim of this work is to explore the potential improvement that the proposed methodology allows us to obtain in terms of the quality of the recovery.

1.2. Prototypes

Due to the potentialities of the microwave imaging modality, several research groups developed and realised different prototypes in recent years. In [20], the authors propose a 2D imaging system which consists of a circular array of 12 printed monopole antennas with some recoveries related to simplified phantoms that mimic the human tissues. Further work more focused on the radar system hardware is reported in [21,22,23,24]. Complementary to the previous articles, the effort of the research community was also focused on software aspects, such as finding hybrid nonlinear strategies to improve the recovery performance [25] or the improvement in sensitivity of radar-based breast imaging [26]. All the previously mentioned work contributed to the development of microwave breast imaging. As a matter of fact, the last years have been characterised by several clinical trials [27,28,29,30,31,32]. Nevertheless, these systems all share the feature of employing more antennas arranged around the object of interest with more (virtual or real) transmitters and receivers with the aim of producing an image of the objects located in the investigation region. These approaches are also called tomographic.
Initial experimental setups consisted of a three-dimensional (3D) prototype for microwave tomography with more transmitters and one receiver which surrounded the imaging domain and were controlled by a motion control system [33]. In this case, the object to be imaged was immersed in deionized water (used as matching medium) in a cylindrical chamber. Other systems followed this setup with an increment in the number of employed antennas [34]. After these preliminary microwave imaging experiences, some clinical prototypes were developed and tested, still preserving the cylindrical arrangement of the sensors around the imaging area [35,36], but employing a different kind of antenna immersed in a saline coupling medium to improve signal penetration inside the tissues. After that, several improvements followed in both hardware and software [37,38,39] which yielded to the latest arrangement of the antennas located in a hemispherical shape to make the exam more comfortable on the patient’s side [40,41].
Nowadays, some of the aforementioned imaging systems have started their active clinical trials, showing the potentialities of this technology for supporting medical decisions [42]. Some results of these trials are summarised in Table 1.

1.3. Algorithms

In order to obtain the permittivity and conductivity maps of the imaged tissues, an inversion algorithm can be adopted for solving the electromagnetic inverse scattering (EIS) problem, which is ill-posed and strongly nonlinear in its general formulation [16]. The non-linearity of the considered inverse problem is related to the multiple scattering interactions between points of different tissues and also inside the same tissue, which makes the problem at hand hard to solve and affected by false solutions. These two challenging issues have stimulated different research groups to develop several methods in order to provide tomographic qualitative as well as quantitative images of the biological tissues under test [16,43]. In this framework, the former approaches aim at identifying the pathology within a region without, necessarily, providing an image of the region under test [10,42] while the latter ones aim at producing full maps of the biological tissues in terms of permittivity and conductivity or in terms of labels which are uniquely related to the considered tissues (i.e., segmentation and classification maps).
The possibility of providing quantitative maps of the electromagnetic properties of tissues under test in a non-invasive way can represent a paramount tool for diagnostic purposes. In this context, the capability of tomographic approaches to provide maps of the investigated biological tissues makes them very attractive for clinical applications. Thus, several imaging strategies have been developed over the past thirty years to solve this problem. With regard to the tomographic approaches, three main categories can be identified: qualitative, approximated and quantitative methods [44].
Qualitative methods aim at solving an inverse obstacle problem [45,46,47] by processing the scattered field and providing an estimation of the total tissue extension in the breast, but not its characterisation (i.e., the type of tissue). Approximated methods exploit some approximations of the scattering phenomena in order to allow easy implementation and to keep the computational complexity low. Although they are quite fast, they suffer from some limitations related to the adopted approximated model, as in the case of the well-known Born and Kirchhoff approximations [48]. Another relevant linearised approach proposed in recent years for biomedical imaging purposes is based on the virtual experiments framework [49], which is based on the combination of real and “virtual” experiments that propose a linear approximation of the EIS problem in the case of non-weak scattering regime. Furthermore, higher-order Born approximations can be exploited to reconstruct the conductivity function of the dielectric tissues under test [50].
In order to overcome the limited retrieving performance related to the aforementioned classes, quantitative approaches can be exploited [51]. With regard to these methods, the class of retrievable objects becomes wider at the expense of higher computational burden and processing time. Different iterative approaches can be employed to solve this problem and to face the issues related to non-linearity which may drive into false solutions due to the presence of local minima. Due to the high computational complexity of these methods, some local minimisation approaches are adopted and thus the choice of the initial step is paramount [52]. However, when the non-linearity of the problem at hand is very strong, then it might be beneficial to use global optimization to avoid local minima [44,53]. Unfortunately, their complexity grows exponentially with the number of unknowns and this makes their use very hard to apply for realistic and/or real-time applications. Moreover, due to non-linearity, the inversion procedure may be more sensitive to modelling errors and uncertainties on the scenario.

1.4. Machine Learning for Quantitative Microwave Imaging

In the framework of inverse problems, fast and reliable non-linear approaches are desirable for addressing the imaging problem in the biomedical area of breast cancer diagnostics. Among the most recent methodologies, artificial neural networks represent a useful and flexible tool for quantitative imaging. As a matter of fact, neural networks and artificial intelligence proved to perform well in the field of computer vision, image processing and classification. First, methodologies based on artificial neural networks were applied to extract some general information about the geometric and electromagnetic properties of the scatterers and tissues at hand [54,55]. Most of these first attempts to face the imaging problem via neural networks used a few spatial as well as electromagnetic parameters to represent the scatterers.
Recently, most of the literature has focused on the use of deep convolutional neural networks (CNNs) for solving the inverse problem [56,57,58,59,60,61,62]. Neural networks with regression features have provided very impressive results on EIS problems. The majority of these articles do not propose a direct inversion scheme, i.e., the approach does not allow us to move directly from the data collected at receivers to an estimate of the profile, but they usually perform a super-resolution of the recovery starting from a raw image obtained via other conventional approaches. One of the most adopted techniques consists of the training of a U-net architecture [63] for obtaining quantitative recoveries via preliminary manipulations, e.g., approximated models and a priori information to move from the data (i.e., the scattered field samples) to contrast/induced currents approximations [57,58]. In order to solve the EIS problem with high contrast, a contrast-source-based neural network combined with a traditional subspace-based optimisation method and CNNs might be employed [59]. Furthermore, neural networks can be also employed as regularisation strategies in conventional inversion approaches [64,65] as well as to obtain super-resolved reconstructions [66].
Contrary to the contemporary scientific literature which focuses on the use of CNNs, in this work we focused on artificial neural networks (ANNs) based on multilayer perceptrons. This kind of network allows us to implement a direct inversion scheme from the scattered field samples to directly retrieve a quantitative map of the dielectric features of breast profiles in a fast, efficient way. Despite the ease in network design, a critical issue lies in the choice of a large enough dataset for training the network, which proves to be of vital importance for the estimation of the links’ strength between nodes [63]. Thus, the main bottleneck of this kind of approaches is related to the computational burden required to train the network. Nevertheless, after the initial training, then a direct mapping between data and unknowns can be obtained, producing reliable images in a considerably fast inversion procedure.
Inspired by the universal approximation theorem (UAT) [67]—which states that any arbitrary non-linear function can be approximated via a proper fully-connected neural network with a large number of neurons in its hidden layers under some mild assumptions—in this manuscript we propose an ANN architecture for the real-time quantitative imaging of female breast dielectric properties. The motivation in choosing such an architecture is supported by the need to define a general approach for the retrieval of whatever non-linear profile, as supported by the UAT. As a matter of fact, in an ANN architecture, all the inputs contribute to every single output, resulting in being more suitable for this kind of application than CNNs, the latter being the default choice when dealing with highly-structured modalities such as images or video.
It is worth noting that the proposed work is focused on the potentialities that deep learning via quantitative tomographic imaging at microwave frequencies can deliver in the framework of biomedical breast imaging. Moreover, the use of an end-to-end network which implements a direct inversion scheme further strengthens the proposed approach, allowing the processing of the scattered field samples to provide a quantitative map of the tissues properties.
Therefore, in short, it is possible to summarize the main element of novelty of the considered work into three main aspects:
  • The use of fully-connected neural networks to perform quantitative imaging of the breast tissues;
  • The use of a direct inversion scheme to obtain the permittivity and conductivity maps of breast tissues;
  • The realistic in-house numerical phantom generator and the corresponding dataset for an overall population of 120,000 elements, which is paramount for a proper training of the neural networks to perform a certain task, and therefore an important element of novelty.
On the other hand, the main limitations of the proposed methodology, as with most deep learning approaches, lie in the initial training time necessary for the correct operation of the neural networks and in the database generation which is task-specific. Nevertheless, it is worth noting that, except for the initial training time, the quantitative imaging of the target tissues is very fast (less than 0.01 s on a Linux Mint machine with an AMD Ryzen Threadripper 3990X processor and 250 GB RAM memory).
The outline of the paper is as follows. In Section 2, an overview of the mathematics involved in electromagnetic inverse scattering (EIS) problem is recapped. In Section 3, an overview of the proposed fully-connected ANN approach is reported with a focus on the dataset generation dealing with a realistic breast-like phantom generator for the database population. Finally, some results on synthetic breast phantoms are shown in Section 4 and a comparison with conventional nonlinear EIS approaches is proposed. Some conclusions are drawn at the end of the article.

2. Problem Statement

In the following, a simplified two-dimensional geometry is considered. The background medium is supposed to be homogeneous with complex permittivity ϵ b and with magnetic permeability μ 0 = 4 π · 10 7 H/m. The reason for such an assumption is related to the fact that biological tissues are characterised by a constant value of the magnetic permeability while a certain variability of the complex permittivity can be observed. The antennas are located along a measurement curve Γ which surrounds the imaging domain Ω . The targets located inside this domain are illuminated via transverse-magnetic electric fields generated by z-oriented current wires located on Γ .
A sketch of the geometry is reported in Figure 1. The scattering phenomena depend on the contrast between the dielectric properties of the background medium and those ones of the target. The aim of the proposed method is to provide the breast tissue properties maps by solving an EIS problem, i.e., their relative complex permittivity ϵ r = ϵ r j σ r ω ϵ 0 , with ϵ r and σ r being the relative permittivity and conductivity maps, respectively. Thus, the electromagnetic scattering equations ruling these phenomena can be written as [16]:
E s = A e ϵ , σ , E t + n ,
E t = E i + A i ϵ , σ , E t ,
in which the dependence on the operating frequency and background dielectric features have been implied. E i and E t represent the incident and total electric field inside the imaging domain Ω , respectively, while E s represents scattered field at receivers locations on Γ . The quantities A i and A e are the radiating operators which depend on the dielectric properties of background medium ϵ b , and n is the noise which affects the collected data, here assumed to be additive, white and Gaussian (AWGN).
As regards the imaging applications, the considered framework can be dealt with as an EIS problem consisting of the retrieval of a quantitative estimate of the unknown relative permittivity ϵ and conductivity σ functions inside Ω from the scattered field samples measured on Γ . As previously stated in Section 1.3, such a problem is both nonlinear and ill-posed [43,68,69], thus finding a solution to this problem is not trivial and requires facing a high computational burden and time-consuming approaches. Nevertheless, reliable and fast algorithms able to provide dielectric properties maps of the tissues under test are desirable for early diagnosis in the biomedical field. In this framework, conventional nonlinear approaches can achieve good recovery performance [70,71,72,73,74] at the expense of high computational burden, which might imply no real time applications. Moreover, considerable a priori information and data pre-processing are required to obtain reliable recoveries.
In this framework, artificial neural networks (ANNs) can be of interest and represent a very attractive alternative to real-time applications with more accurate reconstructions. In this case, the challenge is related to the network design and the dataset generation for the training step which considerably impacts on the recovery performance.

3. Methodology

In order to define the proposed approach to solving the inverse problem at hand, two main aspects have to be described: the characteristics of the database adopted for the ANN training and the layout of the implemented neural network.

3.1. Breast Database

In order to perform quantitative inversion via neural networks in a reliable way, it is fundamental to properly model the specific scenario under test and to have large datasets of realistic breast profiles. Some research groups focused on the generation of these datasets to be used by the scientific community [75,76]. Unfortunately, in most cases the available population is very limited or not flexible enough to be exploited in the framework of neural networks training, since they are limited to a specific measurement configuration and/or in a certain frequency range. In order to overcome these limitations and with the aim of providing a useful simulation tool for the testing of breast microwave imaging algorithms, a numerical realistic breast phantom generator has been developed to populate the database. After that, the corresponding scattering matrices are calculated, exploiting a forward solver based on the method of moments (MoM). These matrices, together with the generated breast profiles, represent the data adopted for the proposed neural network approach. Figure 2 illustrates the main steps of the breasts’ phantom generation, which together with the related scattering matrices define the database. The whole procedure is detailed in the following.
It is worth noting that the aim of the generated dataset is focused on providing a performance assessment of artificial neural networks (ANN) approaches for quantitative microwave breast imaging to be used for clinical applications of breast cancer. Future work will deal with the improvement of recovery performance, tumor detection and characterization at different operating conditions.
The first step of the generator consists of the creation of an elliptical-shape phantom. The dimensions of the ellipse axes vary uniformly in the range [ 6.5 , 12 ] cm. The center of the ellipse is randomly set inside a circle of radius 1 cm and located in the middle of the imaging domain. The orientation of the ellipse and its thickness, emulating the skin layer, are also randomly selected in the range [ 0 , 2 π ] and [ 1.5 , 2.5 ] mm, respectively.
The geometry of breast inner tissues is generated by exploiting a stochastic 2D multi-fractal random field generator [77]. The obtained map is segmented into three tissues, i.e., the fibro-glandular, the transitional and the adipose tissues. Similar to [12], we chose to simulate four breast classes characterized by different tissues percentage, as reported in Table 2, to test the inversion performance.
In order to generate realistic values of the tissues dielectric features, they are generated according to the statistical distributions estimated from the database described in [12]. Finally, spatial correlation among neighbouring pixels has been added to the data. This procedure is repeated to populate the database composed of 120,000 profiles (30,000 per each breast class). For each profile (discretized into 108 × 108 pixels), the scattering matrix has been computed via a fast Fourier transform-conjugate gradient (FFT-CG) forward solver based on the MoM [68] assuming a multiview-multistatic system with transmitters and receivers located in 30 angular equally-spaced locations on a measurement circle of radius 12 cm. The transmitting signal is a line source at a fixed frequency which impinges on an imaging domain of size 15 × 15 cm2 discretized according to Richmond’s rule [78].
The database has been partitioned in order to use the 85% of the profiles for the training phase, the 10% for the validation phase and the remaining 5% for the testing.

3.2. Neural Network Design and Training

Differently from other approaches presented in literature, in this manuscript a fully-connected architecture is considered instead of the commonly adopted CNNs. Some preliminary (and partial) work can be found in [79,80,81].
The key concept of an EIS inversion consists of the fact that a direct mapping between two different spaces, i.e., the scattered field and the dielectric contrast, is established. Most of the work proposed in the scientific literature focuses on the use of CNNs architecture [56,57,58], but these architectures are very useful typically when the input information flow is local, i.e., each output value is related to an input subset, such as denoising or despeckling applications [82,83]. Conversely, in the case of EIS problems, there is not a direct link between the data and unknown spaces, therefore a preliminary extra-mapping is usually required for reaching good accuracy and optimal inversion performance. Practically, an initial raw inversion is mostly performed to move from the data to the unknown space, and then a CNN architecture is usually applied which represents the second part of a two-step inversion procedure [58]. Thus, for imaging purposes, the fully connected ANN’s feature of having global links between all nodes of consecutive layers seems to be the key aspect to us.
A sketch of the adopted neural network architectures is reported in Figure 3. As previously mentioned, a direct mapping between the scattered field samples ( y i , i = { 1 , 2 , , M } , with M number of scattered field samples) and the tissues relative permittivity ϵ r and conductivity σ pixel by pixel ( x j , j = { 1 , 2 , , N } , with N number of pixels in the imaging domain) is established. In the network architectures considered in the following, an adaptive moment estimation method (ADAM) was employed to minimize an l 2 -cost function with an initial learning rate of 5 · 10 5 and 30 epochs per each training phase. It is worth noting that in all the considered numerical tests, a wide enough dataset (composed of 102,000 breast profiles for the training) was adopted to avoid the risk of overfitting, as confirmed by the validation and training loss curves. This aspect, which also drives a heavy computational burden due to the high number of trainable parameters to be learnt, represents one of the main limiting factors compared to convolutional neural networks and other conventional approaches. All the information mentioned previously regarding the adopted neural network architecture and its training is summarized in Table 3 for the sake of clarity.
A proper numerical analysis of the optimal network layout in terms of number of hidden layers and neurons is reported in the performance assessment in Section 4.1. Regarding the computational time, the training of each architecture depends on the number of involved layers and varies from 3 h (1 layer) up to 8 h (5 layers), while the inversion phase is real-time, requiring about 3.5 ms to estimate the permittivity and conductivity maps. These evaluations were performed on a Linux 64 bit workstation with an AMD Ryzen 3990X processor and an NVIDIA Quadro RTX 6000 graphics card.

3.3. Quality Performance Indicator on Testing Population

In order to perform a numerical assessment, some quantitative, fair metrics to determine the quality of the recoveries are required. Conventionally in the image processing community, parameters such as the normalized root mean square error (NRMSE) and the structural similarity index measure (SSIM) are considered. Nevertheless, every metric has its own peculiarities and makes it possible to capture only a limited amount of information regarding the retrieved image. Thus, in order to integrate this information and to evaluate the capability of the approach in terms of details retrieval, another quality measure based on the analysis of the spectral features of the reconstruction profile is here proposed. This approach can be considered as the evaluation of the filtering properties of the inverse radiating operator and of the imaging algorithm.
For every image of the testing set, the 2D power spectrum in polar coordinates ( ν , θ ) was evaluated and averaged with respect to all the images spectra, obtaining a mean-squared spectrum T ( ν , θ ) . By averaging along the angle coordinate, a one-dimensional spectrum S ( ν ) is obtained, i.e.,
T ( ν , θ ) = 1 N n = 1 N | F ( x ˜ n ) | 2 ,
S ( ν ) = 1 2 π π π T ( ν , θ ) d θ ,
in which N = 6000 is the number of testing profiles, F is the Fourier transform operator, x ˜ n is the retrieved estimate of the breast profile obtained via ANN inversion.
The motivation supporting this performance indicator lies in the fact that the closer the recovery is to the true profile, the closer their spectral information will be, providing an indirect measure of the filtering properties of the inverse radiating operator and of its impact at each frequency. This procedure clearly yields a spectral comparison rather than a direct similarity measure. Thus, the radial spectrum indicator of Equation (4) is exploited to further improve the reconstruction performance analysis of the different ANN architectures considered in Section 4.1 and to properly identify a good trade-off between computational burden (for the network training) and high accuracy values.

4. Results

In this section a numerical analysis to properly design the network architecture, i.e., the number of nodes and layers, will be addressed. Moreover, the methodology is tested in the case of different scenarios. The recovery performance was evaluated by comparing the proposed ANN reconstructions with those ones obtained via conventional nonlinear approaches.

4.1. Performance Assessment

In order to find a trade-off between computational complexity and accuracy in the recovery, a performance assessment on different network topologies was carried out. As a first analysis, the impact of the number of nodes was evaluated in the case of a three-layer fully-connected network. As expected, Figure 4 shows that the higher the number of nodes, the better the quality of the recovery since the averaged spectrum S ( ν ) tends to the ideal behavior of the true profile. Nevertheless, the higher the number of nodes, the higher the computational burden, which yields a longer training time. Thus, the 2000-node case was selected since the improvement in the averaged spectrum trend starts to be negligible from this topology and represents a good trade-off for the computational complexity.
Complementary to the previous analysis, a similar study about the number of layers was performed. To this aim, the averaged spectrum S ( ν ) is reported while increasing the number of hidden layers from one to five and adopting 2000 nodes per each layer. The related results are shown in Figure 5. In this case, there is no visible difference in the recovery performance among the networks using more than three layers.
Thus, in order to keep the computational burden reasonably low and, at the same time, a reconstruction quality sufficiently high, a three-layer, 2000-node architecture was selected (except the last regression layer not considered in these analyses). It is worth noting that we selected a number of breast profiles in the dataset high enough for ensuring an effective training for all the architectures considered in this work.

4.2. Comparison with Other Nonlinear Approaches

The performance of the real-time ANN approach proposed has previously been compared with other well-known nonlinear methods, which are the distorted Born iterative method (DBIM) [84,85] and the contrast source inversion (CSI) [86]. In the former case, a solution to the considered inverse problem is sought via a series of linear problems which gradually recover the full nonlinear profile, while in the latter case an iterative inversion scheme based on a functional minimisation involving both Equations (1) and (2) is exploited. These conventional nonlinear inversion approaches are iterative and usually employ local minimisation schemes in order to avoid prohibitive computational time required by global optimisation schemes. Under these assumptions, the choice of the initial guess for starting the minimization procedure becomes paramount for obtaining good recoveries. In more detail, the authors adopted the DBIM implementation proposed in [70], which is a multi-frequency, multi-threshold sparse-based approach named AMTISTA, and, for the CSI implementation, the one proposed in [87] in its cross-correlated formulation (CC-CSI).
Figure 6 and Figure 7 show the retrieved complex permittivity maps (i.e., real and imaginary parts) for four different breast profiles, one per each breast class (A-D) at the operating frequency of 1 GHz and with data corrupted by additive white Gaussian noise (AWGN) with a signal-to-noise ratio (SNR) equal to 30 dB. Further information regarding the considered configuration and the data generation is reported in Section 3.1. More specifically, these figures compare the recoveries obtained by the proposed ANN approach with other conventional nonlinear inversion strategies. It is worth noting that the recoveries obtained via adopting AMTISTA and CC-CSI were carried out with a frequency-hopping process using five frequencies—{200, 400, 600, 800, 1000} MHz—and starting the initial local minimization procedure from the background. Conversely, the proposed neural network approach processes the single-frequency data with no a priori information and in a direct fashion, i.e., no iteration is required, and it is still able to guarantee good accuracy and a better resolution performance compared with the considered conventional nonlinear approaches as reported in Figure 6 and Figure 7.
To realize a quantitative comparison among the considered reconstructions, three quality metrics were employed further to the quality spectral indicator reported in Equation (4): the structural similarity index (SSIM, values in the range [0, 1], higher is better), the normalised root mean square error (NRMSE, values in the range [0, 1], lower is better) and the correlation coefficient (CORR, values in amplitude in the range [0, 1], higher is better). These metrics were evaluated via default Matlab® functions and are reported in Table 4, while the spectral indicators of Equation (4) are shown in Figure 8.
By comparing the recoveries and looking at the quality metrics, which are better in the ANN case than in the AMTISTA and CC-CSI cases, it is clear that the proposed ANN approach outperforms these classical techniques, providing good performance also on the conductivity estimate, which represents the hardest part to be retrieved by conventional methods and also the most important information for diagnostic and therapeutic treatments.
As further proof of the good and stable recovery performance obtained by exploiting the proposed ANN architecture, Figure 9 and Table 5 report the histograms and corresponding means and standard deviations of the considered quality metrics for the whole testing population. It is interesting to note by considering these values that the recovery performance remains quite stable on the whole testing population regardless of the breast type.
Concerning the computational time, the proposed ANN approach results in a real-time procedure, conversely from conventional nonlinear approaches which are time-consuming per frequency and, sometimes, also per iteration (e.g., the AMTISTA case, for which a forward solver has to be run at each inversion iteration, involving at least 1.5-h simulation for 30 iterations for all the five frequencies on the same workstation). It is worth observing that the proposed ANN approach is able to provide information about the shape and location of the breast and, more specifically, about the skin layer, which still represents an issue for its dielectric features and thickness at the considered operating frequency.
Furthermore, the good recovery performance noticeable in the reconstructions related to tissues relative permittivity maps are also evident in the corresponding conductivity maps, which still represents the hardest, most challenging part in nowadays inverse scattering imaging approaches.
In order to test the robustness of the proposed method versus noise, Figure 10 illustrates a comparison among the recoveries of the selected breast phantoms reported in Figure 6 and Figure 7 for a much lower value of the SNR, i.e., 5 dB. The corresponding quality metrics are reported in Table 6. These results prove that the proposed approach is robust against noise and has good inversion performance in very noisy scenarios.
To further prove the robustness of the approach, its recovery performance was tested on realistic numerical breast phantoms generated from magnetic resonance images of female breasts collected by the Cross-disciplinary Electromagnetics Laboratory at the University of Wisconsin-Madison (CEM-UW) [12,75].
Thus, a few slices of these breast phantoms reported in the repository were considered, i.e., breast ID: 012204, slice 54, and breast ID: 062204, slices 47 and 71, which correspond to some examples of classes B and C, respectively. The measurement configuration and frequency were the same as the ones reported in Section 3.1 and an AWGN noise with SNR = 30 dB was applied to the data. As proved by the results shown in Figure 11, the recoveries were still good, especially on the conductivity of the considered breast profiles, and also allows good enough resolution of the skin layer and a fast inversion approach.

5. Conclusions

In this paper, a fully-connected artificial neural network for real-time microwave breast imaging applications has been proposed. The considered network was trained by exploiting an in-house realistic breast-like phantom generator adopted for the dataset population and by in-house forward solver codes for the generation of the scattered field data matrices to be exploited in the training of the network.
The methodology was tested in numerical noisy scenarios with different values of SNR proving good robustness against noise (i.e., 30 and 5 dB) and a performance assessment for the choice of the network architecture was carried out. It is worth noting that the recovery performance depends on the network architecture and each topology provides an improvement in retrieving information in proper bandwidths of the scenario under test.
The results seem good in comparison with conventional nonlinear inverse scattering approaches from a quantitative point of view as confirmed by the metrics reported in Table 4 and Table 5. Indeed, by considering the quality metrics reported in Table 4, it is worth noting that the proposed neural network-based approach allows us to obtain a better recovery performance compared with the considered state-of-the-art approaches in terms of permittivity and conductivity. Furthermore, the values of the metrics reported in Table 5 confirm the general improvement in the quality of the recoveries obtained by the proposed approach on the whole testing population. Lastly, the considered ANN approach proved to be robust versus noise and promising enough to be used on realistic anthropomorphic breast phantoms and, potentially, on experimental data.
The proposed real-time methodology allows us to obtain good recoveries both in terms of permittivity and conductivity maps, the latter usually representing an issue in inverse scattering approaches but also important, useful information for biomedical applications in diagnostics and therapy.
Future work will focus on the design of advanced neural network strategies to improve the quality of the results, enhancing resolution capability, and to process a three-dimensional dataset for the detection and characterization of potential malignant areas, also via exploiting more complex architecture and testing on laboratory experimental data.

Author Contributions

Conceptualization, M.A. and F.B.; resources, M.A. and S.F.; methodology, M.A., S.F. and V.P.; validation, S.F. and M.A.; writing—review and editing, S.F., M.A., F.B. and V.P.; visualization, S.F. and M.A.; supervision, F.B. and V.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data are available prior request to the authors.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
MRIMagnetic Resonance Imaging
SNRSignal to Noise Ratio
EISElectromagnetic Inverse Scattering
CNNConvolutional Neural Network
ANNArtificial Neural Network
UATUniversal Approximation Theorem
AWGNAdditive, White, Gaussian Noise
MoMMethod of Moments
FFT-CGFast Fourier Transform-Conjugate Gradient
NRMSENormalised Root mean Square Error
SSIMStructural Similarity Index Measure
DBIMDistorted Born Iterative Method
CSIContrast Source Inversion
CC-CSICross Correlated-Contrast Source Inversion

References

  1. Fitzmaurice, C.; Dicker, D.; Pain, A.; Hamavid, H.; Moradi-Lakeh, M.; MacIntyre, M.F.; Allen, C.; Hansen, G.; Woodbrook, R.; Wolfe, C.; et al. The global burden of cancer 2013. JAMA Oncol. 2015, 1, 505–527. [Google Scholar] [CrossRef] [PubMed]
  2. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics. CA Cancer J. Clin. 2015, 65, 5. [Google Scholar] [CrossRef]
  3. Zebari, D.; Ibrahim, D.; Zeebaree, D.; Haron, H.; Salih, M.; Damaševičius, R.; Mohammed, M. Systematic review of computing approaches for breast cancer detection based computer aided diagnosis using mammogram images. Appl. Artif. Intell. 2021, 35, 2157–2203. [Google Scholar] [CrossRef]
  4. Jaglan, P.; Dass, R.; Duhan, M. Breast cancer detection techniques: Issues and challenges. J. Inst. Eng. Ser. 2019, 1–8. [Google Scholar] [CrossRef]
  5. Herranz, M.; Ruibal, A. Optical imaging in breast cancer diagnosis: The next evolution. J. Oncol. 2012, 2012. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Sehgal, C.M.; Weinstein, S.P.; Arger, P.H.; Conant, E.F. A review of breast ultrasound. J. Mammary Gland. Biol. Neoplasia 2006, 11, 113–123. [Google Scholar] [CrossRef]
  7. Kapur, A.; Carson, P.L.; Eberhard, J.; Goodsitt, M.M.; Thomenius, K.; Lokh, W.M.; Buckley, D.; Roubidoux, M.A.; Helvie, M.A.; Booi, R.C.; et al. Combination of digital mammography with semi-automated 3D breast ultrasound. Technol. Cancer Res. Treat. 2004, 3, 325–334. [Google Scholar] [CrossRef] [Green Version]
  8. James, A.P.; Dasarathy, B.V. Medical image fusion: A survey of the state of the art. Inf. Fusion 2014, 19, 4–19. [Google Scholar] [CrossRef] [Green Version]
  9. Nikolova, N.K. Microwave imaging for breast cancer. IEEE Microw. Mag. 2011, 12, 78–94. [Google Scholar] [CrossRef]
  10. Aldhaeebi, M.A.; Alzoubi, K.; Almoneef, T.S.; Bamatraf, S.M.; Attia, H.; Ramahi, O.M. Review of microwaves techniques for breast cancer detection. Sensors 2020, 20, 2390. [Google Scholar] [CrossRef]
  11. O’Loughlin, D.; O’Halloran, M.; Moloney, B.M.; Glavin, M.; Jones, E.; Elahi, M.A. Microwave Breast Imaging: Clinical Advances and Remaining Challenges. IEEE Trans. Biomed. Eng. 2018, 65, 2580–2590. [Google Scholar] [CrossRef] [PubMed]
  12. 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. [Google Scholar] [CrossRef] [PubMed]
  13. Fear, E.C.; Hagness, S.C.; Meaney, P.M.; Okoniewski, M.; Stuchly, M.A. Enhancing breast tumor detection with near-field imaging. IEEE Microw. Mag. 2002, 3, 48–56. [Google Scholar] [CrossRef]
  14. Conceição, R.C.; Mohr, J.J.; O’Halloran, M. An Introduction to Microwave Imaging for Breast Cancer Detection; Springer: Berlin/Heidelberg, Germany, 2016. [Google Scholar]
  15. Chandra, R.; Zhou, H.; Balasingham, I.; Narayanan, R.M. On the Opportunities and Challenges in Microwave Medical Sensing and Imaging. IEEE Trans. Biomed. Eng. 2015, 62, 1667–1682. [Google Scholar] [CrossRef] [PubMed]
  16. Colton, D.L.; Kress, R.; Kress, R. Inverse Acoustic and Electromagnetic Scattering Theory; Springer Nature: Berlin/Heidelberg, Germany, 2019; Volume 93. [Google Scholar]
  17. Qin, Y.; Rodet, T.; Lambert, M.; Lesselier, D. Joint Inversion of Electromagnetic and Acoustic Data With Edge-Preserving Regularization for Breast Imaging. IEEE Trans. Comput. Imaging 2021, 7, 349–360. [Google Scholar] [CrossRef]
  18. Qin, Y.; Rodet, T.; Lambert, M.; Lesselier, D. Microwave breast imaging with prior ultrasound information. IEEE Open J. Antennas Propag. 2020, 1, 472–482. [Google Scholar] [CrossRef]
  19. Bevacqua, M.T.; Isernia, T. Shape reconstruction via equivalence principles, constrained inverse source problems and sparsity promotion. Prog. Electromagn. Res. 2017, 158, 37–48. [Google Scholar] [CrossRef] [Green Version]
  20. Tobon Vasquez, J.A.; Scapaticci, R.; Turvani, G.; Bellizzi, G.; Joachimowicz, N.; Duchêne, B.; Tedeschi, E.; Casu, M.R.; Crocco, L.; Vipiana, F. Design and experimental assessment of a 2d microwave imaging system for brain stroke monitoring. Int. J. Antennas Propag. 2019, 2019. [Google Scholar] [CrossRef] [Green Version]
  21. Porter, E.; Kirshin, E.; Santorelli, A.; Coates, M.; Popović, M. Time-domain multistatic radar system for microwave breast screening. IEEE Antennas Wirel. Propag. Lett. 2013, 12, 229–232. [Google Scholar] [CrossRef]
  22. Meaney, P.; Hartov, A.; Bulumulla, S.; Raynolds, T.; Davis, C.; Schoenberger, F.; Richter, S.; Paulsen, K. A 4-channel, vector network analyzer microwave imaging prototype based on software defined radio technology. Rev. Sci. Instrum. 2019, 90, 044708. [Google Scholar] [CrossRef]
  23. Duchesne, L.; Fasoula, A.; Kaverine, E.; Robin, G.; Bernard, J.G. Wavelia microwave breast imaging: Identification and mitigation of possible sources of measurement uncertainty. In Proceedings of the 2019 13th European Conference on Antennas and Propagation (EuCAP), Krakow, Poland, 31 March–5 April 2019; pp. 1–5. [Google Scholar]
  24. Kimura, K.; Kimura, N. Scattering Tomography Method and Scattering Tomography Device. U.S. Patent 10,101,282, 16 October 2018. [Google Scholar]
  25. Fedeli, A.; Maffongelli, M.; Monleone, R.; Pagnamenta, C.; Pastorino, M.; Poretti, S.; Randazzo, A.; Salvadè, A. A tomograph prototype for quantitative microwave imaging: Preliminary experimental results. J. Imaging 2018, 4, 139. [Google Scholar] [CrossRef] [Green Version]
  26. O’Loughlin, D.; Oliveira, B.L.; Glavin, M.; Jones, E.; O’Halloran, M. Advantages and disadvantages of parameter search algorithms for permittivity estimation for microwave breast imaging. In Proceedings of the 2019 13th European Conference on Antennas and Propagation (EuCAP), Krakow, Poland, 31 March–5 April 2019; pp. 1–5. [Google Scholar]
  27. Fear, E.C.; Bourqui, J.; Curtis, C.; Mew, D.; Docktor, B.; Romano, C. Microwave breast imaging with a monostatic radar-based system: A study of application to patients. IEEE Trans. Microw. Theory Tech. 2013, 61, 2119–2128. [Google Scholar] [CrossRef]
  28. Jeon, S.I.; Kim, B.R.; Son, S.H. Clinical trial of microwave tomography imaging. In Proceedings of the 2016 URSI Asia-Pacific Radio Science Conference (URSI AP-RASC), Seoul, Korea, 21–25 August 2016; pp. 1–2. [Google Scholar]
  29. Porter, E.; Coates, M.; Popović, M. An early clinical study of time-domain microwave radar for breast health monitoring. IEEE Trans. Biomed. Eng. 2015, 63, 530–539. [Google Scholar] [CrossRef]
  30. Yang, F.; Sun, L.; Hu, Z.; Wang, H.; Pan, D.; Wu, R.; Zhang, X.; Chen, Y.; Zhang, Q. A large-scale clinical trial of radar-based microwave breast imaging for asian women: Phase I. In Proceedings of the 2017 IEEE International Symposium on Antennas and Propagation & USNC/URSI National Radio Science Meeting, San Diego, CA, USA, 9–14 July 2017; pp. 781–783. [Google Scholar]
  31. Bourqui, J.; Fear, E.C. Systems for ultra-wideband microwave sensing and imaging of biological tissues. In Proceedings of the 2013 7th European Conference on Antennas and Propagation (EuCAP), Gothenburg, Sweden, 8–12 April 2013; pp. 834–835. [Google Scholar]
  32. Sani, L.; Ghavami, N.; Vispa, A.; Paoli, M.; Raspa, G.; Ghavami, M.; Sacchetti, F.; Vannini, E.; Ercolani, S.; Saracini, A.; et al. Novel microwave apparatus for breast lesions detection: Preliminary clinical results. Biomed. Signal Process. Control 2019, 52, 257–263. [Google Scholar] [CrossRef]
  33. Semenov, S.Y.; Svenson, R.H.; Bulyshev, A.E.; Souvorov, A.E.; Nazarov, A.G.; Sizov, Y.E.; Pavlovsky, A.V.; Borisov, V.Y.; Voinov, B.A.; Simonova, G.I.; et al. Three-dimensional microwave tomography: Experimental prototype of the system and vector Born reconstruction method. IEEE Trans. Biomed. Eng. 1999, 46, 937–946. [Google Scholar] [CrossRef]
  34. Joisel, A.; Mallorqui, J.; Broquetas, A.; Geffrin, J.M.; Joachimowicz, N.; Lossera, M.V.; Joire, L.; Bolomey, J.C. Microwave imaging techniques for biomedical applications. In Proceedings of the 16th IEEE Instrumentation and Measurement Technology Conference, Venice, Italy, 24–26 May 1999; Volume 3, pp. 1591–1596. [Google Scholar]
  35. Meaney, P.M.; Fanning, M.W.; Li, D.; Poplack, S.P.; Paulsen, K.D. A clinical prototype for active microwave imaging of the breast. IEEE Trans. Microw. Theory Tech. 2000, 48, 1841–1853. [Google Scholar]
  36. Meaney, P.M.; Fanning, M.W.; Zhou, T.; Golnabi, A.; Geimer, S.D.; Paulsen, K.D. Clinical microwave breast imaging—2D results and the evolution to 3D. In Proceedings of the 2009 International Conference on Electromagnetics in Advanced Applications, Torino, Italy, 14–18 September 2009; pp. 881–884. [Google Scholar]
  37. Grzegorczyk, T.M.; Meaney, P.M.; Kaufman, P.A.; Paulsen, K.D. Fast 3-D tomographic microwave imaging for breast cancer detection. IEEE Trans. Med. Imag. 2012, 31, 1584–1592. [Google Scholar] [CrossRef] [Green Version]
  38. Son, S.H.; Simonov, N.; Kim, H.J.; Lee, J.M.; Jeon, S.I. Preclinical prototype development of a microwave tomography system for breast cancer detection. ETRI J. 2010, 32, 901–910. [Google Scholar] [CrossRef]
  39. Simonov, N.A.; Jeon, S.I.; Son, S.H.; Lee, J.M.; Kim, H.J. 3D microwave breast imaging based on multistatic radar concept system. In Proceedings of the 2011 3rd International Asia-Pacific Conference on Synthetic Aperture Radar (APSAR), Seoul, Korea, 26–30 September 2011; pp. 1–4. [Google Scholar]
  40. Klemm, M.; Craddock, I.; Leendertz, J.; Preece, A.; Benjamin, R. Experimental and clinical results of breast cancer detection using UWB microwave radar. In Proceedings of the 2008 IEEE Antennas and Propag. Society International Symposium, San Diego, CA, USA, 5–11 July 2008; pp. 1–4. [Google Scholar]
  41. Preece, A.W.; Craddock, I.; Shere, M.; Jones, L.; Winton, H.L. Maria m4: Clinical evaluation of a prototype ultrawideband radar scanner for breast cancer detection. J. Med. Imaging 2016, 3, 033502. [Google Scholar] [CrossRef]
  42. Benny, R.; Anjit, T.A.; Mythili, P. An overview of microwave imaging for breast tumor detection. Prog. Electromagn. Res. 2020, 87, 61–91. [Google Scholar] [CrossRef]
  43. Bertero, M.; Boccacci, P. Introduction to Inverse Problems in Imaging; CRC Press: Boca Raton, FL, USA, 1998. [Google Scholar]
  44. Pastorino, M. Microwave Imaging; John Wiley & Sons: Hoboken, NJ, USA, 2010; Volume 208. [Google Scholar]
  45. Salman, A.O.; Gavrilov, S.P.; Vertiy, A.A. Microwave imaging of immersed dielectric bodies: An experimental survey. Electromagnetics 2005, 25, 567–585. [Google Scholar] [CrossRef]
  46. Bolomey, J.C.; Izadnegahdar, A.; Jofre, R.L.; Du, P.; Mezeray, C.; Peronnet, G. Microwave diffraction tomography for biomedical applications. IEEE Trans. Microw. Theory Tech. 1982, 30, 1998–2000. [Google Scholar] [CrossRef] [Green Version]
  47. Colton, D.; Haddar, H.; Piana, M. The linear sampling method in inverse electromagnetic scattering theory. Inverse Probl. 2003, 19, S105. [Google Scholar] [CrossRef]
  48. Cui, T.J.; Chew, W.C.; Aydiner, A.A.; Chen, S. Inverse scattering of two-dimensional dielectric objects buried in a lossy earth using the distorted born iterative method. IEEE Trans. Geosci. Remote Sens. 2001, 39, 339–346. [Google Scholar]
  49. Crocco, L.; Catapano, I.; Di Donato, L.; Isernia, T. The linear sampling method as a way to quantitative inverse scattering. IEEE Trans. Antennas Propag. 2012, 60, 1844–1853. [Google Scholar] [CrossRef]
  50. Cui, T.J.; Qin, Y.; Wang, G.L.; Chew, W.C. Low-frequency detection of two-dimensional buried objects using high-order extended born approximations. Inverse Probl. 2004, 20, S41. [Google Scholar] [CrossRef]
  51. Belkebir, K.; Elissalt, J.M.; Geffrin, J.M.; Pichot, C. Newton-kantorovich and modified gradient-inversion algorithms applied to ipswich data. IEEE Antennas Propag. Mag. 1996, 38, 41. [Google Scholar] [CrossRef] [Green Version]
  52. Nocedal, J.; Wright, S. Numerical Optimization; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  53. Pastorino, M.; Caorsi, S.; Massa, A. A global optimization technique for microwave nondestructive evaluation. IEEE Trans. Instrum. Meas. 2002, 51, 666–673. [Google Scholar] [CrossRef]
  54. Caorsi, S.; Gamba, P. Electromagnetic detection of dielectric cylinders by a neural network approach. IEEE Trans. Geosci. Remote Sens. 1999, 37, 820–827. [Google Scholar] [CrossRef]
  55. Rekanos, I.T. Neural-network-based inverse-scattering technique for online microwave medical imaging. IEEE Trans. Magn. 2002, 38, 1061–1064. [Google Scholar] [CrossRef]
  56. Xu, K.; Wu, L.; Ye, X.; Chen, X. Deep learning-based inversion methods for solving inverse scattering problems with phaseless data. IEEE Trans. Antennas Propag. 2020. [Google Scholar] [CrossRef]
  57. Li, L.; Wang, L.G.; Teixeira, F.L.; Liu, C.; Nehorai, A.; Cui, T.J. Deepnis: Deep neural network for nonlinear electromagnetic inverse scattering. IEEE Trans. Antennas Propag. 2018, 67, 1819–1825. [Google Scholar] [CrossRef]
  58. Wei, Z.; Chen, X. Deep-learning schemes for full-wave nonlinear inverse scattering problems. IEEE Trans. Geosci. Remote Sens. 2018, 57, 1849–1860. [Google Scholar] [CrossRef]
  59. Sanghvi, Y.; Kalepu, Y.; Khankhoje, U.K. Embedding deep learning in inverse scattering problems. IEEE Trans. Comput. Imaging 2019, 6, 46–56. [Google Scholar] [CrossRef]
  60. Guo, R.; Song, X.; Li, M.; Yang, F.; Xu, S.; Abubakar, A. Supervised descent learning technique for 2-D microwave imaging. IEEE Trans. Antennas Propag. 2019, 67, 3550–3554. [Google Scholar] [CrossRef]
  61. Massa, A.; Marcantonio, D.; Chen, X.; Li, M.; Salucci, M. DNNS as applied to electromagnetics, antennas, and propagation—A review. IEEE Antennas Wirel. Propag. Lett. 2019, 18, 2225–2229. [Google Scholar] [CrossRef]
  62. Aggarwal, H.K.; Mani, M.P.; Jacob, M. Modl: Model-based deep learning architecture for inverse problems. IEEE Trans. Med. Imag. 2018, 38, 394–405. [Google Scholar] [CrossRef]
  63. Lucas, A.; Iliadis, M.; Molina, R.; Katsaggelos, A.K. Using deep neural networks for inverse problems in imaging: Beyond analytical methods. IEEE Signal Process. Mag 2018, 35, 20–36. [Google Scholar] [CrossRef]
  64. Ashtari, A.; Noghanian, S.; Sabouni, A.; Aronsson, J.; Thomas, G.; Pistorius, S. Using a priori information for regularization in breast microwave image reconstruction. IEEE Trans. Biomed. Eng. 2010, 57, 2197–2208. [Google Scholar] [CrossRef]
  65. Khoshdel, V.; Asefi, M.; Ashraf, A.; LoVetri, J. Full 3D microwave breast imaging using a deep-learning technique. J. Imaging 2020, 6, 80. [Google Scholar] [CrossRef]
  66. Shah, P.; Moghaddam, M. Super resolution for microwave imaging: A deep learning approach. In Proceedings of the International Symposium on Antennas and Propagation & USNC/URSI National Radio Science Meeting, San Diego, CA, USA, 9–14 July 2017; pp. 849–850. [Google Scholar]
  67. Hornik, K.; Stinchcombe, M.; White, H. Universal approximation of an unknown mapping and its derivatives using multilayer feedforward networks. Neural Netw. 1990, 3, 551–560. [Google Scholar] [CrossRef]
  68. Isernia, T.; Pascazio, V.; Pierri, R. A nonlinear estimation method in tomographic imaging. IEEE Trans. Geosci. Remote Sens. 1997, 35, 910–923. [Google Scholar] [CrossRef]
  69. Bucci, O.M.; Isernia, T. Electromagnetic inverse scattering: Retrievable information and measurement strategies. Radio Sci. 1997, 32, 2123–2137. [Google Scholar] [CrossRef]
  70. Ambrosanio, M.; Kosmas, P.; Pascazio, V. A multithreshold iterative dbim-based algorithm for the imaging of heterogeneous breast tissues. IEEE Trans. Biomed. Eng. 2018, 66, 509–520. [Google Scholar] [CrossRef] [Green Version]
  71. Bevacqua, M.T.; Bellizzi, G.G.; Crocco, L.; Isernia, T. A method for quantitative imaging of electrical properties of human tissues from only amplitude electromagnetic data. Inverse Probl. 2019, 35, 025006. [Google Scholar] [CrossRef]
  72. Mojabi, P.; LoVetri, J. Enhancement of the krylov subspace regularization for microwave biomedical imaging. IEEE Trans. Med. Imag. 2009, 28, 2015–2019. [Google Scholar] [CrossRef]
  73. Shea, J.D.; Kosmas, P.; Van Veen, B.D.; Hagness, S.C. Contrast-enhanced microwave imaging of breast tumors: A computational study using 3D realistic numerical phantoms. Inverse Probl. 2010, 26, 074009. [Google Scholar] [CrossRef] [Green Version]
  74. Rubæk, T.; Meaney, P.M.; Meincke, P.; Paulsen, K.D. Nonlinear microwave imaging for breast-cancer screening using gauss–newton’s method and the cgls inversion algorithm. IEEE Trans. Antennas Propag. 2007, 55, 2320–2331. [Google Scholar] [CrossRef] [Green Version]
  75. Burfeindt, M.J.; Colgan, T.J.; Mays, R.O.; Shea, J.D.; Behdad, N.; Van Veen, B.D.; Hagness, S.C. MRI-derived 3-D-printed breast phantom for microwave breast imaging validation. IEEE Antennas Wirel. Propag. Lett. 2012, 11, 1610–1613. [Google Scholar] [CrossRef] [Green Version]
  76. Reimer, T.; Krenkevich, J.; Pistorius, S. An open-access experimental dataset for breast microwave imaging. In Proceedings of the 2020 14th European Conference on Antennas and Propagation (EuCAP), Copenhagen, Denmark, 15–20 March 2020; pp. 1–5. [Google Scholar]
  77. Schertzer, D.; Lovejoy, S. Nonlinear variability in geophysics: Multifractal simulations and analysis. In Fractals’ Physical Origin and Properties; Springer: Berlin/Heidelberg, Germany, 1989; pp. 49–79. [Google Scholar]
  78. Richmond, J. Scattering by a dielectric cylinder of arbitrary cross section shape. IEEE Trans. Antennas Propag. 1965, 13, 334–341. [Google Scholar] [CrossRef]
  79. Ambrosanio, M.; Franceschini, S.; Baselice, F.; Pascazio, V. Artificial Neural Networks for Quantitative Microwave Breast Imaging. Bioimaging 2020, 204–208. [Google Scholar]
  80. Franceschini, S.; Ambrosanio, M.; Baselice, F.; Pascazio, V. Neural Networks for Inverse Problems: The Microwave Imaging Case. In Proceedings of the 2021 15th European Conference on Antennas and Propagation (EuCAP), Düsseldorf, Germany, 22–26 March 2021; pp. 1–5. [Google Scholar]
  81. Ambrosanio, M.; Franceschini, S.; Baselice, F.; Pascazio, V. Machine learning for microwave imaging. In Proceedings of the 2020 14th European Conference on Antennas and Propagation (EuCAP), Copenhagen, Denmark, 15–20 March 2020; pp. 1–4. [Google Scholar]
  82. Zhang, K.; Zuo, W.; Chen, Y.; Meng, D.; Zhang, L. Beyond a gaussian denoiser: Residual learning of deep cnn for image denoising. IEEE Trans. Image Process. 2017, 26, 3142–3155. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Vitale, S.; Ferraioli, G.; Pascazio, V. Multi-objective cnn-based algorithm for sar despeckling. IEEE Trans. Geosci. Remote Sens. 2020, 1–14. [Google Scholar] [CrossRef]
  84. Chew, W.C.; Wang, Y. Reconstruction of two-dimensional permittivity distribution using the distorted born iterative method. IEEE Trans. Med. Imag. 1990, 9, 218–225. [Google Scholar] [CrossRef] [PubMed]
  85. Chew, W.C. Waves and Fields in Inhomogeneous Media; IEEE Press: Piscataway, NJ, USA, 1995. [Google Scholar]
  86. den Berg, P.M.V.; Abubakar, A. Contrast source inversion method: State of art. Prog. Electromagn. Res. 2001, 34, 189–218. [Google Scholar] [CrossRef] [Green Version]
  87. Sun, S.; Kooij, B.J.; Jin, T.; Yarovoy, A.G. Cross-correlated contrast source inversion. IEEE Trans. Antennas Propag. 2017, 65, 2592–2603. [Google Scholar] [CrossRef]
Figure 1. Sketch of the 2D geometry. The antennas are located on the measurement line Γ and arranged in a multiview-multistatic fashion.
Figure 1. Sketch of the 2D geometry. The antennas are located on the measurement line Γ and arranged in a multiview-multistatic fashion.
Bioengineering 09 00651 g001
Figure 2. Overview of the realistic numerical breast phantom generation procedure. (a) Elliptical-shape phantom generation with corresponding skin layer. (b) Internal-tissue spatial distribution via a random field generator [77] and (c) estimated electromagnetic (i.e., relative permittivity and conductivity) tissues assignment (for the sake of brevity, only the relative permittivity is reported).
Figure 2. Overview of the realistic numerical breast phantom generation procedure. (a) Elliptical-shape phantom generation with corresponding skin layer. (b) Internal-tissue spatial distribution via a random field generator [77] and (c) estimated electromagnetic (i.e., relative permittivity and conductivity) tissues assignment (for the sake of brevity, only the relative permittivity is reported).
Bioengineering 09 00651 g002
Figure 3. Sketch of the inversion procedure. The input of the fully-connected neural network is the scattering matrix E s c t re-arranged in a vector form y i , i = { 1 , 2 , , M } , with M number of scattered field samples. The output of is represented by the tissues relative permittivity ϵ r and conductivity σ pixel by pixel x j , j = { 1 , 2 , , N } , with N number of pixels in the imaging domain.
Figure 3. Sketch of the inversion procedure. The input of the fully-connected neural network is the scattering matrix E s c t re-arranged in a vector form y i , i = { 1 , 2 , , M } , with M number of scattered field samples. The output of is represented by the tissues relative permittivity ϵ r and conductivity σ pixel by pixel x j , j = { 1 , 2 , , N } , with N number of pixels in the imaging domain.
Bioengineering 09 00651 g003
Figure 4. Averaged spectrum S ( ν ) for six different three-layer architectures as a function of nodes number per each hidden layer.
Figure 4. Averaged spectrum S ( ν ) for six different three-layer architectures as a function of nodes number per each hidden layer.
Bioengineering 09 00651 g004
Figure 5. Averaged spectrum S ( ν ) for five different architectures with two-thousand nodes per layer as a function of the number of hidden layers.
Figure 5. Averaged spectrum S ( ν ) for five different architectures with two-thousand nodes per layer as a function of the number of hidden layers.
Bioengineering 09 00651 g005
Figure 6. Comparison between recoveries related to four different breast profiles (one per each class). Retrieved relative permittivity maps (dimensionless quantity) for AMTISTA (DBIM formulation as reported in [70]) (a,e,i,m), cross-correlated CSI (CC-CSI) [87] (b,f,j,n), and proposed ANN approach (c,g,k,o). Reference permittivity maps are reported in (d,h,l,p). Both the x and y axes are in centimeter units.
Figure 6. Comparison between recoveries related to four different breast profiles (one per each class). Retrieved relative permittivity maps (dimensionless quantity) for AMTISTA (DBIM formulation as reported in [70]) (a,e,i,m), cross-correlated CSI (CC-CSI) [87] (b,f,j,n), and proposed ANN approach (c,g,k,o). Reference permittivity maps are reported in (d,h,l,p). Both the x and y axes are in centimeter units.
Bioengineering 09 00651 g006
Figure 7. Comparison between recoveries related to four different breast profiles (one per each class). Retrieved conductivity maps (measure unit: S/m) for AMTISTA (DBIM formulation as reported in [70]) (a,e,i,m), cross-correlated CSI (CC-CSI) [87] (b,f,j,n), and proposed ANN approach (c,g,k,o). Reference conductivity maps are reported in (d,h,l,p). Both the x and y axes are in centimeter units.
Figure 7. Comparison between recoveries related to four different breast profiles (one per each class). Retrieved conductivity maps (measure unit: S/m) for AMTISTA (DBIM formulation as reported in [70]) (a,e,i,m), cross-correlated CSI (CC-CSI) [87] (b,f,j,n), and proposed ANN approach (c,g,k,o). Reference conductivity maps are reported in (d,h,l,p). Both the x and y axes are in centimeter units.
Bioengineering 09 00651 g007
Figure 8. One-dimensional power spectrum indicator of Equation (4) for the retrieved profiles reported in Figure 6 and Figure 7.
Figure 8. One-dimensional power spectrum indicator of Equation (4) for the retrieved profiles reported in Figure 6 and Figure 7.
Bioengineering 09 00651 g008
Figure 9. Histograms evaluated on the whole testing population for the adopted metrics.
Figure 9. Histograms evaluated on the whole testing population for the adopted metrics.
Bioengineering 09 00651 g009
Figure 10. Performance assessment of the proposed ANN approach for an SNR = 5 dB. Retrieved permittivity ((a,e,i,m)—dimensionless quantity) and conductivity ((c,g,k,o)—measure unit: S/m) maps, and corresponding reference profiles (b,f,j,n,d,h,l,p). Both the x and y axes are in centimeter units.
Figure 10. Performance assessment of the proposed ANN approach for an SNR = 5 dB. Retrieved permittivity ((a,e,i,m)—dimensionless quantity) and conductivity ((c,g,k,o)—measure unit: S/m) maps, and corresponding reference profiles (b,f,j,n,d,h,l,p). Both the x and y axes are in centimeter units.
Bioengineering 09 00651 g010
Figure 11. Recoveries obtained via the proposed ANN approach for some realistic breast phantoms proposed by [12,75]. Breast ID: 012204, slice 54: retrieved permittivity and conductivity maps (a,c) and reference (b,d). Breast ID: 062204, slices 47 and 71: retrieved permittivity and conductivity maps (e,i,g,k), and references (f,j,h,l), respectively.
Figure 11. Recoveries obtained via the proposed ANN approach for some realistic breast phantoms proposed by [12,75]. Breast ID: 012204, slice 54: retrieved permittivity and conductivity maps (a,c) and reference (b,d). Breast ID: 062204, slices 47 and 71: retrieved permittivity and conductivity maps (e,i,g,k), and references (f,j,h,l), respectively.
Bioengineering 09 00651 g011
Table 1. Summary of main microwave breast imaging systems with studies on patients [42].
Table 1. Summary of main microwave breast imaging systems with studies on patients [42].
SystemPositionProcessingFrequency [GHz]AntennaScan TimeReference
Bristol Universityproneradar3–8slot30 s[26]
Dartmouth Collegepronetomographic0.7–1.7monopole5 min[22]
ETRI (Korea)pronetomographic3–6monopole15 s/slice[28]
McGill University—tableproneradar2–4TWTLTLA18 min[29]
McGill University—wearablewearableradar2–4microstrip5 min[29]
Southern University of Chinaproneradar4–8horn4 min[30]
Calgary Universityproneradar1.3–7.6vivaldi30 min[31]
Microwave Visionproneradar1–4vivaldi10 min[23]
MammowaveproneHuygens1–9PulsON P20010 min[32]
Kobe Universitysupinetomographic0.05–12UWB30 min[24]
Table 2. Breast phantom tissues composition.
Table 2. Breast phantom tissues composition.
ClassPercentage of Tissue (%)
Fibro-glandularTransitionalAdipose
A5–205–1565–90
B20–3010–2050–70
C30–4015–2040–55
D40–6520–2510–40
Table 3. Details on the implemented neural network architecture.
Table 3. Details on the implemented neural network architecture.
Neural Network Details
Input size1800  30 × 30 × 2
Output size23,328  108 × 108 × 2
Optimization methodADAM, initial learning rate: 5 · 10 5
Training epochs30
Mini-batch size64
Architecture3 layers, 2000 nodes per each
Population120,000 (training: 85%, validation: 10%, testing: 5%)
Table 4. Quality metrics for both permittivity and conductivity retrieved maps of Figure 6 and Figure 7. The bold quantities refer to the best value per row, identifying the approach which provides the best recovery according to the considered metric.
Table 4. Quality metrics for both permittivity and conductivity retrieved maps of Figure 6 and Figure 7. The bold quantities refer to the best value per row, identifying the approach which provides the best recovery according to the considered metric.
Breast
Phantom
Quality
Metric
AMTISTACC-CSIANN
ϵ σ ϵ σ ϵ σ
Breast 1SSIM0.100.020.200.160.450.24
NRMSE0.150.620.357.000.060.22
CORR0.700.560.060.230.890.86
Breast 2SSIM0.050.220.170.240.460.21
NRMSE0.140.360.353.440.100.21
CORR0.780.750.440.180.850.86
Breast 3SSIM0.060.220.190.180.450.27
NRMSE0.140.360.191.700.080.17
CORR0.760.740.660.170.860.89
Breast 4SSIM0.060.040.250.300.400.21
NRMSE0.140.460.100.340.070.17
CORR0.790.790.810.740.870.88
Table 5. Means and standard deviations of the considered quality metrics per each class of the testing population (6000 breast profiles in total). The number of profiles per class is reported in brackets).
Table 5. Means and standard deviations of the considered quality metrics per each class of the testing population (6000 breast profiles in total). The number of profiles per class is reported in brackets).
Breast
Class
Quality
Metrics
Retrieved PermittivityRetrieved Conductivity
MeanStandard DeviationMeanStandard Deviation
A
(1422)
SSIM0.4580.0350.1820.040
NRMSE0.1020.0310.3450.100
CORR0.7950.0660.7790.067
B
(2115)
SSIM0.4530.0370.1790.043
NRMSE0.1020.0300.2530.066
CORR0.8180.0560.8350.045
C
(1837)
SSIM0.4310.0360.1580.037
NRMSE0.1040.0350.2250.062
CORR0.8180.0600.8450.046
D
(626)
SSIM0.4170.0350.1470.035
NRMSE0.0870.0280.1760.043
CORR0.8490.0480.8710.037
Table 6. Quality metrics for retrieved maps of Figure 10 in the case of high and low noise.
Table 6. Quality metrics for retrieved maps of Figure 10 in the case of high and low noise.
SNR30 dB5 dB
Quality MetricSSIMNRMSECORRSSIMNRMSECORR
Breast 1 ϵ 0.500.060.890.460.080.85
σ 0.210.230.860.210.310.81
Breast 2 ϵ 0.460.100.850.440.130.78
σ 0.210.210.860.200.290.80
Breast 3 ϵ 0.450.080.860.360.120.79
σ 0.270.170.890.180.260.82
Breast 4 ϵ 0.450.070.870.380.120.77
σ 0.190.170.880.090.280.80
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ambrosanio, M.; Franceschini, S.; Pascazio, V.; Baselice, F. An End-to-End Deep Learning Approach for Quantitative Microwave Breast Imaging in Real-Time Applications. Bioengineering 2022, 9, 651. https://doi.org/10.3390/bioengineering9110651

AMA Style

Ambrosanio M, Franceschini S, Pascazio V, Baselice F. An End-to-End Deep Learning Approach for Quantitative Microwave Breast Imaging in Real-Time Applications. Bioengineering. 2022; 9(11):651. https://doi.org/10.3390/bioengineering9110651

Chicago/Turabian Style

Ambrosanio, Michele, Stefano Franceschini, Vito Pascazio, and Fabio Baselice. 2022. "An End-to-End Deep Learning Approach for Quantitative Microwave Breast Imaging in Real-Time Applications" Bioengineering 9, no. 11: 651. https://doi.org/10.3390/bioengineering9110651

APA Style

Ambrosanio, M., Franceschini, S., Pascazio, V., & Baselice, F. (2022). An End-to-End Deep Learning Approach for Quantitative Microwave Breast Imaging in Real-Time Applications. Bioengineering, 9(11), 651. https://doi.org/10.3390/bioengineering9110651

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