Improved Metallic Enclosure Electromagnetic Imaging Using Ferrite Loaded Antennas

: Three-dimensional electromagnetic imaging can be used to monitor grain within metallic grain bins. Data acquisition requires multiple antennas surrounding the imaging space, which are used to transmit and receive the electromagnetic energy inside the bin. Due to their presence inside a metallic enclosure and due to very large mechanical forces these antennas are required to be low proﬁle. In addition, since they are part of the imaging domain, they should be simple to model in the imaging software (i.e., using a point source). Existing half-loop magnetic ﬁeld antennas meet these design criteria, but can be improved, particularly with better radiation efﬁciency. Herein, we present an enhanced antenna design: a ferrite-loaded shielded half-loop antenna designed to measure only the tangential component of the magnetic ﬁeld against the metal enclosure wall, while rejecting the normal component of the electric ﬁeld. Experimental results in two bins show that the enhanced design improves the signal level over existing probes by 6–18 dB inside a small-scale enclosure and around 20 dB inside a larger 28 m 3 (800 bushel) bin. Full 3D imaging results of a high-moisture target within a low-moisture grain background inside the test enclosure show that the enhanced antennas improve the quality of the reconstructed results in the smaller bin, particularly where the antenna performance improvements are prominent.


Introduction
Electromagnetic Imaging (EMI) uses electromagnetic waves to qualitatively and/or quantitatively reconstruct the properties of a target or region of interest (ROI). Fields are transmitted into the region and then measured at multiple locations around the ROI with a series of antennas. These measurements are then used as input to an inversion algorithm to reconstruct the electrical properties (i.e., the complex permittivity) of the ROI [1]. EMI has been used in different applications including medicine [2][3][4], geosciences [5], food science [6], stroke imaging [7][8][9][10], breast cancer detection [11][12][13][14][15][16][17][18] and grain bins [19][20][21][22]. Within grain monitoring, the complex permittivity of grain has long been used to estimate the moisture content of the grain [23], and the moisture content of the stored grain is critical to its safe storage.
Antennas are a critical component of EMI systems. They both radiate the EM energy into the grain bin and receive the scattered fields from the ROI that ultimately provide the information about the ROI. Within grain bin EMI, design considerations of these antennas include that (1) they must have a field pattern that can be readily modeled in the inversion algorithm, (2) they must be sensitive enough to detect the fields close to the bin wall, (3) must pick up a well-modeled field polarization (EMI algorithms, e.g., [11,24,25] assume input data of an electric/magnetic field at a point with a known polarization, and antennas that are sensitive to combinations of different fields cannot readily be calibrated in experimental systems) and (4) they must be able to physically survive the rigors of moving grain in a bin. Previous work on designing antennas to satisfy these constraints inside of metallic enclosures has resulted in a half-loop design [20] capable of measuring the H-field perpendicular to the loop's main area, while these antennas work well in grain bins, they suffer from the fact that sub-wavelength loop antennas have a very poor radiation efficiency (i.e., they are very inefficient at converting input power into radiated waves) and most of the signal's incident upon them are reflected. It is thus desirable to have an antenna with better radiation efficiency that still satisfies the above design constraints. Grain bins can have diameters greater than 40 m, and each antenna must be connected to a Vector Network Analyzer (VNA) via a coaxial cable. In such systems, the raw S-parameter signal levels are often in the range of −80 to −100 dB (due to the lossy nature of the grain, the inefficiency of the antennas, and the losses in the long cables). Thus, a more efficient antenna will increase the signal-to-noise ratio in these large bins (as the noise floor of the VNA is often close to −110 dB).
The contribution of this work is to design an improved antenna for imaging in grain bins that satisfies the imaging-system design constraints but can offer an improved signalto-noise level-thus resulting in better final images of spoiled grain. Designing such an antenna in a metallic chamber a few wavelengths across means that many far-field antenna design concepts do not apply to this problem. That is: one cannot take an antenna that works well in the far field in free space and use it in a grain bin (and expect it to work with existing inversion algorithms). Thus, our main contribution is to offer an efficient antenna that works well in metallic bins given the constraints of operating in a grain bin.
Our improved antenna is achieved through designing well-matched antennas that are easy to model using simple point sources and are sensitive to a single polarization of the electromagnetic field at the location of the antenna. In addition, we wish to maintain the low-profile nature of the existing antennas. As a solution to these problems, we present a miniaturized resonant antenna that is used for EMI. This is achieved through the strategic placement of ferrite media around the antenna structure. Like other scaled bin test environments [26], we use a small test bin to experimentally validate our design.
The work herein is an expansion on previously presented preliminary results [27], where we presented initial 2-antenna results that showed the increase in signal strength.
In this work, we show that the full 24 antenna array works with the full 3D imaging procedure in a scaled grain bin (which is not a given strictly from an increase in signal strength). We also present preliminary results with another set of antennas in a much larger (industrial scale) grain bin, thus showing that our proposed antenna can be increased in size to work in these larger bins.

