# Cosine Matching Imaging Method for Rotating Scanning Synthesis Aperture Imaging Radiometer

Authors: Jianfei Chen; Xiaowei Zhu; Sheng Zhang; Yuanjiang Li
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

Millimeter-wave synthesis aperture imaging radiometer (SAIR) can realize high-resolution observations without requiring the real large aperture antenna. Among the SAIRs, the Rotating Scanning SAIR (RS-SAIR) with linear sparse array is a popular system with low redundancy and high reliability. However, due to the lack of matched imaging methods, its imaging precision is usually low. For improving its imaging precision, a novel Cosine Matching Imaging (CMI) method is proposed in this paper. In CMI method, according to the characteristic of rotating projection imaging, the angle-orientation image is constituted by the 1D projection images measured by RS-SAIR in a series of angles firstly. Then, according to the trajectory of the target in angle-orientation image, the pixel values of the brightness temperature image are extracted by cosine matching method from the angle-orientation image one by one. The simulation results demonstrate that the proposed CMI imaging method has higher reconstruction accuracy for the RS-SAIR.


---

## Body

## 1. Introduction

Millimeter-wave synthesis aperture imaging radiometers (SAIR) are powerful sensors for high-resolution observations in various fields [1–3]. By the aperture synthesis technique, SAIR can utilize some small antennas to constitute a large aperture antenna and realize high-resolution imaging observation. Each complex correlation of pairwise antennas is a sample of visibility function, which can be utilized to reconstruct the brightness temperature image of the target scene. Thus, the antenna array configuration is an extremely important component for the SAIR. It directly defines the distribution of the visibility samples and further affects the imaging performance of SAIR system. In fact, how to utilize fewer antennas to achieve the better imaging performance (i.e., array optimization) is always a difficult problem for the SAIR design. Since the first SAIR (ESTAR) developed in 1980s [4, 5], various array configurations have been proposed. Such as 2D rectangular array, sparse “Y” array, “U” array, circular array, and other arrays [6–15]. All of these arrays are looking for achieving the optimal array with low redundancy and larger synthesis aperture constituted by fewer antennas. Unfortunately, there is no effective optimization array which was proposed for 2D SAIR at present. Many antennas are needed for those 2D arrays, which result in high redundancy and system costs.For utilizing antennas effectively and reducing systematic complexity, the scanning SAIRs with sparse array are proposed and applied in some SAIRs [13–17]. Among these SAIRs, the Rotating Scanning SAIR (RS-SAIR) with linear sparse array is the most popular system with low redundancy and high reliability. However, the existing imaging methods may not be accurate enough for some special applications. The sampling points of RS-SAIR is circularly distributed, which is unsuitable for traditional FFT-based methods. In order to deal with these visibility functions with irregular distribution, the Gridding and Nonuniform FFT (NUFFT) methods are proposed for the irregular SAIRs [18–21]. But due to the problem of insertion and resampling, some unnecessary errors will be introduced to the reconstructed images inevitably. Therefore, the effective high-precision imaging method is to find the appropriate inverse method to reconstruct the target images from the irregular visibility function directly. Fortunately, the sampling points of RS-SAIR in radial direction are equispaced and suitable for applying FFT. The 1D radial projection images can be reconstructed by FFT-methods firstly, for example, the back projection (BP) method first to reconstruct the 1D projection images by FFT-methods and then compound the 2D scene images by the back projection operation according to the measuring angle of 1D images [22–24]. But due to the limitations of back projection, BP method cannot effectively deal with the aliasing problem between targets, especially for near-field with wider point spread function (PSF) and serious interference between targets [24–26]. For improving the imaging precision of RS-SAIR further, we study the synthesis aperture projection imaging principle deeply and find that the trajectory of targets in angle-orientation image meets cosine distribution. Based on this fact, a novel Cosine Matching Imaging (CMI) algorithm based on the CLEAN idea is proposed in this paper. In CMI method, the 1D projection images of 2D target scene are measured by the RS-SAIR with linear sparse array in a series of angles firstly. Then, the 1D projection images are arranged to constitute an angle-orientation image according the corresponding array angle. Finally, the pixel values of the brightness temperature image are extracted by the cosine matching method one by one from the angle-orientation image. The numerical simulation results demonstrate that the proposed CMI method can well reconstruct the millimeter images from RS-SAIR with fewer antennas. Its imaging accuracies is better than the BP method under the same conditions.


