Next Article in Journal
Serum Extracellular Vesicle-Derived miRNAs in Patients with Non-Small Cell Lung Cancer—Search for Non-Invasive Diagnostic Biomarkers
Next Article in Special Issue
A Novel Approach on Microwave Hyperthermia
Previous Article in Journal
Immature Plasma Cell Myeloma Mimics Metastatic Renal Cell Carcinoma on 18F-PSMA-1007 PET/CT Due to Endothelial PSMA-Expression
Previous Article in Special Issue
A Simulation-Based Methodology of Developing 3D Printed Anthropomorphic Phantoms for Microwave Imaging Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metasurface-Enhanced Antennas for Microwave Brain Imaging

1
Faculty of Natural and Mathematical Sciences, King’s College London, Strand, London WC2R 2LS, UK
2
Metamaterial Inc., Dartmouth, NS B2Y 4M9, Canada
*
Authors to whom correspondence should be addressed.
Diagnostics 2021, 11(3), 424; https://doi.org/10.3390/diagnostics11030424
Submission received: 2 February 2021 / Revised: 22 February 2021 / Accepted: 28 February 2021 / Published: 3 March 2021
(This article belongs to the Special Issue Electromagnetic Imaging for a Novel Generation of Medical Devices)

Abstract

:
Stroke is a very frequent disorder and one of the major leading causes of death and disability worldwide. Timely detection of stroke is essential in order to select and perform the correct treatment strategy. Thus, the use of an efficient imaging method for an early diagnosis of this syndrome could result in an increased survival’s rate. Nowadays, microwave imaging (MWI) for brain stroke detection and classification has attracted growing interest due to its non-invasive and non-ionising properties. In this paper, we present a feasibility study with the goal of enhancing MWI for stroke detection using metasurface (MTS) loaded antennas. In particular, three MTS-enhanced antennas integrated in different brain scanners are presented. For the first two antennas, which operate in a coupling medium, we show experimental measurements on an elliptical brain-mimicking gel phantom including cylindrical targets representing the bleeding in haemorrhagic stroke (h-stroke) and the not oxygenated tissue in ischaemic stroke (i-stroke). The reconstructed images and transmission and reflection parameter plots show that the MTS loadings improve the performance of our imaging prototype. Specifically, the signal transmitted across our head model is indeed increased by several dB‘s over the desired frequency range of 0.5–2.0 GHz, and an improvement in the quality of the reconstructed images is shown when the MTS is incorporated in the system. We also present a detailed simulation study on the performance of a new printed square monopole antenna (PSMA) operating in air, enhanced by a MTS superstrate loading. In particular, our previous developed brain scanner operating in an infinite lossy matching medium is compared to two tomographic systems operating in air: an 8-PSMA system and an 8-MTS-enhanced PSMA system. Our results show that our MTS superstrate enhances the antennas’ return loss by around 5 dB and increases the signal difference due to the presence of a blood-mimicking target up to 25 dB, which leads to more accurate reconstructions. In conclusion, MTS structures may be a significant hardware advancement towards the development of functional and ergonomic MWI scanners for stroke detection.

1. Introduction

Stroke is a clinical syndrome causing acute dysfunction of a brain area and it is linked to a vascular mechanism resulting in a vessel occlusion or rupture. This definition of stroke refers to ischemic stroke (i-stroke), which involves a central nervous system infarction, but also broadly includes intracerebral hemorrhage (h-stroke) and subarachnoid hemorrhage [1]. Stroke is a very frequent disorder and a major cause of adult disability and death. It is the fourth most common cause of death in the UK, with more than 100,000 cases each year and a cost to society around £26 billion a year [2].
An early diagnosis of stroke is necessary, as the brain loses millions of cells every second after the vessel’s occlusion or rupture and this may result in a permanent damage or even death. Additionally, stroke’s treatments differ based on the stroke type, meaning that an incorrect diagnosis could be lethal for the patient [3]. Thus, a fast access to efficient imaging tools is needed to timely initiate the correct treatment. Also, a strict follow-up in the first days (especially for treated patients) would also be desirable to assess the stroke’s evolution. Currently, stroke detection relies on technologies such as computed tomography (CT), magnetic resonance imaging (MRI) and cerebral angiography. Both CT and MRI are able to confirm the diagnosis and identify the stroke’s site. However, both these tools are not portable and thus cannot be easily used by paramedics at the patient’s bedside or found inside ambulances [4]. Moreover, whereas CT is ionizing and has a limited sensitivity for early ischaemic signs, while MRI exposes the patient to a strong magnetic field and needs more time and patient’s cooperation to be performed [5]. More importantly, MRI is a high cost technology and its availability is limited for emergency imaging.
Microwave imaging (MWI) could be a valid alternative to the current diagnostic methods, as it is not invasive, uses non-ionizing radiations and has a data acquisition time which ranges from milliseconds to a few seconds. In addition, it is a technology which can lead to low-cost and portable brain scanners [6]. This imaging technique is based on the difference in the tissues’ dielectric properties, which results in perturbations of the scattered field. The tissue dielectric contrast can be estimated using radar-based or tomographic reconstruction algorithms that are applied to the acquired data [7]. In particular, in microwave tomography (MWT), it is possible to reconstruct a map of the spatial distribution of the dielectric properties of the region of interest through solving an inverse ill-posed electromagnetic (EM) problem [8,9,10].
A MWI setup consists of several antennas placed around the body to be imaged. Typically, there is a homogeneous matching medium between the antennas and the body. Using a lossy matching medium minimizes reflection and couples the transmitted power to the body, but it also attenuates the useful signal transmitted into the body [7]. As a result, the coupling medium affects the detection of useful scattered “weak” signals. As the antenna array that acquires these signals is confined in a small region surrounding the head, small and compact antennas are commonly used [11,12,13]. In addition to small size, MWT antenna arrays for brain imaging must operate in the 0.5–2.0 GHz frequency range to achieve an optimal trade-off between resolution and penetration depth [14,15,16]. Following these requirements, several MWI head scanners including helmets and headbands [17] have been considered.
To couple the transmitted power into the body and improve the detection of useful scattered “weak” signals, we could take advantage of metasurface (MTS) technology. As demonstrated in our previous studies [18,19], one way to improve the efficacy of a MWT system is incorporating MTS structures into its array, adjacent to the substrate of each antenna. Using this hardware configuration, we have obtained higher quality reconstructed images of a blood-mimicking target placed inside the brain volume of our numerical head model. MTS’s are the two-dimensional equivalent of metamaterials (MM’s), which are engineered materials made of sets of small scatterers or apertures arranged in a regular array throughout a region of space. MTS’s are mainly based on sub-wavelength split-ring resonators (SRRs) which can be designed in order to obtain some desirable EM behaviour which is not found in naturally occurring materials, such as negative refractive index or near-zero constitutive parameters [20]. Thus, we can take advantage of MTS structures to manipulate size, efficiency, bandwidth, and directivity of several EM systems [21,22]. For instance, zero-index MM’s (materials with constitutive parameters ϵ ’ and μ of zero or near-zero values) and negative index MM’s have shown a strong potential in several applications and can be used to fabricate high directivity antennas [23,24,25]. At microwave frequencies, MTS-enhanced antennas are typically modelled tailoring the MTS’s design in order to interact with the emitted wave to form the desired radiation pattern [26]. Recently, layers of MTS have been used in a variety of clinical applications such as biomedical imaging and sensing. For instance, SRRs have been used to enhance the sensing properties of existing biosensors [27,28]. Likewise, MTS’s have been used as flat lenses for near-field imaging [29] and MRI [30].
With the aim of improving the accuracy of MWI for brain stroke detection and localization, we propose three MTS-enhanced antennas. To this end, this work validates our previous designed MTS film by presenting experiments with an elliptical brain phantom including cylindrical stroke-mimicking targets. Our experimental results suggest that the MTS film employed to enhance the sensitivity of our custom-made MWT prototype [31] can have a positive impact when placed on the head, closely fixed to the MWI array’s antennas. Then, we show a detailed simulation study on the performance of a new MTS-enhanced printed square monopole antenna (PSMA) operating in air. Our results show that we can detect a blood-mimicking target placed inside the brain volume of our head model avoiding the use of a thick and bulky matching medium. In addition, the MTS superstrate loading enhances the dielectric contrast between the blood-mimicking target and the surrounding tissue.
The remainder of the paper is organized as follows: Section 2 presents the methodology for the MTS-enhanced antenna designs. It also describes the geometry of the custom-made prototype used to validate our previous MTS designs and the simulation setups used to test the MTS-enhanced PSMA. Section 3 presents reflection and transmission plots from experiments and images reconstructed through our previously developed 2D DBIM-TwIST algorithm. It also includes simulation results which suggest that the PSMA’s MTS loading increases the MWI system’s sensitivity to the signal scattered by a blood-mimicking target. Finally, Section 4 includes a discussion, concluding remarks and future work.

