<figure><img src="image_1.png" /></figure>

HAL
open science

# Application of Spectral Approach Combined with U-NETs for Quantitative Microwave Breast Imaging

Ambroise Diès, Hélène Roussel, Nadine Joachimowicz

► To cite this version:

Ambroise Diès, Hélène Roussel, Nadine Joachimowicz. Application of Spectral Approach Combined with U-NETs for Quantitative Microwave Breast Imaging. Sensors, 2025, 25 (8), pp.2450. 10.3390/s25082450 . hal-05033865

HAL Id: hal-05033865

https://hal.science/hal-05033865v1

Submitted on 14 Apr 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Application of Spectral Approach Combined with U-NETs for Quantitative Microwave Breast Imaging

Ambroise Diès<fnref n="1" /> ID, Hélène Roussel<fnref n="1" /> ID and Nadine Joachimowicz <fnref n="1" />,<fnref n="2" />,* ID

**Abstract:** This study focuses on breast imaging. A spectral approach based on the Fourier diffraction theorem is combined with a pair of U-NETs to perform real-time quantitative human breast imaging. The U-NET pair is trained based on the input of an induced current spectrum and the output of a contrast dielectric spectrum. A spectral database is constructed using combinations of anthropomorphic cavities. The weighted mean absolute percentage error (WMAPE) loss is associated with the Adam optimizer to perform optimization. Numerical results are presented to validate the proposed concept to demonstrate the transformation brought about by the U-NETs.

**Keywords:** microwave imaging; U-NET; deep learning; backpropagation; spectral techniques; anthropomorphic breast model

## 1. Introduction

In inverse scattering for microwave breast imaging, a quantitative description of the breast, corresponding to its dielectric properties, is calculated from the measurement of the scattered field induced by the presence of the breast. From a mathematical viewpoint, this is an ill-posed problem, due in part to the lack of access to evanescent waves, and a non-linear problem because of the multiple scattering effects inside the breast. To solve this issue, iterative optimization methods incorporating prior information have been developed to minimize the deviation between the calculated and measured scattered field distributions. The contrast source inversion method (CSI) [1–3], the distorted born iterative method (DBIM) [4,5], Gauss–Newton (GN) [6], and the Broyden–Fletcher–Goldfarb–Shanno method (L-BFGS) [7] are examples of deterministic optimization algorithms used for medical applications. Recent developments in deep learning methods, such as perceptrons (ANNs) and convolutional neural networks (CNNs), have been under investigation. These are deterministic methods based on a stochastic optimization process that incorporates an entire database as prior information. Numerous studies have shown that it is possible to improve microwave imaging techniques by combining them with deep learning, as proposed in [8–11].

More specifically, Yingying Qin et al. [2] and Mojabi et al. [12] investigated the use of multichannel CNNs and U-NETs, respectively, to combine microwave and ultrasound images using the databases described in [13] and [14]. In [15], Ambrosanio et al. used smart-automated generation of 2D maps to create a large database on which they trained an (ANN) to transform the measured scattered field into dielectric properties.

Flores et al. [16] used a quadratic BIM algorithm, through a complex valued U-NET, with the database in [13]. Khoshdel et al. [17] performed tumor detection using a 3D,

<figure><img src="image_3.png" /></figure>

### check for updates

Academic Editors: Yifan Chen and
Xiong Wang

Received: 19 February 2025

Revised: 14 March 2025

Accepted: 28 March 2025

Published: 13 April 2025

**Citation:** Diès, A.; Roussel, H.; Joachimowicz, N. Application of Spectral Approach Combined with U-NETs for Quantitative Microwave Breast Imaging. *Sensors* **2025**, *25*, 2450. https://doi.org/10.3390/s25082450

**Copyright:** © 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">¹ Sorbonne Université, CNRS, Laboratoire de Génie Electrique et Electronique de Paris, 75252 Paris, France; ambroise.dies@sorbonne-universite.fr (A.D.); helene.roussel@sorbonne-universite.fr (H.R.)</footnote>

<footnote n="2">² Université Paris Cité, F-75006 Paris, France</footnote>

<footnote>* Correspondence: nadine.joachimowicz@geeps.centralesupelec.fr</footnote>
automatically generated random breast tumor database. Marijn Borghouts et al. [18] generated a tumor probability map. Costanzo and Flores used reinforcement learning techniques to perform tumor detection [19].

Our contribution is part of the historical framework of the GeePs-L2S microwave camera [20,21]. We are interested in the possibility of training deep learning models to transform the qualitative images produced by the microwave camera, which represent the distribution of induced currents obtained by backpropagating the measured scattered field into a quantitative image representing the distribution of the object's dielectric properties. Compared with other techniques, the proposed method uses operations in the spectral domain, making the process less sensitive to noise by taking advantage of high-frequency spatial filtering and less time-consuming with respect to computation. Contributions which associate deep learning with spectral imaging exist in the optical domain [22], but to the best of our knowledge, the use of a spectral approach has not been investigated for quantitative microwave imaging of the breast. In addition, most existing papers on deep learning microwave breast imaging are based on realistic databases created from MRI scans performed on a limited number of patients. Our model is trained on a limited three dimensional database, customized (in-house) from anthropomorphic breast phantoms. There are five major reasons for this investigation: to test the validity of the proposed spectral quantitative method, which works with a limited amount of data; to offer low-cost realistic data that can be easily manipulated; to avoid health and personal protection legislation; to provide anthropomorphic printable models inspired by MRI scans to build up a reference database along with physical phantoms for experimentation; and to easily increase and diversify data, by designing a variety of phantoms of different shapes and sizes. Extension of our 3D database, requiring large-scale production of anthropomorphic cavities and associated simulations, is time consuming. Consequently, this will be investigated in the future. This paper focuses on assessment of the original proposed imaging technique, rather than performance metrics relevant for working with an extensive database. The preliminary results are presented to validate the proposed concept.