## 2. The Basic Synthesis Aperture Imaging Theory

Before the description of CMI method for RS-SAIR, the basic synthesis aperture imaging theory is demonstrated here. As Figure1 shows, the antennas are located onOXY, the extend radiation source S is located onoxy, and its brightness temperature image isT(x,y). The radiation source S is dispersed into N small parts. The distances between the i-th radiation source Si and antennas c and l are Ric and Ril, respectively.Figure 1
Geometry diagram of 2D synthetic aperture imaging.The correlation between the received signalsEc and El of pairwise antennas(c, l)is a sample of the so-called visibility function. According to [27], the visibility samples measured by pairwise antennas(c, l) can be expressed as(1)Vc,l=EcRic,t•El∗Ril,t=∑i=1NTxi,yiFcxi,yiFl∗xi,yirc,le-jkRic-Rilwhere < • > denotes time integration operation,(xi,yi) is the coordinate of the point target Si,T is the modified brightness temperature,F(•) is the normalized antenna pattern,k = 2π/λ is angular wavenumber, and rc,l is the so-called fringe-wash function [28].Theexp⁡[-jk(Ric-Ril)] denotes the phase difference of two antennas, which is the key factor for the synthetic aperture imaging. According to Figure 1, the distance Ric and Ril through Taylor expansion can be expressed as (2)Ric=xi-Xc2+yi-Yc2+R2≈R+xi-Xc2+yi-Yc22R(3)Ril≈R+xi-Xl2+yi-Yl22RSubstitute (2) and (3) into (1), and define the spatial coordinatesu=k(Xl-Xc)/R,v=k(Yl-Yc)/R.We can get the relation between visibility function and the brightness temperature image as(4)Vu,v=e-jφu,v∬Tx,yFcFl∗rc,le-jxu+yvdxdywhereφ(u, v)=k(Xc2+Yc2-Xl2-Yl2)/2Ris the phase-modified item, which can improve the imaging accuracy for the near-field SAIR effectively. For the ideal SAIR with identical receivers and antennas, the decorrelation effects rc,l≈1 and antenna patternF(•)can be ignored simply. Then (4) can be rewritten as follows:(5)Vu,v=e-jφu,v∬TOx,ye-jxu+yvdxdy(6)TOx,y=∬Vu,vejφu,vejxu+yvdudvwhereTO(x, y) is the approximate brightness temperature image. For the complete visibility functionV(u, v) measured by the SAIR with “T” array, the brightness temperature imageTO(x,y) can be reconstructed by the FFT-based methods as an approximate solution ofT(x,y).Such as the Modified-FFT (MFFT) imaging method, the brightness temperature images TMF are reconstructed by Fourier transform from the visibility function with phase-modified item directly [26].(7)TMF=FTVu,vejφu,v


## 3. Description of CMI Method for RS-SAIR