2. Materials and Methods

2.1. MTS-Enhanced Antennas Immersed in a Coupling Medium

The triangular patch antenna and the spear patch antenna used in our experiments are shown in Figure 1a,b. These printed monopole antennas were modelled on an FR-4 substrate with a partial ground on the back side. To couple energy more efficiently into the head, both the antennas were designed to operate in a matching medium made of a 90% glycerol-water mixture [17,32,33]. The dimensions of the substrate, patch and transmission line of the antennas are shown in Figure 1a,b.
A 24.75 mm × 29.7 mm MTS superstrate loading was modelled to operate as an enhancer of the antennas described above and glued on the antennas’ radiating elements. This MTS structure is based on the unit cell geometry shown in Figure 2a, which comprises a Jerusalem Cross-shaped copper lattice (thickness = 0.10 mm) embedded between two Rogers 3010 TM substrates (thickness = 1.27 mm, ϵ = 10.2 and tan δ = 0.0022). These two high dielectric substrates are bonded with with Rogers 3001 bonding film ( ϵ = 2.28, tan δ = 0.003). In previous work [18,34], we have performed several simulation studies to optimize this geometry for operating in contact with the human skin when immersed in our coupling medium.

2.2. MTS-Enhanced Antennas: Experimental Validation in a Coupling Medium

To test the MTS-enhanced antennas operating in our coupling medium, we carried out experiments using the custom-made MWT prototype shown in Figure 3a. This setup includes an elliptical array placed inside a 300 mm diameter cylindrical tank and connected to a multiport Keysight M9019A Vector Network Analyser (VNA). Measurements were performed by immersing eight transceivers inside the tank, in a 90% glycerol solution used as matching medium. The antennas were placed as close as possible to the phantom’s external surface. To adjust the antennas’ position, we used horizontal and vertical mounts. By sequentially transmitting from one antenna and receiving by the others, the eight-antenna array produces a scattering matrix to be fed into our DBIM-TwIST algorithm for image reconstructions (see Appendix A).
To examine the performance of the MTS-enhanced triangular patch antenna, two gelatin-oil phantoms mimicking the properties of average brain and blood tissue were fabricated, as described in [35]. A 3D printed elliptic mould made of ABS was used as a holder for our brain liquid phantom that, once solidified, can mimic the brain. We first measured the average brain phantom (“no target” scenario) using both the triangular patch antennas and MTS-enhanced antennas. Then, we inserted a 30 mm diameter cylindrical inclusion resembling blood tissue into the brain-mimicking mixture and performed the “with target” measurements.
To test the MTS-enhanced spear antenna, we fabricated new gelatin-oil phantoms mimicking the properties of average brain, blood tissue and ischaemia. After carrying on the “no target” measurements, we inserted the 30 mm inclusion resembling the bleeding in the brain phantom and we performed the “with target” measurements. Then, we substituted the blood-mimicking inclusion with the ischaemia-mimicking target and we carried out a second round of “with target” measurements.
All the inclusions were positioned at the same coordinates, at an angular position of about 320 degrees relative to the position of the first transmitting antenna, as shown in Figure 3b. The permittivity at 1 GHz of the coupling liquid and each phantom used in these experiments is shown in Table 1.

2.3. MTS-Enhanced Printed Square Monopole Antenna Operating in Air

To investigate the possibility of applying MWT without employing any matching media in between the antennas and the human head, we modelled the printed square monopole antenna (PSMA) shown in Figure 4a using CST Microwave Studio. The proposed PSMA is operating in air and is specifically designed to work in close contact with our numerical head phantom. It is designed on an RT/duroid 5880 LZ low dielectric substrate (thickness = 1.026 mm, ϵ ’ = 1.96, tan δ = 0.0019) and is based on a square patch with a 10 mm partial ground on the back side.
After assessing the PSMA’s performance, we designed the MTS-enhanced PSMA shown in Figure 3b. This MTS-enhanced antenna includes a 10 × 19 unit cell MTS film covering almost the whole antenna’s surface and bonded to the front side of the PSMA through Rogers 3001 bonding film (thickness = 0.03 mm). The MTS unit cell design is based on a 0.01 mm thick Double Square Split-Ring Resonator (DS-SRR) printed on a square RT/duroid 5880 LZ (thickness = 0.635 mm).The dimensions of the DS-SRR are shown in Figure 4b.
The waveguide setup shown in Figure 5a was used to characterize the proposed MTS unit cell structure. The left and right surfaces are assigned as PEC walls, whereas the top and bottom sides (perpendicular to the two waveguide’s ports) are PMC walls. Transmission and reflection were calculated through CST Microwave Studio. The material properties extracted from the S-Parameters are shown in Figure 5.

