A Modified Keystone-Based Forward-Looking Arc Array Synthetic Aperture Radar 3D Imaging Method

An arc array synthetic aperture radar (AA-SAR) is a new type of omnidirectional observation and imaging system. Based on linear array 3D imaging, this paper introduces a keystone algorithm combined with the arc array SAR 2D imaging method and proposes a modified 3D imaging algorithm based on keystone transformation. The first step is to discuss the target azimuth angle, retain the far-field approximation method of the first-order term, analyze the influence of the forward motion of the platform on the along-track position, and realize the two-dimensional focusing of the target slant range–azimuth direction. The second step is to redefine a new azimuth angle variable in the slant-range along-track imaging and use the keystone-based processing algorithm in the range frequency domain to eliminate the coupling term generated by the array angle and the slant-range time. The corrected data are used to perform along-track pulse compression to obtain the focused image of the target and realize the three-dimensional imaging of the target. Finally, in this article, the spatial resolution of the AA-SAR system in the forward-looking state is analyzed in detail, and the change in the spatial resolution of the system and the effectiveness of the algorithm are verified through simulation.


Introduction
Synthetic aperture radar has become one of the most attractive remote sensing detection tools because of its all-day working ability, high resolution, and large viewing angle [1][2][3][4]. Conventional synthetic aperture radars can reconstruct two-dimensional images of the detection area [5][6][7] with all-day operation capability, while avoiding adverse climate effects. In recent years, an arc array SAR imaging system device/method that breaks through the single observation angle of the linear array SAR system and is no longer limited to front-looking, down-looking, or side-looking imaging has been proposed [8][9][10][11]. The existing arc array two-dimensional imaging system is receiving more and more attention. For example, an arc array bistatic SAR with a stationary transmitter has been proposed by relevant scholars. AA-BiSAR has the advantages of both BiSAR and AA-SAR, which can obtain highly efficient data acquisition and reduce the vulnerability of arc array monostatic SAR. In AA-BiSAR, high-range resolution can be achieved by the high-power wide-bandwidth linear frequency modulation signal transmitted by the stationary transmitter far from the passive receiver. In the azimuth direction, multiple antenna array elements are arranged along a cylinder to form an arc-synthetic aperture, and the switched array antenna channels scan quickly over the circular aperture. Based on this, arc array bistatic SAR with a moving transmitter is proposed; it has the advantages of good concealment and the ability to expand the imaging scene, as well as improved flexibility of the system. AA-SAR with one moving transmitter can be applied to improve the ability of a helicopter to perceive and detect the surrounding targets; since the transmitter can move, multiple low-cost receivers can share the same transmitter to reduce the system cost. Because the arc array SAR has a special observation geometric model and azimuth equal-angle sampling method, the traditional imaging algorithm used for the linear array SAR is no longer suitable for the arc array SAR. The imaging processing method of an arc array SAR is necessary. However, in many application scenarios, two-dimensional images cannot meet the requirements. In the existing literature, there is very little research on the 3D imaging of arc arrays. Therefore, how to use the existing system to achieve 3D imaging in low-altitude airspace becomes an urgent problem to be solved. SAR 3D imaging technology is an extension of SAR 2D imaging, which solves problems such as the overlapping and inversion of traditional imaging. It can realize the reconstruction of the spatial position of the target in the observation area and has great application potential in the fields of battlefield reconnaissance, environmental monitoring, terrain mapping, etc. It has become one of the research hotspots in the field of radar imaging. In recent years, researchers at home and abroad have actively explored SAR three-dimensional imaging technology and achieved results in developing imaging models, imaging theory, imaging algorithms, and experimental systems. At present, common three-dimensional imaging technologies can be divided into the following types according to different imaging modes: InSAR and ROSAR. Assuming that only one scatterer in each resolution unit generates echo signals, InSAR 3D imaging obtains the elevation information of the target through the inversion of the phase difference information [12][13][14]. Multi-baseline SAR tomography is realized through multiple altitude baselines [15][16][17]. The altitude baseline can be realized by time-sharing multiple passes or multiple antennas at the same time, but the altitude blur is mainly due to the unavoidable problem of uneven track distribution. The forward-looking 3D imaging of the linear array antenna can realize a certain resolution imaging of the sector area directly in front of the aircraft route [18][19][20]. However, owing to the limitations of aerodynamics and mechanical principles, the array antenna generally cannot be made very large, resulting in poor azimuthal resolution. In ROSAR 3D technology [21,22], the antenna moves in a uniform circular motion along with the wings, while the carrier platform moves forward to form a horizontal composite surface in the air, combined with the wideband signal emitted in the forward direction. ROSAR uses three-dimensional cognition to detect the scene ahead. However, with the large antenna error with mechanical rotation, it is impossible to obtain accurate equivalent sampling point information, which brings great challenges to later imaging processing. The arc array synthetic aperture radar forward-looking 3D imaging technology works in a new imaging mode, which breaks through the shortcomings of the above four geometric configurations and can meet the requirements of specific aircraft platforms for the 3D imaging of forward targets. This paper introduces an accurate and efficient 3D imaging mode for arc array SARthe forward mode. In the imaging algorithm proposed in this paper, the keystone algorithm is introduced based on the three-dimensional imaging of the linear array combined with the two-dimensional imaging method of the arc array to modify the algorithm so that the arc array SAR has three-dimensional imaging capabilities, making the algorithm suitable for arc array. This mode can ensure the safety of the forward flight of the carrier platform and realize the three-dimensional imaging of the scene ahead. This paper simplifies the echo signal and divides the three-dimensional imaging process of the low-altitude groundmounted reconnaissance platform into two processing steps: slant range-azimuth and slant range-along track. However, this processing introduces two problems: On the one hand, the existence of the target azimuth makes the trigonometric function always exist in the target slant-range equation, which is inconvenient for subsequent processing. On the other hand, the range-cell offset is a common problem in various high-resolution SAR systems, and the imaging capability is highly affected by distance migration correction. In view of the above two problems, first, this paper discusses in detail the target azimuth angle, analyzes the influence of the forward motion of the platform on the along-track direction, and uses approximate processing to make subsequent imaging processing more convenient. Second, the trapezoidal transformation is introduced in this paper, correcting the range migration caused by the coupling term between the slant range and azimuth. Keystone transformation is a data scale remapping method [23][24][25], usually used for range offset correction of objects. Because the transformation involves azimuth and trigonometric functions, the uniqueness of the nonlinear transformation is also explained and demonstrated in this paper. In addition, in this mode, the carrier platform moves forward while the arc array antenna adopts the electronic scanning mode. The opening or closing of the antenna element corresponding to the direction of the target position is selected and controlled by the high-speed microwave switch network, thereby forming a horizontal synthetic front in the air, combined with the broadband signal emitted in the oblique direction, so that it can provide three-dimensional imaging of the scene ahead. It does not increase the hardware complexity of the original system but solves the huge challenge caused by the serious coupling of echo data in the three directions of the slant range-azimuth track.
This paper is organized as follows. In Section 2, the imaging geometric model of the arc array synthetic aperture radar and the observation geometric model of the aircraft platform in the low-altitude reconnaissance mode are given. To facilitate subsequent analysis, the corresponding variable transformation between the coordinate systems is given. At the same time, the echo signal model based on imaging geometry is established. Section 3 analyzes the azimuth angle of the target that is not on the track axis and proposes an imaging algorithm based on KT transform. In Section 4, the spatial resolution of arc array synthetic aperture radar in forward mode is analyzed, and the relationship between azimuth and azimuth is given. In Section 5, numerical simulation experiments are carried out on point targets, and the imaging results of the targets verify the proposed imaging method. Finally, conclusions are given in Section 6.