For utilizing the antennas more effectively, RS-SAIR only uses a sparse linear array with fewer antennas to detect the 2D target scene. Its basic imaging theory is the aperture synthesis technology also, but the signal processing and imaging steps are different to other SAIRs. Its measuring principle is illustrated in Figure2. The spares array is placed in theXOY plane and rotated around the center point O. θa is the array rotation angle.Figure 2
Geometry diagram of synthetic aperture projection imaging.For the sparse array with givenθa, its spatial coordinatesu=k(Xl-Xc)/R,v=k(Yl-Yc)/R=0underXrOYr coordinates, and the coordinates of targets should be transformed to the projection coordinates(xr,yr), which are the change of rotation angle θa. Then the visibility function Vθa.(u,v) measured by the linear array with given θa can be rewritten as(8)Vθau,0=e-jφu,0∬Txr,yrejxru+yr0dxrdyrVθauejφu=∫∫Txr,yrdyre-jxrudxr=∫TθaPxre-jxrudxrwhere xr is the projection coordinate on the Xr axis (paralleling to the sparse array at angle θa); yr is the projection coordinate on the Yr axis. TθaP(xr)=∫T(xr,yr)dyr is the 1D projection image, which can be reconstructed by the 1D MFFT method as the following formula: (9)TθaPxr=FT1VθauejφuHowever, the value of 1D projection imageTP is a synthetic value containing all the pixels with the same yr coordinate under angle θa. The targets distributed along the Yr axis will be concealed for each other and hardly to be distinguished. Thus, the large reflectors are applied in the 1D SAIR to compress the antenna beam for achieving the high-resolution on the Yr axis [16]. However, the SAIRs with the large reflector usually have a larger volume and are not suitable for many applications. The BP method uses the back projection method to perform the cumulative imaging with array rotating and can reconstruct the target image effectively. But for the near-field RS-SAIR with wider PSF and serious stacking interference between targets, the back projection method cannot deal with these problems effectively; its recovery precision is not good.In this paper, we seek a better solution from another angle. For the given target Si, when the array is rotating from 1 to 180° (rotating step is 1°), its projection coordinate xri can be expressed as(10)xriθa=Ricos⁡θa+θiRi=xi2+yi2,θi=arctan⁡xiyiFor the given targetSi, itsRi andθi are the constants, which depended on the target coordinates(xi,yi) only. Thus, the projection coordinate xri will change with the angle θa of antenna array, and its curve meets the cosine distribution. As Figure 3 shows, the imagesTP are measured by RS-SAIR at a series of angle θa firstly (Figures 3(a) and 3(b)). Then,TP are arranged to constitute the completed angle-orientation image (Figure 3(c)). Obviously, there are two cosine curves in the angle-orientation image.Figure 3
Demonstrative experiment of RS-SAIR. (a) Target scene. (b) The 1D projection images withθa=1~8°. (c) Completed angle-orientation image. (d) The reconstructed image.
(a) (b) (c) (d)Based on this fact, we can arrange the 1D projection images (with differentθa) to constitute the angle-orientation image according the sequence of angle θa. Then the target will be extracted by the cosine matching method from angle-orientation image one by one. Due to the fact the weak targets are usually hidden in the projection of higher targets, the higher targets will be first s extracted partly for exposing the weak targets. Considering the interference noise between targets is higher for near-field RS-SAIR, the nearby targets of the higher target will be extracted at the same time. In addition, each target Si will be extracted repeatedly in CMI method for maintaining the target continuity. Its detailed process can be summarized as follows:(1)
At the start of the CMI method, initialize the brightness temperature imageTB = 0 and draw the cosine images TS-i of each target Si with same dimension (M, N) of angle-orientation image Tao.(11)TS-im,n=1whenm=Ricos⁡n+θi0forotherm,nwhere (Ri,θi) is the polar orientation of target Si.(2)
At thekth step, calculate the correlation valuesCVS-i between image Ts-i and remaining image Tao-k asCVS-i = sum[Ts-i.∗Tao-k]/N. Then find the highest target Sk0 with maxCVS-k0 and its coordinate (xk0, yk0). Record the highest target Sk0 as TB(xk0,yk0)= TB(xk0, yk0)+CVS-k0. Initial subtraction factor matrix SB=0.(3)
Subtract a fraction of targetSk0 from image Tao-k in advance: Tao-k=Tao-k- (γ∧2CVS-k0)TS-k0, where γ is the minus factor. Then update matrix SB as SB(xk0, yk0)= (1-γ)γCVS-k0.(4)
Search the neighbourhood (contains N pixels) of the targetSk 0. If the value CVS-kn of target Skn is higher than the noise thresholdVnoise, it will be considered as an effective target. Record its value CVS-kn and update TB(xkn,ykn)= TB(xkn,ykn)+ CVS-kn. Then update SB (xkn,ykn)=γCVS-kn.(5)
Perform the subtraction of all selected targets in thekth step to form the new remaining image Tao-(k+1)(12)Tao-k+1=Tao-k-∑j=k0j=knSBxj,yjTS-j(6)
Repeat steps(2)-(5) until all the residualCVS-i are below the noise thresholdVnoise. The final image TB is the brightness temperature image.Thekth flowchart of CMI method is given in Figure 4. The core of this CMI method is the selection of γ,N, andVnoise. As a benefit from the parameter γ, target Si will be extracted repeatedly for reducing the aliasing impact between targets and maintaining the target continuity. We know that a higher value of γ can be further reduced the complexity with fewer steps to reach the noise floor. However, the higher γ is only suitable for simple point targets. For the real-life extended target, the higher γ will lead to the result that some smaller targets will be canceled before be extracted. For our simulations, γ is selected by the following formula adaptively: (13)γ=1-meanTaowheremean(•)represents to calculate the mean of normalized image Tao. Clearly, the value of γ is smaller for real-life target with continuous distribution, and is larger for discrete isolated target. The parameters N andVnoise are the important factor for the search step (4). The larger N and smallerVnoise can reduce the time and complexity of CMI method, but the accuracy will be reduced too. Based on a large number of simulation experiments, we have come to the conclusion that N is usually selected as 2 times of the PSF’s main-lobe width (N =10 in this paper). And the noise thresholdVnoise is usually set from 0.1 to 0.2. Figure 3(d) gives the images reconstructed by CMI method for the above two-point scene. It can be seen that two-point targets are reconstructed clearly by CMI method.Figure 4
The kth flowchart of CMI method.


