Abstract
This article presents a holographic metasurface antenna with stochastically distributed surface impedance, which produces randomly frequency-diverse radiation patterns. Low mutual coherence electric field patterns generated by the holographic metasurface antenna can cover the K-band from 18 to 26Â GHz with 0.1Â GHz intervals. By utilizing the frequency-diverse holographic metasurface (FDHM) antenna, we build a near-field microwave computational imaging system based on reflected signals in the frequency domain. A standard horn antenna is adopted to acquire frequency domain signals radiated from the proposed FDHM antenna. A detail imaging restoration process is presented, and the desired targets are correctly reconstructed using the 81 frequency-diverse patterns through full-wave simulation studies. Compressed sensing technique and iterative shrinkage/thresholding algorithms are applied for the imaging reconstruction. The achieved compressive ratio of this computational imaging system on the physical layer is 30:1.
Introduction
Computational imaging (CI) on the physical layer based on metamaterials is an emerging research hotspot for both metamaterial applications and microwave imaging (). Compared with traditional microwave imaging techniques, for example, range migration, microwave holography, or synthetic aperture radar (; ; ), metamaterial apertures leveraging compressed sensing (CS) for CI possess simple structure, light weight, and low cost. Since major imaging restoration resources are recovery systems using CS theory, we could remarkably reduce the complexity and cost of imaging hardware systems (). The term CI is considered to describe data acquisition on the hardware layer and data post-process (imaging reconstruction) on the CS algorithm layer.
The first metamaterial aperture for CI was reported in the study by ). By utilizing frequency-diverse leaky-wave metamaterial unit cells, randomly distributed radiation patterns, satisfying compressed sampling conditions, at the frequency domain are achieved. Through measuring the reflected frequency-diverse patterns by a horn antenna, imaging of the metallic object was properly reconstructed in the observation region. Further research studies extended this work in various aspects, for instance, high Q-factor metamaterial element (), feed diversity (), dynamic metasurface antenna (; ; ), multi-look imaging (), and polarimetric imaging (). A review of metamaterial antennas for microwave CI is presented in the study by ).
Till now, existing literature mainly focused on leaky-wave metamaterials for microwave CI. We note that such a radiation mechanism only excites parts of the metamaterial elements at the specific frequency, reducing the aperture efficiency. Recent research works indicate that the two-dimensional planar metamaterial, metasurface, exhibits more flexible and powerful electromagnetic (EM) wave tuning abilities (; ; ). To exploit the whole metasurface aperture for radiation, in a recent study, a 1-bit phase-shift programmable reflective metasurface was applied to implement randomly spatial patterns (). However, the spatial feed scheme results in a high-profile structure. In this article, we propose the frequency-diverse holographic metasurface (FDHM) antenna, consisting of planar surface impedance tunable elements (), to realize measurement patterns for CI. The reasons for adopting the FDHM antenna are twofold. First, this type of metasurface antenna guides and leaks surface waves to form spatial patterns, utilizing the whole aperture for radiation. On the other hand, the frequency dispersion nature of the surface impedance element makes FDHM antenna a good candidate to fulfill the random sampling. The presented FDHM antenna shows frequency-diverse patterns within K-band from 18 to 26Â GHz. A near-field microwave CI system is constructed and validated through full-wave simulations. The data acquisition and data post-process for imaging restoration are detailed.
Materials and Methods
Frequency-Diverse Holographic Metasurface Antenna
The concept of holographic metasurface (HM) antenna was first proposed in (), which is a new kind of metasurface antenna applying holographic theory for designing beam patterns. We mentioned that the holographic theory for HM is distinct from holographic imaging. The proposed FDHM antenna is illustrated in Figure 1A, which is composed of randomly distributed square patches printed on a substrate covered by ground plane. The substrate is a Rogers 3006, with a permittivity of 6.15, loss tangent of 0.0025, and thickness of 1.27 mm. The aperture size of the FDHM antenna is 66 mm × 66 mm. The conventional design purposes of HM antennas are beam formation and polarization control (). By contrast, in this work, we use the frequency-diverse characteristics of HM antenna, which is the first attempt to the best of the authors’ knowledge.
FIGURE 1
According to the design principle of holographic metasurface () (), the surface impedance of the element should be determined first. In this design, a commercial eigen solver based on the finite element method is utilized to figure out the operating frequency under a given phase delay. The simulation model is depicted in Figure 1B with the periodic value of L = 2.2Â mm and gap g between the patch and the period. The periodic boundary condition (PBC) is set in the simulator to model infinite planar metasurfaces. Substrate parameters are described as mentioned above. Due to the symmetric property of the square patch, only phase delay Ï• on the x-axis is analyzed. We assign various phase delays in the simulator to solve the eigen frequencies. Moreover, different gaps are also studied to verify that the frequency-diverse surface impedance can be realized as random elements are applied. Once the eigen frequency fe is calculated for a given phase delay Ï• and gap, the effective surface impedance can be expressed as () follows:where Z0 is the free space impedance, c is the light speed, and L is the period of the element.
Near-Field Microwave Computational Imaging System
Based on the proposed FDHM antenna and microwave CI technique, a near-field imaging system is established, as shown in Figure 2. As it can be observed from the scheme, the FDHM antenna as a transmitter (Tx) and a standard horn antenna working at K-band as a receiver (Rx) are the critical components of this system. A basic operation principle is first described, and the microwave CI technique using the frequency-diverse method is then detailed.
FIGURE 2
Generally, microwave CI can be completed in two parts: data acquisition and data post-process. By using the proposed system, data acquisition for post-process is simply measured by a vector network analyzer (VNA). Port one connects the FDHM antenna, while Port two connects the horn antenna. The two antennas are located at the xoy-plane. Here, we assume Port 1 as Tx and Port 2 as Rx. Therefore, the measured S21 contains the required data for image restoration, which is labeled as y. At the inspected plane, the frequency-diverse patterns from Tx sense this area for each inspected voxel.
We detail the positions of the Tx, the Rx, and the inspected plane. As illustrated in Figure 2, the center of the Tx is located at (0, −35, 0) mm with the aperture size of 66 mm × 66 mm, the center of the Rx is located at (0, 10, 0) mm with the aperture size of 16 mm × 22 mm, and the inspected plane is located at (0, 0, 200) mm with the size of 400 mm × 400 mm. The height of the inspected plane is 200 mm, which lies in the Fresnel near field of the Tx/Rx antennas. We emphasize that the FDHM antenna has randomly distributed patterns in the three-dimensional space in the near field. Thus, the inspected plane in Figure 2 located inclined above the proposed antenna also possesses randomly sensing ability.
The reflected signals are captured by Rx and recorded by VNA. When the first born approximation and diffraction-limited imaging, denoted as x, are satisfied, the fields of Tx and Rx form the measurement matrix H. Then, a matrix representation is applied for data post-process, as follows () ():where M is the number of sensing measurements and N is the number of target voxels. As aforementioned, the sensing matrix H is built bywhere m = 1 ... M, and n = 1 ... N. As implied by the CS theory, for a compressible image, the reconstruction can be achieved even though M < N. For this linear inverse problem, iterative shrinkage/thresholding algorithms are effective for realizing imaging restoration () ().
Recognizing materials from background is a characterized capability of microwave imaging which is different from imaging in optics. Thereby, for the proposed imaging system, S-parameter calibration is required, that is (),where is the measured signal with objects and is the measured signal without objects. We note that the S-parameter calibration should be performed within the whole operating frequency band.
Results
Frequency-Diverse Radiation Patterns
Prior to conducting the imaging reconstruction, the randomly distributed patterns in the frequency domain of the FDHM antenna are examined. First, the surface impedance of the square element versus frequency for various gaps is calculated and plotted in Figure 3A. As it can be seen, within the operating band, for each frequency point, distinct gaps are related to different surface impedances. More importantly, for a given gap, the corresponding surface impedance varies with frequency. These characteristics provide the possibility of realizing frequency-diverse patterns through the surface impedance modulation. We mentioned that when the gap is zero, the element shows very high surface impedance because it only supports the guided mode in this case. Therefore, increasing the ratio of this element in an FDHM antenna would increase the surface wave modulation capability accordingly.
FIGURE 3
Without loss of generality, in this design, we use a pseudorandom number generator to construct the element distribution of the FDHM antenna. The FDHM antenna and horn antenna are simulated in Altair FEKO to obtain return losses (|S11|) and near-field radiation patterns. The simulated return losses are plotted in Figure 3B. From the results, the |S11| < −10 dB can be achieved within the working frequency band except for several frequencies at both sides of the band for the FDHM antenna. The results stated that the designed antenna and the standard horn antenna are capable of radiating EM waves effectively. Next, sensing matrix H is built with the near-field electric field (E-field) of the two antennas using (3). The sampling grid number is selected as 51 × 51, that is, N = 2,601. Randomly distributed sensing patterns at 18 GHz, 22 GHz, and 26 GHz are selected and plotted in Figure 3C. Besides, the corresponding sensing frequencies in the measurement matrix are also illustrated for better recognition. The simulated frequency points are M = 81, covering K-band from 18 to 26 GHz with 0.1 intervals. The achieved compressed ratio of this system is about 30:1, which can be calculated by N/M = 51 × 51/81 = 32.11:1 ≈ 30:1.
We further investigate the mutual coherence of the sensing matrix, which is defined aswhere is the sensing matrix with all columns normalized and denotes the transposed sensing matrix. For the proposed FDHM antenna, the calculated mutual coherence is μg = 0.0736, which is sufficient for compressible image restoration ().
Imaging Reconstruction
To demonstrate our near-field imaging system, we conduct theoretical imaging reconstruction to test the sensing matrix and a full-wave imaging reconstruction simulation to verify this system with a realistic metallic object. The two experiments can be applied to confirm each other. The theoretical imaging reconstruction is to build an ideal imaging for vector xideal with zero and one distribution. Element zero stands for free space, while element one stands for the perfect electric conductor (PEC) object. By using this ideal imaging, we could obtain theoretical ytheo for comparison. ytheo is obtained by multiplying the sensing matrix H by ideal imaging xtheo using (2). In the full-wave imaging reconstruction, a real PEC object located at an identical position as ideal imaging is set in the model. The simulated ysimu using (4) should agree well with ytheo. By solving (2), it is found that ytheo and ysimu are corresponding to xtheo and xsimu, respectively. The mean-squared error (MSE) is defined between the ideal image and recovered image, which is expressed aswhere denotes the l2-norm operator and xcalc should be replaced with xtheo and xsimu for calculating MSEs. We mentioned that all the following simulations are performed by using Altair FEKO.
A square object within the inspected plane is first calculated. The designed object has 5 × 5 voxels which are assigned at the center of the inspected plane, as shown in Figure 4A. The magnitude of the calculated ytheo and simulated ysimu are plotted in Figure 4B. As can be seen, the theoretical signals are in line with the simulated ones. The result confirms that the full-wave simulation model is correctly described by the math representation. The size of the square PEC object in the full-wave simulation model is 30 mm × 30 mm. Figure 4C and Figure 4D report the reconstructed imaging of xtheo and xsimu. The calculated MSEs of xtheo and xsimu to xideal are 0.0058 and 0.0075, respectively, which are comparable to the results of the existing literature ().
FIGURE 4
Next, a triangle PEC object is further demonstrated, as shown in Figure 5A. Good agreement of ytheo and ysimu is obtained and plotted in Figure 5B. The reconstructed results of xtheo and xsimu are shown in Figure 5C and Figure 5D, respectively. And the MSEs of xtheo and xsimu are 0.0073 and 0.0085, respectively. From the MSEs of square and triangle objects, it is found that when the shape of the object becomes complex, the recovered errors would increase accordingly.
FIGURE 5
Finally, two square PEC objects are designed to examine the possibility of reconstructing multiple objects, as shown in Figure 6A. Also, we report the comparison of the calculated ytheo and simulated ysimu, as shown in Figure 6B. What stands out in this case is the obvious difference between theoretical and simulated signals. The recovered two objects of xtheo and xsimu are plotted in Figure 6C and Figure 6D, respectively. The corresponding MSEs of xtheo and xsimu are 0.0238 and 0.0535, respectively, which are apparently higher than those in the previous single object cases.
FIGURE 6
Discussion
Observing the above three experiments, basically, the desired near-field CI based on the FDHM antenna is achieved. For single PEC object recovery, excellent results can be observed, and the calculated MSEs also verify the quality of the experiments. As the shape of the weak scatter object becomes complex, it is difficult to clearly reconstruct the boundary of the object. This can be indicated by the slightly increased MSEs.
Nevertheless, in the dual object case, the high MSEs and the discrepancy of theoretical and simulated signals are found, which can be attributed to the increasing diffraction of multiple objects. When the inspected objects are not weak scatters (larger object or multiple objects), artifacts would appear, thus reducing imaging quality. Increasing sensing quantity M would be a possible solution to this issue. In addition, increasing the distance between objects and Tx/Rx antennas might make the scatter field weaker than the near distance. However, as the inspected plane moves forward, the diversity of the sensing patterns also declines. Therefore, a trade-off between compressible sampling and weak scatter condition should be considered in the CS-based microwave CI system.
The main contribution of this article is the demonstration of the feasibility of achieving CI by the FDHM antenna. Even though the demonstration is conducted by numeric analyses, various previous articles indicate that the commercial full-wave simulation software could model practical antennas with a quite high accuracy. In future work, experimental validation will be performed.
Conclusion
This article has successfully demonstrated the feasibility and effectiveness of the FDHM antenna for microwave CI. An FDHM antenna with randomly distributed surface impedance tunable elements is proposed and validated by full-wave simulation. Based on this antenna, a near-field CI system using frequency domain data is presented. The detailed operation principle of this system and the math representation are elaborated. Further, we show the theoretical imaging reconstruction and full-wave simulation imaging reconstruction by using this system. Good agreement is obtained between the two methods, validating our FDHM antenna and the CS theory for the CI system. A single PEC object with a square or triangle shape and two square PEC objects are properly recovered. We could conclude that microwave CI systems centered on the CS technique and the holographic metasurface can challenge conventional microwave imaging systems and produce cost-effective and high-efficient devices.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author contributions
JH and LL proposed the idea. JH, TS, XM, and QF conceived and performed the metasurface antenna design and built the imaging system. HL, YZ, and GS guided the data post-process. All authors analyzed and discussed the results.
Funding
This work was supported in part by the National Natural Science Foundation of China under Grant No. 62001342, in part by the Joint Foundation of Key Laboratory of Shanghai Jiao Tong University-Xidian University, Ministry of Education under Grant No. LHJJ/2020-02, and in part by the Fundamental Research Funds for the Central Universities under Grant No. XJS200207.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AminehR. K.RavanM.KhalatpourA.NikolovaN. K. (2011). Three-Dimensional Near-Field Microwave Holography Using Reflected and Transmitted Signals. IEEE Trans. Antennas Propagat.59 (12), 4777–4789. 10.1109/TAP.2011.2165496
2
Bioucas-DiasJ. M.FigueiredoM. A. T. (2007). A New TwIST: Two-step Iterative Shrinkage/Thresholding Algorithms for Image Restoration. IEEE Trans. Image Process.16 (12), 2992–3004. 10.1109/TIP.2007.909319
3
CossairtO.GuptaM.NayarS. K. (2013). When Does Computational Imaging Improve Performance. IEEE Trans. Image Process.22 (2), 447–458. 10.1109/TIP.2012.2216538
4
DonohoD. L. (2006). Compressed Sensing. IEEE Trans. Inform. Theor.52, 1289–1306. 10.1109/TIT.2006.871582
5
DuarteM. F.DavenportM. A.TakharD.LaskaJ. N.SunT.KellyK. F.et al (2008). Single-pixel Imaging via Compressive Sampling. IEEE Signal. Process. Mag.25 (2), 83–91. 10.1109/MSP.2007.914730
6
FongB. H.ColburnJ. S.OttuschJ. J.VisherJ. L.SievenpiperD. F. (2010). Scalar and Tensor Holographic Artificial Impedance Surfaces. IEEE Trans. Antennas Propagat.58 (10), 3212–3221. 10.1109/TAP.2010.2055812
7
FromentezeT.YurdusevenO.BoyarskyM.GollubJ.MarksD. L.SmithD. R. (2017). Computational Polarimetric Microwave Imaging. Opt. Express25 (22), 27488. 10.1364/OE.25.027488
8
FromentezeT.YurdusevenO.ImaniM. F.GollubJ.DecrozeC.CarsenatD.et al (2015). Computational Imaging Using a Mode-Mixing Cavity at Microwave Frequencies. Appl. Phys. Lett.106 (19), 194104. 10.1063/1.4921081
9
HanJ.LiL.MaX.FengQ.ZhaoY.LiaoG. (2021). A Holographic Metasurface Based on Orthogonally Discrete Unit-Cell for Flexible Beam Formation and Polarization Control. Antennas Wirel. Propag. Lett.20 (10), 1893–1897. 10.1109/LAWP.2021.3099574
10
HanJ.LiL.MaX.GaoX.MuY.LiaoG.et al (2021). Adaptively Smart Wireless Power Transfer Using 2-Bit Programmable Metasurface. IEEE Trans. Ind. Electron., 1. 10.1109/TIE.2021.3105988
11
HanJ.LiL.TianS.LiuG.LiuH.ShiY. (2020). Millimeter-Wave Imaging Using 1-Bit Programmable Metasurface: Simulation Model, Design, and Experiment. IEEE J. Emerg. Sel. Top. Circuits Syst.10 (1), 52–61. 10.1109/JETCAS.2020.2973466
12
HuntJ.DriscollT.MrozackA.LipworthG.ReynoldsM.BradyD.et al (2013). Metamaterial Apertures for Computational Imaging. Science339 (6117), 310–313. 10.1126/science.1230054
13
HuntJ.GollubJ.DriscollT.LipworthG.MrozackA.ReynoldsM. S.et al (2014). Metamaterial Microwave Holographic Imaging System. J. Opt. Soc. Am. A.31 (10), 2109. 10.1364/JOSAA.31.002109
14
ImaniM. F.GollubJ. N.YurdusevenO.DieboldA. V.BoyarskyM.FromentezeT.et al (2020). Review of Metasurface Antennas for Computational Microwave Imaging. IEEE Trans. Antennas Propagat.68 (3), 1860–1875. 10.1109/TAP.2020.2968795
15
LiL.ZhangP.ChengF.ChangM.CuiT. J. (2021). An Optically Transparent Near-Field Focusing Metasurface. IEEE Trans. Microwave Theor. Tech.69 (4). 10.1109/tmtt.2021.3061475
16
LiL.ZhangX.SongC.ZhangW.JiaT.HuangY. (2021). Compact Dual-Band, Wide-Angle, Polarization- Angle -Independent Rectifying Metasurface for Ambient Energy Harvesting and Wireless Power Transfer. IEEE Trans. Microwave Theor. Techn.69 (3), 1518–1528. 10.1109/TMTT.2020.3040962
17
LipworthG.HuntJohn.MrozackAlex.BradyDavid.SmithDavid. R. (2013). Metamaterial Apertures for Coherent Computational Imaging on the Physical Layer. J. Opt. Soc. Amer. A, Opt. Image Sci.30 (8). 10.1109/comcas.2013.6685286
18
Lopez-SahcnezJ. M.Fortuny-GuaschJ. (2000). 3-D Radar Imaging Using Range Migration Techniques. IEEE Trans. Antennas Propagat.48 (5), 728–737. 10.1109/8.855491
19
MoreiraA.Prats-IraolaP.YounisM.KriegerG.HajnsekI.PapathanassiouK. P. (2013). A Tutorial on Synthetic Aperture Radar. IEEE Geosci. Remote Sens. Mag.1 (1), 6–43. 10.1109/MGRS.2013.2248301
20
PandiS.BalanisC. A.BirtcherC. R. (2015). Design of Scalar Impedance Holographic Metasurfaces for Antenna Beam Formation with Desired Polarization. IEEE Trans. Antennas Propagat.63 (7), 3016–3024. 10.1109/TAP.2015.2426832
21
PatelA. M.GrbicA. (2011). A Printed Leaky-Wave Antenna Based on a Sinusoidally-Modulated Reactance Surface. IEEE Trans. Antennas Propagat.59 (6), 2087–2096. 10.1109/TAP.2011.2143668
22
Pulido-ManceraL.FromentezeT.SleasmanT.BoyarskyM.ImaniM. F.ReynoldsM.et al (2016). Application of Range Migration Algorithms to Imaging with a Dynamic Metasurface Antenna. J. Opt. Soc. Am. B33 (10), 2082. 10.1364/JOSAB.33.002082
23
SheenD. M.McMakinD. L.HallT. E. (2001). Three-dimensional Millimeter-Wave Imaging for Concealed Weapon Detection. IEEE Trans. Microwave Theor. Techn.49 (9), 1581–1592. 10.1109/22.942570
24
SleasmanT.BoyarskyM.ImaniM. F.GollubJ. N.SmithD. R. (2016). Design Considerations for a Dynamic Metamaterial Aperture for Computational Imaging at Microwave Frequencies. J. Opt. Soc. Am. B33 (6), 1098. 10.1364/JOSAB.33.001098
25
SleasmanT.F. ImaniM.GollubJ. N.SmithD. R. (2015). Dynamic Metamaterial Aperture for Microwave Imaging. Appl. Phys. Lett.107 (20), 204104. 10.1063/1.4935941
26
YurdusevenO.MarksD. L.FromentezeT.GollubJ. N.SmithD. R. (2017). Millimeter-wave Spotlight Imager Using Dynamic Holographic Metasurface Antennas. Opt. Express25 (15), 18230. 10.1364/oe.25.018230
Summary
Keywords
holographic metasurface, computational imaging, compressed sensing, near-field imaging, frequency-diverse
Citation
Han J, Li L, Tian S, Ma X, Feng Q, Liu H, Zhao Y and Liao G (2021) Frequency-Diverse Holographic Metasurface Antenna for Near-Field Microwave Computational Imaging. Front. Mater. 8:766889. doi: 10.3389/fmats.2021.766889
Received
30 August 2021
Accepted
13 September 2021
Published
19 October 2021
Volume
8 - 2021
Edited by
Wei-Xiang Jiang, Southeast University, China
Updates
Copyright
© 2021 Han, Li, Tian, Ma, Feng, Liu, Zhao and Liao.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Long Li, lilong@mail.xidian.edu.cn
This article was submitted to Metamaterials, a section of the journal Frontiers in Materials
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.