Imaging Geometry
As a new type of airborne SAR observation system, arc array SAR adopts a wide-area circular-beam scanning method, which has a larger observation area and accelerates the speed of radar acquisition of scene information. Among them, the arc array antenna is the most important part of the system, which is related to whether the whole imaging system can realize all-around fast imaging. In the design of the SAR system, the frequencymodulated continuous-wave (FMCW) transceiver is used, the transmitting antenna and the receiving antenna are arranged separately, and the FMCW system is small in size, light in weight, simple in structure, and easy to use. Compared with pulsed radar, the implementation is easy, costs are low, and the burden of the carrier platform is reduced. In terms of the antenna configuration and imaging geometry of the arc array, its sampling model differs from that of the traditional linear array SAR. Linear arrays sample at equal intervals in azimuth, while arc arrays sample at equal angles in azimuth. The sampling models of the two imaging modes are different, so traditional linear array imaging algorithms cannot be directly used for the imaging processing of arc array SAR; thus, the imaging processing method for arc array SAR needs to be studied. Figure 1a shows the structural diagram of the arc antenna. The transmitting array antenna and receiving array antenna are evenly arranged along the arc-angle direction, which is called the array direction here. The angular spacing between adjacent antenna elements is set to the same ∆ϕ Interval , the radius is R arc , ϕ 0 represents the synthetic aperture angle of the arc array antenna, and the azimuth beamwidth of the antenna element is expressed as β. At the same time, ϕ 0 and β jointly determine the observation range of the effective azimuth of the system. The geometric scene of arc array SAR imaging is shown in Figure 1b. The range resolution is realized by transmitting the bandwidth signal. Meanwhile, the azimuth resolution is achieved through the arc synthetic aperture formed by multiple antenna array elements arranged along the circular arc. Meanwhile, the azimuth resolution is achieved through the arc synthetic aperture formed by multiple antenna array elements arranged along the circular arc.
(a) (b)  Figure 2 shows the geometric model of the scene imaging in the forward mode. The forward direction of the helicopter is defined along the X-axis, the intersection direction is defined along the Y-axis, the elevation direction is defined along the Z-axis, and the raxis represents the slant distance. The platform is set to move forward in a straight line at a constant speed v , where  is the time variable along the track, the pitch angle is  , the pitch beam width is  , and the antenna azimuth beam width is  . O moves to a specific position on the X-axis. We converted the data stored in the cylindrical coordinate system to the rectangular coordinate system according to the following formula for the analysis and discussion below.  Figure 2 shows the geometric model of the scene imaging in the forward mode. The forward direction of the helicopter is defined along the X-axis, the intersection direction is defined along the Y-axis, the elevation direction is defined along the Z-axis, and the r-axis represents the slant distance. The platform is set to move forward in a straight line at a constant speed v, where η is the time variable along the track, the pitch angle is θ, the pitch beam width is ε, and the antenna azimuth beam width is β. The origin of the observation scene coordinates is O, and the center point of the arc array antenna element is O . The cylindrical coordinates of the equivalent sampling point P apc of the arc array are expressed as P apc (ϕ, R arc , h arc ), the angle of the horizontal plane of the adjacent equivalent sampling point P apc is expressed as ϕ, the arc array antenna radius is denoted as R arc , and the platform height is denoted as h arc . The position coordinate of the imaging target P n is expressed as P n (ϕ n , R n , h n ), its azimuth angle is expressed as ϕ n , the ground distance is expressed as R n , and the height is expressed as h n . When the platform transmits signals, the helicopter flies forward, and x l indicates that O moves to a specific position on the X-axis. We converted the data stored in the cylindrical coordinate system to the rectangular coordinate system according to the following formula for the analysis and discussion below.