## 4. Simulation Experiments and Results

To verify the validity of the proposed CMI method, many numerical simulation experiments with various scenes have been conducted, and only a few representative ones are demonstrated in this section. The main simulation parameters are listed in Table1. As Table 1 show, the RS-SAIR only needs 16 antennas to constitute the sparse linear array with the larger synthesis aperture. The antennas locations are set based on [29]; its arrangement is shown in Figure 5. The minimum baseline spacing is about 1.1cm, and the synthesis aperture DSA is about 1 m. Then the array rotated around the center to achieve two-dimensional scene detection. That is to say RS-SAIR can achieve the equivalent imaging performance with fewer antennas (about 15% of the “T” array) in theory, and the ratio will be lower for the larger synthesis aperture [29, 30]. In order to simulate RS-SAIR imaging process accurately, the received signals of antennas are obtained by integral operation of the radiation waves generated by all discrete sources, and the visibility samples are calculated by cross-correlated calculation between antennas. The spacing (ΔL) between sources is set to half of the spatial resolution.Table 1
Simulation parameters of RS-SAIR with a sparse array.
Parameters Value Parameters Value Antenna number 16 Rotation angle (θa) 1-180° Wavelength (λ) 3mm Rotating step (Δθ) 2° Synthesis aperture (DSA) 1m Target spacing (ΔL) 4.5mm Baseline spacing (Δu) 1.1cm Imaging distance (R) 3mFigure 5
The arrangement of linear sparse array with 16 antennas (abscissa unit is Δu).In the reconstruction process, the 1D projection images are recovered by the MFFT method firstly. Then the proposed CMI method and BP method are applied to reconstruct the scene images from the 1D projection images. In addition, the NUFFT method [19] is also be applied to reconstruct the images from the irregular visibility function for comparison. All the experiments are performed in MATLAB software. The packaged presentation files used to support the findings of this paper are available from the corresponding author upon request.


### 4.1. Simulation Experiments of the Point Targets

Here, a series of simulation experiments with point targets are carried out to validate the imaging performance of the proposed CMI method for RS-SAIR. It is well known that the point spread function (PSF) is an important parameter for evaluating the resolution of imaging systems. We first make a simulation experiment with a point target. The corresponding results are demonstrated in Figure6.Figure 6
Simulation results of PSF. (a) NUFFT’s PSF. (b) BP’s PSF. (c) CMI’s PSF. (d) Sectional view of three PSFs.
(a) (b) (c) (d)It can be seen from the PSF images that the shape of three PSFs is different. The NUFFT’s PSF is a circular distribution (Figure6(a)) with a fewer noise near the point. The BP’s PSF is similar to the NUFFT’s PSF distribution, but the noise near the point is higher. This is because that the 1D projection images reconstructed by the FFT method usually have higher side-lobe near the targets, then the 1D images are directly used to generate the final image in the BP method by the back projection method. This makes the side-lobe will also focus on the point target just as shown in Figure 6(b). In the proposed CMI method, these side-lobes will be eliminated when the higher target’s cosine images Ts-i are subtracted from image Tao with its nearby targets. As shown in Figure 6(c), the CMI’s PSF is similar to a point distribution compared with the ones of NUFFT and BP, and the noise is eliminated well. For comparing specifically, the 1D sectional view of three PSF images is demonstrated in Figure 6(d). Obviously, CMI’s PSF is closer to a point target. Since the interpolation operation when resampling the visibility function, there is a fewer Gaussian noise near the point target in NUFFT’s PSF. The BP’s PSF is the worst one with the higher cumulative side-lobe noise.For further testing the resolution of the proposed CMI method, a simulation experiment with two-point targets is performed as follows. The space of two-point targets is just the system resolution (i.e., 2ΔL), and the simulation results are shown in Figure7.Figure 7
Simulation results of two-point targets. (a) The angle-orientation image. (b) The result of NUFFT method. (c) The result of BP method. (d) The result of CMI method.
(a) (b) (c) (d)As Figure7 shows, the two-point targets have two cosine curves in the angle-orientation image clearly (Figure 7(a)). However, the images reconstructed by NUFFT and BP method are not very clear; there is a bigger side-lobe noise between two-point targets. It is difficult to distinguish the two points from the Figures 7(b) and 7(c). Actually, as the number of targets increases, this side-lobe noise will increase further. Through the multiple subtraction of higher targets, this side-lobe noise can be eliminated well by the CMI method, the two-point targets can be distinguished from the CMI reconstructed image (Figure 7(d)). This does mean that the imaging accuracy of RS-SAIR can be enhanced by the proposed CMI method effectively.