This paper is organized as follows: Section 2 highlights the database production process; the anthropomorphic breast model files are used to calculate the truncated spectrum of the induced currents in the planar configuration of the microwave camera. The optimization process including the architecture of the U-NETs and the normalization-denormalization process applied to the data is detailed in Section 2.5. Select numerical results that illustrate both the strengths and limitations of the approach are shown in Section 3. Finally, a discussion about the range of this contribution and its future perspectives is given in Section 4.

## 2. Materials and Methods

### 2.1. *Generation of Anthropomorphic Breast Phantoms*

To implement deep learning techniques, we created a database based on three-dimensional anthropomorphic models inspired by the GeePs-L2S breast phantom [23,24]. The breast is divided into three different parts—skin-associated fat, glandular tissues, and tumors—as shown in Figure 1. We developed, in addition to the GeePs-L2S breast model, three models for fat and skin and five for glandular tissue. This enabled us to produce $3 \times 5 - 1 = 14$ models without tumors, as glandular tissue cavity number 5 does not fit inside the smallest fat cavity. This number rose to $14 \times 4 = 56$ without (case 1) or with tumors placed in the adipose (case 2) or glandular (case 3) or both (case 4). In a second step, we considered the mirror of each phantom to double the number of models. Because the phantoms are volumetric, and because of the propagation of the wave, the mirror does not correspond to a 180-degree rotation. Finally, a training set of $56 \times 2 \times 16 = 1792$ images
at the center of the models was obtained by rotating them, considering a full rotation corresponding to 16 different views.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Three-dimensional cavities used to build training database (a); example of how to combine different parts of printable phantom (b).</figcaption></figure>

Different models of tumors are shown in Figure 2a. Their size varies from 0.5 cm to 2 cm and they are located at different places to create various configurations (cases 2, 3, and 4). Figure 2b, shows the distribution of tumor shapes associated with a letter—A, B, C, D, or E—per cavity. The distribution of tumors is chosen so that each tumor is associated with a cavity only once. As a result, there is no correlation between the use of a specific shape in the phantom and the presence of a tumor in the constructed breast.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Tumor dataset (a); repartition of tumors among cavities (b).</figcaption></figure>

However, the training and validation sets are nevertheless correlated, as the current dataset is built from a combination of the various cavities. The data are randomly shuffled and separated, with 80% for training—1433 spectra— and 20% for validation—359 spectra. We set the uncorrelated GeePs-L2S phantom aside for testing. Therefore, the testing set contains 32 spectra (16 rotations and their mirror images). The presented example is used to show that a generalization seems feasible with a wider dataset.

The dielectric properties of the various parts of the phantom are those provided by standard reference databases, e.g., the IFAC [25] or SUPELEC RECIPES [26]. The values of the relative permittivity, $ε_r$, and conductivity, $σ$ [S/m], for each tissue used to perform the simulations are given at 2.45 GHz in Table 1. Because variability exists for human tissues, we added a ±5% variation value on the dielectric properties as a first try. Higher variabilities will be investigated in the future, as studies in the field mention stronger variations of ±10% due to the dehydration of the tissues [27] or based on canine stroke measurements [28,29]. The values are constant within each cavity due to the simulation constraints (the numerical model is based on a method of moments solution to the surface integral equation). The distribution of permittivity variations between the phantoms is uniform.
<table><caption>**Table 1.** $\epsilon_r$ and $\sigma$ for the four tissues of the phantom.</caption><thead><tr><th></th><th>Fat</th><th>Glandular</th><th>Skin (Wet)</th><th>Tumor</th></tr></thead><tbody><tr><td>ε<sub>r</sub>(±5%)</td><td>5</td><td>44</td><td>42</td><td>53</td></tr><tr><td>σ(±5%) [S/m]</td><td>0.1</td><td>1.5</td><td>1.6</td><td>1.8</td></tr></tbody></table>

## 2.2. Experimental Set-Up: GeePs-L2S Planar Camera

The various breast phantoms in the database are studied by placing them inside the GeePs-L2S microwave planar camera as illustrated in Figure 3 [21]. The phantom model is placed in a water tank whose dielectric properties are $\epsilon_r = 73$ and $\sigma = 1$ [S/m]. The tank is illuminated by a monochromatic plane wave polarized along the z-axis at frequency $f = 2.45$ GHz. The scattered field is computed on 64 × 64 points of the retina (representing the location of the dipoles) using the commercial electromagnetic software WIPL-D [30], based on a surface integral formulation solved by method of moments. This constitutes the simulated data for the imaging problem, solved with a backpropagation algorithm. The observation is placed at $d = 10$ cm from the retina, corresponding to the location of a tumor, inside the fatty region of the GeePs-L2S testing phantom, shown in pink in Figure 3. This location choice helps us in investigating other three-dimensional algorithms such as diffraction tomography or deep learning models for volumetric reconstruction. We then compute the filtered spectrum of the induced currents at the tumor position using the Fourier diffraction theorem.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Configuration of GeePs-L2S microwave camera.</figcaption></figure>

## 2.3. Backpropagation Projection Algorithm

The forward electromagnetic problem states that the induced currents $j$ can be computed by solving the integral equation given in (1), where $G$ is Green's function; $E_t$ is the total electric field; $k(x,y,z)$ is the wavenumber at the point $(x,y,z)$, $k_{water}$ is the wavenumber of the coupling medium (water), $\nabla$ is the operator Nabla, and $I$ is the identity matrix. A single bar below quantities indicates that they are considered a vector and a double bar stands for matrices.

$$ \underline{E}_t(x, y, z) = \underline{E}_i(x, y, z) + \left[ \underline{I} + \frac{\nabla \cdot \nabla}{k_{\text{water}}^2} \right] \iiint_V j(x', y', z') G(x, y, z; x', y', z') dx' dy' dz' \quad (1) $$