EMI as a Grain Permittivity Measurement Instrument
The use of electromagnetic instrumentation and measurement of grain parameters, particularly measurement of permittivity, has a long history (see, e.g., [28][29][30][31][32][33][34][35]). The moisture content of grain has long been measured with electromagnetic devices that measure the permittivity of small samples of the grain [23]. The EMI system we use in this paper [19][20][21][22] can be considered an extension of such permittivity measurement instruments that operate in two modes, both of which have some improvements on previously existing grain permittivity measurement instrumentation. The first mode [36] treats the entire grain bin as a single permittivity measurement instrument and outputs the average permittivity of the bin contents, and the second mode is a system for providing a 3D map of the permittivity of each small region in the grain (in our case, each small volume is represented by a tetrahedral). In this sense, by providing an improved EMI antenna, we are developing an improved instrument for measuring the permittivity of grain in a wide variety of grain storage bins, essentially transforming any metallic grain bin into a bulk and 3D permittivity measurement instrument.

Paper Organization
The paper is organized as follows: the imaging algorithm is briefly described in Section 2. In addition, cost functions for optimization algorithms are presented. Section 3 covers antenna design, simulation, fabrication, and measurement results. The important step of data collection is presented in Section 4. Section 5 represents the inversion process steps. In Section 6, parametric inversion results have been presented. In Section 7, we present images produced by the designed antennas. In Section 8, we have presented our designed large FLA and implemented it inside a large grain bin for EMI. Finally Section 9 concludes.

Problem Definition and Imaging Algorithm
A basic schematic of the EMI problem is presented in Figure 1. In this figure, the yellow region represents the background media with relative permittivity b (r) (note that this background permittivity can be heterogeneous). In addition, inside the enclosure, there is a region of unknown permittivity r (r). In our application, we are concerned with changes in permittivity in the media. Equation (1) represents complex permittivity in media, where 0 is the permittivty of free space, r (r) is the complex relative permittivity with r (r) and r (r)) as it relative real and imaginary part of relative permittivity, respectively. = 0 r (r) = 0 ( r (r) + j r (r)) (1) The antennas under consideration in this paper are sensitive to the magnetic field, H. We break up the total field into incident and scattered fields associated with antenna j as H tot The incident fields are associated with the set of background media b (r), and we are free to define the background media to best solve the inverse problem. The scattered fields are associated with the contrast between the true media and our selected background; this contrast is defined as: For example, in the case of the grain bin used in this work, we choose to define the background media b (r) to be a mixture of air and a relatively homogeneous mass of grain with some height and cone angle. The contrast (and thus scattered fields) is thus associated with deviations from the homogeneous grain (e.g., higher moisture regions cause the scattered fields). In general, non-linear inverse problems are more successful when the contrast is as small as possible, so choosing, e.g., the background to be entirely free-space leads to much worse inversion results.
The core of the inversion solver in this paper uses the Finite Element Method (FEM) to simulate the electric field equation [24]: is the ( inhomogeneous) background wave-number and w(r) is called contrast source and defined as follows: While the FEM code directly simulates only the E-field, the magnetic fields, H, are generated via a transform on the FEM mesh [20] (and all inversion happens with the magnetic fields).
In the FEM code, the source term at each antenna position is a magnetic dipole (i.e., a magnetic current) wherev is the unit-vector in the direction of the antenna orientation (i.e., perpendicular to the loop). Finally, to compute the magnetic fields at the antenna positions (which is required at the receiver positions) r j , the FEM code applies Faraday's law to the FEM-electric field basis functions [20].