Echo Signal Model
The arc array SAR imaging adopts the frequency-modulated continuous-wave (FMCW) system, and the transmitted signal is expressed as where 0 f refers to the carrier frequency, r K is the frequency modulation rate,  is the distance-to-fast time variable, and p T is the pulse width of the transmitted signal. The echo signal scattered by the imaging target is received by the equivalent sampling point and expressed as

Echo Signal Model
The arc array SAR imaging adopts the frequency-modulated continuous-wave (FMCW) system, and the transmitted signal is expressed as where f 0 refers to the carrier frequency, K r is the frequency modulation rate, τ is the distance-to-fast time variable, and T p is the pulse width of the transmitted signal. The echo signal scattered by the imaging target is received by the equivalent sampling point and expressed as In Formula (3), c represents the speed of light, σ n (ϕ n ; R n ; h n ) is the backscattering coefficient of the target, and ϕ s represents the beam width of the antenna array element along the horizontal array direction. Mixing the echo signal with the transmitted signal, which is De-chirp processing, the resulting demodulated signal can be expressed as The second exponential term jπK r (τ − 2R m /c) 2 in Equation (4) is the residual video phase term (RVP), and the De-chirp signal after compensation by the RVP term can be given by

Target Azimuth Analysis
According to the azimuth section view shown in Figure 3, the instantaneous slope distance R m from any target in the imaging area to the equivalent sampling point can be written as where ( ) represents the ground distance of the point target when the axis position is at mid x . Define  as the shortest distance from the array center ' O to the target. For the convenience of expression in the following text, substitute (7) into the replacement and perform a Taylor-series expansion, which is expressed as follows: Here, when the platform position moves to mid x , the shortest distance between ' O and the target is expressed as When the platform center O moves to x l , and the carrier platform detects the target in the front area, L = h arc [tan(θ + ε/2) − tan(θ − ε/2)] is defined as the detection range along the track, which is the synthetic aperture length along the track. Here, for the convenience of analysis, the platform moves to the center point x mid of the detectable range L in an ideal situation, x l = x mid + ∆x l is re-expressed, and ϕ mid = a tan[y 0 /(x 0 − x mid )] represents the where R mid = (x 0 − x mid ) 2 + y 0 2 represents the ground distance of the point target when the axis position is at x mid .
Define Ω as the shortest distance from the array center O to the target. For the convenience of expression in the following text, substitute (7) into the replacement and perform a Taylor-series expansion, which is expressed as follows: Here, when the platform position moves to x mid , the shortest distance between O and the target is expressed as At this time, the instantaneous slope distance can be rewritten as