### 4.2. Simulation Experiments of the Extended Targets

For further demonstrating the imaging performance of the proposed CMI method, the 2D simulation experiments with extended targets are performed here. The target scene radiation intensity (brightness temperature) images are shown in Figures8(a) and 9(a), with their gray value as the amplitude of the discrete point targets, and the distance between point targets is set as ΔL. The corresponding reconstructed images are shown in Figures 8(b)–8(d) and 9(b)–9(d).Figure 8
Simulation results of Tank&Car. (a) Original image. (b) The NUFFT’s reconstructed image. (c) The BP’s reconstructed image. (d) The CMI’s reconstructed image.
(a) (b) (c) (d)Figure 9
Simulation results of Boats. (a) Original image. (b) The NUFFT’s reconstructed image. (c) The BP’s reconstructed image. (d) The CMI’s reconstructed image.
(a) (b) (c) (d)It can be seen from the simulation results that the images reconstructed by the NUFFT and BP methods have serious noise pollution near the target caused by the serious side-lobe and interference between targets (Figures8(b), 8(c), 9(b), and 9(c)). Since the BP method reconstructs the images from 1D projection images with higher noise directly, its reconstructed images are fuzzy. Some weaker targets are even covered by higher noises. Compared to BP’s images, more target information are reconstructed from the resampling visibility function as shown in Figures 8(b) and 9(b). However, due to the influence of resampling error and lack of corresponding noise filtering means, the serious noise pollution still exists in the NUFFT’s images. As a benefit from the multiple extraction operations of block target, these interference noises can be eliminated by the proposed CMI method effectively. As shown in Figures 8(d) and 9(d), the images reconstructed by CMI method are more clearer; the noise is well eliminated while the target information is well preserved. For comparing specifically, the sectional drawing (the middle row) of the reconstructed images (Figures 8(b)–8(d) and 9(b)–9(d)) is demonstrated in Figure 10. Obviously, the noise pollution around target of CMI method is lower than the NUFFT and BP method, and its target information is closer to the real target.Figure 10
1D sectional drawing of two-target scene. (a) Sectional drawing of Tank&Car. (b) Sectional drawing of Boats.
(a) (b)In order to compare the simulation results objectively, the structural similarity (SSIM), peak signal-to-noise ratio (PSNR) of images (Figures8 and 9) and the root-mean square errors (RMSE) of the sectional drawing (Figure 10) are calculated by the following formula:(14)RMSET′,T=∑i,jT′i,j-Ti,j2∑i,jTi,j2(15)SSIMT′,T=2μT′μT+C1σT′T+C2μT′2+μT2+C1σT′2+σT2+C2(16)PSNRT′,T=10log10⁡max⁡T2∑i=1.j=1i=M,j=NT′i,j-Ti,j2/M×NwhereT’ is the reconstructed image and T is the original one. μT​’ and μT are the mean of imageT’ and T separately. σT​’ and σT are the standard deviation of imageT’ and T, respectively.σT’T is the covariance of imageT’ andT. C1 and C2 are the smallest positive constants. max⁡(T) denotes to find the max value of the original imageT.The corresponding calculation results are presented in Table2. Clearly, the images of CMI have the highest SSIM and PSNR, and the RMSE of CMI method is also the lowest one. This indicates that the proposed CMI method can realize a higher imaging accuracy for RS-SAIR and reconstruct the target information accurately while reducing the noise.Table 2
Comparison of objective data between NUFFT, BP, and CMI methods.
Target scene Tank&car Boats RMSE SSIM PSNR RMSE SSIM PSNR NUFFT 0.1683 0.2733 10.6947 0.1288 0.2205 13.7518 BP 0.1730 0.2568 9.7559 0.1345 0.2102 13.0721 CMI 0.1449 0.7256 16.9164 0.1162 0.6893 18.9312


