Next Article in Journal
Measurements of Night Sky Brightness in the Veneto Region of Italy: Sky Quality Meter Network Results and Differential Photometry by Digital Single Lens Reflex
Next Article in Special Issue
Image-Based RCS Estimation from Near-Field Data
Previous Article in Journal
Methods for Assessment and Monitoring of Light Pollution around Ecologically Sensitive Sites
Previous Article in Special Issue
A Minimum Rank Approach for Reduction of Environmental Noise in Near-Field Array Antenna Diagnosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Stopping Criteria for Optimization-Based Microwave Breast Imaging Algorithms

Electrical and Computer Engineering, University of Manitoba, Winnipeg, MB R3T 5V6, Canada
*
Author to whom correspondence should be addressed.
Submission received: 25 April 2019 / Accepted: 14 May 2019 / Published: 22 May 2019
(This article belongs to the Special Issue Microwave Imaging and Electromagnetic Inverse Scattering Problems)

Abstract

:
A discontinuous Galerkin formulation of the Contrast Source Inversion algorithm (DGM-CSI) for microwave breast imaging employing a frequency-cycling reconstruction technique has been modified here to include a set of automated stopping criteria that determine a suitable time to shift imaging frequencies and to globally terminate the reconstruction. Recent studies have explored the use of tissue-dependent geometrical mapping of the well-reconstructed real part to its imaginary part as initial guesses during consecutive frequency hops. This practice was shown to improve resulting 2D images of the dielectric properties of synthetic breast models, but a fixed number of iterations was used to halt DGM-CSI inversions arbitrarily. Herein, a new set of stopping conditions is introduced based on an intelligent statistical analysis of a window of past iterations of data error using the two-sample Kolmogorov-Smirnov (K-S) test. This non-parametric goodness-of-fit test establishes a pattern in the data error distribution, indicating an appropriate time to shift frequencies, or terminate the algorithm. The proposed stopping criteria are shown to improve the efficiency of DGM-CSI while yielding images of equivalent quality to assigning an often liberally overestimated number of iterations per reconstruction.

Graphical Abstract

1. Introduction

Microwave imaging (MWI) has made steady progress towards widespread clinical application in breast cancer detection and monitoring over the past decade, offering a safety advantage over established x-ray-based modalities due to its use of non-ionizing radiation and a significant cost benefit over magnetic resonance imaging (MRI). While its currently attainable spatial resolution does not match that of cancer screening tools like mammography, concurrent use of MWI of the breast in its current state could lead to decreased false positive rates among certain population groups, particularly pregnant women or those with radiographically dense (“Category C” and “Category D”) breasts, as classified by the American College of Radiology’s Breast Imaging Reporting and Data System (BI-RADS) [1,2].
Methods of utilizing MWI for biomedical applications have been investigated for almost half a century [3,4]. Some recent representative research can be found in [5]. MWI algorithms can be split into those producing qualitative images, where the reconstructed value at each pixel is a relative quantity, and those producing quantitative reconstructions of a physical property, typically the complex-valued permittivity. Examples of qualitative algorithms are those utilized in so-called radar imaging techniques [6], and those retrieving only qualitative aspects of the target, such as its support [7,8]. The interest in this paper are quantitative imaging methods that solve the full electromagnetic inverse problem.
It is well known that the electromagnetic inverse problem associated with MWI, where one attempts to reconstruct the complex-valued permittivity of an object-of-interest (OI), is ill-posed [9]. The ill-posedness stems from the non-uniqueness of the inverse-source problem where one attempts to reconstruct the electromagnetic sources responsible for a remotely measured field. In the electromagnetic inverse problem, these sources are the contrast-sources that are illuminated by the interrogating field. Various regularization techniques have been developed over several decades of research to deal with the ill-posedness of the problem [10]. The electromagnetic inverse problem is not only ill-posed but also non-linear with repect to the two unknowns within the inaccessible imaging domain, the electromagnetic field and the permittivity. The Contrast Source Inversion (CSI) technique effectively linearizes the problem by casting the data-error norm in terms of contrast-source variables, which vary with the interrogating field, and regularization is performed by introducing the so-called Maxwellian regularizer which is written in terms of the contrast-sources and the contrast [11]. Much research has been performed on advancing the CSI algorithm since it was first reported, including the addition of a multiplicative regularizer [12].
For biomedical imaging applications, which are of interest in the context of the present work, the incorporation of discretized numerical inversion models into CSI that are based on either the finite-element method (FEM) or the discontinuous Galerkin method (DGM) forward solvers have provided much flexiblity and several advantages for the overall inversion process [13,14,15]. Most importantly, these forward solvers have allowed one to incorporate prior information into the inversion model that has the effect of regularizing the inversion, enabling high-contrast reconstructions such as are required for breast imaging [16,17,18,19]. The flexibility of these partial differential equation (PDE) based solvers has also allowed advancements in data-acquisition systems utilized to acquire scattered-field data for MWI [20,21]. Other potential advancements in MWI for breast tumor detection include the use of contrast agents [22].
Recent work in multi-frequency imaging has shown promise in increasing the spatial resolution and correspondingly, the amount of distinguishable anatomic detail in MWI reconstructions. For instance, while the use of “frequency hopping” has long been shown to be an effective means of obtaining images from high-frequency data by feeding low-frequency solutions as “initial guesses” into inversion algorithms [23], a so-called “frequency cycling” and “tissue-dependent mapping” framework has been adopted in recent studies to perform multi-frequency imaging demonstrating significant enhancement in the quality of resulting 2D images [24,25]. In those works, an arbitrary fixed number of iterations was used as the stopping condition for the imaging algorithm, but it was noted that further exploration into more intelligent stopping criteria was justified, particularly based on the observation that running the algorithm longer than necessary could deteriorate the imaging results. In fact, to date, there has been little exploration into optimizing the conditions to which a microwave imaging algorithm should adhere when determining the appropriate point to move on from intermediate inversions to the next frequency in a multi-frequency sequence, or when to terminate reconstructions entirely.
While monitoring the convergence of data error has been the most common objective approach to halting microwave imaging algorithms, as it is often difficult to ascertain whether a chosen error limit will undershoot, overshoot, or consistently meet the attainable image reconstruction quality for a given dataset, several early studies focusing on synthetic and experimental imaging performance fell into the same habit of simply pre-assigning an arbitrary number of iterations for the algorithm to complete, chosen for convenience or through trial and error [11,26,27,28].
Although there have been previous attempts to develop more intelligent stopping criteria for imaging algorithms in recent literature [17], a more thorough statistical analysis of the data error is undertaken here and employed as part of a logical multi-variable framework of stopping criteria for a frequency-cycling tissue-dependent mapping formulation of Contrast Source Inversion (CSI). For comparison purposes, some 2D imaging results from a recent study are included in this work to demonstrate the significant gains in efficiency and image quality granted through the incorporation of the described stopping criteria [25]. Reconstruction results are subjected to quantitative analysis using common image error metrics to provide an objective measurement of these improvements.

2. Materials and Methods

2.1. DGM-CSI Algorithm

