PhaselessMicrowaveTomographyAssessment forBreast Imaging: Preliminary Results

In the present work, a phaseless approach for microwave imaging applications is presented. (e proposed solution strategy is based on the formulation of the scattering phenomena in terms of contrast source, while no phase-recovery stage is involved into the numerical procedure, thus providing a phaseless single-step resolution method. (e image recovering potentialities of the discussed method are numerically validated by successfully distinguishing different tissues of a slice breast model, with a tumor located wherein. (e above preliminary assessment encourages the adoption of the proposed solution in the framework of biomedical imaging.


Introduction
Recent trends in biomedical imaging have raised some concerns about the safety of the imaging apparatus, as well as about the overall system costs. Microwave tomography (MWT) can be regarded as a safe, cost-effective, noninvasive supplement to the widely adopted imaging techniques. However, limited penetration depth and relatively low resolution put some constraints on a large-scale application of MWT. A major reason to move towards the microwave spectrum is given by the possibility to identify possible pathological alterations on biological tissues according to their water content. Furthermore, the use of nonionizing radiation in MWT techniques represents one of the main advantages as compared to computed tomography (CT) and positron emission tomography (PET). In this work, according to the preliminary assessment presented in [1], a Phaseless Contrast Source Inversion (P-CSI) method is successfully implemented for breast tissue reconstruction. Specifically, starting from the contrast source (CS) formulation described in [2,3], the inverse scattering problem is solved with no linearization procedure, but recasting it into an iterative optimization problem, where the two problems unknowns, namely, the contrast source and the dielectric contrast, are alternatively updated according to a conjugate gradient scheme. e inversion procedure is performed by exploiting the amplitude-only data of the measured total field, locally defined as the sum of the incident field and the scattered field, the former obtained as a base-line measurement in the absence of the object under test (OUT), and the latter due to the interaction of the incident field with the OUT. e full-data information of the incident field inside the imaging domain, obtained through a wave expansion starting from the base-line measurement, is also required into the reconstruction process.

Phaseless Inverse Scattering Formulation
Let us consider a 2D tomography problem, aiming at retrieving the shape, the location and the dielectric properties of a generic OUT, hereby denoted as B. A TM-polarized incident field is assumed, and cylindrical targets are investigated.
In the whole scenario, a magnetic permittivity equal to that of free space is assumed, such as in a typical biomedical scenario, due to the nonmagnetic properties of the biological tissues. erefore, the electromagnetic properties are fully described uniquely by the dielectric permittivity values. A circular acquisition setup is considered, where the measurement points are displaced on the acquisition curve S, as shown in the general scheme of Figure 1. It describes a multistatic and multiview setup, where the transmitter location is alternatively changed, thus determining N TX angles of incidence and resulting into a number of N R ·N TX measurements. e imaging domain, hereby indicated as D, fully contains the unknown OUT.
e inverse scattering problem is governed by the electrical field integral equations (EFIEs), known in the literature as data equation and state equation, and obtained as solution of the Helmholtz equation applied to the aforementioned setup [4]: In the above expressions, E s represents the scattered field, while E t , E i give the total and the incident field, respectively; G denotes the Green's function of the background medium, while parameter χ gives the contrast function, normalized with respect to the permittivity of the background medium: e considered approach is mathematically arranged into a local optimization problem, wherein a proper cost function is to be defined. e solution of the scattering problem is obtained as an iterative interplay of the EFIEs expressed by equations (1) and (2), which is formulated by the following cost function: with In the above expressions, parameters α S-D denote proper normalization factors. e expressions shown in equation (5) represent the discrepancy between the measured total field in S, indicated as f, and the relative value according to the model defined into equation (1), while equation (6) considers the error affected by equation (2). erefore, the cost function F combines the cumulative difference between the analytical data and the measured field for every illumination scenario v. e incident field in the OUT is determined starting from the incident electric field measurements performed along the acquisition curve S, by expressing the incident field as the sum of a limited number N of Hankel functions of the second kind [5,6], so that the following relation holds true: where ρ and ϑ represent the polar coordinates. e truncation number N is determined by evaluating the variation of the Hankel function with respect to N, for a fixed parameter ρ equal to the radius R of the sampling curve, as shown in Figure 2. e computed c ] coefficients are subsequently implemented in the evaluation of the incident field inside the OUT [7][8][9]. e iterative process starts by updating the contrast source ω v,n , through the adoption of a nonlinear conjugate gradient scheme [10,11], while the contrast function χ n , after the evaluation of the current value of the total field at the n-th step, E t v,n � E i v + G S ω v,n , is obtained as a minimizer of the state function, indicated as F D (E t v,n , ω v,n ). Since an iterative local-based optimization problem is involved, an initial guess for the contrast source ω v,0 occurs. erefore, a proper initial value for the contrast source is considered, able to exploit all the available information provided by the measurement step. e initial guess is obtained by applying the steepest descent method to equation (5), thus resulting into the following expression: where β is the step size, while the symbol ( * ) indicates the adjoint operator. According to equation (8), all the available data are fully exploited in the initialization of the contrast source, with no arbitrary assumption on the unknown phase of the total field, which is instead proposed in [10,11]. Furthermore, no a priori information about the dielectric contrast inside D is required, as compared to the approach outlined in [12].

Human Breast Modelling and Imaging Results
To perform a numerical validation of the proposed inverse strategy, a proper microwave imaging setup is simulated on COMSOL Multiphysics ® [13], according to the proposed reconstruction method. e simulated setup is shown in Figure 3. It consists of a multistatic measurement setup, where N TX point current source is uniformly arranged in a circular fashion, by surrounding the OUT. e 2-GHz operating frequency is considered, which is the most suitable one for microwave breast imaging, as it offers a considerable trade-off between spatial resolution and penetration depth into human breast tissues [14]. A perfectly matched layer (PML) surrounding the background medium is assumed as boundary condition. e forward data computation is performed by alternatively evaluating the total field and the incident field for a fixed position of the transmitting antenna. In order to speed up the simulation process, out-ofplane line currents are implemented as TM sources. e dielectric profile of the breast is constructed similarly to the approach outlined in [14]. In particular, a simplified 3-tissue breast model is considered, whose dielectric properties are listed in Table 1. e implemented human breast model consists of a cylinder with a 45.5 mm radius, a fibroglandular region having a radius equal to

Results and Discussion
Firstly, the cancerous portion is placed inside the adipose tissue (Figure 4), where a higher dielectric contrast can   be experienced. A second case is then considered, with the tumor located inside the fibroglandular tissue ( Figure 5), thus resulting into a lower dielectric contrast. As clearly shown in the contrast function map of Figures 6 and 7, for the two scenarios, the tumor is properly localized in both cases.

Conclusions
A single-step phaseless microwave tomography method for biomedical imaging has been proposed in this work. Firstly, the analytical formulation of the inverse scattering procedure has been discussed, by introducing a novel initialization for the algorithm. Secondly, some reconstruction results of a slice breast model have been presented, where the reconstruction capabilities of the proposed methodology have been confirmed. e discussed results represent a preliminary assessment for the use of MWT as a potential low-cost tool in the field of breast cancer detection and breast imaging applications. Further developments are currently being performed in order to achieve a quantitative reconstruction in terms of dielectric properties for the under test scenario.

Data Availability
e simulation data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.