Inversion Algorithm
The process of solving for the total permittivity from the measured scattered fields is known as the inverse scattering problem. This process assumes knowledge of both the incident field in the background permittivity, as well as measurements of the total fields at the receiver points. We use the Contrast Source Inversion algorithm with the forward solver implemented in the Finite Element Method (FEM) algorithm, known as FEM-CSI [24]. This algorithm optimizes the following cost functional: where L represents the inverse of finite element discrete matrix which depends on the boundary condition of the enclosure, and M stands for interpolation matrix operator to convert fields values along edges of meshes to vector fields at locations of transceivers or points of a domain. Moreover, represents the element-wise multiplication of matrices of the same size. The input data to the optimization Equation (6) are the scattered magnetic fields measured by the antennas (see below for antenna design for discussion as to why the antennas are magnetic field detectors). The CSI algorithm optimizes this functional through by applying the conjugate gradient minimization algorithm to the two unknowns, χ and w [24].

Antenna Design
Against a metallic boundary, only tangential magnetic fields and normal electric fields exist. Due to the modeling complications and limitations of the inversion algorithms, we design an antenna that measures a single field component [19,20] in this case, the tangential magnetic fields. A half loop magnetic field probe, depicted in Figure 2, can be used to measure a single polarization of magnetic fields. This was previously achieved using a co-axial cable with a small cut on the outer conductor at the center of the half-loop (see Figure 2). The inner conductor serves as the loop that is sensitive to the magnetic field, while the outer conductor shields the arms of the loop parallel to the normal electric field [20]. The such antenna measures the magnetic field passing perpendicularly through its surface [20]. As described previously, this type of antenna suffers from poor radiation efficiency, resulting in a high magnitude of the reflection coefficient. To improve the antenna performance, we propose the use of ferrite material in the antenna design. We refer to the ferrite type antennas as Ferrite-Loaded Antennas (FLA), and the original antennas [20] as Non-Ferrite Loaded Antennas (NFLA).
In this work, we design the antennas around a small-scale grain bin, previously used in [20,21]. Due to the size of this small-scale grain bin, a frequency range of 200 MHz to 600 MHz is selected for design and imaging purposes. This range is selected such that the number of wavelengths across the bin radius is approximately 1-3 wavelengths, which is a result of a trade-off between computation time and resolution, while we focus on this scaled bin, the principals involved will also scale to larger antennas (required for larger bins).

Ferrite Media
The general categorization of magnetic media can be found in [37]. At microwave frequencies, ferrites are dispersive solid ceramics with low conductivity and relative permittivity between 10 and 15 [37,38].
Unlike non-magnetic media, both permeability (µ), as well as permittivity ( ), must be considered when using ferrites in antenna design. In general, because of their dispersion, for a wide-band antenna, the dispersion needs to be modeled. A 2nd order Debye model [38] is used for both permittivity and permeability (µ = µ 0 µ r (r) = µ 0 (µ r (r) + jµ r (r))) where µ 0 is the permeability of free space, i.e., The parameters for the second order Debye model for permittivity and permeability of the ferrite material in the antenna design are presented in Table 1.
Using the parameters in Table 1 results in permittivity and permeability over the desired maximum design frequency range are shown in Figure 3.

Ferrite Loaded Antenna Simulation
Through a design process that carefully considered cost, mechanical stability (ferrite is very brittle, and grain bins are harsh environments), and ferrite availability, we attempted multiple configurations of ferrite media in the half-loop antenna. These possible configurations were simulated, but we finally settled on the design shown in Figure 4 (where the ferrite material is shown in aqua color). As this antenna must be mounted in a bin containing grain, and the movement of grain can produce large forces (and rip antennas off the bin wall), we have mounted ferrite plates directly to the bin walls to prevent extra protrusions into the bin (the protrusions increase the forces on the antenna). With this design, the antenna remains electrically small [39], with a resonance occurring at approximately 260 MHz.
The proposed design was simulated with a 50 Ω excitation port using ANSYS Electromagnetic Suite 2021 (HFSS). The plot of the magnitude of the reflection coefficient of the antenna (i.e., S 11 ) from the simulation is shown in Figure 5. It is worth noting that the simulation model included a metallic backplane to simulate the enclosure's wall, instead of the full resonant metallic enclosure (to reduce the computation cost of the simulations), and the antenna was surrounded by a homogeneous region of the air. For a Voltage Standing Wave Ratio (VSWR) with criteria VSWR < 2:1, the design has a bandwidth of 60 MHz, centered at 260 MHz. At this frequency band, the inclusion of the ferrite has created a well-matched antenna.