An implementation of CSI using a high-order frequency-domain formulation of Maxwell’s curl equations has been used for all 2D breast images presented here, employing the discontinuous Galerkin method (DGM) as a forward solver. Nodal coefficients in high-order polynomial expansions represent unknown field and contrast quantities in the resulting DGM-CSI algorithm. Independent from this choice of implementation, among the benefits of CSI is its use of operators that are functions of frequency and the material properties of the background medium only, which do not change from iteration to iteration [11]. As aforementioned, further improvement in performance has been demonstrated through the use of a tissue-dependent mapping process and the practice of cycling back to low-frequency reconstructions [25], and a brief overview of this process is provided in Section 2.2. Further details specific to DGM-CSI can be found in [14,29], but a summary of the relevant quantities and definitions involved in the CSI cost functional is included here, since it pertains to the discussion of data and domain error that follows.
The microwave imaging problem in Figure 1 depicts an object of interest (OI) with an unknown dielectric contrast χ ( r ) , a function of position vector r , typically defined as
χ ( r ) = ε ^ r ( r ) ε ^ b ε ^ b ,
where ε ^ r is the complex relative permittivity and ε ^ b is the complex background permittivity. In breast imaging applications, the background permittivity is usually represented as the homogeneous immersion medium in which the breast is submerged and the imaging domain D is an area typically chosen to be contained by the outer skin layer of the breast. An array of t = 1 , 2 , , T transmitters is used to generate probing incident fields; the resulting scattered fields are measured by an array of R receivers that make up a discrete measurement surface S .
It is assumed that the reader is familiar with the standard CSI algorithm, and only the aspects important to the proposed stopping criteria will be reviewed. The CSI algorithm solves for the contrast χ ( r ) of breast tissue within D from this scattered field data by minimizing the following two-part cost functional:
F CSI ( χ , w t ) = F S ( w t ) + F D ( χ , w t ) ,
where both the contrast χ ( r ) and the transmitter-dependent contrast sources w t ( r ) = χ ( r ) u t tot (a product of the contrast and the total field u t tot ) are unknown and confined to the domain D .
The first term of Equation (2) represents the “data error”, given by
F S ( w t ) = t u t sct L S { w t } S 2 t u t sct S 2 ,
and correspondingly, the second term represents the “domain error”,
F D ( χ , w t ) = t χ u t inc w t + χ L D { w t } D 2 t χ u t inc D 2 ,
with · D and · S denoting the L 2 norms on D and S . The measured scattered field data is represented above by u t sct at R receiver locations per transmitter t, and the incident field inside D is u t inc . The forward operator L S converts contrast source estimates w t in D to scattered field values at R receiver points per transmitter t on the measurement surface S ; L D performs a similar function, but transforms contrast sources to scattered field values within D . For further simplicity, the “data error” per iteration i can be defined from the numerator of Equation (3) as
ρ t ( i ) = d t L S { w t ( i ) }
and correspondingly, the “domain error” from the numerator of Equation (4) as
r t ( i ) = χ ( i ) u t inc w t ( i ) + χ ( i ) L D { w t ( i ) } .
These error values, particularly the data error of Equation (5), are subject to analysis for evaluating convergence behaviour of DGM-CSI in this study, described in Section 2.3.

2.2. Frequency-Cycling Tissue-Dependent Mapping Technique

The well-documented dielectric properties of breast tissues lead to the observation that the geometry of the real and imaginary parts of tissues’ permittivity profiles should be qualitatively similar [24,25]. The quantitative permittivity values of different types of breast tissue also form reasonably discrete ranges that are highly correlated between real and imaginary components, making it easy to stratify these values into expected tissue types (i.e., fat, transitional, fibroglandular, and cancer) [30]. Through a tissue-dependent mapping of the imaginary part based on geometric and quantitative reconstructions of the dielectric constant, an “artificial” imaginary component was proposed as part of the initial guess for the next frequency dataset in the imaging sequence, to good effect [24,25].
In addition to tissue-dependent mapping, this previous work invoked “frequency cycling” as opposed to conventional frequency hopping, where frequency cycling simply continues the reconstruction process past the highest frequency in a multi-frequency sequence by returning to the lowest frequency dataset (with an initial guess based on the highest frequency solution) [25]. The initial studies using DGM-CSI have shown that this practice preserves the fine detail acquired from high-frequency data, and produces better results than a single incremental frequency-hopping sequence operating for the same total number of iterations. In a practical microwave imaging scenario without the benefit of foreknowledge of this total number of iterations, however, a potential drawback to frequency cycling would be the increased number of iterations required to produce images when arbitrary or overly-stringent stopping conditions are applied.

2.3. Stopping Criteria for Single-Frequency Reconstructions

As described in Section 2.1, iterative inversion methods such as DGM-CSI attempt to solve the inverse scattering problem by minimizing the cost functional of Equation (2), which includes a data error ρ t ( i ) (Equation (5)) represented by the difference between the measured and computed fields on S for a given transmitter t. If U is the total number of data points across all receivers R and sources T, then the U × 1 data-error vector ρ ¯ ( i ) can be represented as a concatenation of every transmitter’s data errors, such that
ρ ¯ ( i ) = [ ρ 1 ( i ) ; ρ 2 ( i ) ; ; ρ T ( i ) ]
After all contrast sources w t ( i ) are updated for each transmitter t, the inversion process by convention performs a convergence check that calculates the new data-error vector ρ ¯ ( i + 1 ) . The goal is to reduce this new data error as much as possible, and aside from simply assigning a maximum number of iterations, one could monitor the data error until it falls below an arbitrarily-chosen tolerance level δ such that some form of a termination condition examining its sum, maximum value, or norm applies, commonly implemented as
if ( ρ ¯ ( i + 1 ) < δ ) stop .
Alternatively, the difference between the normalized data-error vector at two successive iterations, ρ ¯ ( i ) and ρ ¯ ( i + 1 ) , can be calculated until it falls within a prescribed tolerance to check for convergence [26].
Such conditions would seem to be reasonable criteria by which termination of the inversion algorithm should take place. However, the total data error often has little bearing on the accuracy of the resulting reconstruction, despite properly converging to a preset minimum. While liberally large values of error tolerance will logically terminate the imaging algorithm prematurely, it has also been observed that conservatively small choices of error tolerance will result in a poorer quality image. This outcome is likely due to the supposition that the minimization of a cost functional towards zero arguably means accounting for any and all error in modeling, calibration, and measurements (i.e., noise) through non-physical changes in the reconstructions. The apparent inadequacy of having the data error converge to within an arbitrary tolerance alone served as motivation to develop a multifaceted approach to stopping criteria that improves the reconstructions of a frequency-cycling reconstruction technique through statistical analysis of data and domain error of each iteration.
This proposed framework stops the imaging algorithm based primarily on the statistical distribution of the contributing terms to the data-error vector ρ ¯ . Referring the reader back to Equation (7), assuming a best-case scenario where the forward operator L S is able to perfectly model the field behaviour of the imaging system and the contrast sources w t had been iteratively solved to the exact solution, the data error ρ t would be reduced to the difference between the field values d t experimentally measured at the receivers for all given transmitters and the ideal noise-free field values generated by L S { w t ( i ) } . This difference would amount to the noise of the experimental system, which would presumably converge to a known pattern, such as a Gaussian distribution.
However, it is clear that even if one could acquire completely noise-free measurement data, it is not practical, or even possible, to find a numerical model that could exactly match it. This so-called modeling error, attributed to mismatches between the approximate representation of the imaging environment L S (which varies by formulation, mesh granularity, order number, etc.) and the true, unknowable functions that perfectly predict field behaviours within the system, invalidates the assumption that the remaining data-error vector ρ ¯ (upon convergence of the imaging algorithm to a solution) consists simply of noise. Modeling error is the primary reason why data calibration must be performed on experimental data, and it is not a trivial issue to address in MWI [31].
Correspondingly, since the contributions to the data error are multifactorial in nature and vary with each experimental set-up and possibly with each target, it will not necessarily converge to a known family of statistical distributions. Use of any of the myriad of statistical normality tests available (such as the Anderson–Darling or Shapiro–Wilk test) on the data error as part of the proposed algorithmic stopping criteria would therefore be unreliable. However, the profile of the data error distribution (regardless of its final form) should reach a steady state when the algorithm has sufficiently converged to a particular solution. To confirm that this trend occurs during a reconstruction requires comparison of the current iteration’s data error distribution to those of multiple previous iterations, which is accomplished in this work with an implementation of the two-sample Kolmogorov-Smirnov (K-S) test.
The two-sample K-S test is a non-parametric goodness-of-fit hypothesis test that evaluates the maximum absolute difference D between the cumulative distribution functions (CDF) of two sample data vectors over a range of x in each data set. Recall that the CDF ( F X ) of a continuous random variable X can be expressed as the integral of its probability density function (PDF) f x (that is, a function whose value at any given sample can be interpreted as providing a relative likelihood that the value of the random variable would be equal to that sample), such that
F X ( x ) = x f X ( t ) d t .
In the discrete case, F X ( x ) would be equivalent to the proportion of values in f X ( x ) less than or equal to a specified value of x.
Suppose the data-error vector PDFs from two consecutive iterations of the DGM-CSI imaging algorithm, ρ ¯ ( i ) and ρ ¯ ( i + 1 ) , each have some observed CDF, say P ^ ( i ) ( x ) and P ^ ( i + 1 ) ( x ) ; that is, for any specified value of x, the value of P ^ ( x ) is the proportion of error values less than or equal to x in the corresponding dataset. Consider the quantity
D = max x ( | P ^ ( i + 1 ) ( x ) P ^ ( i ) ( x ) | ) ,
that represents the maximum error between two CDFs. The K-S test makes use of this D value to return a probability, in the form of a scalar asymptotic p-value in the range of [ 0 , 1 ] , which represents the likelihood that the two samples were drawn from the same distribution. More precisely, the p-value is defined as the probability of observing a test statistic (D) as extreme or more extreme than the observed value, for the null hypothesis that the data in these two vectors are from the same continuous distribution [32,33,34,35]. In the implementation of the test, the decision to reject the null hypothesis is based on comparing the returned p-value with a preset significance level (commonly 5% for many applications). The test becomes very accurate for large sample sizes, and is believed to be reasonably accurate for sample sizes U 1 and U 2 , such that ( U 1 × U 2 ) / ( U 1 + U 2 ) 4 , which is the case for a collection of data error values in a system of 24 transmitters and receivers operating at a single frequency (representing a total of 576 data points).
While the strength of the two-sample K-S test lies in determining when two samples are from differing distributions, the nature of the test guarantees that the smaller the value of D, the higher the returned p-value will be. This feature allows the test’s p-value to be used as a numerical indicator for how likely the two samples were drawn from the same distribution (i.e., how closely their CDFs match), as opposed to its primary use in rejecting the null hypothesis.
Empirical studies of the real and imaginary parts of the data error of DGM-CSI reconstructions indeed revealed that the shape of the probability distribution curves fluctuate significantly with large variances in early iterations (Figure 2) but gradually reach a steady state in later iterations as expected, once all remarkable subjective change in the resulting images had ceased (Figure 3). Since any two neighboring iterations’ data error may only change slightly and be interpreted as being drawn from the same distribution by the two-sample K-S test, it was decided that a sliding window of past iterations’ data errors should be examined and compared to that of the current iteration. Moreover, the threshold value that the returned p-value should surpass for any given K-S test within this sliding window needed to be high enough to similarly reject subtler variations in two iterations of data error that would otherwise pass the goodness-of-fit test. This wider breadth of analysis ensured that more gradual changes were detected so that any corresponding stopping criteria would not halt the reconstruction prematurely.
Therefore, the window size of iterations under scrutiny, the p-value threshold (assumed constant over the window), and the percentage of that window’s iterations reaching or surpassing that p-value threshold, are all parameters to the stopping criteria for the algorithm based on the K-S test, which is applied separately for the real and imaginary parts of each iteration’s data error. A preliminary study attempting to establish the effects of these parameters was undertaken for the frequency-cycled reconstruction of a 2D breast model, described in Section 3.1.
It should be noted that the statistical analysis governing the aforementioned stopping criteria only makes use of the data error, which may be sufficient for certain MWI algorithms (such as the distorted Born iterative method) whose minimized cost functional consists solely of a data error term. However, as shown in Section 2.1, CSI-based iterative inversion algorithms make use of cost functionals containing a second domain error term dealing with total field values inside the imaging domain [11,26]. Therefore, another ad hoc stopping condition based on the domain error is also prudently employed in the algorithm, as described in the next subsection.