Arc Array SAR 3D Imaging Algorithm
This section analyzes the azimuth angle of the target, which is very important for the approximate processing of the distance between the platform and the target and gives the specific imaging algorithm flow. The complete imaging flow chart of the arc array synthetic aperture radar low-flying reconnaissance mode is shown in Figure 4. After demodulating the echo signal, perform RVP compensation, and focus the processed echo signal in the azimuth and slant directions to obtain a single SAR image. Then, the coupling between the range frequency domain variable and the azimuth angle is successfully eliminated by keystone transformation, and the focused image is outputted after eliminating the two-dimensional coupling term and along-track matched filtering.

Arc Array SAR 3D Imaging Algorithm
This section analyzes the azimuth angle of the target, which is very important for the approximate processing of the distance between the platform and the target and gives the specific imaging algorithm flow. The complete imaging flow chart of the arc array synthetic aperture radar low-flying reconnaissance mode is shown in Figure 4. After demodulating the echo signal, perform RVP compensation, and focus the processed echo signal in the azimuth and slant directions to obtain a single SAR image. Then, the coupling between the range frequency domain variable and the azimuth angle is successfully eliminated by keystone transformation, and the focused image is outputted after eliminating the two-dimensional coupling term and along-track matched filtering.

Azimuth-Slant Range Processing
In Equation (5), let ; the echo signal demodulated to the baseband signal can be written as

Azimuth-Slant Range Processing
In Equation (5), let f = f 0 + K r τ; the echo signal demodulated to the baseband signal can be written as ss(τ; ϕ; R m ) = σ n (ϕ n ; r n ; h n )rect We can easily find that when the carrier platform moves forward at a certain sampling point along the heading, a synthetic aperture is obtained in the azimuth direction so that the target in the front scene can be distinguished. After the above echo signals are postprocessed by focusing on the azimuth and oblique directions, the obtained single SAR image can be expressed as where A 0 is the amplitude of the target signal after azimuth-slant range focusing, the signal distance to bandwidth is expressed as B r , and L eff is the effective synthetic aperture of the target azimuth [26,27]. By intercepting any azimuth direction ϕ n = ϕ x l , the two-dimensional data expression of the slant-range direction and the along-track direction can be obtained as where A 1 = σ n (ϕ n ; r n ; h n )sin c L eff λR m ϕ − ϕ x l . We ignore it in subsequent imaging processing to simplify the process.