2.4. PSMA and MTS-Enhanced PSMA: Comparison between Different MWT Scanners

To test the performance of our new PSMA and MTS-enhanced PSMA operating in air, we compared our previous developed brain scanner operating in an infinite lossy matching medium to two tomographic systems: an 8-PSMA system and an 8-MTS-enhanced PSMA system. Our previous developed MWT scanner for brain imaging consists of 12 spear patch antennas immersed in a 90% glycerol-water mixture and placed around EN 50361 Specific Anthropomorphic Mannequin (SAM) head model [18,19]. A similar setup (“System 1”), including 8 antennas placed elliptically around our numerical head model, was modelled in CST Microwave Studio. The antenna array was immersed in an infinite mixture of 90% glycerol-water matching medium. Using the same head model, which is made of a nylon mould ( ϵ ’ = 3.2, tan δ = 0.013) containing an average brain numerical phantom inside ( ϵ ’ = 45.8 and σ = 0.76 S/m), we studied other two MWT scanners operating in air: “System 2”, which comprises 8 PSMA, and “System 3”, which includes 8 MTS-enhanced PSMA. In these MWT systems, the antennas were closely fixed to the head’s surface. The MWT antenna arrays described above are shown in Figure 6.
For each antenna array, the S-Parameters were measured over the 0.5–2.0 GHz frequency range for the “no target” configuration (head model including average brain only) and the “with target” configuration. This last configuration includes a cylindrical blood-mimicking target (30 mm diameter and 35 mm height) inserted inside the brain volume and placed in the back side of the head model, close to antenna 1. Then, the signal difference “with target–no target” (dB) at relevant frequencies was calculated. Finally, images were reconstructed at the single frequencies 0.9 GHz, 1.0 GHz and 1.1 GHz, by applying our 2D DBIM-TwIST algorithm [36] to the simulated data.

3. Results

3.1. Experimental Results for the MTS-Enhanced Antennas Operating in a Coupling Medium

In our previous study [19], we have already shown the benefits of using MTS-enhanced antennas. In this section, we present experimental results for both the arrays described in Section 2.2. The reconstructed images and transmission and reflection parameter plots suggest that the MTS superstrate loading illustrated in Figure 2 has the potential to improve the performance of our custom-made tomographic system. The S-Parameter plots are shown for the “no target” and “with target” configuration, with and without MTS superstrate loading. Moreover, to provide more complete information on the enhancement of the system response in the presence of the MTS, we also present field amplitude distributions from accurate numerical simulations of our system interacting with the brain phantom in Section 3.3.

3.1.1. 8-MTS-Enhanced Triangular Patch Antenna Array

Figure 7a shows the reflection coefficient for one of the triangular patch antennas (number 8 in Figure 6b), which is placed in front of the blood-mimicking target, while Figure 7b shows an example of transmission signal levels (in dB) across the brain phantom, along the direction containing the target. The antenna’s reflection coefficient in the “no target” configuration is reduced by around 3 dB, suggesting that the antenna’s matching is enhanced at the operating frequency when the MTS is present. Moreover, the plotted transmission parameter is improved of around 8 dB over all the considered frequency range.
Two-dimensional (2D) multi-frequency (frequency hopping) reconstructions were performed through our DBIM-TwIST algorithm [36], which is discussed in detail in Appendix A. The reconstruction results of the dielectric contrast between the blood-mimicking target and the surrounding brain tissue are shown in Figure 8, where three frequencies (0.7 GHz, 1.1 GHz and 1.3 GHz) were used. An improvement in the quality of the reconstructed images is shown when the MTS layers are added to each antenna element. In particular, a reduction of artefacts and a better localization of the blood target are observed.

3.1.2. 8-MTS-Enhanced Spear Patch Antenna Array

Figure 9a plots the reflection coefficient vs. frequency for the spear patch antenna number 8, while Figure 9b plots the transmission coefficient vs. frequency for a pair of antennas opposite the brain phantom along the direction containing the haemorrhage-mimicking or ischaemic-mimicking target. These plots suggest that similar improvement in performance is achieved if the MTS is used with the spear antennas, similar to the case of the triangular antennas.
The blood-mimicking target and the ischaemia-mimicking target were reconstructed via our 2D DBIM-TWiST algorithm. Multi-frequency reconstructions were performed using frequencies 0.7 GHz, 1.1 GHz and 1.3 GHz. The reconstructed blood phantom’s permittivity and ischaemia-mimicking tissue’s permittivity are shown in Figure 10 and Figure 11, respectively. A better localization of the targets and a reduction of the artefacts are shown when the MTS films are integrated in the setup.

3.2. PSMA and MTS-Enhanced PSMA: Simulation Results with Different Brain Scanners

For each of the systems described in Section 2.4, the S-Parameters were calculated through CST Microwave Studio over the 0.5–2.0 GHz frequency range. Then, the signal difference “with target - no target” (dB) at relevant frequencies was plotted. Figure 12 shows the reflection parameter for each of the antennas of “System 1”, “System 2” and “System 3”. The plots show that presence of the MTS superstrate loading improves the reflection coefficient of the PSMA by 5 dB.
Figure 13 shows the signal difference “with target – no target” (dB) as a function of receiver location at frequencies 1 GHz and 1.1 GHz for “System 1”, “System 2” and “System 3”. As shown in the graphs, when the PSMA is loaded with our MTS superstrate, it leads to an overall improvement of the signal difference (up to 25 dB) compared to “System 2”. Also, the coupling liquid affects the “weak” signal scattered from the blood-mimicking inclusion, which falls below the noise level.
Finally, to test the MTS’s impact on tomographic reconstructions, we have applied our 2D DBIM-TwIST algorithm to the simulated data plotted in Figure 13. We carried out single-frequency reconstructions, assuming approximate knowledge of the dielectric properties of plastic and average brain tissue as initial guess for our algorithm. The estimated dielectric properties at frequencies 0.9 GHz, 1 GHz and 1.1 GHz is shown in Figure 14. When employing our new PSMA operating in air (“System 2”), the blood-mimicking target is detected correctly. However, a higher contrast between the haemorrhage-mimicking target and the surrounding brain tissue is observed when the MTS is integrated in the scanner (“System 3”).

3.3. Near-Field Analysis of the MTS-Enhanced Antennas

In near-field MWI, spatial distributions of the electric field magnitude are important measures of antenna performance, which is affected by the impact of evanescent k-space contributions [37]. To evaluate performance in the near field, we have used our simulation setup to calculate examples of spatial distributions due to the transmitting antenna with and without the MTS loading. Examples of these distributions for the real part of the electric field calculated at the antennas’ resonant frequencies are shown in Figure 15 and Figure 16 for the triangular patch antenna and the spear patch antenna, respectively. These results confirm that the antenna-radiated probing fields are significantly stronger in the presence of the MTS superstrate loading, leading to an increased sensitivity of our MWI system.