Measurement
The antenna was fabricated and installed inside the bin ( Figure 5). Note that there is a white plastic cover on the coaxial cable that protects the cable from mechanical damage but does not affect the electromagnetic performance. After installing the antenna in the metallic enclosure, the magnitude of the reflection coefficient of the antenna was measured using a VNA as shown in Figure 5. For comparison purposes, the magnitude of the reflection coefficient for the NFLA is also graphed on the same plot.
The differences in S 11 are obvious between the fabricated FLA and NFLA antennas: most of the incoming energy is reflected from the NFLA for all the design frequencies (approx. (200-600) MHz). For the FLA, there are small differences between simulated and measured values. The differences can be attributed to the fact that the fabricated antenna was placed inside the open bin (where the simulation was in free space), the physical antenna has imperfections in its shape, as well as the fact that the simulation uses a Perfect Electric Conductor (PEC) while the actual antenna is made of copper. Given these differences, there is a good match between simulation and measurement.

Experimental Comparison of NFLA and FLA Antennas
The main goal of the proposed FLA antennas is to increase the signal strength-in practice the S 21 -above those of the existing NFLA antennas. To show the higher S 21 signal levels with the proposed design, we measured S NFLA 21 for two NFLA (original) antennas with the test bin partially full. We next added the ferrite media to the proposed antennas (making them FLA) and re-took the measurement, resulting in S FLA 21 . The grain and rest of the measurement apparatus remained untouched between the measurements. Figure 6 shows the |S 21 | comparisons of FLA with NFLA over the designed frequency ranges. This plot shows that the new antenna receives up to 19 dB more S 21 signal than the NFLAs with the strongest signal improvement occurring at the resonance frequency of the antenna.

Imaging Tests
While the higher-signal (and thus lower-noise) advantages of the FLA antennas are clear from Figure 6, there are other design parameters that must be satisfied by these antennas. Ultimately, they must be able to provide useful images inside of the resonant enclosure. For example, it is possible to design a well-matched antenna that does not meet other design criteria, such as producing an incident field well modeled by a magnetic dipole, as well as primarily receiving the H-field polarization perpendicular to the antenna loop. To show that these design parameters are met, we have performed a series of experiments in the small-scale grain bin with both FLA and the original NFLA antennas. These experiments are described below.

Small Scaled Grain Bin Measurement Setup
The scaled grain bin has been depicted in Figure 7. The bin has a depth of approximately 70 cm and a radius of 28.5 cm, with a total volume of approximately 190 L. Previous results with this scaled grain bin have been presented in [19,21]. We installed either 24 FLA or NFLA antennas, as shown in Figure 7. Each antenna was connected to the GrainViz imaging system that contains a Vector Network Analyzer (VNA), amplifiers, control computers, and a 24 × 2 port switch. Antennas were connected to the GrainViz system through 1.5 m long coaxial cables. Antennas are staggered (with a 15 • shift per layer) in four layers. This arrangement of antennas was chosen to maximize the spread of antennas on the inner wall of the enclosure and ensure proper field sampling. The VNA was used to transmit and receive data through the antennas via the switching network, resulting in 24 × 23 (552) data points per frequency. The measured data was used as input to imaging algorithms.