## 4.1. Simulation Experiments of the Point Targets

Here, a series of simulation experiments with point targets are carried out to validate the imaging performance of the proposed CMI method for RS-SAIR. It is well known that the point spread function (PSF) is an important parameter for evaluating the resolution of imaging systems. We first make a simulation experiment with a point target. The corresponding results are demonstrated in Figure6.Figure 6
Simulation results of PSF. (a) NUFFT’s PSF. (b) BP’s PSF. (c) CMI’s PSF. (d) Sectional view of three PSFs.
(a) (b) (c) (d)It can be seen from the PSF images that the shape of three PSFs is different. The NUFFT’s PSF is a circular distribution (Figure6(a)) with a fewer noise near the point. The BP’s PSF is similar to the NUFFT’s PSF distribution, but the noise near the point is higher. This is because that the 1D projection images reconstructed by the FFT method usually have higher side-lobe near the targets, then the 1D images are directly used to generate the final image in the BP method by the back projection method. This makes the side-lobe will also focus on the point target just as shown in Figure 6(b). In the proposed CMI method, these side-lobes will be eliminated when the higher target’s cosine images Ts-i are subtracted from image Tao with its nearby targets. As shown in Figure 6(c), the CMI’s PSF is similar to a point distribution compared with the ones of NUFFT and BP, and the noise is eliminated well. For comparing specifically, the 1D sectional view of three PSF images is demonstrated in Figure 6(d). Obviously, CMI’s PSF is closer to a point target. Since the interpolation operation when resampling the visibility function, there is a fewer Gaussian noise near the point target in NUFFT’s PSF. The BP’s PSF is the worst one with the higher cumulative side-lobe noise.For further testing the resolution of the proposed CMI method, a simulation experiment with two-point targets is performed as follows. The space of two-point targets is just the system resolution (i.e., 2ΔL), and the simulation results are shown in Figure7.Figure 7
Simulation results of two-point targets. (a) The angle-orientation image. (b) The result of NUFFT method. (c) The result of BP method. (d) The result of CMI method.
(a) (b) (c) (d)As Figure7 shows, the two-point targets have two cosine curves in the angle-orientation image clearly (Figure 7(a)). However, the images reconstructed by NUFFT and BP method are not very clear; there is a bigger side-lobe noise between two-point targets. It is difficult to distinguish the two points from the Figures 7(b) and 7(c). Actually, as the number of targets increases, this side-lobe noise will increase further. Through the multiple subtraction of higher targets, this side-lobe noise can be eliminated well by the CMI method, the two-point targets can be distinguished from the CMI reconstructed image (Figure 7(d)). This does mean that the imaging accuracy of RS-SAIR can be enhanced by the proposed CMI method effectively.


## 4.2. Simulation Experiments of the Extended Targets