$$ j(x, y, z) = \left[ k^2(x, y, z) - k_{\text{water}}^2 \right] \underline{E}_t(x, y, z) \quad (2) $$

The dipoles placed on the retina of the camera measure the vertical component of the total electric field. We deduce the field scattered by the phantom $E_{s,z}$ by subtracting the
incident field along the z-axis. Assuming that the depolarization inside the phantom is negligible, it is known that Equation (1) can be reduced to a scalar equation. In that way, the scattered field along the z-component on the retina, $E_{s,z}$ at a distance $d$ from the object, is then written as follows:

$$E_{s,z}(x, d, z) = \iiint_V j_z(x', y', z') G(x, d, z; x', y', z') dx' dy' dz' \quad (3)$$

Because of Weyl's angular spectrum expansion, Green's function can be expressed as follows:

$$G(x, d, z; x', y', z') = - \iint_{-\infty}^{+\infty} \frac{i}{8\pi^2\gamma} e^{-i(k_x(x-x') + \gamma|d-y'| + k_z(z-z'))} dk_x dk_z \quad (4)$$

$$\gamma = \sqrt{k_{water}^2 - k_x^2 - k_z^2} \quad (5)$$

If we only consider the propagative waves, Equation (5) enforces that $\gamma \in \mathbb{R}$. This condition can only be verified in the so-called visible range, a region bounded by a circle of radius $k_{water} = 438 \text{ m}^{-1}$, according to

$$k_x^2 + k_z^2 \le k_{water}^2 \quad (6)$$

Some supplementary calculations lead to the Fourier diffraction theorem (7) [31,32]. It stipulates that the Fourier transform of the induced currents $\hat{j}_z$ is linked to the filtered spectrum, limited to the visible domain, of the scattered field $\hat{E}_{s,z}$.

$$\hat{E}_{s,z}(k_x, k_z) = -\frac{i}{2\gamma}\hat{j}_z(k_x, \gamma - k_{water}, k_z)e^{-i\gamma d} \quad (7)$$

Then, we can easily deduce the backpropagation projection equation given in (8) [20], which gives an approximation of the induced currents inside the phantom.

$$j_z(x, z) = FT_{2D}^{-1} [2i\gamma \hat{E}_{s,z}(k_x, k_z)e^{i\gamma d}] \quad (8)$$

This formula integrates the induced currents of the entire 3D phantom into a single projective plane, making $|j_z|$ given in Figure 4a difficult to interpret. It also means that the projective plane includes information from the different parts of the phantom, making the tumor localization difficult. For assessing the possibility of tumor detection, we calculate the difference between the induced currents with and without the tumor Figure 4b. It appears that tumor information is present in the differential image, demonstrating the sensitivity of the method to the presence of the tumor.

<figure><img src="image_7.png" /><figcaption>Figure 4. Magnitude of the reconstructed induced current distribution (a); differential image (b). The white rectangle shows the tumor location.</figcaption></figure>
## 2.4. Spectral Filtering

Once the distribution of the scattered field on the retina of the camera has been obtained, it is possible to deduce its spectrum using fast Fourier transform, followed by the spectrum of its currents using (7). Let us introduce the Normalized Spectral Density of $\hat{j}_z$, noted NSD, given by the following:

$$NSD = 10 \log_{10} \left( \frac{|\hat{j}_z(k_x, \gamma - k_{water}, k_z)|}{\max|\hat{j}_z(k_x, \gamma - k_{water}, k_z)|} \right)^2 \quad (9)$$

The measurement is conducted in water with 64x64 dipoles spaced $\frac{\lambda_{water}}{4}$. The spectrum $\hat{E}_{s,z}$ is therefore composed of 4096 pixels, and its extent is $1751 \text{ m}^{-1}$. The Fourier diffraction theorem naturally filters out the high spatial frequencies of the current spectrum $\hat{j}_z$. The range of the spectrum, according to Equation (6), is bounded by a circle of diameter $876 \text{ m}^{-1}$, referred to as the visible range, which reduces the number of pixels to 805 pixels.

In Figure 5, the modulus of the $\hat{j}_z$ distribution of the test phantom and its NSD in decibel [dB] are represented. The two white circles represent the limit of the visible range (outer circle) and this same limit divided by two.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Spectrum (a) and NSD [dB] (b) of the induced currents. The white circles represent the visible and half-visible domains.</figcaption></figure>

The figure shows that the spectrum is almost evanescent outside the visible half-domain. Filtering spatial frequencies outside this circle should therefore not significantly degrade the image. This operation reduces the spectrum to 197 pixels.

Hence, spectral filtering offers several advantages. Firstly, spectrum truncation reduces the number of unknowns to be determined, thereby reducing the complexity of the problem posed, from 4096 to 197 pixels. Secondly, filtering the invisible spectrum and beyond removes the information contained in evanescent waves, which are not accessible due to measurement noises. This natural truncation in the back-propagation algorithm can be assimilated to a regularization technique, making the spectral technique less sensitive to noise addition.

The real and imaginary parts of the induced current distributions are reconstructed for different signal-to-noise ratios and shown in Figure 6. Additive White Gaussian Noise (AWGN) is added on $j_z$ before filtering. As expected, the method is very robust, and a signal-to-noise ratio (SNR) of 10 dB has a limited effect on the reconstruction. The spectral version of the induced currents—in the visible domain divided by two—will be taken in the following as the input of the U-NETs.
<figure><img src="image_9.png" /><figcaption>Figure 6. Effect of noise on filtered $j_z$ for various SNRs. Top: real part. Bottom: imaginary part. Noiseless (a); SNR$_{20dB}$ (b); SNR$_{10dB}$ (c).</figcaption></figure>

## 2.5. From Qualitative to Quantitative Imaging