Grain Description
We used two treatments of hard red spring wheat as our test grain. The bulk of the wheat was very dry (7.5% moisture content, wet basis). This grain was harvested in approximately 2018 and had been stored in a stable indoor environment, and thus was assumed to be homogeneous. From the available literature, the approximate relative permittivity of the hard red spring wheat at ≈260 MHz is r = 3 − j0.3 [40]. It is worth noting that this relative permittivity is approximate as the density of the grain can have significant effects on the permittivity (we did not directly measure the permittivity of the grain) (Temperature also affects the permittivity of grain, but [40] was at 25 • C and our grain was approximately at room temperature (≈21 • C), which has relatively small effects on permittivity).
We next created a wet-grain target by placing 1410 g of dry wheat mixed with 210 g of tap water inside two plastic bags. A photo of the wet-grain target in position is shown in Figure 8. This target was allowed to sit at standard refrigeration temperatures for at least one week to allow the water to be absorbed by the grain. Assuming all of the water was absorbed over this time, the wet-grain target had a moisture content of 19.5% and would have an approximate permittivity of r = 4.5 − j0.6 at ≈260 MHz [40]. At room temperatures, such high-moisture wheat could only be stored for a couple of days without spoilage [41], while still malleable (and thus not purely in the shape of a rectangular cuboid) the approximate dimensions of the wet-grain target were 7.5 cm × 11 cm × 21 cm.

Data Collection
We began the data collection process by filling the grain bin with the dry hard red spring wheat. To achieve a flat grain surface, the initial fill of grain was smoothed by hand. The grain leveling was done without the aid of leveling tools, so the accuracy of the flattened surface was unknown. We had enough grain to fill in slightly less than half full. This was followed by placing a plastic bucket with a diameter of 0.1 m and a height of 0.25 m inside the bin. The bucket is needed to accurately place the wet-grain target as it is very difficult to dig inside of the free-moving grain. The thin bucket is made of plastic and has a permittivity like the bulk grain, so it is ignored in the inversion process. Figure 8 shows the bucket before it is covered with wheat (but with the wet target), near antenna #15. We then re-filled the bucket with dry wheat back to the original level. Once the bucket was buried in dry grain, we then took an incident S-parameter measurement herein referred to as S inc . This data set is then subsequently used to estimate the background permittivity, b , using a parametric imaging algorithm [36].
Next, we removed the dry grain inside of the bucket and placed the wet-grain target inside (see Figure 8). After re-filling, the bucket back to level with the dry grain (not pictured), We then took a total S-parameter measurement, which has been labeled S tot . Finally, note that we repeated each of the steps above for both the proposed FLA and NFLA.

Inversion Process
Standard inversion work-flow for imaging inside of grain bins has been previously outlined in [20,22] and is as follows: 1.
Obtain b (r): Using raw S inc (the measurements without the wet grain target), estimate the background permittivity, b (r) through phaseless parametric inversion [36]. This algorithm outputs the grain height and cone angle, as well as an average permittivity. The background estimation algorithm uses the magnitude data and a gradient-free optimization to find the best average permittivity, grain height, and cone-angle.

2.
Generate Incident Fields: Using the inhomogeneous air/grain background permittivity, b (r), as the input to a FEM forward solver, the incident magnetic field for the numerical model, called H bkg , is obtained.

3.
Calibrate Scattered Fields: We generate calibrated scattered fields from the incident magnetic fields, and the total/incident S-parameter measurements via the equation: where t x , r x denote the transmitter and the receiver index.

4.
Full Inversion: Finally, we pass the calibrated scattered fields from Equation (8), along with the background permittivity b (r), to the FEM-CSI algorithm described in Section 2.
For this work, we choose to run the FEM-CSI independently at 10 frequencies in two ranges: 240, 244, 250, 260, 270, 280 and 540, 550, 560, 570 (all in MHz). Frequencies as the ≈250 MHz group was close to the FLA resonance, while the higher frequencies are chosen as this is where the NFLAs have been previously used successfully [20,21].

Parametric Inversion Results
We first present phaseless parametric inversion for both the FLA and NFLA. The objective function for the parametric inversion function, Equation (6), is L 2 norm based on the magnitudes of the data. As a way of determining how successful the parametric inversion was at computing the correct result, we use the Pearson correlation coefficient between the magnitude of the raw data, |S inc | and the fields generated from the computational model |H bkg |.