2.4. Global Termination of Multi-Frequency Reconstructions

The stopping conditions described in Section 2.3 are appropriate for halting single-frequency inversions, but multiple-frequency reconstructions have thus far not been addressed. In a frequency-hopping scenario, the stopping criteria based on K-S tests alone can be easily employed to determine when each individual inversion should terminate and move on to the next frequency, with global termination occurring following reconstruction of the highest frequency.
However, in a frequency-cycling reconstruction scheme, the cycle of reconstructions could theoretically continue ad infinitum without another explicit set of rules in place governing the global termination of the imaging process. A separate, secondary ad hoc global termination criterion was therefore implemented, based on a previous study [17]; once each frequency’s dataset in the cycle is subjected to the K-S test stopping criteria at least once, if the percentage of relative change in domain error between two successive iterations falls below 0.1%, then it is deemed appropriate to globally terminate the frequency cycle and end the reconstruction. This provision could also serve as a fallback contingency in the event that the data error-based conditions turn out to be too strict to satisfy for a given individual frequency.

2.5. Full Description of Multi-Frequency Imaging Procedure

The key concepts involved in the proposed multi-frequency imaging procedure, while implemented herein for the DGM-CSI algorithm, can easily be applied (with minor modifications) to other iterative algorithms widely used in MWI, including Gauss-Newton Inversion (GNI) and distorted Born Iterative Method (DBIM). As such, a descriptive approach to the technique is undertaken here in order to keep its applicability as general as possible. Taking into account the stopping criteria described throughout Section 2.3 and using the procedural description of the frequency-cycling reconstruction approach with tissue-dependent mapping described in [25], the new consolidated multi-frequency imaging procedure is proposed as follows:
  • The real and imaginary parts of the complex permittivity are reconstructed using the lowest frequency data available (e.g., 1.0 GHz). The termination point of this reconstruction is dictated by the results of successive two-sample K-S tests performed on both the real and imaginary parts of the data error separately, comparing the current iteration to those of a sliding window of past iterations, governed by a choice of parameters for p-value, window size, and percentage of windowed iterations reaching this p-value threshold. For robustness, a back-up termination condition may be implemented, either related to the relative change in domain error for CSI-based algorithms as described in Section 2.4, or a maximum number of iterations.
  • A point-by-point search through the reconstructed real part of each nodal basis coefficient in the DGM-CSI mesh (or more generally, each mesh element or pixel of the reconstructed image) classifies the type of breast tissue. This classification is based solely on the range of expected values of dielectric constant at that frequency, as outlined in [25].
  • An initial guess for the next imaging frequency (e.g., 2.0 GHz) is generated using the tissue-dependent mapping process [24,25]. It consists of the unmodified real parts of the reconstructed ε r at the mesh nodal points, and a new imaginary part created from a simple linear interpolation of the expected range of dielectric loss values, based on the appropriate Cole-Cole models of tissues classified in Step 2. This technique preserves the geometry of the real and imaginary parts of the solution.
  • This new initial guess for the complex permittivity is used to run the inversion algorithm at the next frequency (e.g., 2.0 GHz). As per the procedure outlined in [25], the user may choose to keep the imaginary part constant during this inversion and update only the real part to converge to a new solution. This “anchoring” process has been shown to improve overall imaging results due to the tendency of CSI-based inversion algorithms to cause significant deterioration of the imaginary part at high-frequency reconstructions. Again, the aforementioned parameterized stopping criteria would be primarily employed to determine the appropriate point to halt this reconstruction.
  • If more than two frequencies are used in the frequency hop, steps 2–4 are repeated as necessary until the reconstruction of the final frequency of the succession is complete (e.g., 3.0 GHz). This succession may include “frequency cycling”; that is, returning the inversion algorithm to the lowest frequency data and incrementally stepping through each frequency again [25].