Deep learning techniques are online optimization processes. To optimize the spectrum of induced currents, a reference must be chosen. The simplest and most logical approach would be to conduct a qualitative-to-qualitative optimization to optimize the approximate induced currents with respect to the theoretical currents calculated using Equation (2). However, from a computational viewpoint, this choice is not satisfactory. Considering the variations in the field inside the phantom complicates the process of optimizing the induced currents map. The use of the complex contrast C, given in Equation (10), seems to be a better choice, as it is strongly linked to different parts of the breast and more challenging.

$$C(x, z) = (\epsilon'(x, z) - \epsilon'_{water}) - i(\epsilon''(x, z) - \epsilon''_{water}) \quad (10)$$

$$\epsilon''_{water} = \frac{\sigma_{water}}{\omega\epsilon_0} \quad (11)$$

In a first step, the data must be standardized and normalized, before training the U-NETs. The choice of scaling method is linked to the values of the parameters to be optimized. Here, the imaginary part of the dielectric contrast and $j_z$ currents is small compared to their real part, which can lead to balance problems when the error is backpropagated to the neural network. Min-max normalization ensures that real and imaginary parts follow equivalent dynamics. For this reason, $j_z$ and the reference contrast C are rescaled to the interval [0, 1] using the standard min-max scaling formula on the real and imaginary parts.

$$j_{min-max} = \frac{\Re(j_z) - \min(\Re(j_z))}{\max(\Re(j_z)) - \min(\Re(j_z))} + i \frac{\Im(j_z) - \min(\Im(j_z))}{\max(\Im(j_z)) - \min(\Im(j_z))} \quad (12)$$

$$C_{min-max} = \frac{\Re(C) - \min(\Re(C))}{\max(\Re(C)) - \min(\Re(C))} + i \frac{\Im(C) - \min(\Im(C))}{\max(\Im(C)) - \min(\Im(C))} \quad (13)$$

To avoid the exploding gradient issue during the U-NET backpropagation process and because deep learning models are pseudo-probabilistic models, the input data should be close to 0. It is therefore desirable for the real and imaginary parts of the spectrum to lie within the interval [-1, 1]. Hence, a simple normalization is applied to generate the input $\hat{j}_{norm}$ and reference spectra $\hat{C}_{norm}$ used in the optimization process. The applied operation
preserves the Fourier relationship as it only modifies the magnitude and not the phase of the quantities:

$$ \hat{j}_{norm}(k_x, k_z) = \frac{\hat{j}_{min-max}(k_x, k_z)}{\max(|\hat{j}_{min-max}(k_x, k_z)|)} \qquad (14) $$

$$ \hat{C}_{norm}(k_x, k_z) = \frac{\hat{C}_{min-max}(k_x, k_z)}{\max(|\hat{C}_{min-max}(k_x, k_z)|)} \qquad (15) $$

The aim of the U-NET pair is to perform a full transformation of the input $\hat{j}_{norm}$ into the output contrast $\hat{C}_{net}$, using $\hat{C}_{norm}$ as a reference. The input, reference, and output of the U-NETs for the test example are shown in Figure 7.

<figure><img src="image_10.png" /><figcaption>Figure 7. Input, reference, and output of the U-NETs for the testing spectrum. The white circle radius is half the visible domain. Top: real part. Bottom: imaginary part. $\hat{j}_{norm}$ (a); $\hat{C}_{norm}$ (b); $\hat{C}_{net}$ (c).</figcaption></figure>

In a second step, the denormalization and destandardization process required to compute $\epsilon_{rnet}$ and $\sigma_{net}$ values from the $C_{net}$ contrast reconstructed by the U-NETs is conducted as follows. We introduce $M_{\Re}$ and $M_{\Im}$, the mean of the real or imaginary value of all the pixels outside the phantom on the reconstructed image (in water), respectively. As the dielectric properties of water are known, the operations given in (16) and (17) calibrate the reconstructed values, $\epsilon_{cal}$ and $\sigma_{cal}$, with respect to the water, making the process less sensitive to oscillations caused by spectral reconstruction errors.

$$ \epsilon_{cal}(x, z) = \left( \frac{\Re(C_{net}(x, z))}{M_{\Re}} - 1 \right) \epsilon_{water} \qquad (16) $$

$$ \sigma_{cal}(x, z) = \left( \frac{\Im(C_{net}(x, z))}{M_{\Im}} - 1 \right) \sigma_{water} \qquad (17) $$

By using min-max scaling a second time, we can rescale $\epsilon_{cal}$ and $\sigma_{cal}$ in the interval $[min(|\epsilon_{cal}|), max(|\epsilon_{cal}|)]$ and $[min(|\sigma_{cal}|), max(|\sigma_{cal}|)]$ to enforce that the maximum and the minimum values reconstructed by the U-NETs are positive. That way, Equations (18) and (19) give us the dielectric properties $\epsilon_{rnet}$ and $\sigma_{net}$ reconstructed by the U-NETs.

$$ \epsilon_{rnet} = min(|\epsilon_{cal}|) + (max(|\epsilon_{cal}|) - min(|\epsilon_{cal}|)). \left[ \frac{\epsilon_{cal} - min(\epsilon_{cal})}{max(\epsilon_{cal}) - min(\epsilon_{cal})} \right] \qquad (18) $$

$$ \sigma_{net} = max(|\sigma_{cal}|) + (min(|\sigma_{cal}|) - max(|\sigma_{cal}|)). \left[ \frac{\sigma_{cal} - min(\sigma_{cal})}{max(\sigma_{cal}) - min(\sigma_{cal})} \right] \qquad (19) $$
Finally, the application of the contour adds prior information by enforcing the constant values of $\epsilon_{rwater}$ and $\sigma_{water}$ outside the phantom.

## 2.6. U-NETs for Spectral Reconstruction of Dielectric Contrast

