# Spatial Spectrum-Based Imaging for UWB Through-the-Wall MIMO Arrays

Authors: Biying Lu; Xin Sun; Yang Zhao; Zhimin Zhou
Journal: International Journal of Antennas and Propagation
 (2014)

---

## Abstract

To keep the system complexity at a reasonable level and conform to the propagation demands, MIMO arrays are usually sparse in through-the-wall applications, which results in corrupted and gapped data. The corresponding imaging results are seriously affected by the high-level sidelobes. To solve this problem, a new imaging model for ultra-wideband (UWB) MIMO arrays is constructed via spatial spectrum theory in this paper. Based on the model, the characteristics of the spatial spectrum for the MIMO array and its effects on imaging are discussed. To improve the imaging quality, a through-the-wall imaging enhancement method is proposed via spatial spectrum estimation. Synthetic and experimental results show that, unlike the conventional amplitude weighting methods and nonlinear techniques, the proposed method can efficiently suppress sidelobes in the imagery, especially for the sparse MIMO array, and consequently improve the target image quality without degrading the mainlobe resolution. The proposed method has been successfully used in our real through-the-wall radar system.


---

## Body

## 1. Introduction

Ultra-wideband (UWB) through-the-wall imaging (TWI) approaches that can detect objects through obstacles, such as walls, doors, and other opaque materials, are considered powerful tools for a variety of civilian and military applications [1–5].In TWI applications, the imaging component of the application is considered the most important because it is usually the first step for the subsequent processes, such as detection, identification, and wall parameters estimation [6–10]. Currently, to obtain a satisfying target image, two types of radars are widely used: synthetic aperture radar (SAR) and multiple input multiple output (MIMO) radar. Although SAR has better resolution, it has a heavy time cost. By using the high-speed electronic switch, the time to acquire a dataset in a MIMO system is greatly reduced, compared to SAR systems. Therefore, MIMO radar is preferred over SAR in real applications, especially for moving target imaging.By using the proper array design method, we can obtain an optimal array configuration. However, in certain real cases, the equipment complexity and the shape may be our first consideration. Therefore, we make the tradeoff between size and performance [11]. For example, to achieve the Nyquist sampling criterion, the interelement space (d) must be kept below half of the wavelength (λ) for the MIMO array [12]. However, conforming to this criterion will lead to a large number of array elements, even for a small aperture. Usually, when a MIMO array is used in TWI applications, the element spacing is made significantly higher than λ/2 to keep the system complexity at reasonable levels and to increase the element size to achieve an acceptable SNR. Furthermore, for a typical TWI radar system, the most commonly used frequency range is from 1 GHz to 3 GHz to support the range resolution and wall propagation ability. Therefore, for the ultra-wideband signal, even if more elements can be placed in the equipment, the elements are usually dense in the low frequency band but sparse in the high frequency band. In such a case, the MIMO array will not be optimal but it will be sparse with gapped virtual elements, which would otherwise diminish the array imaging performance. As a result, the image quality of the TWI results, in real applications, is significantly limited by the ratio of the main to sidelobe amplitude.To suppress the sidelobes and improve the image quality, many imaging methods for through-the-wall imaging, including the back projection (BP) method [13, 14], the beamforming method [15, 16], and the tomography method [17, 18], are presented in recent years. In these methods, the sidelobes are reduced by applying an amplitude weighting function to the data prior to the final IFFT. However, the sidelobes have been reduced at the expense of the main lobe width, which determines the ultimate resolution of the imagery [19]. For example, the Hanning main lobe is twice as wide (null-to-null) as the sinc function. These methods are consequently a compromise between a narrow main lobe (high resolution) and low sidelobes.To retain the main lobe resolution while reducing the sidelobes, several nonlinear signal processing methods are introduced into radar imaging. Typical methods include spatially variant apodization (SVA), super-SVA, and the CLEAN technique [19–22]. By using interpolation or extrapolation operations, these methods are successfully used in SAR signal data processing to minimize the effects of corrupted and gapped data. However, for MIMO radar, because of the more complicated signal channels, the distribution of the received data is significantly different from that in SAR. In this situation, the performance of these methods is seriously affected.Based on the rigorous derivation of the UWB MIMO array and experimental validation via real TWI radar systems, we proposed in this paper a through-the-wall imaging enhancement method via spatial spectrum theory. Unlike the conventional amplitude weighting methods and nonlinear techniques, the proposed method can effectively suppress the sidelobes from imagery, especially for the UWB sparse MIMO array, and consequently enhance the target image quality without degrading the main lobe resolution.This paper is organized as follows. In Section2, the imaging model for the MIMO array is constructed via spatial spectrum theory. Then, the spatial spectrum of the UWB MIMO array is deeply analyzed. The effects of the spatial spectrum distribution on PSF are discussed, and the spatial spectrum characteristics for the typical TWI UWB MIMO array are obtained. In Section 4, to improve the image quality, an imaging enhancement method by spatial spectrum estimation is proposed. Synthetic and experimental processing results are given in Section 5. Conclusions end this paper.