4. Discussion

In this paper, we have presented experiments including brain tissue-mimicking phantoms and cylindrical targets mimicking strokes as well as numerical simulations with several brain scanners. Our study demonstrates that MTS structures have the potential of improving brain imaging when integrated in MWT scanners. The MTS unit cell shown in Figure 2a and the MTS superstrate loadings based on this design were studied experimentally using the custom-made prototype described in Section 2.2. The MTS loadings were shown to enhance transmission and improve the matching of our in-house fabricated antennas in the 0.5–2.0 GHz frequency range. In particular, the signal transmitted across the head phantom, along the direction containing the blood-mimicking target, is overall increased over all the frequency range. Moreover, the antenna reflection coefficient is reduced by several dB when the MTS films are incorporated into the system. This improvement in the near-field antennas performance translates into higher-quality reconstruction images.
Furthermore, an innovative MTS-enhanced antenna design was proposed to test our 2D DBMI-TwIST algorithm in air. To this end, three MWT scanners were tested and compared using a numerical head model made of average brain tissue and including a blood-mimicking target. For each system, the S-Parameters were measured over the 0.5–2.0 GHz frequency range and the signal difference “with target–no target” (dB) at relevant frequencies was calculated. Single frequency reconstructions were also performed. Our results indicate that, when our PSMA is loaded with the MTS superstrate, the signal difference due to the presence of the target is increased up to 25 dB. This translates into a higher contrast between the target and the surrounding tissue. Another interesting point is to note that the signal scattered by the blood-mimicking inclusion is severely affected by the coupling medium. Thus, it is possible to detect a haemorrhage-mimicking target placed inside the brain volume of our numerical head model without using a bulky matching liquid.
To the best of our knowledge, this is the first study presenting experimental and simulation results which demonstrate the possibility of enhancing the detection of stroke-mimicking targets by incorporating MTS antenna enhancers for different MWT arrays. While previous work by various research groups has argued that MWT is applicable for brain imaging (e.g., [38]), improving the sensitivity and hence the quality of the data is critical to achieve clinical accuracy. This feasibility study suggests that MTS technology might be an important step towards the goal of developing functional, portable, and ergonomic MWI scanners with the desired clinical accuracy.
It is of course important to note that most clinical MWI applications involve a three-dimensional (3D) non-linear inverse problem, which is complicated due to the complex structure of the brain and the high dielectric contrast between its different tissues [39]. Thus, in addition to optimising the data acquisition system, the ability of the DBIM-TwIST algorithm to tackle this highly non-linear problem within the distorted Born approximation must be demonstrated. Our previous work has shown that the algorithm can recover complex structures of high dielectric contrast in realistic numerical breast phantoms [36], and some initial simulation results with complex brain phantoms have also shown promise [40]. These studies have shown the importance of the initial guess and prior information as well as frequency hopping in getting the DBIM-TwIST algorithm to converge to meaningful reconstructions of these complex structures.
Our vision of developing a portable microwave stroke detection scanner must of course tackle several additional challenges. These include the requirement to place antenna sensors conformally to the head surface, uncertainty and variability in the properties of the skin-hair-scalp interface, and inhomogeneities and variability in the dielectric properties of brain tissues (including the presence of blood vessels), which are lossy and result in compromising imaging resolution to ensure sufficient penetration. These result in a non-linear, very challenging 3D EM inverse scattering problem, which will certainly require some prior information on the brain tissues distribution and properties. Regardless of the challenges in microwave stroke detection, we must emphasise, however, that our results suggest that MTS-enhanced arrays can be advantageous for various other medical MWI applications (e.g., breast cancer detection) and algorithms (e.g., radar-based methods).
Based on these considerations, our future work will focus on conducting more experimental studies using more complex setups. Regarding the antenna arrays immersed in the coupling liquid, we aim to validate our MTS design using the 3D version of our tomographic algorithm, by employing 24-antenna arrays distributed in two or three rows around the head phantom. Then, we will work towards optimizing the MTS-enhanced PSMA and carrying out experimental studies to assess the feasibility of applying our algorithm without employing any matching media. Before evaluating our MTS-enhanced system with clinical data, we also plan to fabricate and test our system with more realistic inhomogeneous brain phantoms. In particular, we will use digital volumetric brain phantoms with complex structure and different tissues (e.g., grey matter, white matter, bone, muscle and skin) to test the capabilities and requirements of our DBIM-TwIST approach.

Author Contributions

Conceptualization, E.R. and P.K.; methodology, E.R., P.L., W.G., O.K., I.S. and N.G.; validation, E.R., P.L., O.K., I.S. and N.G.; resources, P.K., E.K. and G.P.; writing—original draft preparation, E.R.; supervision, P.K.; project administration, P.K., E.K. and G.P.; funding acquisition, P.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the EMERALD project funded from the European Union’s Horizon 2020 under the Marie Skłodowska-Curie grant agreement No. 764479. The research was also funded in part by Innovate UK grant number 103920 and in part by the Engineering and Physical Sciences Research Council grant number EP/R013918/1.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; and in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
MWIMicrowave Imaging
MTSMetasurface
PSMAPrinted Square Monopole Antenna
CTComputed Tomography
MRIMagnetic Resonance Imaging
MWTMicrowave Tomography
EMElectromagnetic
SRRSplit-Ring Resonator
VNAVector Network Analyzer
DS-SRRDouble-Square-Split-Ring Resonator
PECPerfect Electric Conductor
PMCPerfect Magnetic Conductor
SAMSpecific Anthropomorphic Mannequin

Appendix A. The DBIM-TwIST Algorithm