The model was fully developed hand-made in Python, using functions from the open-source PyTorch library [33]. For spectrum enhancement, we chose to use a pair of identical U-NETs [34], with the architecture illustrated in Figure 8. One U-NET is designed to optimize the real part of the spectrum $\hat{C}_{net}$, while the second U-NET focuses on the imaginary part. The spectral filter is applied both on the spectral input and output of the U-NET pair. That way, the U-NETs exclusively optimizes the filtered part and performs a proper pixel-to-pixel transformation of the spectrum.

<figure><img src="image_11.png" /><figcaption>Figure 8. U-NET for image enhancement through spectrum transformation.</figcaption></figure>

To perform the reconstructions, the spectra are grouped into batches of three, containing a total of $n = 17 \times 17 \times 3$ parameters. In the compression phase, we implement a series of filters with decreasing sizes, using zero-padding adjustments to ensure gradual compression of the spectrum without employing pooling operations. This allows the information to be compressed through convolutions, reducing the dimensions of each feature map from $17 \times 17$ to $7 \times 7$. The various applied filters gradually generate an increasing number of feature maps, gathering the statistical information needed to modify the spectrum.

During the decompression phase, the filter dimensions are unchanged, and the padding is gradually increased to the initial spectrum size at the end of the network. The features of the U-NETs architecture are presented in Table 2.

<table><caption>Table 2. Features of U-NETs.</caption><thead><tr><th rowspan="2">Step</th><th colspan="4">Compression</th><th colspan="2">Linear Layer</th><th colspan="3">Decompression</th></tr><tr><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th><th>6</th><th>7</th><th>8</th><th>9</th></tr></thead><tbody><tr><td>Feature map size</td><td>13 × 13</td><td>11 × 11</td><td>9 × 9</td><td>7 × 7</td><td>7 × 7</td><td>7 × 7</td><td>9 × 9</td><td>11 × 11</td><td>13 × 13</td></tr><tr><td>Feature map number</td><td>32</td><td>64</td><td>128</td><td>256</td><td>-</td><td>256 + 256</td><td>128 + 128</td><td>64 + 64</td><td>32 + 32</td></tr><tr><td>Filter size</td><td>9 × 9</td><td>7 × 7</td><td>5 × 5</td><td>3 × 3</td><td>-</td><td>3 × 3</td><td>3 × 3</td><td>3 × 3</td><td>3 × 3</td></tr><tr><td>Padding</td><td>2</td><td>2</td><td>1</td><td>0</td><td>-</td><td>2</td><td>2</td><td>2</td><td>3</td></tr></tbody></table>

Before each convolution, batch normalization is applied to speed up the learning process, and the ReLu function is introduced to add non-linearity to the U-NETs. The transition from compression to decompression is a linear layer.

The Adam optimizer [35] minimizes the Weighted Mean Absolute Percentage Error (WMAPE) loss for the real and imaginary components of the spectrum, given by Equa-
tions (20) and (21), respectively. The optimizer parameters are set to their default values: learning rate ($lr$) = 0.001 and inertial coefficients ($\beta1, \beta2$) = (0.9, 0.999).

$$ loss_{\mathcal{R}} = \frac{\sum_{i=1}^{n} |\Re(\hat{C}_{net}) - \Re(\hat{C}_{norm})|}{\sum_{i=1}^{n} |\Re(\hat{C}_{norm})|} \qquad (20) $$

$$ loss_{\mathcal{F}} = \frac{\sum_{i=1}^{n} |\Im(\hat{C}_{net}) - \Im(\hat{C}_{norm})|}{\sum_{i=1}^{n} |\Im(\hat{C}_{norm})|} \qquad (21) $$

The U-NETs were trained over 30 epochs. The computation time for the learning phase was less than an hour and a half, on an NVIDIA RTX A2000 8GB. The computed loss in Figure 9 shows the WMAPE loss function as a function of the number of epochs. The loss values displayed are average values, taking into account all the batches processed. A similar convergence is observed for the learning and validation process, for both real and imaginary parts.

<figure><img src="image_12.png" /><figcaption>Figure 9. U-NETs: WMAPE loss on validation and training sets versus number of epochs.</figcaption></figure>

# 3. Results

## 3.1. Validation Results

The validation set of 359 spectra contains 20% of the data. The U-NETs were trained with noise-free data in order to assess their effectiveness on well-controlled data.

As the models are built from various combinations of fat, skin, glandular tissues, and tumors, the validation dataset is correlated with the training dataset, which generates overfitting issues in our model. It should also be noted that, since the database models are three-dimensional, tumors are not necessarily present in the central plane of observation. The results are shown in Figure 10. Six models are considered, including three with tumors in the central slice. Four images are presented for each model; the first column corresponds to the exact profile of the dielectric constant (top) and conductivity (bottom), and the second column gives the corresponding profiles reconstructed by the U-NETs.

The network generally produces a high-quality image of permittivity and conductivity distributions. The shapes of the various tissues—skin, fat, and glandular—are well reproduced, with values close to the reference ones. The result is more questionable when the model contains tumors, or more generally for “unusual” situations, leading to failure cases. In the examples shown in Figure 10, we have intentionally selected failure situations. As we can see, the tumor appears in the right place in tumor Models Figure 10a,b, but not in Model Figure 10f. The reason is probably linked to the training dataset, since 90% of it is generated from tumor-free breast models in the central breast slice. Similarly, as shown in Figure 10e, the neural network faces issues when reconstructing glandular-free models, creating artifacts inside the breast, as only 5% of the generated slices are devoid of glandular tissues.
<figure><img src="image_13.png" /><figcaption>**Figure 10.** Validation examples: phantoms with tumor (a,b); phantoms without tumor (c,d); failure cases (e,f). Color scale is fixed to [1,73] on $\epsilon_r$ and [0,2] [S/m] on $\sigma$ to fit all values distant of ±5% from standard ones.</figcaption></figure>