Correlation Coefficient
The correlation coefficient is insensitive to scaling differences between the fields and the measured S-parameters and can be used as an intuitive metric to compare how well the computational model fits the data: 100% correlation is a perfect match between the real world and the model. Any deviation from 100% means the computational model is not a perfect fit. A low correlation coefficient for a given antenna would mean that it is not a good fit for the computational model (e.g., if the antenna did not act like a magnetic dipole, the correlation coefficient would be low). For S inc data sets taken with both FLAs and NFLAs, Figure 9a displays the correlation coefficients at each frequency between |S inc | and |H bkg |. In every case, the correlation coefficient for FLAs is higher than for the NFLAs. The FLA correlations are also higher near the lower, resonant frequency band.

Height, Angle, and Background Permittivity Detection
The height, cone angle, and real permittivity estimated by the parametric inversion for both FLA and NFLA S inc data sets are shown in Figure 9b Errors appear larger in the angle of repose of the grain, shown in Figure 9c. We have assumed the true repose angle is 0 • , but the repeated positive repose angle for both data sets likely means that our manual flattening process was leaving a small peak of grain. These errors are larger than the discretization errors inherent in the super-position of the cone-angle onto our unstructured mesh (see Figure 10 for a visualization of how the cone angle is superimposed on the mesh). For example, considering the FLA data set, the average angle is 10 • (rounded to the nearest degree). This represents a peak height of 5 cm in the center of the bin (or about 1.5 tetrahedral edge lengths). Thus, unlike the height error, the error in the angle of repose is slightly larger than the discretization error. The values for the real permittivity at the chosen frequencies have been shown in Figure 9d. The real permittivity is ultimately used to detect the moisture of grain [23], and the imaginary permittivity is much more variable on grain variety/density, etc. [23]. Therefore, we present only the real permittivity results (although the imaginary part of the permittivity is reconstructed in the inversion algorithm), while many factors (most importantly density) can affect the permittivity of the wheat, the obtained values are comparable to the literature values. Using the permittivity values in [23] for wheat at 7.5% at 25 • C, we have added Nelson's model points to Figure 9d. It is worth noting that [23] does not have data above 300 MHz. Comparing Re( b ) in plots of Figure 9d, the results from the FLAs follow the expected value in our lower frequency band, while at the upper-frequency band they are very similar. Overall, the FLAs are likely to provide a more accurate estimate of real permittivity than the NFLAs.

Full 3D Imaging Results
We followed the imaging procedure outlined in Section 6 and collected the total-field data S tot with the wet-grain target buried in the grain mass for both FLA and NFLA. Then the data as per Equation (8), is calibrated and the FEM-CSI inversion algorithm is run.
For each data set, we present full 3D results for two different frequencies: 260 MHz and 570 MHz (Note that other frequencies in the two ranges provided similar results, and the results only at these two frequencies are presented in order to save space). As mentioned earlier, only real permittivity results are used to eventually measure the moisture content of the grain. All inversion results are presented on the same color scale, and this color scale is shown in Figure 10. This common color scale allows the direct comparison of imaging results between the FLAs and NFLAs. For the 3D threshold images, we use a threshold of Re( r ) > 3.6, which is 90% of the maximum color scale of 4.
We first present the 3D Re( r (r)) reconstructions of the grain with the wet target for the FLAs in Figure 10, along with some metrics in Table 2. Reconstructions of the same target with the NFLAs have been shown in Figure 11 and Table 3. Switching between FLA and NFLA requires fully emptying the grain bin and removing the target, so the target is not in the same position between tests (Ideally, we would have placed the target in as close a position for both the FLA and NFLA tests, but experimental error led us to place the target in different positions each time. However, we do note that the target is at the same depth in each test and that its x/y position is rotated approximately 180 • . Since the bin is symmetric to 180-degree rotations, we do not expect the differing target position to have a major effect on image quality).
For the FLA results in Figure 10, the target center was measured about (x, y, z) = (0, −14, −57) cm. Based on Figure 10, we can estimate the size and position of the target; these values have been presented in Table 2. In this table, dimensions of the target at each frequency are d x , d y , and d z where |e x |, |e y |, and |e z | are corresponding relative errors to true dimensions in each direction and have been represented in percentage format. Moreover, p x , p y , and p z are the coordinates for the center of the target. We note that all errors must be taken as approximate as the actual target was not a rectangular cuboid. To calculate the dimensions d i and center, p i , of the target we visually marked the red zones in each direction (x, y, z). We then calculated the distance between each mark and averaged them over each plane. For instance, from Figure 10, we obtain 12 cm and 13 cm for length of target along x-direction from xy and xz side views at 260 MHz, respectively. We then present the obtained dimension as d x = 12.5 cm. Absolute percent errors, %e i in each dimension are based on the approximate target dimensions of 7.5 × 11 × 21 cm 3 . We again emphasize that at 250 MHz, 1 cm is about ( 1 70 ) th of a wavelength in grain.  From the NFLA results shown in Figure 11 and Table 3, it is clear that the NFLAs do not perform as well as the FLAs at the lower frequency (260 MHz) (i.e., The peak Re( r ) is not as high, and the reconstructed target is much larger than the true target, particularly in the y-direction). At the higher frequencies, the NFLA images appear to have a more well-delineated target than the FLA images, especially in the z direction. That said, the NFLA-570 MHz results give an incorrectly low estimate of the z-height of the target.