Improved Transformation Based on Keystone
Fourier transform is performed on the signal in Equation (12) along the slant range; the result can be computed as follows: where f r represents the range frequency domain variable. It can be obtained from the above analysis that in the actual situation R arc Ω, the Taylor expansion of the slant range is performed, and the first-order term far-field approximation method is retained, so R m can be approximately expressed as Substituting the approximate slope distance into Equation (13), the signal can be written as For the convenience of analysis, Formula (15) can be rewritten as According to the above formula, in the second exponential term exp j4π f r c R arc cos θ cos ϕ − ϕ x l , the range frequency domain variable f r is coupled with the azimuth angle ϕ. The distance cell migration can be corrected by removing the coupling before along-track focusing is performed. In this paper, the method of keystone transformation is used to write the transformation relationship between the target azimuth angle ϕ and the virtual sample α as Equation (17), and the distance travel term is eliminated by redefining a virtual azimuth sample α for conversion.
Since the transformation involves trigonometric functions, we next analyze the uniqueness of the result as follows. In the actual situation, because of the difference in the distance of the observation target, the relationship between target ground distance and array distribution affects the antenna synthetic aperture angle ϕ 0 . In general, it is best to select the synthetic aperture angle of the array antenna to be smaller than the beam width of the antenna array, which is ϕ 0 β. However, here, we consider taking ϕ 0 = β in the ideal state to obtain a better resolution in the array direction because the arc array antenna selects and controls the opening or closing of the corresponding azimuth antenna element according to the target position through the high-speed microwave switch. Therefore, in the actual case ϕ − ϕ x l ≤ β/2 ≤ 90 • , so in the angular range of [0 • , 90 • ], cos ϕ − ϕ x l is a monotonic function, so Formula (17) can guarantee the uniqueness after transformation. Therefore, after substituting (17) into (16) and undergoing keystone transformation, the equation can be re-expressed as Equation (18): It can be seen that after the keystone transformation, the coupling term between the range frequency domain variable f r and the azimuth angle ϕ is successfully eliminated, and the range cell migration is corrected.
For the AA-SAR in the forward mode, its main task is to image the point target in the front or oblique front area. Within the detectable range of the radar, the scene's width is relatively small, which means the azimuth angle of the target is small. Therefore, within an appropriate range, for the same point target, we ignore the azimuth difference between the different axis x l positions when the aircraft moves forward; in the following processing, it can be approximated as ϕ x l = ϕ mid . According to Expansion (10), the expression of Equation (18) in the along-track frequency domain can be written as where K w represents the distance wavenumber domain variable and K x represents the along-track wavenumber domain variable. It can be seen from the above formula that the fifth exponential term is a two-dimensional coupling term, and Equation (20) gives the decoupling function to correct the coupling term.
The signal after the correction of the two-dimensional coupling term is subjected to the inverse fast Fourier transform of the slant range, and the signal in the frequency domain along the track in the time domain of the slant range is obtained as follows: The fourth exponential term in Equation (21) is the along-track phase modulation information, and the along-track matching function is constructed as After eliminating the phase information of the along-track direction, an inverse Fourier transform of the along-track direction is performed to complete the along-track imaging; the signal is compressed as It can be obtained from the above formula that after along-track matching imaging, the actual position of the target at Ω mid is shifted by x mid cos θ cos(α − ϕ mid ) along the slant range, and a corresponding geometric correction is required for this. The slice data corresponding to different azimuths in the oblique range and along the track are selected in turn for two-dimensional processing and rearranging the acquired two-dimensional data according to the above steps to obtain the three-dimensional imaging of the front target. The expression is written as where A 6 represents the following: In Equation (24), B r represents the range bandwidth, B a = L eff λR m is the azimuth bandwidth, and B x = 2ε sin θ cos 2 (α−ϕ mid ) λ is the signal along-track bandwidth. So far, the three-dimensional imaging of the front target of the entire arc array carrier platform has been completed.