## 2. Spatial Spectrum-Based Imaging Model for a MIMO Array

We assume a MIMO array has M transmitters and N receivers, as shown in Figure1. For convenience, in the following derivation we choose one transmitter and one receiver, and the polar coordinate origin is assumed at the target center. We assume a stepped frequency waveform as the transmitted signal, which has a frequency range of [f0,fP-1] and a frequency increment of Δf:
(1)s(t)=∑p=0P-1ej2πfpt,fp=(f0+pΔf),p=0,1,…,P-1,
where f0 is the start frequency. The variable P denotes the frequency numbers.Figure 1
The geometry of the MIMO array imaging scene.We consider a general bistatic radar scene with a target as shown in Figure1. The transmitter located at rt=(Rt,θt) transmits the EM wave. The wave arrives at the reflection point located at rs=(ρ,θ) and is backscattered. Under the assumption of free space propagation, the total received backscattered field of the distributed target by the receiver at rr=(Rr,θr) can be expressed as
(2)G(fp,rt,rr)=∫ρ∫θσ(ρ,θ)ej2π(fp/c)(|rt-rs|+|rr-rs|)dρdθ,
where σ(ρ,θ) is the reflection function of the target. Assuming ρ≪ρt and ρ≪ρr, the following approximation holds:
(3a)|rt-rs|≈ρt-ρcos⁡(θ-θt),(3b)|rr-rs|≈ρr-ρcos⁡(θ-θr).Then, (2) becomes
(4)G(fp,rt,rr)=ej2π(fp/c)(ρt+ρr)×∫ρ∫θσ(ρ,θ)e-j2π(fp/c)[ρcos⁡(θ-θt)+ρcos⁡(θ-θr)]dρdθ.The fixed valueej2π(fp/c)(ρt+ρr) outside the integrals can be neglected. We define x=ρcos⁡φ, y=ρsinφ in the equivalent Descartes coordinates and the wavenumber k=2πfp/c in the expression. Then, we rewrite (2) in the Descartes coordinate system:
(5)G(k,rt′,rr′)=∫x∫yσ(x,y)e-jk[(cos⁡θt+cos⁡θr)x+(sinθt+sinθr)y]dxdy=∫x∫yσ(x,y)e-j[(kt,x+kr,x)x+(kt,y+kr,y)y]dxdy,
where rt′ and r′ denote the positions of the transmitter and receiver, respectively, in Descartes coordinates. kt,x=kcos⁡θt and kt,y=ksinθt are the transmitter wavenumbers in the directions of the x-axis and y-axis. kr,x=kcos⁡θr and kr,y=ksinθr are the receiver wavenumbers in the directions of the x-axis and y-axis.We further define the virtual wavenumberskv,x=kt,x+kr,x and kv,y=kt,y+kr,y. Then, the spatial spectrum G(kv,x,kv,y) can be expressed as the 2D Fourier transform of the target reflection function, σ(x,y)(6)G(kv,x,kv,y)=∫x∫yσ(x,y)e-j(kv,xx+kv,yy)dxdy.If the wavenumbers are expressed in the form of vectors askv=[kvx,kv,y]T, kt=[kt,xkt,y]T, and kr=[kr,xkr,y]T, then
(7)kv=kt+kr=2kcos⁡θt-θr2[cos⁡θt+θr2sinθt+θr2]T.
Here, the vector kv is the wavenumber for the virtual element, which, as shown in Figure 2, is the vector sum of the transmitter wavenumber kt and the receiver wavenumber kr.Figure 2
The relationship between the virtual element wavenumber, the transmitter wavenumber, and the receive wavenumber.Therefore, by taking the inverse Fourier transform of the spatial spectrum of the target scattering, that is,GM(kx,ky), the target image can be obtained:
(8)I(x,y)=∫kx∫kyGM(kx,ky)exp⁡[j(kxx+kyy)]dkxdky.As a result, for the MIMO array constructed by M transmitters and N receivers, the spatial spectrum of the received signal is determined by the array structure. By using an ideal point target located at the origin, that is,σ(x,y)=δ(x)δ(y), the point spread function (PSF) of the array system can be demonstrated by (9) in the spatial spectrum domain:
(9)PM(kx,ky)={1,(kx,ky)∈K0,(kx,ky)∉K,
where K is the support area of the spatial spectrum. After the inverse Fourier transform is performed, the spatial expression of the PSF is
(10)PSFM(x,y)=∫kx∫kyPM(kx,ky)kkkkk×exp⁡[j2π(kxx+kyy)]dkxdky.According to (9) and (10), for a MIMO imaging system, the PSF is equivalent to the support area of the spatial spectrum, and the performance of the imaging system can be depicted by the support area of the spatial spectrum.


## 3. Spatial Spectrum Analysis for the UWB MIMO Array

As we know, an imaging system can be fully characterized by the point spread function (PSF) defined as the response of the imaging system to an ideal point source. Equation (10) presents the relationship between the spatial spectrum distribution and the PSF. In this section, the effects of the spatial spectrum distribution on PSF are first discussed. Then, the spatial spectrum of the typical through-the-wall UWB MIMO array is analyzed.


### 3.1. Effects of the Spatial Spectrum Distribution on the PSF

For an imaging system, the ideal support area of the spatial spectrum is an evenly sampled rectangle, as shown in Figure3(a). The corresponding PSF is
(11)PSFr(x,y)=∫-kx0kx0exp⁡(j2πxkx)dkx·∫ky0ky1exp⁡(j2πyky)dky=ϕsin⁡(πΔkxx)πxsin⁡(πΔkyy)πy,
where Δkx=2kx0, Δky=ky1-ky0, and ϕ=ejπ(ky1+ky0)y are the phase term. The resolution of the image is(12a)ρx=1Δkx,(12b)ρy=1Δky.In this situation, the highest level of the sidelobe for the sinc function is −13.3 dB, and the sidelobes are orthogonal and spread along the x- and y-axes (see Figure 4(d)).Geometric shapes of the spatial spectrum support area.
(a)
Rectangle
(b)
Isosceles trapezoid
(c)
Gapped rectangleThree different support areas of the spatial spectrum and the corresponding target images: (a) rectangle, (b) isosceles trapezoid, (c) gapped rectangle, (d) target image for Figure4(a), (e) target image for Figure 4(b), and (f) target image for Figure 4(c).
(a)
(b)
(c)
(d)
(e)
(f)Unfortunately, the support area of a wideband wide beam imaging system is an annulus sector. In this situation, a trapezoid is used instead of a rectangle to approximate the actual spectral support area to increase the resolution [23, 24]. Furthermore, the support area of the acquired spectrum is often not continuous due to incomplete data. In radar imaging, a target may be illuminated by the radar from several separate angular regions to improve cross-range resolution, and, consequently, the sampled data will have missing columns [25, 26]. For the imaging system using fixed aperture, sparse arrays rather than filled arrays are preferred to reduce the number of array elements, which means that spatial samples can more easily be missed [22, 27, 28].For the isosceles trapezoid shown in Figure3(b), the PSF is given as
(13)PSFt(x,y)=∫ky0ky1∫-kytanθ/2kytanθ/2exp⁡(j2πxkx)exp⁡(j2πyky)dkxdky=(cos⁡⁡[π(ky1+ky0)(y+xtanθ2)]kk×sin⁡[π(ky1-ky0)(y+xtanθ2)])×(πx(y+xtanθ2))-1+(cos⁡⁡[π(ky1+ky0)(y-xtanθ2)]kkll×sin⁡[π(ky1-ky0)(y-xtanθ2)])×(πx(y-xtanθ2))-1,
where θ is the processing angle shown in Figure 3(b). Equation (13) indicates that the sidelobes are nonorthogonal and spread along the linear axes with y=±xtanθ/2 (see Figure 4(e)).For the gapped spatial spectrum shown in Figure3(c), the PSF is given as
(14)PSFg(x,y)=(∫-kx0kx0exp⁡⁡(j2πkxx)dkx-∫-kg0kg0exp⁡⁡(j2πkxx)dkx)×∫ky0ky1exp⁡(j2πky)dky=sin(2πkx0x)-sin(2πkg0x)πx·sin(2π(ky1-ky0)y)πy.To analyze the effects of the gapped spatial spectrum on the imaging result, Figure4 shows three support areas of the spatial spectrum, that is, rectangle, isosceles trapezoid and gapped rectangle, and the corresponding target images. For the rectangle and gapped rectangle support area, the sidelobes are orthogonal and spread along the x- and y-axes. When the support area is an isosceles trapezoid, its sidelobes are no longer orthogonal. Additionally, by comparing Figure 4(d) and Figure 4(f), we can find that when the support area is gapped, the sidelobes are significantly enhanced. As a result, the target imaging quality in Figure 4(f) is seriously affected. The PSF figures in the cross-range direction for the ideal support area and the gapped support area are given in Figure 5. The red lines and the blue lines correspond to Figure 4(d) and Figure 4(f), respectively. The main lobe is approximately unchanged, but the sidelobes are significantly different. The highest level of the sidelobes is −13.3 dB in the case of ideal support area and −7.0 dB when the support area is gapped, as in Figure 4(c).Figure 5
The PSF figures in the cross-range direction. In the figure, the blue lines denote the PSF for the ideal support area, that is, Figure4(a). The red lines denote the PSF for the gapped support area, that is, Figure 4(c).According to the derivation and figures, the PSF performance is determined by the spatial spectrum distribution. Furthermore, the distribution and levels of the sidelobes are determined by the shape and density of the spatial spectrum, respectively.


### 3.2. Spatial Spectrum of the Typical TWI UWB MIMO Array

Three characteristics of the spatial spectrum for the UWB MIMO array can be determined. First, a large processing angle is needed in the low frequency band to obtain a satisfying azimuth resolution. However, in real practice, the processing angle is usually limited. Therefore, the support area is no longer nearly rectangular but is described by an annulus sector. In this situation, the trapezoid is usually used instead of the rectangle to approximate the actual spectral support area. Second, to keep the system complexity reasonable, the elements are usually sparse and the element spacing is significantly higher thanλ/2. For the ultra-wideband signal, the spatial spectrum will be dense in the low frequency band and sparse in the high frequency band. Third, for some MIMO arrays with special structures, the corresponding spatial spectrums may be gapped.In order to explain the characteristics, here we take a spare array, for example. The configuration of the array is shown in Figure6(a), and the receiver array consists of 15 elements with an interelement spacing of 0.2 m. Two transmit elements are located at (−2 m, 0 m) and (2 m, 0 m), at the ends of the array. The physical aperture length of the array is 4 m. Using the concept of a virtual phase center, a virtual aperture with 30 virtual transmit/receive (T/R) elements is synthesized, as shown in Figure 6(b). The interelement spacing of the synthesized virtual aperture is 0.1 m except for a gap of 0.6 m at the midpoint, synthesizing a total aperture length of 3.4 m. The gap at the midpoint is due to the separation of the transmit elements.Split transmit virtual aperture array and its virtual elements.
(a)
Split transmit aperture array
(b)
Virtual elementsThe transmitted signal is a stepped frequency waveform, with a range from 0.5 GHz to 1 GHz. The increment frequency is 2 MHz. For a point target located at (0 m, 5 m), we obtain the spatial spectrum support area shown in Figure7. There are 30 dotted lines in the figure, corresponding to the 30 virtual elements, and the envelope of the support area is an isosceles trapezoid. For each dotted line, the points are dense when ky is small (corresponding to low frequency) and sparse when ky is large (corresponding to high frequency). Unfortunately, the support area is missing in the middle region, which is caused by the gap in the virtual elements.Figure 7
Support area of the spatial spectrum for this array.To analyze the effects of the gapped spatial spectrum on target imaging, the result and profile in the cross-range direction are presented in Figures8 and 9, respectively. The whole target image is severely affected, and the highest level of the sidelobes is −10.67 dB.Figure 8
The target imaging result.Figure 9
The profile in the cross-range direction of the target image.


## 3.1. Effects of the Spatial Spectrum Distribution on the PSF

For an imaging system, the ideal support area of the spatial spectrum is an evenly sampled rectangle, as shown in Figure3(a). The corresponding PSF is
(11)PSFr(x,y)=∫-kx0kx0exp⁡(j2πxkx)dkx·∫ky0ky1exp⁡(j2πyky)dky=ϕsin⁡(πΔkxx)πxsin⁡(πΔkyy)πy,
where Δkx=2kx0, Δky=ky1-ky0, and ϕ=ejπ(ky1+ky0)y are the phase term. The resolution of the image is(12a)ρx=1Δkx,(12b)ρy=1Δky.In this situation, the highest level of the sidelobe for the sinc function is −13.3 dB, and the sidelobes are orthogonal and spread along the x- and y-axes (see Figure 4(d)).Geometric shapes of the spatial spectrum support area.
(a)
Rectangle
(b)
Isosceles trapezoid
(c)
Gapped rectangleThree different support areas of the spatial spectrum and the corresponding target images: (a) rectangle, (b) isosceles trapezoid, (c) gapped rectangle, (d) target image for Figure4(a), (e) target image for Figure 4(b), and (f) target image for Figure 4(c).
(a)
(b)
(c)
(d)
(e)
(f)Unfortunately, the support area of a wideband wide beam imaging system is an annulus sector. In this situation, a trapezoid is used instead of a rectangle to approximate the actual spectral support area to increase the resolution [23, 24]. Furthermore, the support area of the acquired spectrum is often not continuous due to incomplete data. In radar imaging, a target may be illuminated by the radar from several separate angular regions to improve cross-range resolution, and, consequently, the sampled data will have missing columns [25, 26]. For the imaging system using fixed aperture, sparse arrays rather than filled arrays are preferred to reduce the number of array elements, which means that spatial samples can more easily be missed [22, 27, 28].For the isosceles trapezoid shown in Figure3(b), the PSF is given as
(13)PSFt(x,y)=∫ky0ky1∫-kytanθ/2kytanθ/2exp⁡(j2πxkx)exp⁡(j2πyky)dkxdky=(cos⁡⁡[π(ky1+ky0)(y+xtanθ2)]kk×sin⁡[π(ky1-ky0)(y+xtanθ2)])×(πx(y+xtanθ2))-1+(cos⁡⁡[π(ky1+ky0)(y-xtanθ2)]kkll×sin⁡[π(ky1-ky0)(y-xtanθ2)])×(πx(y-xtanθ2))-1,
where θ is the processing angle shown in Figure 3(b). Equation (13) indicates that the sidelobes are nonorthogonal and spread along the linear axes with y=±xtanθ/2 (see Figure 4(e)).For the gapped spatial spectrum shown in Figure3(c), the PSF is given as
(14)PSFg(x,y)=(∫-kx0kx0exp⁡⁡(j2πkxx)dkx-∫-kg0kg0exp⁡⁡(j2πkxx)dkx)×∫ky0ky1exp⁡(j2πky)dky=sin(2πkx0x)-sin(2πkg0x)πx·sin(2π(ky1-ky0)y)πy.To analyze the effects of the gapped spatial spectrum on the imaging result, Figure4 shows three support areas of the spatial spectrum, that is, rectangle, isosceles trapezoid and gapped rectangle, and the corresponding target images. For the rectangle and gapped rectangle support area, the sidelobes are orthogonal and spread along the x- and y-axes. When the support area is an isosceles trapezoid, its sidelobes are no longer orthogonal. Additionally, by comparing Figure 4(d) and Figure 4(f), we can find that when the support area is gapped, the sidelobes are significantly enhanced. As a result, the target imaging quality in Figure 4(f) is seriously affected. The PSF figures in the cross-range direction for the ideal support area and the gapped support area are given in Figure 5. The red lines and the blue lines correspond to Figure 4(d) and Figure 4(f), respectively. The main lobe is approximately unchanged, but the sidelobes are significantly different. The highest level of the sidelobes is −13.3 dB in the case of ideal support area and −7.0 dB when the support area is gapped, as in Figure 4(c).Figure 5
The PSF figures in the cross-range direction. In the figure, the blue lines denote the PSF for the ideal support area, that is, Figure4(a). The red lines denote the PSF for the gapped support area, that is, Figure 4(c).According to the derivation and figures, the PSF performance is determined by the spatial spectrum distribution. Furthermore, the distribution and levels of the sidelobes are determined by the shape and density of the spatial spectrum, respectively.


## 3.2. Spatial Spectrum of the Typical TWI UWB MIMO Array

Three characteristics of the spatial spectrum for the UWB MIMO array can be determined. First, a large processing angle is needed in the low frequency band to obtain a satisfying azimuth resolution. However, in real practice, the processing angle is usually limited. Therefore, the support area is no longer nearly rectangular but is described by an annulus sector. In this situation, the trapezoid is usually used instead of the rectangle to approximate the actual spectral support area. Second, to keep the system complexity reasonable, the elements are usually sparse and the element spacing is significantly higher thanλ/2. For the ultra-wideband signal, the spatial spectrum will be dense in the low frequency band and sparse in the high frequency band. Third, for some MIMO arrays with special structures, the corresponding spatial spectrums may be gapped.In order to explain the characteristics, here we take a spare array, for example. The configuration of the array is shown in Figure6(a), and the receiver array consists of 15 elements with an interelement spacing of 0.2 m. Two transmit elements are located at (−2 m, 0 m) and (2 m, 0 m), at the ends of the array. The physical aperture length of the array is 4 m. Using the concept of a virtual phase center, a virtual aperture with 30 virtual transmit/receive (T/R) elements is synthesized, as shown in Figure 6(b). The interelement spacing of the synthesized virtual aperture is 0.1 m except for a gap of 0.6 m at the midpoint, synthesizing a total aperture length of 3.4 m. The gap at the midpoint is due to the separation of the transmit elements.Split transmit virtual aperture array and its virtual elements.
(a)
Split transmit aperture array
(b)
Virtual elementsThe transmitted signal is a stepped frequency waveform, with a range from 0.5 GHz to 1 GHz. The increment frequency is 2 MHz. For a point target located at (0 m, 5 m), we obtain the spatial spectrum support area shown in Figure7. There are 30 dotted lines in the figure, corresponding to the 30 virtual elements, and the envelope of the support area is an isosceles trapezoid. For each dotted line, the points are dense when ky is small (corresponding to low frequency) and sparse when ky is large (corresponding to high frequency). Unfortunately, the support area is missing in the middle region, which is caused by the gap in the virtual elements.Figure 7
Support area of the spatial spectrum for this array.To analyze the effects of the gapped spatial spectrum on target imaging, the result and profile in the cross-range direction are presented in Figures8 and 9, respectively. The whole target image is severely affected, and the highest level of the sidelobes is −10.67 dB.Figure 8
The target imaging result.Figure 9
The profile in the cross-range direction of the target image.


## 4. Through-the-Wall Imaging Enhancement via Spatial Spectrum Estimation

As shown in the analysis above, the serious sidelobes problem in the STVA system is caused by the missing spatial spectrum. To suppress the sidelobes, conventional methods are applied using a weighting function, such as the Hanning, Hamming, or Blackman functions. However, these weighting methods suppress the sidelobes at the expense of the main lobe resolution. A through-the-wall imaging enhancement method via a spatial spectrum estimation for suppressing the sidelobes without degrading the main lobe resolution is presented in this section.The principle of our method is to extrapolate the missing support area using the existing spatial spectrum. Then, according to the filled support area, the sidelobes-suppressed image can be obtained by applying the two-dimensional inverse Fourier transform.The common method for obtaining the support area of the spatial spectrum for the target image is the Stolt interpolation, which is widely used in far-field-based imaging. However, in TWI MIMO radar applications, the distribution of the spatial spectrum is so complicated that the interpolation processing is inaccurate. According to the derived relationship between the image and the spatial spectrum, in our method, the support area is obtained by taking the 2D Fourier transform of the image. Thus, the steps of our method are as follows.Step
1.Perform the imaging processing for the received echo. Here, we take the BP imaging method as an example. In the imaging operation, the region of interest is divided into a finite number of pixels in the range and cross-range directions. The complex amplitude image value for the pixel located at x=(x,y) is obtained by applying frequency-dependent phases and weights to all the received data:
(15)I(x,y)=∭w·S(fm,xT,xR)kkkkk·exp⁡(j2πfmr(xT,x,xR)c)dxTdxRdfm,
where w is the weighting function to shape the beam. S(·) is the received echo. xT and xR denote the transmitter and receiver positions. r(xT,x,xR) is the compensation signal distance for pixel x, which is computed as follows:
(16)r(xT,xq,xR)=(xT-xq)2+(yT-yq)2+(xq-xR)2+(yq-yR)2.Step
2.Take the 2D Fourier transform of the imaging result, I(x,y). Then, the corresponding support area of the spatial spectrum is obtained:
(17)E(kx,ky)=∬I(x,y)exp⁡[-j(kxx+kyy)]dxdy.
Here, E(kx,ky) is the obtained support area of the spatial spectrum.Step
3. Extrapolate the missing spatial spectrum. As shown in Figure 10, the black points denote the sample data in the original spatial spectrum and the white points are the missing spatial spectrum. Therefore, the missing spatial spectrum can be forecasted by the existing data. To obtain an accurate result, the left and right data are used to make forward and back forecasting, respectively. The forward forecasting and the back forecasting AR models are given in (18) and (19), respectively:
(18)EL(kx,p,ky,q0)=-∑i=1LbL(i)E(kx,p-i,ky,q0),(19)EH(kx,p,ky,q0)=-∑i=1HbH(i)E(kx,p+i,ky,q0),
where the variables L and H denote the order of each model. The variables bL(i) and bH(i) denote the coefficient in the models. Here, the model order is evaluated by the forecasting error, and the model coefficient is computed by the Burg algorithm.The missing sample data of the spatial spectrum can be estimated by the existing data. (a) The target is in the middle line of the array. (b) The target is not in the middle line of the array. In this situation, a transfer matrix is needed.
(a)
(b)Step
4. Update the ky,q0 to estimate the missing spatial spectrum. Then, the filled spatial spectrum, denoted by E′(kx,ky), is obtained.Step
5. Take the 2D inverse Fourier transform of E′(kx,ky). Then, the enhanced image, I′(x,y), is obtained:
(20)I′(x,y)=∬E′(kx,ky)exp⁡[j(kxx+kyy)]dkxdky.When the target is not in the middle line of the array, the spatial spectrum will not be symmetrical alongky but will be symmetrical along kv (see Figure 2). In this situation, (18) and (19) need to be reconstructed in a new coordinate system. By using the transfer matrix, the variables kx and ky in these expressions are replaced by kx′ and ky′ in the new coordinate system. The transfer matrix is
(21)[kx′ky′]=[sinθt+θr2-cos⁡θt+θr2-cos⁡θt+θr2sinθt+θr2][kxky],
where (θt+θr)/2 is the angle of vector kv, which is marked in Figure 2.Additionally, in TWI MIMO radar applications, different targets have different incident angles and received angles. Therefore, the spatial spectrums of the targets will be located at different places. To minimize the effects of the overlapped spectrum, the whole image area can be divided into several small subregions. Then, the extrapolation can be applied in each subregion.


## 5. Simulated and Experimental Results

The simulation and experiments used to validate the proposed method are described in this section. In the simulation, the abovementioned array, which is shown in Figure6, is adopted. In the coordinate system, the x-axis is along the array and the origin is the middle point of the array. Three point targets, located at (0 m, 10 m), (−2 m, 9 m), and (2 m, 9 m), are in free space. The stepped frequency signal, with a range from 1 GHz to 2 GHz, is used. The increment frequency is 5 MHz.Using the back projection (BP) imaging method, the original image is given in Figure11. Figure 12 is the processed result made by our proposed method. The sidelobes of the three targets are shown to be efficiently suppressed. To analyze the algorithm performance quantitatively, the peak sidelobe ratio (PSLR) for each target is computed, which are suppressed from −10.0 dB to −15.9 dB, −17.1 dB, and −15.9 dB (see Figure 13). The main lobes for the three targets are unchanged.Figure 11
The original imaging result by the back projection imaging method.Figure 12
The processed result by the proposed method.Figure 13
The sidelobes are suppressed by the proposed method. The PSLR for each target decreased from −10.0 dB to −15.9 dB, −17.1 dB, and −15.9 dB, respectively. The main lobes in this figure are all unchanged.Because the spatial spectrums of the targets overlap (see Figure14), in our processing, we divide the imaging area into several subregions. Therefore, the spatial spectrum for each target is obtained by the 2D Fourier transform, which is shown in Figure 15. From Figures 15(a)–15(c), the spatial spectrums for these three targets are shown to be gapped. After spatial spectrum estimation, the missing spectrums are filled in (see Figures 15(d)–15(f)).Figure 14
The spatial spectrums for the three targets overlap largely when we take the 2D Fourier transform of the whole image area.By dividing the imaging area into several subregions, the spatial spectrum for each target is obtained. (a–c) The gapped spatial spectrum for three targets, namely, P1, P2, and P3, respectively. (e-f) The extrapolated spatial spectrum for each target, that is, P1, P2, and P3.
(a)
(b)
(c)
(d)
(e)
(f)To validate the performance of the proposed method, the through-the-wall imaging experiments are processed in a real environment. In this experiment, a sparse STVA array, which has two transmitters and six receivers, is used. The length of the array is 4.1 m, and two transmitters are placed at the two ends of the array. The height of the array is 1.5 m and the interelement space of the receivers is 0.25 m. The antennas used in the system are Archimedes antennas. A transreceiver module is designed to transmit and receive the EM wave. The waveform used in the system is stepped frequency signal. Its frequency range is from 1 GHz to 2 GHz, with the increment step of 2 MHz. The principle of the radar system is given in Figure16.Figure 16
The principle diagram of through-the-wall radar system.Imaging data were collected by the radar system. Radar system is placed at the left side of a cinderblock building at a distance of 23.7 m. The antennas are parallel to the side wall. As shown in Figure17, the thickness of the left-side cinderblock wall is 30 cm. The width of the left wall is 10 m and its extent is 4 m of the first floor. During the acquisition of the data, a 172-cm-tall man stands behind the wall and moves slightly (see Figure 17).Imaging of a person inside a building: (a) the geometry of the imaging and (b) a man stands behind the left-side cinderblock wall.
(a)
(b)A standard differential back projection (BP) imaging algorithm is adopted to process the acquired data. It is noted that the wall parameters are estimated and compensated by using the image-domain method (see [14] for details). Figure 18 is the original imaging result. The human imaging has strong sidelobes and the image quality is severely affected. Using the proposed method, the processed result is shown in Figure 19. Compared with Figure 18, the sidelobes are efficiently suppressed and the image quality is significantly enhanced.Figure 18
The original imaging result by the standard differential BP imaging algorithm.Figure 19
The processing result by the method proposed in this paper.To validate the algorithm performance for the static target behind the wall, another experiment is processed. In the experiment, we use a 30 cm trihedral as the target. The trihedral is placed in the abovementioned building, with 2 m behind the side cinderblock wall (see Figure20). The radar system is placed at the left side of the building, with a distance of 9 m.Figure 20
In the static target experiment, a trihedral is used as the target and placed 2 m behind the wall.By using the BP imaging method and background subtract technique, the original imaging result is obtained (see Figure21). As shown in the figure, the trihedral is present in the image, but it has strong sidelobes. Then we use the proposed method to extrapolate its spatial spectrum, and consequently, the processing result is obtained (see Figure 22). It is shown that, after processed by the proposed method, the target imaging is significantly enhanced and its sidelobes are suppressed.Figure 21
The original imaging result of the trihedral. In the figure, target has strong sidelobes.Figure 22
By using the proposed method, the sidelobes are suppressed.


## 6. Conclusion

In this paper, we construct an imaging model for an UWB MIMO radar via the spatial spectrum. The rigorous derivation of the model shows that the more spatial spectrum is used, the better imaging performance will be obtained. Therefore, when designing a MIMO array, the best solution is to make full use of the spatial spectrum. Unfortunately, to keep the system complexity at a reasonable level and conform to the propagation demands, MIMO arrays are usually sparse in through-the-wall applications, which results in corrupted and gapped data. The corresponding imaging results are seriously affected by the high-level sidelobes.Aiming at this problem, we proposed a spatial spectrum-based imaging enhancement method in this paper. By estimating the missing spatial spectrum, the effects of the gapped virtual elements can be significantly minimized. The processing results of the synthetic and experimental data show that the proposed method can efficiently improve the imaging quality for both the moving target and static target in through-the-wall applications. Unlike the conventional amplitude weighting methods and nonlinear techniques, the proposed method does not degrade the main lobe resolution when suppressing the sidelobes. At present, the proposed method has been successfully applied to our real through-the-wall radar system.



---

*Source: 825403-2014-07-21.xml*