Research Article  Open Access
Hui Liu, Qing Xu, Guowang Jin, Hongmin Zhang, "MIMO Antenna Polynomial Weighted Average Design Method of DownwardLooking Array SAR", International Journal of Antennas and Propagation, vol. 2017, Article ID 3029847, 18 pages, 2017. https://doi.org/10.1155/2017/3029847
MIMO Antenna Polynomial Weighted Average Design Method of DownwardLooking Array SAR
Abstract
MIMO antenna polynomial weighted average design method of downwardlooking array SAR was proposed from the angle of surveying and mapping in this paper, in order to solve the illposed problem that an equivalent virtual array can be implemented by a variety of physical transmitterreceiver arrays for bistatic MIMO linear array. For wave band, resolution, elevation precision, and working height concerned by the applications of surveying and mapping, the length of equivalent virtual array and actual physical array meeting the needs of large scale topographical mapping was solved. Then array numbers and position vectors of MIMO downwardlooking array SAR of real aerial mapping platform were optimized. According to this design, some simulation experiments and comparisons were processed. The results proved the rationality and effectiveness of this array configuration by comparing the differences of 3D imaging results and the original simulation scene, counting mean and standard deviation of elevation reconstruction error eliminating the influence of shadow areas, and counting the probability of elevation reconstruction error within half a resolution of the whole scene and individual building area.
1. Introduction
As a new type of threedimensional synthetic aperture radar (3DSAR) system, downwardlooking array SAR has been paid more and more attention by scholars both at home and abroad. In this system, a linear array antenna is distributed along the crosstrack direction, working at downwardlooking model. And the threedimensional resolution in azimuth direction, range direction, and crosstrack direction is, respectively, achieved by synthetic aperture technique, pulse compression technique, and beamforming technique [1, 2]. The new system effectively resolves the problems about geometric distortion, the blind of the bottom, and the left and right fuzzy of the traditional SAR system. It plays an important role in threedimensional reconstruction of urban areas, topographic survey of complex terrain, disaster prediction and assessment, aircraft navigation, and blind landing.
The crosstrack resolution of downwardlooking array SAR depends on the length of the array antenna. The longer the length of the array is, the higher the resolution of this dimension will be. Then, as the number of required array elements increases more, the design cost and the implementation complexity of the linear array antenna will double, which will make it impossible for us to physically implement such a large array. The introduction of multiple input multiple output (MIMO) technology [3] solves this problem. If we want to obtain the equivalent array with virtual array elements, we can reduce the number of real array elements to by reasonable arrangement of the array elements’ positions, thus greatly reducing the hardware cost and system complexity.
At present, MIMO radar array design is mainly realized by the principle of spatial convolution [4, 5] and equivalent phase center (EPC) [6]. In 2006, the airborne radar for threedimensional imaging and nadir observation (ARTINO) system was developed by the FGANFHR, which used the EPC principle to form 1440 virtual array elements with a spacing of 0.0028 m [7–11]. In 2010, through the derivation of the spatial convolution principle, Su et al. proposed the method using selfconvolutional inverse operation to solve the position of transmitterreceiver common array element [12]. Mitra et al. introduced the array pattern and performed variable substitution and factorization of its analytic expression, so as to realize the design of uniform superposition equivalent array [13]. Chen et al. realized the design of broadband MIMO radar array by combining two uniform linear arrays into an equivalent uniform array [14]. Du et al. proposed an optimized array antenna disposal way for MIMO downwardlooking 3DSAR based on the EPC principle [1]. Zhang et al. proposed a novel MIMO radar array structure optimization algorithm with the highest array element utilization [15]. In 2011, Wang Q. studied the method of MIMO radar array design using polynomial theory [16]. Ye et al. optimized the antenna configuration of a “oneactive” sparse array 3DSAR imaging system with twoterminaltransmission and midreception mode by the EPC principle [17]. In 2012, Peng et al. constructed array antenna configuration of MIMO downwardlooking airborne 3DSAR [18]. Li et al. proposed an improved ARTINO antenna structure, compact midtransmission model, and compact midreception model [19]. In 2013, Wang et al. designed the MIMO antenna array with the polynomial factorization method [20]. Wu et al. solved the multigroup MIMO linear array design scheme by the spatial convolution inverse operation and analyzed the advantages and disadvantages of each group’s design scheme, which provided technical support for the actual construction of airborne linear array 3DSAR imaging system [21]. In 2014, Yu et al. proposed a MIMOSAR array configuration method with midtransmission and twoterminalreception mode and proved that this model can get better imaging performance in downwardlooking 3D imaging [22]. Li et al. gave a sparse array optimization design method based on simulated annealing algorithm [23]. However, for the array configuration of aerial photogrammetric aircraft platform commonly used in surveying and mapping, there is no specific research in the existing literatures.
It is an illposed problem to obtain the actual physical transmitterreceiver array by a given equivalent virtual array for MIMO linear array [21]. If the number of the transmitting array elements and the number of the receiving array elements are factors of the number of equivalent array elements without 1 and itself, the same equivalent virtual array can be constructed by exchange multiplication, which brings the difficulties and opportunities for the optimal design of the array. In this paper, MIMO antenna polynomial weighted average design method of downwardlooking array SAR is introduced from the angle of surveying and mapping. Wave band, resolution, elevation precision, and working height concerned by surveying and mapping application are introduced; the array number and position vectors of MIMO downwardlooking array SAR of real mapping aerial platform are optimized. According to this configuration, some simulation experiments are carried out, proving that the array configuration is rationality and effectiveness.
2. Relationship between Convolution and Polynomial Weighted Average
An example of polynomial multiplication is firstly introduced:
As shown in Figure 1, the coefficients of polynomial are recorded as , which are arranged in descending order. The coefficients of polynomial are recorded as , which are arranged in ascending order. Firstly, let the highest power of the two polynomials be vertically aligned and multiply. Then whenever the ascending polynomial shifts one term to the right, the vertically aligned terms of the two polynomials are multiplied one time, until the lowest power of the two polynomials is aligned. Finally, the result of the multiplication is added to obtain the result polynomial, just as , which is the coefficients of polynomial .
Through the calculation method in Figure 1, it is easy to get
Without loss of generality, we suppose the polynomial coefficients are and , and the result polynomial coefficients can be obtained by multiplying the two polynomials:
It is not difficult to find that (3) is also a convolution expression for and . Therefore, we associate the convolution with the polynomial weighted average. We can take as the weighted average for function to function , where is for the weight role and is for the averaged role.
Literatures [4, 5, 12] have pointed out that the positional relationship of the equivalent virtual array is determined by the spatial position convolution of the transmitting array element and receiving array element under certain conditions. Therefore, from the above analysis, we can see that the equivalent virtual array positional relationship can also be obtained by polynomial weighted average of the spatial position of the transmitting array element and receiving array element under certain conditions.
3. Design Principle of Polynomial Weighted Average for CrossTrack MIMO Antenna Array
The imaging geometry relationship of the MIMO downwardlooking array SAR system is shown in Figure 2 (4 transmitting antennas and 4 receiving antennas). The sparse array antennas distributed along the wing are divided into two parts: the transmitting antenna subarray and the receiving antenna subarray , in which contains transmitting array elements and contains receiving array elements.
From the above analysis, we can see that both the transmitting array and the receiving array can be represented by a polynomial. The position of each element has a certain relation with the power of the polynomial. As shown in (4), we suppose that an array corresponds to a unique polynomial :where is a power corresponding to the position of array elements, is the number of array elements, and is the coefficient of the th term. For a physical array, marks whether there are array elements at the position indicated by the power, just as (5). For a virtual array, it marks the overlapping times of the virtual array element at the position indicated by the power, which represents the weight of the virtual array element:
MIMO antenna design of array SAR is the process of solving transmitting array polynomial and receiving array polynomial by virtual array polynomial. However, a polynomial may not be decomposed into products of two polynomials whose coefficients are positive integers. Therefore, the solution is not unique; that is to say, the transmitting array and the receiving array are not unique. So it needs to add certain constraint conditions to obtain the optimal solution to meet the requirements.
The conventional MIMO antenna design of array SAR used to adopt the twoterminaltransmission and midreception mode. That is, the transmitting antennas are evenly distributed at both ends of the whole array, and the receiving antennas are evenly distributed in between. However, the length of the actual physical array corresponding to this distribution is always longer than the length of the equivalent virtual array, which is intolerable to very limited crosstrack spatial resources. Therefore, according to the symmetrical distribution characteristics of array SAR, this paper seeks an array distribution of equivalent virtual array and real physical array with the same length.
In order to avoid the problem that the nonuniform virtual array can degrade the spatial resolution and image quality of array SAR system, it is desirable that the coefficients of the virtual array polynomial are all 1, and the power of the independent variables is continuous. We suppose that the required uniform virtual array is polynomial :where is the number of virtual array elements and is even due to the symmetrical distribution characteristics. If can be expressed as the product of prime numbers,
Then can be expressed as [16, 20]withwhere .
Through combination multiplication of , we can get multiple sets of transmitting arrays and receiving arrays. However, it should be noted that the changes of the order of will affect the value of .
MIMO antenna design of array SAR requires transmitting antennas to be distributed evenly at both ends, so the spatial position convolution of MIMO antenna array can be equivalent as follows: when the number of transmitting antennas is even, half of them can be placed at the other end of the array. Figure 3 shows an equivalent receiving array based on the principle of spatial position convolution, where 1, 2, and 3 represent the number of times that the virtual array element is overlapped at this position. Figure 4 shows the results when half of the transmitting antennas are placed at the opposite end of the array. The equivalent receiving array is exactly the same to the former.
Therefore, (8) should be added to some constraint conditions to obtain the antenna configuration formula of array SAR:(1)Since the array elements need to be distributed symmetrically, the first prime number of the prime factorization of the equivalent virtual array must be 2.(2)Polynomial, whose all powers are even numbers, is the receiving array, and the power is the position of the receiving array element.(3)Polynomial with uninterrupted arrangement power from 0 to is the transmitting array.(4)According to the symmetrical distribution of the array, transmitting array polynomial is arranged in ascending order and the number of array elements is divided into two parts. The first array elements are located at one end of the physical array, and the power is the position of the transmitting antenna. The left array elements are located at the other end of the physical array, and is the position of the transmitting antenna.
So (8) can be changed to
As an example, we suppose and substitute it into (7). Under the premise of meeting the first constraint condition, 16 can be decomposed into 4 prime factors . Therefore, we can get , substitute it into (11), and further obtain the virtual array polynomial.
Combined with the constraint conditions in (10), we can get the transmitting array polynomial and receiving array polynomial meeting the requirements. The original twelve transmitterreceiver arrays available through any combination of factors have only left the following three:
Combination 1
Combination 2
Combination 3
Combination 1. We consider the polynomial as averaged factor and the polynomial as weight. The polynomial with uninterrupted arrangement power is the transmitting array, and the polynomial is the receiving array, whose all powers are even numbers after expansion. This combination corresponds to the actual physical array with 2 transmitting antennas and 8 receiving antennas. As shown in Figure 5, the positions of the transmitting array and receiving array are, respectively, and .
Combination 2. We consider the polynomial as averaged factor and the polynomial as weight. The polynomial with uninterrupted arrangement power is the transmitting array, and the polynomial is the receiving array, whose all powers are even numbers after expansion. This combination corresponds to the actual physical array with 4 transmitting antennas and 4 receiving antennas. As shown in Figure 6, the positions of the transmitting array and receiving array are, respectively, and .
Combination 3. We consider the polynomial as averaged factor and the polynomial as weight. The polynomial with uninterrupted arrangement power is the transmitting array, and the polynomial is the receiving array, whose all powers are even numbers after expansion. This combination corresponds to the actual physical array with 8 transmitting antennas and 2 receiving antennas. As shown in Figure 7, the positions of the transmitting array and receiving array are, respectively, and .
The length of virtual array obtained by EPC principle is always shorter than the length of the physical array, but the physical array designed by this method is equal in length to the virtual array. Therefore, in the case of limited crosstrack resources, this method has more advantages than EPC principle.
4. MIMO Antenna Design of DownwardLooking Array SAR
MIMO antenna design is the key to the physical realization of downwardlooking array SAR airborne system. As shown in Figure 8, it is necessary to consider the problems about topographic mapping requirements, aircraft platform loads, radar device dimensions, crosstrack sampling intervals (the spatial distribution of equivalent virtual elements should satisfy the spatial sampling law), and so on:(1)Firstly, it should select the highfrequency band with more application potential and determine the interval of the equivalent virtual array element based on the sampling law (in order to avoid the influence of the grating lobe, the distance between two adjacent array elements should be less than or equal to half wavelength [24]).(2)Secondly, according to the requirements about flight height and resolution of topographic mapping and the determined wavelength, the length of equivalent virtual array in corresponding scale mapping is solved.(3)Then, the total number of the array elements is obtained by dividing the length of equivalent virtual array by the interval of equivalent virtual array element. The number of the transmitting arrays and receiving arrays is determined by analyzing the decomposable factors of total number.(4)Finally, the position of the array element in the actual airborne platform is determined by the corresponding method.
In this paper, MIMO antenna array configuration of the real mapping aerial platform is studied in an alternate transmitting mode. In which, only one transmitting array element transmits signal for each pulse repetition time, and all receiving array elements receive the echo of the scene synchronously. After all transmitting array elements sequentially transmit signal, a uniform virtual array with equal spacing can be obtained.
4.1. Topographic Accuracy Requirements
As shown in Table 1, the following requirements for topographic accuracy are listed in the specification named “Basic Geographic Information 3D Model Product Specification” [25]. In this paper, we study MIMO antenna configuration of downwardlooking array SAR based on above specification. In which, the index of restraining the crosstrack resolution of array SAR is “grid spacing,” and the mean square error (MSE) of height in the last four columns restrains the range resolution of array SAR.