Target Tracking
As a final experiment in the small bin, we performed a series of target tracking experiments with the FLAs where we placed the targets in four different positions. For each position, we performed the full imaging experiment process (i.e., incident field measurement, bury the target, total field measurement, run the parametric optimization on the incident field, calibrate, then run FEM-CSI). The goal of these experiments is to show that the target is properly located in the inversion images and to determine if the position of the target caused large differences in the inversion images (e.g., perhaps the target is reconstructed more accurately depending on its position). We moved the target to positions close to antennas 7, 15, 3, and 19. Figure 12 presents the threshold images for this experiment. Note that, to save space, we present two different positions of the target in the same image (but would like to be clear that each image represents two different experiments of imaging a single target). The target tracking experiments support the claim that inversion results are representative of the true target position, and the reconstructed target size/position is not highly dependent on the target position.

A Preliminary Experiment in 28 m 3 (800 bu) Grain Bins
To verify that newly designed FLAs work inside larger grain bins (i.e., not a small-sized experiment), we present preliminary results from larger-sized FLAs in a 28 m 3 (800 bushel) bin filled with canola at the University of Manitoba's Grain Storage Research facility. Unlike the scaled test, we are unable to directly compare FLA and NFLA antennas of the same size design (i.e., with and without the ferrite) as the physical size of the bins precludes changing the antenna structure in situ. However, we do have a comparison bin with FLA antennas of a much larger size as we have an installation of NFLAs in the bin next to the FLA bin. A photo of the two bins is shown in Figure 13a. The two bins have a flat bottom with a diameter of 2.73 m and a height of 4.58 m. We again used 24 antennas in each bin. In bins this size, we have typically used inversion frequencies of (120-180) MHz. The design and construction of our larger FLA has been shown in Figures 13b and 14a, and is made of semi-rigid coaxial cable RF195. Unlike the small FLA in previous sections, this design includes a metallic plate as its ground plane. This plate is used to attach the antenna to the bin's wall using magnets (see Figure 13c). Two ferrite plates close to each port of the antenna are in the antenna structure. If we cover the entire metallic plate with ferrite material, the antenna's S 11 will improve by less than 2 dB, therefore to reduce the net cost of the antenna and also consider mechanical issues in large bins, we decided to reduce the number of ferrite plates in the antenna structure. To protect the larger FLAs from the moving grain in the bin, the antenna was enclosed in a hard plastic box (Figure 13c) with dimensions 20 × 12 × 7.5 cm 3 .
(a) (b) Figure 14. (a) Dimensions of the large FLA in mm are: w g = 56; l g = 160; h g = 10; w f = 56; l f = 100; l l = 140; w l = 30; g = 1. Note that aqua colors, two regions with dimensions l f × w f (= w g ), are Ferrite Plates (FP1 and FP2). They both have a thickness of 2.2 mm which has not been shown in the figure. (b) Measured magnitude of the reflection of large FLA and large NFLA.
Critically important to our results in comparing the two bins, the nominal area of the loop for the FLA antenna is 14 × 3 = 42 cm 2 , while the nominal area of the large NFLA is 47.1 × 11.2 = 527.5 cm 2 . Given that the signal in the loop antenna is properly improved by the area of the antenna, we would expect the larger NFLA to have 12.55 times (or approx 22 dB) more signal than the smaller antenna structure (without the use of the ferrite).
We measured the S 11 of the large FLA and the large NFLA and these results are shown in Figure 14b. These plots show that larger FLA resonance occurs at 120 MHz.