Figure 11 shows the structural similarity (SSIM) [36,37] and the peak signal-to-noise ratio (PSNR) that [37] computed with torchmetrics default parameters, over the entire validation dataset of 359 spectra. As simulations are run in batches of 3, the metrics are plotted versus the 119 batches. The yellow and blue curves correspond to an evaluation of image quality at the input and output of the U-NETs, respectively. As can be seen, the overall improvement brought about by the U-NETs is major, with an average SSIM of 0.97 and a gain of over 20 dB according to the PSNR metric, and on both the real and imaginary parts of the spectrum. In comparison with Figure 11d,c, the imaginary part of the spectrum appears slightly less well optimized than the real part of the spectrum.

To further improve these results, the training and validation datasets need to be expanded. This can be achieved easily and significantly by adding tumors in the center and diversifying the geometric shapes of the breast models, as this will have a large impact on the total number of possible combinations. Also, with the enlargement of the database, a better decorrelation between the training and the testing data could be achieved. These results, obtained with a limited database and a qualitative input image of very low resolution, are very encouraging as preliminary results.
<figure><img src="image_14.png" /><figcaption>**Figure 11.** SSIM (a,b) and PSNR (c,d) metrics on $\hat{j}_{norm}$ input, shown in orange, and $\hat{C}_{net}$ output shown in blue. Dotted lines represent average metric over validation dataset. Real part (a,c); imaginary part (b,d).</figcaption></figure>

## 3.2. Test Results

For now, the testing set is limited to the 32 slices of the GeePs-L2S breast phantom. Figure 12 shows the dielectric constant and conductivity profiles provided by the network for the slice passing through the tumor and for different signal-to-noise ratios. The effect of noise is almost negligible, confirming the method's low sensitivity to noise for the reasons outlined in Section 3.

<figure><img src="image_15.png" /><figcaption>**Figure 12.** Test example. **Top:** $\epsilon_r$. **Bottom:** $\sigma$ [S/m]. From left to right: reference; noiseless reconstruction; SNR 20 dB reconstruction; SNR 10 dB reconstruction. The white rectangle shows the tumor location.</figcaption></figure>

An improvement can be seen both in terms of spatial resolution and of what the image represents. The U-NET improves spectrum quality by transforming the low-spatial-resolution qualitative image provided by the backpropagation algorithm into a higher-spatial-resolution quantitative image, providing access to the phantom's dielectric properties, but not those of the tumor.

As before, SSIM and PSNR values were calculated on the 32 spectra of the test database and depicted in Figure 13. As expected, given the size of the database, the network output
(blue curves) shows a lower SSIM than in the previous case, a value between 0.8 and 0.85. This result aligns with the lower-quality image reconstruction depicted in Figure 12. This reflects the neural network's difficulty in reconstructing a novel tumor in a new location. However, if we compare these values with those obtained with the network's input quantities (yellow curves in the figures), the improvement remains significant for both real and imaginary parts. U-NETs bring an average improvement of 10 dB on the PSNR of the reconstructed spectrum.

<figure><img src="image_16.png" /><figcaption>**Figure 13.** SSIM (a,b) and PSNR (c,d) metrics on $\hat{j}_{norm}$ input, shown in orange, and $\hat{C}_{net}$ output shown in blue. Dotted lines represent average metric over test dataset. Real part (a,c); imaginary part (b,d).</figcaption></figure>

## 4. Discussion

The main objective of this paper was to develop a combination of U-NETs with a spectral method to produce real-time quantitative breast imaging. The preliminary results presented in the previous section validate the method and confirm its potential in terms of efficiency with an SSIM of 0.97 in validation and 0.83 in testing and robustness to a low SNR of 10 dB. A second achievement was the development of an in-house, realistic, diverse, low-cost, and easy-to-use 3D database. To our knowledge, both developments are new.

However, this work is a first step. Indeed, as research on model architectures such as CNN, U-NET, ANN, and CVCNN remain important, the most essential aspect of deep learning remains the quality and quantity of the database and data pre/post-processing, as the neural network parameters are optimized from the available dataset. Unfortunately, the number of breast phantoms currently available worldwide remains limited. In the literature, studies that clearly avoid data leakage by using uncorrelated testing and training sets, achieving deep learning generalization results during image enhancement, dispose either of a large amount of automatically generated data [15] or of information from other modalities [12]. This is not the case here. Consequently, overfitting is mandatory with such a small database, and the correlation between training and validation sets is ensured by the combinatorial nature of the training/validation set. This is clearly an obstacle to
generalization, and techniques such as regularization or dropout are useless when facing a
lack of data.

To solve this issue, we propose in this study to begin with a small in-house database, with the idea of optimizing and extending it in the future. The idea is to study both the reconstruction process and the data themselves in an optimized way. This would enable us to better understand the behavior of the neural network as a function of the database, i.e., its size, processing, and potential biases. The database would become a customizable input, and the optimization process would be seen as a unified optimization of the model and data.

This database must be realistic, diversified, and easy to develop. This is why we propose a combination of anthropomorphic breast cavities, derived from modifications made to a realistic phantom, generated from an MRI scanner [13]. Working with phantoms not only frees us from legal and ethical issues, but also allows us to control the entire database development program. Moreover, the existence of printable twins of constructed cavities can be used for experiments or to calibrate measurements.

However, conducting research in this way comes at a cost both in terms of development time and effort. This is why this study focuses primarily on assessing the validity of the method rather than its performance.

Does the combination of U-NETs with a spectral method of very low computational cost, compared to classical optimization methods—DBIM, CSI, GN, etc.—yield promising results? The answer is positive, if we refer to the images and associated metrics presented in the previous section.

**5. Conclusions**

The deep learning approach implemented in this paper seems well suited to the complex problem of retrieving the dielectric properties of an anthropomorphic breast model from the spectrum of induced currents approximated by a sufficiently well-trained neural network. Preliminary results are impressive, with significant improvements in both image contrast and spatial resolution. However, this efficiency would be even more effective with a well-dimensioned database whose development considers the behavior of the neural network in the face of statistical biases and correlations that the database may introduce. Although this database is not currently available, the choice of working in the spectral domain with the backpropagation algorithm seems attractive. The proposed method is a real-time process, robust to noise, and easy to implement for the purpose of evaluating the effectiveness of deep learning. It is potentially a good candidate for future imaging systems.