If a simple frequency-hopping scheme is employed, the reconstruction will terminate after the highest-frequency inversion, again as decided by the parameterized stopping criteria based on K-S tests of the data error. However, if a frequency-cycling scheme is used, and if the imaginary part of the solution has been “anchored” in place following the first inversion in Step 1, there are further considerations that come into play for global termination of the imaging process, and to address concerns of full CSI optimization:
  • When each available dataset in the frequency cycle has been used at least once to contribute to the overall image reconstruction, a global termination criterion will become active, which will monitor the relative change in the domain error between successive iterations (Section 2.4). If this relative change falls below 0.1% at any point, the current reconstruction is halted and the frequency cycle is broken.
  • Regardless of the frequency at which the algorithm was halted by this relative domain error threshold, if the imaginary part of the solution has been continuously held constant during the frequency cycle after Step 1, one last initial guess is generated as in Step 3 and a final reconstruction is run at the lowest frequency available (e.g., 1.0 GHz) with both the real and imaginary parts allowed to converge to a solution (i.e., the imaginary part is no longer “anchored”). This final inversion is terminated by the parameterized stopping criteria or a relative change of domain error between successive iterations falling below 0.1%, whichever occurs first. The purpose of this final run is to demonstrate the stability of the final solution and ensure that its imaginary part, despite being originally based on the geometry and tissue properties of the real part, does indeed satisfy full CSI optimization.

2.6. Synthetic Breast Models

Initial testing of the described stopping criteria and comparison to fixed-iteration scenarios were carried out on synthetic transverse magnetic (TM) data collected from a two-dimenstional MRI-derived BI-RADS Category C (heterogeneously dense) 2D breast model supplied by the University of Calgary, derived from an MRI slice of a cancer patient with a breast tumour visible at the “3 o’clock” position. To demonstrate the robustness of the technique across different breast phantoms, an additional two-dimensional MRI-derived BI-RADS Category D (“extremely dense”) healthy breast model supplied by the University of Wisconsin’s public database (ID: 070604PA2) was employed. Both models are depicted at three different frequencies (1.0 GHz, 2.0 GHz, and 3.0 GHz) in Figure 4 and Figure 5, illustrating the similarities between breast tissues’ dielectric permittivity at the low-gigahertz band. The values of the tissue-dependent complex permittivity in all cases were calculated for every frequency using an appropriate fitted single-pole Cole-Cole model (or equivalently, a single-pole Debye model as the exponent parameter α = 0 ), and for the University of Calgary model, then subjected to random perturbations of ±10% [36].
Data was collected for both models in a low-loss background of ϵ r = 23 1.13 i with 24 transmitters and 24 receivers evenly distributed at a radius of 10 cm, employing a finite-element method (FEM) forward solver with a finely-discretized mesh independent from that used for the DGM-CSI inversion. The synthetic data was collected for open boundary conditions at 1.0 GHz, 2.0 GHz, and 3.0 GHz. To test the robustness of the stopping criteria to problems with PEC boundaries, and to better reflect the capabilities of a recent University of Manitoba imaging system prototype [37], data was collected a second time with the Category C model surrounded by a circular PEC boundary with a radius of 15 cm at a reduced bandwidth. The frequencies used for the PEC-bounded case were 1.0 GHz, 1.25 GHz, and 1.5 GHz. Uniformly-distributed noise at 5% of the maximum field magnitude was used to corrupt the E z scattered electric field data for the Category C model, and several noise levels (3%, 5%, 7.5%, and 10%) were tested for the Category D model to evaluate the effect of noise on the stopping criteria’s performance. The skin thickness and inner skin boundary both remained unknown; the only prior information used during inversions was the outer skin boundary.

2.7. Error Calculation

An objective method of evaluating the DGM-CSI algorithm’s performance for differing imaging scenarios was accomplished by comparing the reconstruction results to the original Category C and Category D synthetic models of the current frequency being used in the reconstruction cycle, unless otherwise specified. The models and reconstructions each employed different underlying meshes, so their complex permittivity values were interpolated onto a common 250 × 250-pixel square grid (N = 62,500) that encompassed an area covering ( x 1 , y 1 ) = ( 0.07 , 0.07 ) [m] to ( x 2 , y 2 ) = ( 0.07 , 0.07 ) [m]. The relative error of L 1 and L 2 norms of the difference between the model (actual) complex permittivity ε ^ r act and the reconstructed complex permittivity ε ^ r rec over each pixel k were used as primitives to evaluate the performance of the imaging algorithm. For convenience, these will hereafter be referred to as the relative error norms REN 1 and REN 2 , and are calculated and expressed as percentages, where
REN 1 ( ε ^ r act , ε ^ r rec ) = ε ^ r act ε ^ r rec 1 ε ^ r act 1 × 100 % = k = 1 N | ε ^ r act ( k ) ε ^ r rec ( k ) | k = 1 N | ε ^ r act ( k ) | × 100 % ,
REN 2 ( ε ^ r act , ε ^ r rec ) = ε ^ r act ε ^ r rec 2 ε ^ r act 2 × 100 % = k = 1 N [ ε ^ r act ( k ) ε ^ r rec ( k ) ] 2 k = 1 N [ ε ^ r act ( k ) ] 2 × 100 % .
The raw numerical value of the norms will be slightly reduced since the square grid covers an area outside of the imaging domain of the problem, which will consist of common unaltered background permittivity in both the model and reconstruction and thus represent pixels of zero error. Since the REN values serve primarily as a means of monitoring iterative trends in solution convergence and to demonstrate comparative improvement between imaging scenarios with equally-sized interpolated grids and imaging domains, the impact of this artificial error reduction is of little significance.

3. Results and Discussion

3.1. Imaging with Open Boundaries

As mentioned in the introduction of the stopping criteria in Section 2.3, the approach had three variables: the p-value threshold that any given two-sample K-S test would use to conclude the two data error samples were drawn from the same distribution, the window size of past iterations’ data errors compared to that of the current iteration, and the percentage of that window’s iterations that would need to to reach or surpass the chosen p-value threshold in order to halt the reconstruction at the current frequency, deemed the “pass percentage”. For the same Category C breast model in Figure 4, DGM-CSI was run in a frequency-cycled configuration using the rules described in Section 2.5 for every combination of p-value thresholds of 0.90, 0.95 and 0.99, window sizes of 10, 30, and 50, and pass percentages of 80%, 90%, and 100%. There were 27 simulations whose total number of iterations would vary based on the laxity or stringency of stopping criteria associated with low and high variable values respectively. Error norms ( REN 1 and REN 2 ) for each scenario were plotted for every frequency jump (Figure 6), and the best variable combination (resulting in the lowest error metric magnitudes upon termination) was found to be a p-value threshold of 0.99, a window size of 30, and a pass percentage of 80%. This combination of values was subsequently used for all future simulations employing the stopping criteria. Note that in Figure 6 and in all subsequent plots of REN values, the x-axis’ number of iterations represents the cumulative total throughout all frequency-hopping or frequency-cycling steps of the reconstruction.
Having established an appropriate set of variables to govern the stopping criteria, a series of simulations adhering to the full description of the technique in Section 2.5 could be carried out to demonstrate its benefit in efficiency over conventional methods. The REN values of three scenarios (A, B, and C from Table 1) comparing fixed-iteration reconstructions (with and without tissue-dependent mapping) to the use of both tissue-dependent mapping and the stopping criteria, are shown in Figure 7. The final imaging reconstructions of these three cases are also depicted in Figure 8. As the proposed stopping conditions introduce an unpredictable halt to the reconstruction cycle that could occur at any frequency and also includes a mandatory return to the lowest frequency (in this case, 1.0 GHz) for the final optimization, each REN data point may not represent consecutive frequencies of the cycle. For consistency, the REN values shown in Figure 7 are therefore all calculated based on the same 2.0 GHz model from Figure 4. For clarity, the frequency associated with each data point plotted in Figure 7, along with the triggering condition that caused a change in frequency or termination in the reconstruction (when relevant), is shown in Table 1.
Although DGM-CSI, like other iterative microwave imaging algorithms, performs reasonably well in recovering the real component of breast models it reconstructs, it has difficulty producing an accurate profile of the imaginary part, most notably at higher frequencies. This deficiency is illustrated from the Scenario A inversion of the Category C breast model in the top image of Figure 8. Subjectively, the benefit of employing the tissue-dependent mapping in the fidelity of the imaginary component is obvious in both the middle and bottom images. It has been already shown elsewhere that the described mapping technique and “anchoring” of the imaginary component at each subsequent frequency stage, along with the practice of frequency cycling, imparts a benefit to the quality of the real part of the reconstructions following the same number of iterations [24,25]. However, these other studies only tested reconstructions that used fixed numbers of iterations as termination conditions; the addition of the stopping criteria for the latter case here (Scenario C) does not appreciably degrade the quality of the final image, while reducing the total number of iterations by over 35%.