The distorted Born iterative method (DBIM) is an iterative algorithm for solving the electromagnetic (EM) inverse scattering problem, which is used for for medical imaging applications to estimate the spatial distribution of the tissue’s dielectric properties inside a region V in the human body (see, for example [41]). The nonlinear inverse scattering problems for each transmitter-receiver (TR) pair is linearized and approximated by the following equation,
E s ( r n , r m ) = E t ( r n , r m ) E b ( r n , r m ) = ω 2 μ 0 ϵ 0 V G b ( r n , r ) E b ( r , r m ) δ ϵ ( r ) d r ,
where E t , E s , E b are the total, scattered and background electric fields respectively, r n and r m are the transmitting and receiving antenna locations, ω is the angular frequency, μ 0 is the permeability of free space, ϵ 0 is the permittivity of free space, and G b is the Green’s function for the background medium. As we only consider point sources, the Green’s functions can be calculated from the electric field as, G b ( r n , r ) = j ω μ 0 J E b ( r , r n ) . The difference δ ϵ between the relative complex permittivity of the reconstructed region, ϵ r ( r ) , and the background medium, ϵ b ( r ) , is defined as δ ϵ ( r ) = ϵ r ( r ) ϵ b ( r ) .We note that the scalar integral equation above assumes two-dimensional (2D) TM propagation, and is only an approximation of the three-dimensional (3D) inverse problem at hand. Despite this loss in information, the 2D approximation can produce images of acceptable quality in many microwave imaging problems arising in medical applications.
At each DBIM iteration i, the integral equation can be discretized for each TR pair as,
E s ( r n , r m ) = j ω ϵ 0 δ ϵ ( r ) E b ( r n , r ) E b ( r , r m ) ,
leading to an ill-posed linear system as,
A δ ϵ = b
where A is an M × N matrix ( M N ) , with M transmit-receive pairs and N voxels of the reconstruction region V, b is a M × 1 vector of the scattered fields. The A is calculated at each DBIM iteration by the forward solver which yields E b for a known background ϵ b . The background field is used to build the linear system above, which is then solved by an inverse solver. Finally, the background profile is updated by ϵ b i + 1 ( r ) = ϵ b i ( r ) + δ ϵ ( r ) and the DBIM continues to next iteration i + 1 .
The forward solver of our DBIM-TwIST algorithm uses the finite difference time domain (FDTD) method. The FDTD method simulates the EM wave propagation of the direct, “forward” problem based on Maxwell’s equations. As mentioned above, this work considers only 2D transverse magnetic (TM) waves to reconstruct 2D geometries. Furthermore, our FDTD implementation uses a single-pole Debye model to model frequency-dependent materials such as brain tissues as, ε ε 0 = ε + Δ ε 1 + j ω τ + σ s j ω ε 0 .
We employ the two-step iterative shrinkage/thresholding (TwIST) algorithm as the solver of the ill-posed linear inverse problem at each DBIM iteration. Thresholding algorithms solve the ill-conditioned linear system A x = b , by finding a solution x which minimizes the least squares error function as F ( x ) = 1 2 A x b 2 + λ x 1 with a regularization term λ x 1 to stabilize the solution by limiting its l 1 -norm. The general structure of the TwIST algorithm for solving this minimization problem is given by [42],
x t + 1 = ( 1 α ) x t 1 + ( α β ) x t + β C λ ( x t ) C λ ( x ) = Ψ λ ( x + A T ( y A x ) )
The parameters for the TwIST algorithm are calculated as κ = ξ 1 ξ m , ρ = 1 κ 1 + κ , α = ρ 2 + 1 , β = 2 α ξ 1 + ξ m , where ξ 1 and ξ m are the smallest and largest eigenvalues of A T A respectively. The shrinkage/thresholding operation is a soft-thresholding function, calculated as Ψ λ ( x ) = s i g n ( x ) max { 0 , | x | λ } . At each TwIST step, the new solution is updated based on two previous solutions and the soft-thresholding function. As the linear system can be extremely ill-conditioned for microwave imaging problems, the TwIST parameters above must be optimised specifically for the problem at hand [43].
The stopping criterion of the TwIST algorithm can be set based on a tolerance value, which is the normalised difference between the previous and current values of F ( x ) and is defined as t o l = F ( x k + 1 ) F ( x k ) F ( x k ) . The TwIST algorithm stops when t o l is smaller than a preset value, usually in range between 10 4 and 10 1 . This early termination of the iterative algorithm serves an additional regulariser, with a similar effect to a Tikhonov approach [41]. Our l 1 -problem can be regularised further by employing the Pareto curve method, which which yields the optimal tradeoff between the residual error A x b 2 and the norm of the solution x 1 (similar to the L-curve method for l 2 -problems) [36]. To reduce the computational cost, a practical approach for selecting λ is based on the form λ = δ A T b , where δ is factor with 0 < δ < 1 [36].