**Author Contributions:** Conceptualization, H.R. and N.J.; Software, A.D.; Validation, A.D.; Writing—original draft, A.D.; Writing—review & editing, H.R. and N.J.; Supervision, H.R. and N.J. All authors have read and agreed to the published version of the manuscript.

**Funding:** The authors would like to thank the Institut Universitaire d'Ingénierie en Santé (IUIS) of Sorbonne Université for financing their doctoral studies.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The original contributions presented in this study are included in the article. Further inquiries can be directed to the corresponding author.

**Conflicts of Interest:** The authors declare no conflicts of interest.
## References

1. Berg, P.M.V.D.; Kleinman, R.E. A Contrast Source Inversion Method. *Inverse Probl.* **1997**, *13*, 1607–1620. https://doi.org/10.1088/0266-5611/13/6/013.

2. Qin, Y. Early Breast Anomalies Detection with Microwave and Ultrasound Modalities. Ph.D. Thesis, Université Paris-Saclay, Gif-sur-Yvette, France, 2021.

3. Mariano, V.; Vasquez, J.A.T.; Rodriguez-Duarte, D.O.; Vipiana, F. Field-Based Discretization of the 3-D Contrast Source Inversion Method Applied to Brain Stroke Microwave Imaging. *IEEE J. Electromagn. RF Microwaves Med. Biol.* **2024**, *8*, 290–297. https://doi.org/10.1109/JERM.2024.3414196.

4. Chew, W.; Wang, Y. Reconstruction of Two-Dimensional Permittivity Distribution Using the Distorted Born Iterative Method. *IEEE Trans. Med. Imaging* **1990**, *9*, 218–225. https://doi.org/10.1109/42.56334.

5. Neira, L.M.; Van Veen, B.D.; Hagness, S.C. High-Resolution Microwave Breast Imaging Using a 3-D Inverse Scattering Algorithm with a Variable-Strength Spatial Prior Constraint. *IEEE Trans. Antennas Propag.* **2017**, *65*, 6002–6014. https://doi.org/10.1109/TAP.2017.2751668.

6. Joachimowicz, N.; Pichot, C.; Hugonin, J. Inverse Scattering: An Iterative Numerical Method for Electromagnetic Imaging. *IEEE Trans. Antennas Propag.* **1991**, *39*, 1742–1753. https://doi.org/10.1109/8.121595.

7. Borzooei, S.; Tournier, P.H.; Dolean, V.; Pichot, C.; Joachimowicz, N.; Roussel, H.; Migliaccio, C. Numerical Modeling for Shoulder Injury Detection Using Microwave Imaging. *IEEE J. Electromagn. RF Microwaves Med. Biol.* **2024**, *8*, 282–289. https://doi.org/10.1109/JERM.2024.3411799.

8. 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.* **2019**, *67*, 1819–1825. https://doi.org/10.1109/TAP.2018.2885437.

9. Huang, Y.; Song, R.; Xu, K.; Ye, X.; Li, C.; Chen, X. Deep Learning-Based Inverse Scattering With Structural Similarity Loss Functions. *IEEE Sens. J.* **2021**, *21*, 4900–4907. https://doi.org/10.1109/JSEN.2020.3030321.

10. Yago Ruiz, Á.; Cavagnaro, M.; Crocco, L. A Physics-Assisted Deep Learning Microwave Imaging Framework for Real-Time Shape Reconstruction of Unknown Targets. *IEEE Trans. Antennas Propag.* **2022**, *70*, 6184–6194. https://doi.org/10.1109/TAP.2022.3162320.

11. Zhang, Y.; Lambert, M.; Fraysse, A.; Lesselier, D. Unrolled Convolutional Neural Network for Full-Wave Inverse Scattering. *IEEE Trans. Antennas Propag.* **2023**, *71*, 947. https://doi.org/10.1109/TAP.2022.3216999.

12. Mojabi, P.; Hughson, M.; Khoshdel, V.; Jeffrey, I.; LoVetri, J. CNN for Compressibility to Permittivity Mapping for Combined Ultrasound Microwave Breast Imaging. *IEEE J. Multiscale Multiphys. Comput. Tech.* **2021**, *6*, 62–72. https://doi.org/10.1109/JMMCT.2021.3076827.

13. 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. https://doi.org/10.1109/LAWP.2012.2236293.

14. Omer, M.; Fear, E. Anthropomorphic Breast Model Repository for Research and Development of Microwave Breast Imaging Technologies. *Sci. Data* **2018**, *5*, 180257. https://doi.org/10.1038/sdata.2018.257.

15. 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.

16. Costanzo, S.; Flores, A.; Buonanno, G. Fast and Accurate CNN-Based Machine Learning Approach for Microwave Medical Imaging in Cancer Detection. *IEEE Access* **2023**, *11*, 66063–66075. https://doi.org/10.1109/ACCESS.2023.3291076.

17. Khoshdel, V.; Asefi, M.; Ashraf, A.; LoVetri, J. Full 3D Microwave Breast Imaging Using a Deep-Learning Technique. *J. Imaging* **2020**, *6*, 80. https://doi.org/10.3390/jimaging6080080.

18. Borghouts, M.; Ambrosanio, M.; Franceschini, S.; Autorino, M.M.; Pascazio, V.; Baselice, F. Microwave Breast Sensing via Deep Learning for Tumor Spatial Localization by Probability Maps. *Bioengineering* **2023**, *10*, 1153. https://doi.org/10.3390/bioengineering10101153.