4.2. CrossTrack Resolution of MIMO DownwardLooking Array SAR
Literature [1] (p. 32) pointed out that the crosstrack resolution of the MIMO array SAR iswhere is the wavelength, is the platform flight height, and is the synthetic aperture length of the virtual array. Crosstrack resolution is proportional to the wavelength and the flight height and inversely proportional to the synthetic aperture length. In order to obtain higher crosstrack resolution, and should be minimized, and should be larger. Therefore, the Kaband that could obtain wide bandwidth and narrow beam and the lowaltitude platform can be used to achieve the physical realization. Because this technology is limited to lowaltitude mapping, it is not suitable for small and medium scale mapping but suitable for large scale mapping such as 1 : 500, 1 : 1000, and 1 : 2000.
4.3. TransmitterReceiver Antenna Configuration
Airborne platform of array SAR generally includes two types: large aircraft and unmanned aerial vehicles. This paper only studies the former. In the surveying and mapping production, common aerial photogrammetric aircrafts include transport 5, transport 12, PC6, and the air king, whose aircraft wingspan and maximum load are shown in Table 2. The wingspan directly limits the maximum length of the equivalent virtual array, and the maximum load determines the maximum number of array elements that the aircraft can tolerate.

Learning from the experience of optical aerial photography, flight height depends on the pixel size, focal length, and resolution. For different cameras, flight height is roughly distributed in the 580–8000 m range. When the mapping scales are, respectively, 1 : 500, 1 : 1000, and 1 : 2000, relative flying heights are generally around 600 m, 1200 m, and 2500 m. Thus, according to these height constraints of aerial photography, relative flying height of the corresponding scale can be used as a reference to analyzing the length of equivalent virtual array. Furthermore, the “LowAltitude Digital Aerial Photography Specification” [26] also gives other height constraint conditions. The relative flying height of lowaltitude platform is generally not more than 1500 m, and the maximum does not exceed 2000 m. So the reference flying height of 1 : 500, 1 : 1000, and 1 : 2000 and the corresponding 3D topographic accuracy in Table 1 can be, respectively, substituted into (15):
Because MIMO antenna design of downwardlooking array SAR should be applicable to all large scale mapping tasks, we should take the maximum length of virtual array . Combined with the array element spacing the number of required array elements is . However, 1201 is a prime number and can not effectively represent the idea of array configuration. 1200, its approximate value, is used to discuss the polynomial weighted average design method. As shown in Table 3, 1200 equivalent virtual array elements include the following 28 kinds of transmitterreceiver combinations. In order to minimize the number of array elements and the length of actual array and reduce the complexity of industrial manufacturing, it is recommended to select the transmitterreceiver array with 30 transmitting antennas and 40 receiving antennas.
For the array with 1200 equivalent virtual array elements, 1200 can be decomposed into 7 prime factors, when is equal to 6:
Substituting into (10), we can get the virtual array polynomial:
According to (10), the transmitting array of physical transmitterreceiver array with 30 transmitting antennas and 40 receiving antennas is
The receiving array is
5. Experiment and Analysis
In this paper, the MIMO antenna array configuration with 30 transmitting antennas and 40 receiving antennas is used to simulate the echo data of downwardlooking 3D building scene, which is equivalent to 1200 virtual array elements. The simulation parameters are shown in Table 4. The original simulation scene is shown in Figure 9. There exists an Ltype building with 30 meters height in the range of 200 m 200 m. According to the geometric model of downwardlooking array SAR, the height transition between the building and the ground will inevitably produce shadow areas. The existence of shadows is determined by the radar incident angle calculated from the 3D position of the simulation building, and the backscattering coefficient of the shadow areas is set to zero. Figures 10 and 11, respectively, show the top and shadow information of the simulated building.