Large Bin Results
We first show a representative S 12,2 result for the FLAs and NFLAs by frequency in Figure 15a. This is followed up by the average of the S x,y signal level over all possible antenna pairs in Figure 15b. For both the FLA and NFLA large-bin results, we ran the parametric inversion solver to obtain the height/cone angle and average permittivity of the grain at frequencies of 120, 125, 130, 135, 140, 150, 160, and 170 MHz, while not shown in figure form, the estimated heights, cone angles, and average permittivities were very similar between the FLA and NFLA results. For example, FLA heights were in the range of 5.23-5.46 m, NFLA heights were 5.23-5.45 m, FLA real permittivity was between 2.35 and 2.47, and NFLA real perm. was 2.37-2.45, the FLA cone angle was between 1 and 25 • , and the NFLA cone angle was between 4 and 25 • . Overall, these were very similar to each other. Of more interest were the correlation coefficients: we again computed the correlation coefficients between the computed fields H and the raw data S. The correlation coefficients for the NFLA and FLA cases are shown in Figure 16.

Large Bin Results Discussion
In turn, each of the bins was filled with approximately 24.5 m 3 (700 bushels) of canola. When the bin was full, we took a full set of electromagnetic measurements. This canola had been stored at room temperature for months and was at a moisture content of 7%. Given the large difference in the size of the loops in the large bin FLA and NFLA antennas, we have shown that the introduction of the ferrite does again increase the raw signal strength using the ferrite. Given that we can expect 12.5 times as much signal from the NFLAs based on loop size, the fact that the average signal strength is within 3-4 dB across the imaging frequency range (see Figure 15b), we can conclude that ferrite antennas have a stronger signal (about 20 dB, similar to the smaller FLAs). Alternatively, one can view that with the ferrite we can reduce the area of the antennas by more than 10 times while maintaining roughly the same signal strength. This has implications for the manufacturing costs and physical robustness of the antennas.
The correlation coefficient results are more puzzling. Unlike the small-bin case, the FLA antennas show a reduced correlation coefficient for all tested frequencies. This means that the bin with the FLA antennas does not match the computational model in the optimization algorithm as well as the bin with the NFLA antennas. We do note that this loss of correlation coefficient does not seem to cause large differences in the result in terms of estimating the height/cone angle/permittivity of the stored grain.
We see several possibilities for the difference in correlation coefficients: (a) due to the enclosure, there is a permanent air gap around the FLA, even when the canola is present, (b) the two bins' dimensions could be different and the FLA bin may have differences in its size/height or deformities in the bin that are not modeled, (c) the fact that the ferrite is not attached directly to the bin wall but is on the steel mounting plate may make a difference, (d) the larger FLAs do not model a magnetic dipole point source as well as the smaller (barrel) FLA antennas, and (e) the FLAs being installed vertically instead of horizontally is making a difference. The cause of the lower correlation coefficients could be a combination of any of these causes. Sorting out these reasons remains the subject of future work.

Discussion and Conclusions
We have designed an improved ferrite-loaded antenna for imaging applications inside metallic enclosures (i.e., grain storage bins). For two experiments of very different bins, the FLAs provide on the order 18-22 dB more S xy signal level than the NFLAs. Inside a smaller experiment, this signal improvement occurs while maintaining (or even improving on) the important characteristics of the NFLAs. Inside the larger bin, the match between the computational model and the antennas appears to be worse for the FLAs. In both large and small experiments, the FLA and NFLA antennas provide similar estimates of the grain height/permittivity and cone angle. Such FLAs remain a promising avenue for improving signal levels in much larger bins and reducing the antenna size in smaller bins where the signal level is already at an acceptable level.
For the small-bin test, we were able to take several measurements of high-moisture targets in the bin which showed that the FLA antennas provided improvements in the detection of such high-moisture targets. Future work will focus on re-designing this antenna concept for use in larger industrial bins and determining methods of improving the computational model of the antenna in the imaging algorithms.