Spatial Resolution Analysis
In this paper, a forward-looking AA-SAR 3D imaging model is constructed, and the corresponding imaging algorithm is proposed. In this chapter, the spatial resolution of AA-SAR is analyzed in detail when moving forward, and the simulation experiments are carried out on point targets to verify the validity of the model and algorithm. Spatial resolution, one of the important indicators used to evaluate the effectiveness of radar imaging algorithms, reflects the ability of radar to distinguish targets in different positions. The system parameters of the carrier platform are shown in Table 1, simulating the lowflying situation of the carrier platform and selecting a smaller pitch angle for the observation to reach a farther scene. To ensure that there is no ambiguity problem along the track and, at the same time, the low-altitude aircraft is safe, the aircraft is usually required to fly at a slow speed. Since the azimuth angle information corresponds to the azimuth direction information, the relationship between the target azimuth angle and the azimuth is y 0 = R x l sin ϕ x l , so we derive the azimuth resolution according to the angular resolution.
Using the geometric model, we know that the arc array SAR antenna array is arranged in turn along the angle, so the angular resolution is the array resolution. Because the instantaneous frequency in the azimuth direction is always tangent to the arc where the antenna element is located, the special antenna structure makes the direction of the instantaneous frequency change with the sampling point. The azimuth instantaneous phase of the point target can be expressed as ; from this, the azimuth instantaneous frequency formula can be deduced as When the observation target is very far away from the radar antenna, it is known that the azimuth half-wave width will not exceed under normal circumstances, so the value range of ϕ − ϕ x l in Formula (26) is less than or equal to half of the azimuth beam width. Within this range, the azimuth instantaneous frequency is a monotonic function, and the angular resolution can be approximately expressed as Formula (27) shows that the azimuth resolution is related to the antenna radius, target slant range, and platform height and changes with the target ground distance. According to the relationship between azimuth and azimuth, to reduce the influence of azimuth changes, the scope of the detection scene is limited, and the azimuth angle is taken as a smaller value. Thus, we can approximate the relationship between the azimuth angle and azimuth as y 0 = R x l ϕ x l and obtain the azimuth information, and the azimuth resolution can be expressed as To observe the effect of target height, ground distance, and point azimuth angle on the azimuth resolution, a comparative experiment was carried out. Figure 5 shows that the azimuth resolution does not depend on the azimuth angle of the target or whether the target is at the edge of the beam. The azimuth resolution of the target at the same distance is the same, but it has a certain relationship with the antenna radius and the azimuth beam width.
To observe the effect of target height, ground distance, and point azimuth angle on the azimuth resolution, a comparative experiment was carried out. Figure 5 shows that the azimuth resolution does not depend on the azimuth angle of the target or whether the target is at the edge of the beam. The azimuth resolution of the target at the same distance is the same, but it has a certain relationship with the antenna radius and the azimuth beam width.

Along-Track Resolution
After the above-mentioned spatial position correction and 3D imaging processing, the 3D imaging expression of the target can be obtained as The comparison experiment of the along-track resolution is shown in Figure 6. It can be concluded that the along-track resolution will be affected by the target height, ground distance, and azimuth. In addition, it will also be affected by the forward-looking SAR beam front view angle and the elevation beam width.

Along-Track Resolution
After the above-mentioned spatial position correction and 3D imaging processing, the 3D imaging expression of the target can be obtained as According to the along-track information contained in Equation (29), the along-track resolution can be expressed as The comparison experiment of the along-track resolution is shown in Figure 6. It can be concluded that the along-track resolution will be affected by the target height, ground distance, and azimuth. In addition, it will also be affected by the forward-looking SAR beam front view angle and the elevation beam width.

Height-Wise Resolution
Using the target slope distance Ω = R x l − R arc 2 + (h arc − z 0 ) 2 , it can be known that the height information is hidden in the slope distance information. To obtain the height resolution, we introduce the definition of the fuzzy function [28,29] and construct the fuzzy function of the target P n , which is different from the point target P n in the height direction by ∆h as follows: According to (31), the height resolution can be expressed as In (32), the height resolution of the forward-looking AA-SAR imaging area depends on the slant-range bandwidth of the radar transmit signal and the imaging geometry of the radar; its value is equal to the slant-range resolution of a typical SAR multiplied by the reciprocal of the sine value of the viewing angle of a beam front. This shows that the height resolution mainly comes from the effective bandwidth or effective time width of the radar transmit signal in the height direction. Figure 7 is a comparison experiment of different target height resolutions. It can be seen that the height resolution is affected by the target slant range and target height in the imaging geometry, but the target azimuth does not affect the height resolution.