References

  1. Sacco, R.L.; Kasner, S.E.; Broderick, J.P.; Caplan, L.R.; Connors, J.J.; Culebras, A.; Elkind, M.S.; George, M.G.; Hamdan, A.D.; Higashida, R.T.; et al. An updated definition of stroke for the 21st century: A statement for healthcare professionals from the American heart association/American stroke association. Stroke 2013, 44, 2064–2089. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Stroke Association. Available online: https://www.stroke.org.uk/what-is-stroke/stroke-statistics (accessed on 18 December 2020).
  3. Powers, W.J.; Rabinstein, A.A.; Ackerson, T.; Adeoye, O.M.; Bambakidis, N.C.; Becker, K.; Biller, J.; Brown, M.; Demaerschalk, B.M.; Hoh, B.; et al. Guidelines for the Early Management of Patients With Acute Ischemic Stroke: A Guideline for Healthcare Professionals From the American Heart Association/American Stroke Association. Stroke 2018, 49, 46–99. [Google Scholar] [CrossRef] [PubMed]
  4. Mobashsher, A.T.; Nguyen, P.T.; Abbosh, A. Detection and localization of brain strokes in realistic 3-D human head phantom. In Proceedings of the 2013 IEEE MTT-S International Microwave Workshop Series on RF and Wireless Technologies for Biomedical and Healthcare Applications, IMWS-BIO, Singapore, 9–11 December 2013. [Google Scholar]
  5. Vymazal, J.; Rulseh, A.M.; Keller, J.; Janouskova, L. Comparison of CT and MR imaging in ischemic stroke. Insights Imaging 2012, 3, 619–627. [Google Scholar] [CrossRef] [Green Version]
  6. Hopfer, M.; Planas, R.; Hamidipour, A.; Henriksson, T.; Semenov, S. Electromagnetic tomography for detection, differentiation, and monitoring of brain stroke: A virtual data and human head phantom study. IEEE Trans. Antennas Propag. 2017, 59, 86–97. [Google Scholar] [CrossRef]
  7. De Azevedo-Marques Mazzoncini, P.; Mencattini, A.; Salmeri, M.; Rangayyan, R.M. Medical Image Analysis and Informatics: Computer-Aided Diagnosis and Therapy; CRC Press: Boca Raton, FL, USA, 2017; pp. 451–466. [Google Scholar]
  8. Wei, H.Y.; Soleimani, M. Electromagnetic tomography for medical and industrial applications: Challenges and opportunities. Proc. IEEE 2013, 101, 559–565. [Google Scholar] [CrossRef]
  9. Scapaticci, R.; Bucci, O.M.; Catapano, I.; Crocco, L. Differential microwave imaging for brain stroke followup. Int. J. Antennas Propag. 2014, 2014, 1–11. [Google Scholar] [CrossRef] [Green Version]
  10. Aldhaeebi, M.A.; Alzoubi, K.; Almoneef, T.S.; Bamatra, S.M.; Attia, H.; Ramahi, O.M. Review of microwaves techniques for breast cancer detection. Sensors 2020, 20, 2390. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Alqadami, A.S.M.; Nguyen-Trong, N.; Stancombe, A.E.; Bialkowski, K.; Abbosh, A. Compact Flexible Wideband Antenna for On-Body Electromagnetic Medical Diagnostic Systems. IEEE Trans. Antennas Propag. 2020, 68, 8180–8185. [Google Scholar] [CrossRef]
  12. Abbosh, A.M.; Kan, H.K.; Bialkowski, M.E. Compact ultra-wideband planar tapered slot antenna for use in a microwave imaging system. Microw. Opt. Technol. Lett. 2013, 55, 2562–2568. [Google Scholar] [CrossRef]
  13. Sohani, B.; Khalesi, B.; Ghavami, N.; Ghavami, M.; Dudley, S.; Rahmani, A.; Tiberi, G. Detection of haemorrhagic stroke in simulation and realistic 3-D human head phantom using microwave imaging. Biomed. Signal Proces. 2020, 61, 3–10. [Google Scholar] [CrossRef]
  14. Scapaticci, R.; Di Donato, L.; Catapano, I.; Crocco, L. A Feasibility Study on Microwave Imaging for Brain Stroke Monitoring. Prog. Electromagn. Res. 2014, 40, 305–324. [Google Scholar] [CrossRef] [Green Version]
  15. Rodriguez-Duarte, D.O.; Vasquez, J.A.T.; Scapaticci, R.; Crocco, L.; Vipiana, F. Brick-Shaped Antenna Module for Microwave Brain Imaging Systems. IEEE Antennas Wirel. Propag. Lett. 2020, 19, 2057–2061. [Google Scholar] [CrossRef]
  16. Zhurbenko, V. Challenges in the design of microwave imaging systems for breast cancer detection. Adv. Electr. Comp. 2011, 11, 91–96. [Google Scholar] [CrossRef] [Green Version]
  17. Wei, G.; Ahsan, S.; Kosmas, P. Portable Microwave Imaging Head Scanners for Intracranial Haemorrhagic Detection. In Proceedings of the Asia-Pacific Microwave Conference (APMC), Singapore, 10–13 December 2019. [Google Scholar]
  18. Razzicchia, E.; Sotiriou, I.; Cano-Garcia, H.; Kallos, E.; Palikaras, G.; Kosmas, P. Feasibility study of enhancing microwave brain imaging using metamaterials. Sensors 2019, 19, 5472. [Google Scholar] [CrossRef] [Green Version]
  19. Karadima, O.; Razzicchia, E.; Kosmas, P. Image Improvement Through Metamaterial Technology for Brain Stroke Detection. In Proceedings of the 14th European Conference on Antennas and Propagation (EuCAP 2020), Copenaghen, Denrmark, 15–20 March 2020. [Google Scholar]
  20. Holloway, C.L.; Kuester, E.F.; Gordon, J.A.; O’Hara, J.; Booth, J.; Smith, D.R. An overview of the theory and applications of metasurfaces: The two-dimensional equivalents of metamaterials. IEEE Trans. Antennas Propag. 2012, 54, 10–35. [Google Scholar] [CrossRef]
  21. Ziolkowski, R.W. Metamaterial-based antennas: Research and developments. IEICE Trans. Electron. 2006, E89-C, 1267–1275. [Google Scholar] [CrossRef]
  22. Ziolkowski, R.W.; Erentok, A. Metamaterial-based efficient electrically small antennas. IEEE Trans. Antennas Propag. 2006, 54, 2113–2130. [Google Scholar] [CrossRef] [Green Version]
  23. Ziolkowski, R.W. Propagation in and scattering from a matched metamaterial having a zero index of refraction. Phys. Rev. E 2004, 70, 1–15. [Google Scholar] [CrossRef] [PubMed]
  24. Li, L.W.; Li, Y.N.; Yeo, T.S.; Mosig, J.R.; Martin, O.J. A broadband and high-gain metamaterial microstrip antenna. Appl. Phys. Lett. 2010, 96, 5–8. [Google Scholar] [CrossRef]
  25. Mahmud, Z.; Islam, M.T. A Negative Index Metamaterial to Enhance the Performance of Miniaturized UWB Antenna for Microwave Imaging Applications. Appl. Sci. 2017, 7, 1149. [Google Scholar] [CrossRef] [Green Version]
  26. Bulu, I.; Caglayan, H.; Aydin, K.; Ozbay, E. Compact size highly directive antennas based on the SRR metamaterial medium. New J. Phys. 2005, 7, 1–4. [Google Scholar] [CrossRef] [Green Version]
  27. Jakšić, Z.; Vuković, S.; Matovic, J.; Tanasković, D. Negative refractive index metasurfaces for enhanced biosensing. Materials 2010, 4, 1–36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Puentes, M.; Maasch, M.; Schubler, M.; Jakoby, R. Frequency multiplexed 2-dimensional sensor array based on split-ring resonators for organic tissue analysis. IEEE Trans. Microw. Theory Tech. 2012, 60, 1720–1727. [Google Scholar] [CrossRef]
  29. Wang, G.; Fang, J.; Dong, X. Resolution of near-field microwave target detection and imaging by using flat LHM lens. IEEE Trans. Antennas Propag. 2007, 55, 3534–3541. [Google Scholar] [CrossRef]
  30. Slobozhanyuk, A.P.; Poddubny, A.N.; Raaijmakers, A.J.; Van Den Berg, C.A.; Kozachenko, A.V.; Dubrovina, I.A.; Melchakova, I.V.; Kivshar, Y.S.; Belov, P.A. Enhancement of Magnetic Resonance Imaging with Metasurfaces. Adv. Mater. 2016, 28, 1832–1838. [Google Scholar] [CrossRef]
  31. Karadima, O.; Ghavami, N.; Sotiriou, I.; Kosmas, P. Performance assessment of microwave tomography and radar imaging using an anthropomorphic brain phantom. In Proceedings of the 2020 XXXIIIrd General Assembly and Scientific Symposium of the International Union of Radio Science, Rome, Italy, 29 August–5 September 2020. [Google Scholar]
  32. Ahsan, S.; Guo, Z.; Miao, Z.; Sotiriou, I.; Koutsoupidou, M.; Kallos, E.; Palikaras, G.; Kosmas, P. Multiple-Frequency Microwave Tomography System. Sensors 2018, 18, 3491. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Ahsan, S.; Koutsoupidou, M.; Razzicchia, E.; Sotiriou, I.; Kosmas, P. Advances towards the Development of a Brain Microwave Imaging Scanner. In Proceedings of the 13th European Conference on Antennas and Propagation (EuCAP 2019), Krakow, Poland, 31 March–5 April 2019. [Google Scholar]
  34. Razzicchia, E.; Koutsoupidou, M.; Cano-garcia, H.; Sotiriou, I.; Kallos, E.; Palikaras, G. Metamaterial Designs to Enhance Microwave Imaging Applications. In Proceedings of the ICEAA 2019 International Conference on Electromagnetics in Advanced Applications, Granada, Spain, 9–13 September 2019. [Google Scholar]
  35. Lazebnik, M.; Madsen, E.L.; Frank, G.R.; Hagness, S.C. Tissue-mimicking phantom materials for narrowband and ultrawideband microwave applications. Phys. Med. Biol. 2005, 50, 4245–4258. [Google Scholar] [CrossRef]
  36. Miao, Z.; Kosmas, P. Multiple-Frequency DBIM-TwIST Algorithm for Microwave Breast Imaging. IEEE Trans. Antennas Propag. 2017, 65, 2507–2516. [Google Scholar] [CrossRef] [Green Version]
  37. Gaikovich, K.P.; Maksimovitch, Y.S.; Sumin, M.I. Inverse scattering problems of near-field subsurface pulse diagnostics. Inverse Probll. Sci. Eng. 2018, 26, 1590–1611. [Google Scholar] [CrossRef]
  38. Semenov, S.; Planas, R.; Hopfer, M.; Hamidipour, A.; Vasilenko, A.; Stoegmann, E.; Auff, E. Electromagnetic tomography for brain imaging: Initial assessment for stroke detection. In Proceedings of the IEEE Biomedical Circuits and Systems Conference: Engineering for Healthy Minds and Able Bodies, BioCAS, Atlanta, GA, USA, 22–24 October 2015. [Google Scholar]
  39. Semenov, S. Microwave tomography: Review of the progress towards clinical applications. Philos. Transact. A Math. Phys. Eng. Sci. 2009, 367, 3021–3042. [Google Scholar] [CrossRef] [Green Version]
  40. Lu, P.; Ahsan, S.; Kosmas, P. Preliminary Study on the Feasibility of Reconstructing Anatomically Complex Numerical Brain Phantoms with Limited Prior Information. In Proceedings of the 2020 33rd General Assembly and Scientific Symposium of the International Union of Radio Science, URSI GASS 2020, Rome, Italy, 29 August–5 September 2020. [Google Scholar]
  41. Shea, J.D.; Kosmas, P.; Hagness, S.C.; Van Veen, B.D. Three-dimensional microwave imaging of realistic numerical breast phantoms via a multiple-frequency inverse scattering technique. Med. Phys. 2010, 37, 4210–4226. [Google Scholar] [CrossRef]
  42. Bioucas-Dias, J.M.; Figueiredo, M.A. A new TwIST: Two-step iterative shrinkage/thresholding algorithms for image restoration. IEEE Trans. Image Process. 2007, 16, 2992–3004. [Google Scholar] [CrossRef] [Green Version]
  43. Miao, Z.; Kosmas, P. Microwave breast imaging based on an optimized two-step iterative shrinkage/thresholding method. In Proceedings of the 9th European Conference on Antennas and Propagation (EuCAP 2015), Lisbon, Portugal, 13–17 April 2015. [Google Scholar]