For further demonstrating the imaging performance of the proposed CMI method, the 2D simulation experiments with extended targets are performed here. The target scene radiation intensity (brightness temperature) images are shown in Figures8(a) and 9(a), with their gray value as the amplitude of the discrete point targets, and the distance between point targets is set as ΔL. The corresponding reconstructed images are shown in Figures 8(b)–8(d) and 9(b)–9(d).Figure 8
Simulation results of Tank&Car. (a) Original image. (b) The NUFFT’s reconstructed image. (c) The BP’s reconstructed image. (d) The CMI’s reconstructed image.
(a) (b) (c) (d)Figure 9
Simulation results of Boats. (a) Original image. (b) The NUFFT’s reconstructed image. (c) The BP’s reconstructed image. (d) The CMI’s reconstructed image.
(a) (b) (c) (d)It can be seen from the simulation results that the images reconstructed by the NUFFT and BP methods have serious noise pollution near the target caused by the serious side-lobe and interference between targets (Figures8(b), 8(c), 9(b), and 9(c)). Since the BP method reconstructs the images from 1D projection images with higher noise directly, its reconstructed images are fuzzy. Some weaker targets are even covered by higher noises. Compared to BP’s images, more target information are reconstructed from the resampling visibility function as shown in Figures 8(b) and 9(b). However, due to the influence of resampling error and lack of corresponding noise filtering means, the serious noise pollution still exists in the NUFFT’s images. As a benefit from the multiple extraction operations of block target, these interference noises can be eliminated by the proposed CMI method effectively. As shown in Figures 8(d) and 9(d), the images reconstructed by CMI method are more clearer; the noise is well eliminated while the target information is well preserved. For comparing specifically, the sectional drawing (the middle row) of the reconstructed images (Figures 8(b)–8(d) and 9(b)–9(d)) is demonstrated in Figure 10. Obviously, the noise pollution around target of CMI method is lower than the NUFFT and BP method, and its target information is closer to the real target.Figure 10
1D sectional drawing of two-target scene. (a) Sectional drawing of Tank&Car. (b) Sectional drawing of Boats.
(a) (b)In order to compare the simulation results objectively, the structural similarity (SSIM), peak signal-to-noise ratio (PSNR) of images (Figures8 and 9) and the root-mean square errors (RMSE) of the sectional drawing (Figure 10) are calculated by the following formula:(14)RMSET′,T=∑i,jT′i,j-Ti,j2∑i,jTi,j2(15)SSIMT′,T=2μT′μT+C1σT′T+C2μT′2+μT2+C1σT′2+σT2+C2(16)PSNRT′,T=10log10⁡max⁡T2∑i=1.j=1i=M,j=NT′i,j-Ti,j2/M×NwhereT’ is the reconstructed image and T is the original one. μT​’ and μT are the mean of imageT’ and T separately. σT​’ and σT are the standard deviation of imageT’ and T, respectively.σT’T is the covariance of imageT’ andT. C1 and C2 are the smallest positive constants. max⁡(T) denotes to find the max value of the original imageT.The corresponding calculation results are presented in Table2. Clearly, the images of CMI have the highest SSIM and PSNR, and the RMSE of CMI method is also the lowest one. This indicates that the proposed CMI method can realize a higher imaging accuracy for RS-SAIR and reconstruct the target information accurately while reducing the noise.Table 2
Comparison of objective data between NUFFT, BP, and CMI methods.
Target scene Tank&car Boats RMSE SSIM PSNR RMSE SSIM PSNR NUFFT 0.1683 0.2733 10.6947 0.1288 0.2205 13.7518 BP 0.1730 0.2568 9.7559 0.1345 0.2102 13.0721 CMI 0.1449 0.7256 16.9164 0.1162 0.6893 18.9312


## 5. Conclusion

RS-SAIR with linear sparse array can utilize some small aperture antennas to realize high-resolution 2D imaging observation. However, due to lack of matched imaging method, its imaging precision is usually low. For improving the imaging accuracy of RS-SAIR further, a novel imaging method (named Cosine Matching Imaging (CMI)) is proposed in this paper. In CMI method, the 1D projection images are reconstructed by the MFFT method at every angle firstly. Then, the 1D projection images are arranged to constitute an angle-orientation image according the measured angle. The brightness temperature image will be extracted by the CMI method from the angle-orientation image in the end. For reducing interference noise of RS-SAIR and maintaining the target continuity, each target is extracted with its nearby targets repeatedly in CMI method. The effectiveness of the proposed CMI is validated by the simulation experiments. The simulation results demonstrate that the CMI method can realize a higher imaging accuracy for RS-SAIR; the images reconstructed by CMI method are clearer than that by NUFFT and BP methods. The imaging error caused by the interference between targets is eliminated well by the proposed CMI method. In addition, although the proposed CMI method is analyzed and verified in near-field condition in this paper, it is also suitable for the far-field RS-SAIR.



---

*Source: 7482753-2018-09-26.xml*