3.2. Imaging with PEC Boundaries

A set of simulations similar to those carried out in the previous section were performed using PEC boundaries and data from frequencies of 1.0 GHz, 1.25 GHz and 1.50 GHz, to demonstrate this method’s robustness to smaller frequency bandwidths and different boundary conditions, as mentioned in Section 2.6. Scenarios D and E are imaging simulations recreated from an earlier study that employed fixed-iteration reconstructions exclusively [25]; they are included in order to emphasize the benefit of adding stopping criteria in Scenario F. The imaging results are shown in Figure 9 with a breakdown of the reconstruction progression documented in Table 2. A plot of the relative error norms again similar to the open boundary cases, employing REN calculations based on the original 1.25 GHz model, is also provided in Figure 10.
As explained in [25], it was expected that the results would not be as impressive as open-boundary scenarios since PEC-bounded problems are generally more difficult to solve, and data was taken at lower frequencies. However, it is again clear from Figure 9 that there is an improvement in the recovery of the imaginary component when tissue-dependent mapping is used (middle and bottom images), more so for Scenario F. A modest improvement in the real part of ε r for Scenarios E and F is also noted. Without the use of stopping criteria, however, Scenario E overshoots the expected values of dielectric constant in several areas. Its REN values depicted in Figure 10 correspondingly suffer increases relating to the fact that the reconstruction cycle has perhaps been allowed to run too long, as pointed out in [25]. The introduction of the stopping criteria has effectively prevented this degradation by terminating the reconstruction cycle close to the observed nadir of these error values from Scenario E, while once again demonstrating a sizable increase in computational efficiency by reducing the number of iterations by an impressive 66%.

3.3. Effect of Noise Levels

To validate that the stopping criteria developed and tested on the BI-RADS Category C 2D model would perform similarly on different breast models and be robust to varying levels of noise, another full set of simulations were carried out on the aforementioned Category D model from the University of Wisconsin database using open boundary data at 1.0 GHz, 2.0 GHz, and 3.0 GHz (Figure 5). The uniformly-distributed noise used to corrupt the E z scattered electric field data was generated at 3%, 5%, 7.5%, and 10% of the maximum field magnitude. The imaging results are shown in Figure 11 with a breakdown of the reconstruction progression documented in Table 3. A plot of the relative error norms again similar to the Category C cases, employing REN calculations based on the original 2.0 GHz Category D model, is also provided in Figure 12.
As shown in the final images of Figure 11, there is an obvious improvement using tissue-dependent mapping and stopping criteria over the simple frequency-hopping approach (Scenario G), especially in the imaginary part, whose geometry is not reconstructed properly and contains a large overshooting boundary artifact. However, remarkably, there is very little subjective difference in the images reconstructed from data corrupted with different levels of noise (from 3% to 10% for Scenarios H through K), with subtle changes in the real and imaginary parts visible only with careful scrutiny.
The relative error norms of all scenarios depicted in Figure 12 and the reconstruction break-down in Table 3 offer more objective measurements of the effect of noise. Specifically, lower final REN values are observed for images inverted from data containing lower noise levels, as would be expected. Also, the stopping criteria transitions frequencies and halts the entire reconstruction cycle sooner for higher noise levels. These trends are consistent for all reconstructions employing tissue-dependent mapping and stopping criteria in each level of noise used in the simulations. For instance, data containing 3% noise (Scenario H) runs longest, for a total of 759 iterations; this total is observed to gradually decrease at the 5% and 7.5% noise levels, and data containing 10% noise (Scenario K) runs for only 448 iterations.
These observations are congruous with the theory introduced in Section 2.3 that data containing higher magnitudes of noise would have less useful information for the optimization algorithm to extract before the distribution of data error reached a steady state (i.e., before the algorithm began to reconstruct noise) and continuing the inversion would be superfluous, or perhaps even detrimental to the final image. Overall, while these additional 2D simulations do not necessarily represent a fully comprehensive investigation across all possible breast models and imaging scenarios, the performance of the proposed stopping criteria appears satisfactory on the Category D model, and behaves as expected on noisier datasets, even with parameter choices derived from the Category C model reconstructions.

4. Conclusions

It has been demonstrated that a stopping condition based on the similarities of the statistical distribution of data error at successive iterations has the capability of significantly speeding up the DGM-CSI algorithm while having no appreciable detrimental effect on the final imaging results, even yielding superior solutions than fixed-iteration scenarios. This boost of efficiency is granted through novel use of sequential two-sample K-S tests on the data error to determine appropriate times to halt the imaging algorithm, truncating likely unnecessary (and in some cases counterproductive) reconstruction iterations.
Whether or not the best set of parameters used in this analysis represents the true optimal values for the stopping criteria for other breast models, or if these same values are applicable to 3D inversions, is open to debate. However, given the promising imaging results obtained from this initial set of tests, and the gains in efficiency provided through significant reduction of the number of iterations needed to converge to these solutions, it was judged appropriate to keep the best parameters from this introductory exploration of the method and leave a more comprehensive examination of parameter optimization for future investigations.
The most interesting form of the proposed multi-frequency imaging procedure employs tissue-dependent mapping and “anchoring” of the imaginary part along with the new stopping criteria. This ad hoc procedure also returns the reconstruction cycle to stable low-frequency data after recovering anatomical detail from higher frequencies, allowing both real and imaginary parts to properly converge before global termination. The advantages gained by this use of both tissue-dependent mapping and novel stopping criteria have also been demonstrated to be robust for different breast models, frequency bandwidths, boundary conditions, and noise levels.
In fact, the stopping criteria have a more dramatic effect on the imaging results obtained from the more complex PEC-bounded problem and data with higher noise levels, suggesting this technique may be especially useful for avoiding image degradation and reconstruction artifacts in difficult cases such as full 3D reconstructions in semi-resonant enclosures, including experimental measurements that may be excessively noisy or contain poor frequency data unsuitable for inversion. To explore this possibility, preliminary extension of this technique to the 3D DGM-CSI algorithm developed at the University of Manitoba is a logical next step, along with parameter optimization studies of the stopping criteria for data from experimental systems.

Author Contributions

Conceptualization, C.K., I.J. and J.L.; Methodology, C.K., I.J. and J.L.; Resources, C.K., I.J. and J.L.; Software, C.K. and I.J.; Supervision, J.L.; Validation, C.K. and I.J.; Writing—original draft, C.K.; Writing—review & editing, C.K., I.J. and J.L.

Funding

This research is partially supported by the Natural Sciences and Engineering Research Council of Canada (NSERC).

Acknowledgments