19. Costanzo, S.; Flores, A. Reinforcement Learning to Enhanced Microwave Imaging for Accurate Tumor Detection in Breast Images. In Proceedings of the 2024 IEEE International Conference on Metrology for eXtended Reality, Artificial Intelligence and Neural Engineering (MetroXRAINE), St Albans, UK, 21–23 October 2024; pp. 101–106. https://doi.org/10.1109/MetroXRAINE62247.2024.10796628.

20. Bolomey, J.C.; Pichot, C. Microwave Tomography: From Theory to Practical Imaging Systems. *Int. J. Imaging Syst. Technol.* **1990**, *2*, 144–156. https://doi.org/10.1002/ima.1850020210.

21. Henriksson, T. Contribution to Quantitative Microwave Imaging Techniques for Biomedical Applications. Ph.D. Thesis, Mälardalen University, Paris, France, 2009.

22. Chen, T.; Su, X.; Li, H.; Li, S.; Liu, J.; Zhang, G.; Feng, X.; Wang, S.; Liu, X.; Wang, Y.; et al. Learning a Fully Connected U-Net for Spectrum Reconstruction of Fourier Transform Imaging Spectrometers. *Remote Sens.* **2022**, *14*, 900. https://doi.org/10.3390/rs14040900.
23. Joachimowicz, N.; Duchêne, B.; Conessa, C.; Meyer, O. Anthropomorphic Breast and Head Phantoms for Microwave Imaging. *Diagnostics* **2018**, *8*, 85. https://doi.org/10.3390/diagnostics8040085.

24. Abedi, S. Standard Phantoms for EM Device Testing. Ph.D. Thesis, Sorbonne Université, Paris, France, 2022.

25. Andreuccetti, D.; Fossi, R.; Petrucci, C. An Internet Resource for the Calculation of the Dielectric Properties of Body Tissues in the Frequency Range 10 Hz–100 GHz. 1997. Available online: http://niremf.ifac.cnr.it/tissprop/ (accessed on 11 April 2025).

26. Joachimowicz, N.; Roussel, H. Standard Uwb Phantom ELECTromagnetic RECIPES for Microwave Imaging. In Proceedings of the 2023 IEEE Conference on Antenna Measurements and Applications (CAMA), Genoa, Italy, 15–17 November 2023; pp. 227–229, Available online: http://applis.iut.univ-paris-diderot.fr/CNRS/ (accessed on 11 April 2025). https://doi.org/10.1109/CAMA57522.2023.10352905.

27. Borzooei, S.; Tournier, P.H.; Dolean, V.; Migliaccio, C. Microwave Digital Twin Prototype for Shoulder Injury Detection. *Sensors* **2024**, *24*, 6663. https://doi.org/10.3390/s24206663.

28. Semenov, S.Y.; Svenson, R.H.; Posukh, V.G.; Nazarov, A.G.; Sizov, Y.E.; Bulyshev, A.E.; Souvorov, A.E.; Chen, W.; Kasell, J.; Tatsis, G.P. Dielectrical Spectroscopy of Canine Myocardium during Acute Ischemia and Hypoxia at Frequency Spectrum from 100 kHz to 6 GHz. *IEEE Trans. Med. Imaging* **2002**, *21*, 703–707. https://doi.org/10.1109/TMI.2002.800590.

29. Tournier, P.H.; Bonazzoli, M.; Dolean, V.; Rapetti, F.; Hecht, F.; Nataf, F.; Aliferis, I.; El Kanfoud, I.; Migliaccio, C.; de Buhan, M.; et al. Numerical Modeling and High-Speed Parallel Computing: New Perspectives on Tomographic Microwave Imaging for Brain Stroke Detection and Monitoring. *IEEE Antennas Propag. Mag.* **2017**, *59*, 98–110. https://doi.org/10.1109/MAP.2017.2731199.

30. WIPL-D Pro. 3D Electromagnetic Solver; WIPL-D d.o.o.: Belgrade, Serbia, 2011.

31. Devaney, A.J. A Filtered Backpropagation Algorithm for Diffraction Tomography. *Ultrason. Imaging* **1982**, *4*, 336–350. https://doi.org/10.1016/0161-7346(82)90017-7.

32. Kirisits, C.; Quellmalz, M.; Setterqvist, E. Generalized Fourier Diffraction Theorem and Filtered Backpropagation for Tomographic Reconstruction. *arXiv* **2024**, arXiv:2407.0179. https://doi.org/10.48550/arXiv.2407.01793.

33. Paszke, A.; Gross, S.; Massa, F.; Lerer, A.; Bradbury, J.; Chanan, G.; Killeen, T.; Lin, Z.; Gimelshein, N.; Antiga, L.; et al. PyTorch: An Imperative Style, High-Performance Deep Learning Library. *arXiv* **2019**, arXiv:1912.01703. https://doi.org/10.48550/arXiv.1912.01703.

34. Ronneberger, O.; Fischer, P.; Brox, T. U-Net: Convolutional Networks for Biomedical Image Segmentation *arXiv* **2015**, arXiv:1505.04597. https://doi.org/arXiv:1505.04597.

35. Kingma, D.P.; Ba, J. Adam: A Method for Stochastic Optimization. *arXiv* **2017**, arXiv:1412.6980. https://doi.org/10.48550/arXiv.1412.6980.

36. Zhou Wang,; Bovik, A.; Sheikh, H.; Simoncelli, E. Image Quality Assessment: From Error Visibility to Structural Similarity. *IEEE Trans. Image Process.* **2004**, *13*, 600–612. https://doi.org/10.1109/TIP.2003.819861.

37. Detlefsen, N.S.; Borovec, J.; Schock, J.; Jha, A.H.; Koker, T.; Di Liello, L.; Stancl, D.; Quan, C.; Grechkin, M.; Falcon, W. TorchMetrics—Measuring Reproducibility in PyTorch. *J. Open Source Softw.* **2022**, *7*, 4101. https://doi.org/10.21105/joss.04101.

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.