The original echo data is a chirp signal along range direction as shown in Figure 12. After threedimensional (3D) rangeDoppler (RD) imaging processing, Figures 13, 14, and 15 give the section views of the azimuthangle plane, the anglerange plane, and the azimuthrange plane, respectively. But these imaging results and the original terrain are not in the same coordinate system, so they can not be comparatively analyzed. The imaging results need to be interpolated with the same grid spacing as the original terrain according to the relationship of angle, slope range and height, and slant range. After being mapped in the same coordinate system, the simulation results can be obtained as shown in Figure 16. The original appearance of the Ltype building is basically restored with only some glitches existing in the marginal area, which proves the effectiveness of the method. After the final imaging results are differentiated from the original terrain, the differential contrast image can be obtained. As shown in Figure 17, the red part of the figure shows that the original scene is higher than the reconstruction scene, and the blue part shows the original scene is lower than the reconstruction scene. In this figure, the edge area of the building is mostly red or blue, and the interior of the building is basically green, which indicates that the reconstruction height is more accurate in the interior of the building and the difference is larger in the edge area of the building. The reason is that low crosstrack resolution, inaccurate shaded area heights, and inaccurate phase compensation will directly affect the imaging accuracy.
In order to further verify the effectiveness of the array design, we compare the imaging results of the method in this paper with the results obtained by spatial convolution principle and EPC principle. In order to ensure that the contrast reference is consistent, all three methods should adopt the array configuration with 30 transmitting antennas and 40 receiving antennas, which can be equivalent to 1200 virtual array elements. Different principles can lead to different distances between the transmitting and receiving array elements relative to the same target area, so the compensation phase is also different. In order to simplify the difficulty of imaging, image processing of each principle will do the appropriate approximation.
The equivalent virtual array obtained by spatial convolution principle is equal in length to the actual physical array. And the transmitting array is
The receiving array is
The imaging results corresponding to this proposed method are shown in Figures 18–23.
The equivalent virtual array obtained by EPC principle is shorter than the actual physical array by unit array element distance, where is the number of transmitting antennas. So the transmitting array is
The receiving array is
The imaging results corresponding to this proposed method are shown in Figures 24–29.
From the visual point, the imaging results using above two principles are very similar to those using the method in this paper. And these three methods can well restore the original appearance of the building. In order to quantitatively compare the imaging results of all three methods, Table 5 summarizes the mean and standard deviation of height reconstruction error eliminating the influence of shadow areas and the probability of height reconstruction error within half a resolution of the whole scene and the individual building area. The smaller the mean and standard deviation are and the greater the probability of the error within half a resolution is, the more accurate the elevation will be after imaging. In this paper, the theoretical resolution in the range direction is 1m calculated by [1] (p. 30). From the statistical values in Table 5, the results indicators of EPC principle are higher than spatial convolution principle and this method, which are very limited. And the latter two are equivalent in all indicators. The main reason is that this method is evolved based on spatial convolution principle. When the equivalent uniform array is obtained using the two methods, the different transmitting and receiving array elements correspond to their specific equivalent position, resulting in the compensation phase of each target point being more complicated and more approximate than EPC principle. However, from both the visual point and the statistical values, this method is still able to get good imaging results, which proves that the array design is reasonable and effective. This proposed method can provide another effective strategy for MIMO array design of array SAR.