The authors would like to thank the University of Calgary for providing the 2D Category C synthetic breast model used for the reconstructions presented in this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Trentham-Dietz, A.; Kerlikowske, K.; Stout, N.K.; Miglioretti, D.L.; Schechter, C.B.; Ergun, M.A.; van den Broek, J.J.; Alagoz, O.; Sprague, B.L.; van Ravesteyn, N.T.; et al. Tailoring breast cancer screening intervals by breast density and risk for women aged 50 years or older: Collaborative modeling of screening outcomes. Ann. Intern. Med. 2016, 165, 700–712. [Google Scholar] [CrossRef] [PubMed]
  2. D’Orsi, C.J.; Mendelson, E.B.; Morris, E.A. (Eds.) ACR BI-RADS Atlas: Breast Imaging Reporting and Data System, 5th ed.; American College of Radiology: Reston, VA, USA, 2012. [Google Scholar]
  3. Boerner, W.M.; Brand, H.; Cram, L.A.; Gjessing, D.T.; Jordan, A.K.; Keydel, W.; Schwierz, G.; Vogel, M. (Eds.) Inverse Methods in Electromagnetic Imaging; D. Reidel Publishing: Dordrecht, The Netherlands, 1985. [Google Scholar]
  4. Larsen, L.E.; Jacobi, J.H. (Eds.) Medical Applications of Microwave Imaging; IEEE Press: New York, NY, USA, 1986. [Google Scholar]
  5. Kosmas, P.; Crocco, L. Introduction to Special Issue on “Electromagnetic Technologies for Medical Diagnostics: Fundamental Issues, Clinical Applications and Perspectives”. Diagnostics 2019, 9, 19. [Google Scholar] [CrossRef]
  6. Sill, J.M.; Fear, E.C. Tissue sensing adaptive radar for breast cancer detection—Experimental investigation of simple tumor models. IEEE Trans. Microw. Theory Tech. 2005, 53, 3312–3319. [Google Scholar] [CrossRef]
  7. Delbary, F.; Brignone, M.; Bozza, G.; Aramini, R.; Piana, M. A visualization method for breast cancer detection using microwaves. SIAM J. Appl. Math. 2010, 70, 2509–2533. [Google Scholar] [CrossRef]
  8. Cakoni, F.; Colton, D.; Monk, P. Qualitative Methods in Inverse Electromagnetic Scattering Theory: Inverse Scattering for Anisotropic Media. IEEE Antennas Propag. Mag. 2017, 59, 24–33. [Google Scholar] [CrossRef]
  9. Colton, D.; Kress, R. Inverse Acoustic and Electromagnetic Scattering Theory; Springer Science & Business Media: Berlin, Germany, 2012; Volume 93. [Google Scholar]
  10. Engl, H.W.; Hanke, M.; Neubauer, A. Regularization of Inverse Problems; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1996. [Google Scholar]
  11. Van den Berg, P.M.; Kleinman, R.E. A contrast source inversion method. Inverse Probl. 1997, 13, 1607–1620. [Google Scholar] [CrossRef]
  12. Abubakar, A.; Van den Berg, P.M.; Mallorqui, J.J. Imaging of biomedical data using a multiplicative regularized contrast source inversion method. IEEE Trans. Microw. Theory Tech. 2002, 50, 1761–1771. [Google Scholar] [CrossRef] [Green Version]
  13. Zakaria, A.; Gilmore, C.; LoVetri, J. Finite-element contrast source inversion method for microwave imaging. Inverse Probl. 2010, 26, 115010. [Google Scholar] [CrossRef]
  14. Jeffrey, I.; Geddert, N.; Brown, K.; LoVetri, J. The Time-Harmonic Discontinuous Galerkin Method as a Robust Forward Solver for Microwave Imaging Applications. Prog. Electromagn. Res. 2015, 154, 1–21. [Google Scholar] [CrossRef]
  15. Brown, K.G.; Geddert, N.; Asefi, M.; LoVetri, J.; Jeffrey, I. Hybridizable Discontinuous Galerkin Method Contrast Source Inversion of 2-D and 3-D Dielectric and Magnetic Targets. IEEE Trans. Microw. Theory Tech. 2019, 67, 1766–1777. [Google Scholar] [CrossRef]
  16. Baran, A.; Kurrant, D.; Zakaria, A.; Fear, E.; LoVetri, J. Breast Imaging Using Microwave Tomography with Radar-Based Tissue-Regions Estimation. Prog. Electromagn. Res. 2014, 149, 161–171. [Google Scholar] [CrossRef]
  17. Kurrant, D.J.; Baran, A.; Fear, E.C.; LoVetri, J. Integrating prior information into microwave tomography Part 1: Impact of detail on image quality. Med. Phys. 2017, 44, 6461–6481. [Google Scholar] [CrossRef]
  18. Kurrant, D.; Baran, A.; LoVetri, J.; Fear, E. Integrating prior information into microwave tomography part 2: Impact of errors in prior information on microwave tomography image quality. Med. Phys. 2017, 44, 6482–6503. [Google Scholar] [CrossRef]
  19. Abdollahi, N.; Kurrant, D.; Mojabi, P.; Omer, M.; Fear, E.; LoVetri, J. Incorporation of Ultrasonic Prior Information for Improving Quantitative Microwave Imaging of Breast. IEEE J. Multiscale Multiphys. Comput. Tech. 2019, 4, 98–110. [Google Scholar] [CrossRef]
  20. Nemez, K.; Baran, A.; Asefi, M.; LoVetri, J. Modeling Error and Calibration Techniques for a Faceted Metallic Chamber for Magnetic Field Microwave Imaging. IEEE Trans. Microw. Theory Tech. 2007, 65, 4347–4356. [Google Scholar] [CrossRef]
  21. Asefi, M.; Baran, A.; LoVetri, J. An Experimental Phantom Study for Air-Based Quasi-Resonant Microwave Breast Imaging. IEEE Trans. Microw. Theory Tech. 2019. [Google Scholar] [CrossRef]
  22. 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]
  23. Chew, W.; Lin, J. A frequency-hopping approach for microwave imaging of large inhomogeneous bodies. IEEE Microw. Guided Wave Lett. 1995, 5, 439–441. [Google Scholar] [CrossRef]
  24. Kaye, C.; Jeffrey, I.; LoVetri, J. Enhancement of multi-frequency microwave breast images using a tissue-dependent mapping technique with discontinuous Galerkin contrast source inversion. In Proceedings of the AES 2016—4th Advanced Electromagnetics Symposium, Malaga, Spain, 26–28 July 2016; pp. 37–39. [Google Scholar]
  25. Kaye, C.; Jeffrey, I.; LoVetri, J. Improvement of Multi-Frequency Microwave Breast Imaging through Frequency Cycling and Tissue-Dependent Mapping. IEEE Trans. Antennas Propag. 2019. submitted. [Google Scholar]
  26. Abubakar, A.; van den Berg, P.M.; Habashy, T.M. Application of the multiplicative regularized contrast source inversion method on TM- and TE-polarized experimental Fresnel data. Inverse Probl. 2005, 21, S5–S13. [Google Scholar] [CrossRef] [Green Version]
  27. Van den Berg, P.M.; van Broekhoven, A.L.; Abubakar, A. Extended contrast source inversion. Inverse Probl. 1999, 15, 1325–1344. [Google Scholar] [CrossRef]
  28. Bloemenkamp, R.F.; Abubakar, A.; van den Berg, P.M. Inversion of experimental multi-frequency data using the contrast source inversion method. Inverse Probl. 2001, 17, 1611–1622. [Google Scholar] [CrossRef]
  29. Jeffrey, I.; Zakaria, A.; LoVetri, J. Microwave Imaging by Mixed-Order Discontinuous Galerkin Contrast Source Inversion. In Proceedings of the 2014 XXXIst URSI General Assembly and Scientific Symposium (URSI GASS 2014), Beijing, China, 16–23 August 2014; pp. 255–258. [Google Scholar]
  30. Lazebnik, M.; Popovic, D.; McCartney, L.; Watkins, C.B.; Lindstrom, M.J.; Harter, J.; Sewall, S.; Ogilvie, T.; Magliocc, M.A.; Breslin, T.M.; et al. A large-scale study of the ultrawideband microwave dielectric properties of normal, benign and malignant breast tissues obtained from cancer surgeries. Phys. Med. Biol. 2007, 52, 6093–6115. [Google Scholar] [CrossRef] [PubMed]
  31. Kaye, C. Development and Calibration of Microwave Tomography Imaging Systems for Biomedical Applications Using Computational Electromagnetics. Master’s Thesis, University of Manitoba, Winnipeg, MB, Canada, 2009. [Google Scholar]
  32. Massey, F.J. The Kolmogorov-Smirnov Test for Goodness of Fit. J. Am. Stat. Assoc. 1951, 46, 68–78. [Google Scholar] [CrossRef]
  33. Miller, L.H. Table of Percentage Points of Kolmogorov Statistics. J. Am. Stat. Assoc. 1956, 51, 111–121. [Google Scholar] [CrossRef]
  34. Marsaglia, G.; Tsang, W.; Wang, J. Evaluating Kolmogorov’s Distribution. J. Stat. Softw. 2003, 8, 1–4. [Google Scholar] [CrossRef]
  35. MathWorks. Two-Sample Kolmogorov-Smirnov Test (R2018a). Available online: www.mathworks.com/help/stats/kstest2.html (accessed on 16 March 2018).
  36. Kurrant, D.J.; Fear, E.C. Regional estimation of the dielectric properties of inhomogeneous objects using near-field reflection data. Inverse Probl. 2012, 28, 075001. [Google Scholar] [CrossRef]
  37. Nemez, K.; Asefi, M.; Baran, A.; LoVetri, J. A faceted magnetic field probe resonant chamber for 3D breast MWI: A synthetic study. In Proceedings of the 2016 17th International Symposium on Antenna Technology and Applied Electromagnetics (ANTEM), Montreal, QC, Canada, 10–13 July 2016; pp. 322–324. [Google Scholar]