Point Target Imaging Analysis
To verify the effectiveness of the model and algorithm, we conducted simulation experiments on point targets. The simulation results were analyzed as follows. Figure 8a is the distribution map of the real position of the target in the imaging scene, and Figure 8b is the imaging result of the target in the observation scene.

Point Target Imaging Analysis
To verify the effectiveness of the model and algorithm, we conducted simulation experiments on point targets. The simulation results were analyzed as follows. Figure 8a is the distribution map of the real position of the target in the imaging scene, and Figure 8b is the imaging result of the target in the observation scene.

Point Target Imaging Analysis
To verify the effectiveness of the model and algorithm, we conducted simulation experiments on point targets. The simulation results were analyzed as follows. Figure 8a is the distribution map of the real position of the target in the imaging scene, and Figure 8b is the imaging result of the target in the observation scene. Point targets P5 and P9 were selected for measurement and analysis to further verify the effectiveness of the proposed imaging method. Figure 9 shows the single-point 3D imaging of point targets P5 and P9, which shows that the imaging point target has good focusing performance. Point targets P5 and P9 were selected for measurement and analysis to further verify the effectiveness of the proposed imaging method. Figure 9 shows the single-point 3D imaging of point targets P5 and P9, which shows that the imaging point target has good focusing performance.  Figure 10 shows the cross-sectional view and response function of point target P5 in three dimensions, and Figure 11 shows the cross-sectional view and response function of point target P9 in three dimensions, which shows that the focusing effect of the two-point targets in the along-track, azimuth, and height directions is good.  Figure 10 shows the cross-sectional view and response function of point target P5 in three dimensions, and Figure 11 shows the cross-sectional view and response function of point target P9 in three dimensions, which shows that the focusing effect of the two-point targets in the along-track, azimuth, and height directions is good.  Figure 10 shows the cross-sectional view and response function of point target P5 in three dimensions, and Figure 11 shows the cross-sectional view and response function of point target P9 in three dimensions, which shows that the focusing effect of the two-point targets in the along-track, azimuth, and height directions is good.  The two-dimensional slice diagram of the point target in Figure 10 shows that threedimensional imaging can be well realized because there is no azimuth difference at this position. The point targets in Figure 11 have an azimuth difference, so there is a partial mismatch in the two-dimensional slice image. In addition, the three-dimensional response functions of point targets are all sinc functions, which shows that the resolution of the The two-dimensional slice diagram of the point target in Figure 10 shows that threedimensional imaging can be well realized because there is no azimuth difference at this position. The point targets in Figure 11 have an azimuth difference, so there is a partial mismatch in the two-dimensional slice image. In addition, the three-dimensional response functions of point targets are all sinc functions, which shows that the resolution of the airborne forward AA-SAR three-dimensional mode meets the imaging requirements, which further verifies the effectiveness of the arc array SAR three-dimensional imaging algorithm based on keystone transform studied in this paper.

Conclusions
Arc array 3D imaging will have greater advantages in the future, with a wider market and a wider imaging field of view. To quickly obtain omnidirectional imaging and overcome the shortcomings of traditional array SAR 3D imaging, this paper proposes a 3D scene imaging mode suitable for arc arrays in a low-altitude ground-based reconnaissance mode. Based on the premise of not increasing the complexity of the original system, the AASAR platform is used to move forward to obtain the second synthetic aperture, and a corresponding imaging algorithm is proposed for the above modes. Because of the special geometric configuration of AASAR, the external loader platform moves forward at a constant speed, resulting in complex antenna motion trajectories, which brings difficulties to subsequent imaging. This paper proposes a corresponding solution, which simplifies the AASAR forward mode into two steps: First, the target distance equation is simplified, the target azimuth angle is discussed, and the influence of the forward motion of the platform on the along-track position is analyzed to further realize the two-dimensional focusing of the target slant range-azimuth direction. Second, a new azimuth angle variable is redefined in the slant-range on-track imaging, and a 3D imaging algorithm based on keystone transformation is proposed to eliminate the coupling term generated by the array angle and slant range time to realize the 3D imaging of the target. In addition, the spatial resolution of arc array SAR is analyzed in detail, focusing on the relationship between the target azimuth and azimuth. Finally, the algorithm's effectiveness is verified through numerical simulation experiments involving point targets.