Figure 1. Dimensions of the substrate, patch and transmission line of (a) the spear patch antenna and (b) the triangular patch antenna.
Figure 1. Dimensions of the substrate, patch and transmission line of (a) the spear patch antenna and (b) the triangular patch antenna.
Diagnostics 11 00424 g001
Figure 2. (a) Dimensions of the MTS unit cell (b) Photograph of the triangular patch antenna and its MTS-enhanced version (c) Photograph of the spear patch antenna and its MTS-enhanced version.
Figure 2. (a) Dimensions of the MTS unit cell (b) Photograph of the triangular patch antenna and its MTS-enhanced version (c) Photograph of the spear patch antenna and its MTS-enhanced version.
Diagnostics 11 00424 g002
Figure 3. (a) Custom-made tomographic experimental setup (b) Elliptical mould containing the average brain phantom with the blood-mimicking target.
Figure 3. (a) Custom-made tomographic experimental setup (b) Elliptical mould containing the average brain phantom with the blood-mimicking target.
Diagnostics 11 00424 g003
Figure 4. (a) Dimensions of PSMA (b) PSMA with MTS superstrate loading and square DS-SRR based unit cell.
Figure 4. (a) Dimensions of PSMA (b) PSMA with MTS superstrate loading and square DS-SRR based unit cell.
Diagnostics 11 00424 g004
Figure 5. (a) Waveguide setup, with an incident wave propagating along Z-direction (b) Real and imaginary values of permittivity and permeability for the MTS.
Figure 5. (a) Waveguide setup, with an incident wave propagating along Z-direction (b) Real and imaginary values of permittivity and permeability for the MTS.
Diagnostics 11 00424 g005
Figure 6. (a) 8-PSMA array (“System 1”) (b) 8-MTS-Enhanced PSMA array (“System 3”) and its transversal view. “System 2” is equal to “System 2” except for including non-enhanced PSMAs.
Figure 6. (a) 8-PSMA array (“System 1”) (b) 8-MTS-Enhanced PSMA array (“System 3”) and its transversal view. “System 2” is equal to “System 2” except for including non-enhanced PSMAs.
Diagnostics 11 00424 g006
Figure 7. (a) Measured S 88 Parameter for the triangular patch antenna array in the “with target” (red lines) and “no target” (blue lines) configuration, with (solid lines) and without MTS (dashed lines) (b) Measured S 71 Parameter for the triangular patch antenna array in the “with target” (red lines) and “no target” (blue lines) configuration, with (solid lines) and without MTS (dashed lines).
Figure 7. (a) Measured S 88 Parameter for the triangular patch antenna array in the “with target” (red lines) and “no target” (blue lines) configuration, with (solid lines) and without MTS (dashed lines) (b) Measured S 71 Parameter for the triangular patch antenna array in the “with target” (red lines) and “no target” (blue lines) configuration, with (solid lines) and without MTS (dashed lines).
Diagnostics 11 00424 g007
Figure 8. Frequency hopping reconstructions of the blood-mimicking target obtained from experimental data for the triangular patch antenna array (a) without MTS and (b) with MTS.
Figure 8. Frequency hopping reconstructions of the blood-mimicking target obtained from experimental data for the triangular patch antenna array (a) without MTS and (b) with MTS.
Diagnostics 11 00424 g008
Figure 9. (a) Measured S 88 Parameter for the spear patch antenna array in the “with target” (red lines) and “no target” (blue lines) configuration, with (solid lines) and without MTS (dashed lines) (b) Measured S 71 Parameter for the spear patch antenna array in the “with target” (red lines) and “no target” (blue lines) configuration, with (solid lines) and without MTS (dashed lines).
Figure 9. (a) Measured S 88 Parameter for the spear patch antenna array in the “with target” (red lines) and “no target” (blue lines) configuration, with (solid lines) and without MTS (dashed lines) (b) Measured S 71 Parameter for the spear patch antenna array in the “with target” (red lines) and “no target” (blue lines) configuration, with (solid lines) and without MTS (dashed lines).
Diagnostics 11 00424 g009
Figure 10. Frequency hopping reconstructions of the blood-mimicking target obtained from experimental data for the spear patch antenna array (a) without MTS and (b) with MTS.
Figure 10. Frequency hopping reconstructions of the blood-mimicking target obtained from experimental data for the spear patch antenna array (a) without MTS and (b) with MTS.
Diagnostics 11 00424 g010
Figure 11. Frequency hopping reconstructions of the ischaemia-mimicking target obtained from experimental data for the spear patch antenna array (a) without MTS and (b) with MTS.
Figure 11. Frequency hopping reconstructions of the ischaemia-mimicking target obtained from experimental data for the spear patch antenna array (a) without MTS and (b) with MTS.
Diagnostics 11 00424 g011
Figure 12. Calculated reflection coefficients in the “no target” scenario for (a) the spear patch antennas, (b) the PSMAs and (c) the MTS-enhanced PSMAs.
Figure 12. Calculated reflection coefficients in the “no target” scenario for (a) the spear patch antennas, (b) the PSMAs and (c) the MTS-enhanced PSMAs.
Diagnostics 11 00424 g012
Figure 13. Calculated signal difference “with target - no target” (dB) for the 8-spear patch antenna array operating in infinite matching medium (“System 1”), the 8-PSMA array (“System 2”) and the 8-MTS-enhanced PSMA array (“System 3”) at (a) 1 GHz and (b) 1.1 GHz.
Figure 13. Calculated signal difference “with target - no target” (dB) for the 8-spear patch antenna array operating in infinite matching medium (“System 1”), the 8-PSMA array (“System 2”) and the 8-MTS-enhanced PSMA array (“System 3”) at (a) 1 GHz and (b) 1.1 GHz.
Diagnostics 11 00424 g013
Figure 14. Real part of permittivity reconstructed from simulation data at frequencies 0.9 GHz, 1 GHz and 1.1 GHz for (a) the 8-spear patch antenna array operating in infinite matching medium (“System 1”) (b) the 8-PSMA array (“System 2”) and (c) the 8-MTS-enhanced PSMA array (“System 3”) operating in air.
Figure 14. Real part of permittivity reconstructed from simulation data at frequencies 0.9 GHz, 1 GHz and 1.1 GHz for (a) the 8-spear patch antenna array operating in infinite matching medium (“System 1”) (b) the 8-PSMA array (“System 2”) and (c) the 8-MTS-enhanced PSMA array (“System 3”) operating in air.
Diagnostics 11 00424 g014
Figure 15. Real part of the E-field calculated at 0.9 GHz and plotted on a transverse section of the simulated brain phantom including a blood-mimicking target. A higher field intensity in the target region is shown when the MTS is placed adjacent on the triangular antenna’s substrate (right plot). In particular, the E-Field values calculated at the centre of the target are 23.14 dB (V/m) and 27.60 dB (V/m) for the triangular patch antenna and the MTS-enhanced triangular patch antenna, respectively.
Figure 15. Real part of the E-field calculated at 0.9 GHz and plotted on a transverse section of the simulated brain phantom including a blood-mimicking target. A higher field intensity in the target region is shown when the MTS is placed adjacent on the triangular antenna’s substrate (right plot). In particular, the E-Field values calculated at the centre of the target are 23.14 dB (V/m) and 27.60 dB (V/m) for the triangular patch antenna and the MTS-enhanced triangular patch antenna, respectively.
Diagnostics 11 00424 g015
Figure 16. Real part of the E-field calculated at 1.55 GHz and plotted on a transverse section of the simulated brain phantom including a blood-mimicking target. A higher field intensity in the target region is shown when the MTS is placed adjacent on the spear antenna’s substrate (right plot). In particular, the E-Field values calculated at the centre of the target are 20.90 dB (V/m) and 24.90 dB (V/m) for the spear patch antenna and the MTS-enhanced spear patch antenna, respectively.
Figure 16. Real part of the E-field calculated at 1.55 GHz and plotted on a transverse section of the simulated brain phantom including a blood-mimicking target. A higher field intensity in the target region is shown when the MTS is placed adjacent on the spear antenna’s substrate (right plot). In particular, the E-Field values calculated at the centre of the target are 20.90 dB (V/m) and 24.90 dB (V/m) for the spear patch antenna and the MTS-enhanced spear patch antenna, respectively.
Diagnostics 11 00424 g016
Table 1. Permittivity of glycerol solution, brain phantom, blood-mimicking target and ischaemia-mimicking target measured at 1 GHz.
Table 1. Permittivity of glycerol solution, brain phantom, blood-mimicking target and ischaemia-mimicking target measured at 1 GHz.
Type of MaterialPermittivity at 1 GHz
90% Glycerol-water Mixture15
Average Brain Phantom38
Blood Phantom65
Ischaemia Phantom22
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Razzicchia, E.; Lu, P.; Guo, W.; Karadima, O.; Sotiriou, I.; Ghavami, N.; Kallos, E.; Palikaras, G.; Kosmas, P. Metasurface-Enhanced Antennas for Microwave Brain Imaging. Diagnostics 2021, 11, 424. https://doi.org/10.3390/diagnostics11030424

AMA Style

Razzicchia E, Lu P, Guo W, Karadima O, Sotiriou I, Ghavami N, Kallos E, Palikaras G, Kosmas P. Metasurface-Enhanced Antennas for Microwave Brain Imaging. Diagnostics. 2021; 11(3):424. https://doi.org/10.3390/diagnostics11030424

Chicago/Turabian Style

Razzicchia, Eleonora, Pan Lu, Wei Guo, Olympia Karadima, Ioannis Sotiriou, Navid Ghavami, Efthymios Kallos, George Palikaras, and Panagiotis Kosmas. 2021. "Metasurface-Enhanced Antennas for Microwave Brain Imaging" Diagnostics 11, no. 3: 424. https://doi.org/10.3390/diagnostics11030424

APA Style

Razzicchia, E., Lu, P., Guo, W., Karadima, O., Sotiriou, I., Ghavami, N., Kallos, E., Palikaras, G., & Kosmas, P. (2021). Metasurface-Enhanced Antennas for Microwave Brain Imaging. Diagnostics, 11(3), 424. https://doi.org/10.3390/diagnostics11030424

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