Figure 1. Two-dimensional representation of the imaging problem. The object of interest (OI) has an unknown contrast χ , where D is the imaging domain and S is the surface containing the transmitters (Tx) and receivers (Rx).
Figure 1. Two-dimensional representation of the imaging problem. The object of interest (OI) has an unknown contrast χ , where D is the imaging domain and S is the surface containing the transmitters (Tx) and receivers (Rx).
Jimaging 05 00055 g001
Figure 2. Examples of histograms demonstrating estimates of the probability distribution of the real and imaginary parts of the data error ( ρ E ) during a DGM-CSI inversion of an arbitrary 2D synthetic breast model early in the reconstruction process (iteration 20).
Figure 2. Examples of histograms demonstrating estimates of the probability distribution of the real and imaginary parts of the data error ( ρ E ) during a DGM-CSI inversion of an arbitrary 2D synthetic breast model early in the reconstruction process (iteration 20).
Jimaging 05 00055 g002
Figure 3. Examples of histograms demonstrating estimates of the probability distribution of the real and imaginary parts of the data error ( ρ E ) during the same DGM-CSI inversion as Figure 2, later in the reconstruction process (iteration 190).
Figure 3. Examples of histograms demonstrating estimates of the probability distribution of the real and imaginary parts of the data error ( ρ E ) during the same DGM-CSI inversion as Figure 2, later in the reconstruction process (iteration 190).
Jimaging 05 00055 g003
Figure 4. Complex dielectric properties of University of Calgary 2D synthetic breast model at 1.0 GHz (top), 2.0 GHz (middle) and 3.0 GHz (bottom).
Figure 4. Complex dielectric properties of University of Calgary 2D synthetic breast model at 1.0 GHz (top), 2.0 GHz (middle) and 3.0 GHz (bottom).
Jimaging 05 00055 g004
Figure 5. Complex dielectric properties of University of Wisconsin 2D synthetic breast model at 1.0 GHz (top), 2.0 GHz (middle) and 3.0 GHz (bottom).
Figure 5. Complex dielectric properties of University of Wisconsin 2D synthetic breast model at 1.0 GHz (top), 2.0 GHz (middle) and 3.0 GHz (bottom).
Jimaging 05 00055 g005
Figure 6. Relative error norms of open-boundary 2D DGM-CSI reconstructions across several choices of stopping criteria parameter values. Each curve in the figure corresponds to a frequency-cycled inversion for a particular choice of the three parameters (Section 3.1), with the window sizes coded by color for convenience. Data points on each curve correspond to REN values for the reconstructed model at each frequency change. Arrows point to the final REN values of simulations that terminated with the lowest relative error norms, indicating the best combination of parameters among those tested for this imaging scenario.
Figure 6. Relative error norms of open-boundary 2D DGM-CSI reconstructions across several choices of stopping criteria parameter values. Each curve in the figure corresponds to a frequency-cycled inversion for a particular choice of the three parameters (Section 3.1), with the window sizes coded by color for convenience. Data points on each curve correspond to REN values for the reconstructed model at each frequency change. Arrows point to the final REN values of simulations that terminated with the lowest relative error norms, indicating the best combination of parameters among those tested for this imaging scenario.
Jimaging 05 00055 g006
Figure 7. Relative error norms of frequency-hopping and frequency-cycled reconstructions: Scenario A (solid line)—without tissue-dependent mapping at 400 iterations per frequency terminating after first inversion of 3.0 GHz, Scenario B (dashed line)—with tissue-dependent mapping at 200 iterations per frequency, cycling through reconstruction frequencies once (with imaginary component “anchored” following initial 1.0 GHz inversion), Scenario C (dotted line)—with tissue-dependent mapping and stopping criteria in place, terminating after two consecutive inversions of 1.0 GHz data (one with the imaginary component “anchored” and the final run with the imaginary component freely optimized according to the guidelines of Section 2.5). See Table 1 for further details.
Figure 7. Relative error norms of frequency-hopping and frequency-cycled reconstructions: Scenario A (solid line)—without tissue-dependent mapping at 400 iterations per frequency terminating after first inversion of 3.0 GHz, Scenario B (dashed line)—with tissue-dependent mapping at 200 iterations per frequency, cycling through reconstruction frequencies once (with imaginary component “anchored” following initial 1.0 GHz inversion), Scenario C (dotted line)—with tissue-dependent mapping and stopping criteria in place, terminating after two consecutive inversions of 1.0 GHz data (one with the imaginary component “anchored” and the final run with the imaginary component freely optimized according to the guidelines of Section 2.5). See Table 1 for further details.
Jimaging 05 00055 g007
Figure 8. Final results of DGM-CSI frequency-hopping and frequency-cycled complex dielectric property reconstruction of synthetic breast model using 1.0–3.0 GHz data, without modification of intermediate initial guesses (Scenario A—top), using tissue-dependent mapping at fixed 200 iterations per frequency (Scenario B—middle), and employing stopping criteria and tissue-dependent mapping (Scenario C—bottom).
Figure 8. Final results of DGM-CSI frequency-hopping and frequency-cycled complex dielectric property reconstruction of synthetic breast model using 1.0–3.0 GHz data, without modification of intermediate initial guesses (Scenario A—top), using tissue-dependent mapping at fixed 200 iterations per frequency (Scenario B—middle), and employing stopping criteria and tissue-dependent mapping (Scenario C—bottom).
Jimaging 05 00055 g008
Figure 9. Final results of DGM-CSI frequency-hopping and frequency-cycled complex dielectric property reconstruction of PEC-bounded synthetic breast model using 1.0–1.5 GHz data, without modification of intermediate initial guesses (Scenario D—top), using tissue-dependent mapping at fixed 200 iterations per frequency (Scenario E—middle), and employing stopping criteria and tissue-dependent mapping (Scenario F—bottom).
Figure 9. Final results of DGM-CSI frequency-hopping and frequency-cycled complex dielectric property reconstruction of PEC-bounded synthetic breast model using 1.0–1.5 GHz data, without modification of intermediate initial guesses (Scenario D—top), using tissue-dependent mapping at fixed 200 iterations per frequency (Scenario E—middle), and employing stopping criteria and tissue-dependent mapping (Scenario F—bottom).
Jimaging 05 00055 g009
Figure 10. Relative error norms of reconstructions with PEC boundaries: Scenario D (solid line)—without tissue-dependent mapping at 400 iterations per frequency terminating after first inversion of 1.5 GHz data, Scenario E (dashed line)—with tissue-dependent mapping at 200 iterations per frequency, cycling through reconstruction frequencies once (with imaginary component “anchored” following initial 1.0 GHz inversion), Scenario F (dotted line)—with tissue-dependent mapping and stopping criteria in place, terminating after two consecutive inversions of 1.0 GHz data (one with the imaginary component “anchored” and the final run with the imaginary component freely optimized according to the guidelines of Section 2.5). See Table 2 for further details.
Figure 10. Relative error norms of reconstructions with PEC boundaries: Scenario D (solid line)—without tissue-dependent mapping at 400 iterations per frequency terminating after first inversion of 1.5 GHz data, Scenario E (dashed line)—with tissue-dependent mapping at 200 iterations per frequency, cycling through reconstruction frequencies once (with imaginary component “anchored” following initial 1.0 GHz inversion), Scenario F (dotted line)—with tissue-dependent mapping and stopping criteria in place, terminating after two consecutive inversions of 1.0 GHz data (one with the imaginary component “anchored” and the final run with the imaginary component freely optimized according to the guidelines of Section 2.5). See Table 2 for further details.
Jimaging 05 00055 g010
Figure 11. Final results of DGM-CSI frequency-hopping and frequency-cycled complex dielectric property reconstruction of open-boundary Category D synthetic breast model using 1.0–3.0 GHz data: (a) Scenario G—without tissue-dependent mapping at 400 iterations per frequency terminating after first inversion of 3.0 GHz (5% noise), (b) Scenario H—with tissue-dependent mapping and stopping criteria in place at 3% noise, (c) Scenario I—with tissue-dependent mapping and stopping criteria in place at 5% noise, (d) Scenario J—with tissue-dependent mapping and stopping criteria in place at 7.5% noise, (e) Scenario K—with tissue-dependent mapping and stopping criteria in place at 10% noise. See Table 3 for further details.
Figure 11. Final results of DGM-CSI frequency-hopping and frequency-cycled complex dielectric property reconstruction of open-boundary Category D synthetic breast model using 1.0–3.0 GHz data: (a) Scenario G—without tissue-dependent mapping at 400 iterations per frequency terminating after first inversion of 3.0 GHz (5% noise), (b) Scenario H—with tissue-dependent mapping and stopping criteria in place at 3% noise, (c) Scenario I—with tissue-dependent mapping and stopping criteria in place at 5% noise, (d) Scenario J—with tissue-dependent mapping and stopping criteria in place at 7.5% noise, (e) Scenario K—with tissue-dependent mapping and stopping criteria in place at 10% noise. See Table 3 for further details.
Jimaging 05 00055 g011
Figure 12. Relative error norms of frequency-hopping and frequency-cycled reconstructions of Category D model: Scenario G (blue solid line)—without tissue-dependent mapping at 400 iterations per frequency terminating after first inversion of 3.0 GHz (5% noise), Scenario H (red dashed line)—with tissue-dependent mapping and stopping criteria in place at 3% noise, Scenario I (black dotted line)—with tissue-dependent mapping and stopping criteria in place at 5% noise, Scenario J (magenta solid line)—with tissue-dependent mapping and stopping criteria in place at 7.5% noise, Scenario K (green dashed line)—with tissue-dependent mapping and stopping criteria in place at 10% noise. See Table 3 for further details.
Figure 12. Relative error norms of frequency-hopping and frequency-cycled reconstructions of Category D model: Scenario G (blue solid line)—without tissue-dependent mapping at 400 iterations per frequency terminating after first inversion of 3.0 GHz (5% noise), Scenario H (red dashed line)—with tissue-dependent mapping and stopping criteria in place at 3% noise, Scenario I (black dotted line)—with tissue-dependent mapping and stopping criteria in place at 5% noise, Scenario J (magenta solid line)—with tissue-dependent mapping and stopping criteria in place at 7.5% noise, Scenario K (green dashed line)—with tissue-dependent mapping and stopping criteria in place at 10% noise. See Table 3 for further details.
Jimaging 05 00055 g012
Table 1. Reconstruction progression of open boundary scenarios (Figure 7).
Table 1. Reconstruction progression of open boundary scenarios (Figure 7).
Frequency:
ScenarioTMSCNo. of Iterations (Stopping Condition)Total
Components Reconstructed
1.0 GHz:2.0 GHz:3.0 GHz:
ANoNo400 (F)400 (F)400 (F) 1200
Re, ImRe, ImRe, Im
1.0 GHz:2.0 GHz:3.0 GHz:1.0 GHz:2.0 GHz:3.0 GHz:
BYesNo200 (F)200 (F)200 (F)200 (F)200 (F)200 (F)1200
Re, ImReReReReRe
1.0 GHz:2.0 GHz:3.0 GHz:1.0 GHz:2.0 GHz:1.0 GHz:
CYesYes *140 (KS)190 (KS)255 (KS)53 (KS)89 (DE)41 (KS)768
Re, ImReReReReRe, Im
TM = tissue mapping; SC = stopping criteria; F = fixed; KS = Kolmogorov-Smirnov test window; DE = Domain error; Re = Real part; Im = Imaginary part. * 80% of a 30-iteration data error window must reach or surpass a K-S test p-value threshold of 0.99.
Table 2. Reconstruction progression of PEC-bounded scenarios (Figure 10).
Table 2. Reconstruction progression of PEC-bounded scenarios (Figure 10).
Frequency:
ScenarioTMSCNo. of Iterations (Stopping Condition)Total
Components Reconstructed
1.0 GHz:1.25 GHz:1.5 GHz:
DNoNo400 (F)400 (F)400 (F) 1200
Re, ImRe, ImRe, Im
1.0 GHz:1.25 GHz:1.5 GHz:1.0 GHz:1.25 GHz:1.5 GHz:
EYesNo200 (F)200 (F)200 (F)200 (F)200 (F)200 (F)1200
Re, ImReReReReRe
1.0 GHz:1.25 GHz:1.5 GHz:1.0 GHz:1.0 GHz:
FYesYes *123 (KS)124 (KS)97 (KS)14 (DE)48 (DE) 406
Re, ImReReReRe, Im
TM = tissue mapping; SC = stopping criteria; F = fixed; KS = Kolmogorov-Smirnov test window; DE = Domain error; Re = Real part; Im = Imaginary part. * 80% of a 30-iteration data error window must reach or surpass a K-S test p-value threshold of 0.99.
Table 3. Reconstruction progression of Category D model noise-variant scenarios (Figure 12).
Table 3. Reconstruction progression of Category D model noise-variant scenarios (Figure 12).
Frequency:
Scenario (Noise %)TMSCNo. of Iterations (Stopping Condition)Total
Components Reconstructed
1.0 GHz:2.0 GHz:3.0 GHz:
G (5%)NoNo400 (F)400 (F)400 (F) 1200
Re, ImRe, ImRe, Im
1.0 GHz:2.0 GHz:3.0 GHz:1.0 GHz:2.0 GHz:1.0 GHz:
H (3%)YesYes *155 (KS)140 (KS)265 (KS)50 (KS)52 (DE)43 (DE)759
Re, ImReReReReRe, Im
1.0 GHz:2.0 GHz:3.0 GHz:1.0 GHz:1.0 GHz:
I (5%)YesYes *139 (KS)140 (KS)265 (KS)34 (DE)8 (DE) 586
Re, ImReReReRe, Im
1.0 GHz:2.0 GHz:3.0 GHz:1.0 GHz:1.0 GHz:
J (7.5%)YesYes *91 (KS)112 (KS)234 (KS)37 (DE)13 (DE) 487
Re, ImReReReRe, Im
1.0 GHz:2.0 GHz:3.0 GHz:1.0 GHz:1.0 GHz:
K (10%)YesYes *87 (KS)90 (KS)228 (KS)11 (DE)32 (DE) 448
Re, ImReReReRe, Im
TM = tissue mapping; SC = stopping criteria; F = fixed; KS = Kolmogorov-Smirnov test window; DE = Domain error; Re = Real part; Im = Imaginary part. * 80% of a 30-iteration data error window must reach or surpass a K-S test p-value threshold of 0.99.

Share and Cite

MDPI and ACS Style

Kaye, C.; Jeffrey, I.; LoVetri, J. Novel Stopping Criteria for Optimization-Based Microwave Breast Imaging Algorithms. J. Imaging 2019, 5, 55. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging5050055

AMA Style

Kaye C, Jeffrey I, LoVetri J. Novel Stopping Criteria for Optimization-Based Microwave Breast Imaging Algorithms. Journal of Imaging. 2019; 5(5):55. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging5050055

Chicago/Turabian Style

Kaye, Cameron, Ian Jeffrey, and Joe LoVetri. 2019. "Novel Stopping Criteria for Optimization-Based Microwave Breast Imaging Algorithms" Journal of Imaging 5, no. 5: 55. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging5050055

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