6. Conclusions
The resolution of the crosstrack for downwardlooking array SAR depends on the array antenna length. The longer the length of the array is, the higher the resolution of this dimension will be. Then, as the number of required array elements increases more, the design cost and the implementation complexity of the linear array antenna will double. In this paper, MIMO antenna polynomial weighted average design method of downwardlooking array SAR is proposed by combining MIMO technology with array SAR. Furthermore, wave band, resolution, elevation precision, and working height concerned by surveying and mapping application are introduced, and the array number and position vectors of MIMO downwardlooking array SAR of real mapping aerial platform are optimized. The simulation experiment proves the rationality and validity of the array design.
In order to promote the practical design of the array antenna, we need to further consider the influence of beam azimuth width and azimuth resolution on the size of subarray, the size of array, and the number of elements and further verify the elevation inversion accuracy, providing a solution to the illposed problem of complex array antenna design with MIMO mode.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China under Grant no. 41071296, no. 41474010, and no. 61401509. In addition, the authors would like to gratefully acknowledge Professor Liangsheng Lou for his array SAR technical support.
References
 L. Du, Study on Model, Algorithm and Experiment for DownwardLooking Synthetic Aperture Radar ThreeDimensional Imaging Based on Linear Array Antennas, Institute of Electronics, Chinese Academy of Sciences, Beijing, China, 2010.
 Y. B. Wang, Research on Key Technologies of New ThreeDimensional Array SAR, University of Electronic Science and Technology of China, Chengdu, China, 2009.
 J. H. G. Ender, “Alongtrack array processing for MIMOSAR/MTI,” in Proceedings of the European Conference on Synthetic Aperture Radar (EUSAR '08), pp. 4032–4036, Friedrichshafen, Germany, 2008. View at: Google Scholar
 L. B. White and P. S. Ray, “Signal design for MIMO diversity systems,” in Proceedings of the Conference Record of the 38th Asilomar Conference on Signals, Systems and Computers, pp. 973–977, Pacific Grove, Calif, USA, November 2004. View at: Publisher Site  Google Scholar
 F. C. Robey, S. Coutts, D. Weikle, J. C. McHarg, and K. Cuomo, “MIMO radar theory and experimental results,” in Proceedings of the Conference Record of the 38th Asilomar Conference on Signals, Systems and Computers, pp. 300–304, Lexinton, Ky, USA, November 2004. View at: Google Scholar
 A. Bellettini and M. A. Pinto, “Theoretical accuracy of synthetic aperture sonar micronavigation using a displaced phasecenter antenna,” IEEE Journal of Oceanic Engineering, vol. 27, no. 4, pp. 780–789, 2002. View at: Publisher Site  Google Scholar
 M. Weiß and J. H. G. Ender, “A 3D imaging radar for small unmanned airplanes—ARTINO,” in Proceedings of the 2nd European Radar Conference ( EURAD '05), pp. 209–212, IEEE, Paris, France, October 2005. View at: Publisher Site  Google Scholar
 J. Klare, “A new airborne radar for 3D imagingsimulation study of ARTINO,” in Proceedings of the European Conference on Synthetic Aperture Radar (EUSAR '06), pp. 3842–3845, Dresden, Germany, 2006. View at: Google Scholar
 J. Klare, D. CeruttiMaori, A. Brenner, and J. Ender, “Image quality analysis of the vibrating sparse MIMO antenna array of the airborne 3D imaging radar ARTINO,” in Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS '07), pp. 5310–5314, Barcelona, Spain, June 2007. View at: Publisher Site  Google Scholar
 M. Weiß, O. Peters, and J. Ender, “First flight trials with ARTINO,” in Proceedings of the 7th European Conference on Synthetic Aperture Radar (EUSAR '08), pp. 1–4, Friedrichshafen, Germany, June 2008. View at: Google Scholar
 M. Weiß and M. Gilles, “Initial ARTINO radar experiments,” in Proceedings of the 8th European Conference on European Conference on Synthetic Aperture Radar (EUSAR '10), pp. 1–4, Aachen, Germany, 2010. View at: Google Scholar
 Y. Su, Y. T. Zhu, W. X. Yu et al., “Design theory and algorithm of multichannel radar antenna arrays,” Chinese Science: Information Science, vol. 40, pp. 1372–1383, 2010. View at: Google Scholar
 S. K. Mitra, K. Mondal, M. K. Tchobanou, and G. J. Dolecek, “General polynomial factorizationbased design of sparse periodic linear arrays,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 57, no. 9, pp. 1952–1966, 2010. View at: Publisher Site  Google Scholar
 A.L. Chen, D.W. Wang, X.Y. Ma, and Y. Su, “Array design and imaging method for wideband MIMO radar,” Systems Engineering and Electronics, vol. 32, no. 12, pp. 2541–2546, 2010. View at: Publisher Site  Google Scholar
 J. Zhang, L.R. Zhang, and N. Liu, “Array optimization algorithm for MIMO radar with a maximum physical sensors efficiency,” Journal of Xidian University, vol. 37, no. 1, pp. 86–90, 2010. View at: Publisher Site  Google Scholar
 Q. Wang, Study on Approaches of MIMO Radar Array Design,, National University of Defense Technology, Changsha, China, 2011.
 Y. Ye, G. Y. Liu, and Z. Meng, “Airborne downwardlooking 3DSAR systems and its imaging using sparse array antennas,” Journal of CAEIT, vol. 6, pp. 96–100, 2011. View at: Google Scholar
 X.M. Peng, Y.P. Wang, W.X. Tan, W. Hong, and Y.R. Wu, “Airborne downwardlooking MIMO 3DSAR imaging algorithm based on crosstrack thinned array,” Journal of Electronics and Information Technology, vol. 34, no. 4, pp. 943–949, 2012. View at: Publisher Site  Google Scholar
 Y. T. Li, Y. G. Chen, Y. K. Deng, and C. B. Yin, “A study on airborne MIMOSAR linear array antenna configuration,” Modern Radar, vol. 34, pp. 4–7, 2012. View at: Google Scholar
 W.Q. Wang, H. Shao, and J. Cai, “MIMO antenna array design with polynomial factorization,” International Journal of Antennas and Propagation, vol. 2013, Article ID 358413, 9 pages, 2013. View at: Publisher Site  Google Scholar
 Z.B. Wu, Y.T. Zhu, Y. Su, Y. Li, and X.J. Song, “MIMO array design for airborne linear array 3D SAR imaging,” Journal of Electronics and Information Technology, vol. 35, no. 11, pp. 2672–2677, 2013. View at: Publisher Site  Google Scholar
 C. Y. Yu, Y. N. Xie, H. L. Zhao et al., “Airborne downwardlooking MIMO 3DSAR imaging properties based on the new array configuration,” Electronic Measurement Technology, vol. 37, pp. 64–69, 2014. View at: Google Scholar
 D. J. Li, Y. N. Hou, X. M. Teng et al., Radar Technology and Application of Sparse Array Antenna, Science Press, Beijing, China, 1st edition, 2014.
 D. W. Bliss and K. W. Forsythe, “Multipleinput multipleoutput (MIMO) radar and imaging: degrees of freedom and resolution,” in Proceedings of the 37th Asilomar Conference on Signals, Systems and Computers, pp. 54–59, Pacific Grove, Calif, USA, November 2003. View at: Publisher Site  Google Scholar
 National Administration of Surveying and Mapping, “CH/Z, basic geographic information threedimensional model product specification (draft),” Standard, Surveying and Mapping Press, Beijing, China, 2012. View at: Google Scholar
 National Administration of Surveying and Mapping, “Lowaltitude digital aerial photography specification,” Standard CH/Z 30052010, Surveying and Mapping Press, Beijing, China, 2010. View at: Google Scholar
Copyright
Copyright © 2017 Hui Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.