Position and Orientation System Error Analysis and Motion Compensation Method Based on Acceleration Information for Circular Synthetic Aperture Radar

: Circular synthetic aperture radar (CSAR) possesses the capability of multi-angle observation, breaking through the geometric observation constraints of traditional strip SAR and holding the potential for three-dimensional imaging. Its sub-wavelength level of planar resolution, resulting from a long synthetic aperture, makes CSAR highly valuable in the field of high-precision mapping. However, the motion geometry of CSAR is more intricate compared to traditional strip SAR, demanding high precision from navigation systems. The accumulation of errors over the long synthetic aperture time cannot be overlooked. CSAR exhibits significant coupling between the range and azimuth directions, making traditional motion compensation methods based on linear SAR unsuitable for direct application in CSAR. The dynamic nature of flight, with its continuous changes in attitude, introduces a significant deformation error between the non-rigidly connected Inertial Measurement Unit (IMU) and the Global Positioning System (GPS). This deformation error makes it difficult to accurately obtain radar position information, resulting in imaging defocus. The research in this article uncovers a correlation between the deformation error and radial acceleration. Leveraging this insight, we propose utilizing radial acceleration to estimate residual motion errors. This paper delves into the analysis of Position and Orientation System (POS) errors, presenting a novel high-resolution CSAR motion compensation method based on airborne platform acceleration information. Once the system deformation parameters are calibrated using point targets, the deformation error can be directly calculated and compensated based on the acceleration information, ultimately resulting in the generation of a high-resolution image. In this paper, the effectiveness of the method is verified with airborne flight test data. This method can compensate for the deformation error and effectively improve the peak sidelobe ratio and integral sidelobe ratio of the target, thus improving image quality. The introduction of acceleration information provides new means and methods for high-resolution CSAR imaging.


Introduction
Synthetic aperture radar (SAR) is widely applied in many fields due to its all-weather and all-day capabilities [1][2][3].Circular synthetic aperture radar (CSAR) has attracted widespread attention due to its unique geometric observation perspective.CSAR is a SAR operating mode that was proposed in the 1990s [4,5].CSAR makes a circular trajectory in the air through the radar platform, and the radar beam always illuminates the scene to achieve all-round observation of the target.The long-time circular synthetic aperture allows CSAR to achieve plane resolution at a sub-wavelength level, eliminating the dependence on system bandwidth [6].Additionally, CSAR has the potential for three-dimensional imaging [7], breaking through the two-dimensional limitations of traditional SAR and mitigating or even eliminating issues such as shadows, geometric deformation, and overlap.With these unique advantages, CSAR holds significant potential for applications in high-precision mapping, ground object classification, target recognition, and more [8][9][10][11][12][13].As a result, CSAR has emerged as a prominent area of research within the SAR field.
The advantages of CSAR are closely tied to its 360-degree observation capability.However, the limited spotlight region is one of the factors that hamper the development of CSAR.The spotlight region is affected by many factors, such as pulse repetition frequency, beam width, and range ambiguity [14].In 2014, Ponce et al. [15] conducted theoretical and technical research on Multiple-Input Multiple-Output (MIMO) CSAR, applying MIMO technology to CSAR to achieve high-resolution imaging of large scenes.In 2021, the AIRCAS conducted research on multi-beam CSAR technology to expand the spotlight area of CSAR.
Compared to traditional linear SAR, CSAR offers unparalleled advantages.However, owing to the intricate motion trajectory of CSAR, traditional imaging algorithms designed for linear trajectories are no longer applicable.Currently, CSAR imaging algorithms can be broadly categorized into two types: time-domain imaging algorithms and frequency-domain imaging algorithms.Time-domain imaging algorithms, such as the crosscorrelation algorithm [16], confocal projection algorithm [6,17], and backprojection (BP) algorithm [18,19], find application in CSAR.Notably, the BP algorithm stands out due to its versatility in handling various imaging geometries.Nonetheless, it is crucial to acknowledge that the BP algorithm is hampered by high computational complexity and low efficiency.On the other hand, frequency-domain imaging algorithms achieve focused imaging through Fourier transforms in the frequency domain.Examples include the polar format algorithm (PFA) [20] and the wavefront reconstruction algorithm proposed by Soumekh [5].Nevertheless, these algorithms encounter challenges associated with their limited adaptability to diverse scenes and complex implementation processes [21].Furthermore, the inherent approximation procedures can affect the imaging accuracy of high-resolution CSAR.
For airborne SAR platforms, imaging errors can be introduced due to airflow disturbances and carrier mechanical vibrations [22].The trajectory of CSAR is intricate, and the synthetic aperture time is prolonged, making the impact of motion errors significant.Unlike traditional linear SAR, CSAR motion compensation does not require aligning with an ideal trajectory.It only needs to correct slant-range errors and phase errors arising from trajectory measurement inaccuracies and other factors [23].The accuracy of the navigation system plays a crucial role in obtaining high-resolution CSAR images.To achieve superior imaging quality, the trajectory measurement errors need to be controlled within 1/16 of the wavelength [24].Additionally, CSAR faces challenges in motion compensation due to the strong coupling between the range and azimuth, along with the complexity of slant-range expressions.These factors contribute to the increased difficulty in motion compensation for CSAR.
In 2011, Lin et al. [25] introduced a CSAR focusing algorithm based on echo generation.This method performs reverse echo generation by extracting the phase error of point targets after coarse imaging to achieve motion compensation.Zhou et al. [22] conducted a detailed analysis in 2013 of the impact of three-dimensional position errors of the radar on the envelope and phase of the echo.In 2015, Zeng et al. [26] proposed a high-resolution SAR autofocus algorithm in polar format to address the challenge of insufficient inertial navigation accuracy for spotlight SAR, thereby achieving high-resolution imaging.In the same year, Zhang et al. [27] presented a CSAR autofocus algorithm based on sub-aperture synthesis to address the issue of diminished imaging resolution resulting from variations in target scattering coefficients with observation angles.In 2016, Luo et al. [28] introduced an extended factor geometric autofocus algorithm for timedomain CSAR.This innovative approach combines factorized geometrical autofocus (FGA) with fast factorized backprojection (FFBP) processing, relying on gradually changing trajectory parameters to obtain a clear image.In 2017, Wang et al. [29] introduced a subbandwidth phase-gradient autofocus method to accomplish motion compensation for high resolution.In 2018, Kan et al. [30] quantitatively analyzed the peak-drop coefficient of the target in the presence of motion errors and assessed the impact of motion errors on imaging quality based on the peak-drop coefficient.In 2023, Li et al. [31] proposed an autofocus method based on the Prewitt operator and particle swarm optimization (PSO).This method aims to enhance imaging quality by optimizing the autofocus process using the Prewitt operator and the PSO algorithm.In general, current motion compensation methods for CSAR can be broadly categorized into two types: one that relies on a high-precision Position and Orientation System (POS), and another based on echo data.However, relying on the motion compensation of a POS is generally challenging, as it is difficult to meet the accuracy requirements of high-resolution CSAR [24,32,33].On the other hand, motion compensation based on echo data involves estimating the motion-induced errors in the echo signal, which can be a complex process and is susceptible to interference from clutter signals [7].The long synthetic aperture in CSAR introduces more complex motion errors compared to linear SAR, making it difficult to directly apply classical autofocus algorithms to CSAR [22].Therefore, specialized motion correction algorithms are needed to address the unique challenges of CSAR.
Currently, motion compensation in CSAR primarily investigates the impact of radar positioning errors on imaging, with a notable scarcity of research on compensating for motion errors using acceleration data.In fact, the degree of deformation experienced by the arm connected to the Inertial Measurement Unit (IMU) and the Global Positioning System (GPS) is proportional to the acceleration.Errors caused by deformation can be determined by analyzing the acceleration information.The innovation of this paper is to analyze the relationship between the acceleration information of the airborne platform and the motion error and propose a high-resolution CSAR motion compensation technology based on the platform's acceleration information.After calibrating system parameters, it uses acceleration data to directly calculate the amount of motion error caused by deformation.In contrast to the traditional method, there is no need to employ point targets for calibrating the phase error of each scene.Point targets only need to be introduced after calibrating the system's deformation error parameters.Then, the phase errors of different scenes or different sub-apertures can be directly calculated according to the acceleration.This method reveals the existence of a correlation between the acceleration and residual phase errors and reduces the need for point targets.Finally, the effectiveness of this method is verified through airborne flight experiments.The results demonstrate that the algorithm significantly improves the peak sidelobe ratio (PSLR) and integral sidelobe ratio (ISLR) of the target, thereby achieving high-quality imaging for CSAR.
This article is structured as follows.Section 2 introduces the signal model of multibeam CSAR and the motion error analysis of CSAR.In Section 3, a high-resolution CSAR motion compensation algorithm is proposed, utilizing the acceleration information of the carrier platform.In Section 4, we introduce the airborne flight experiments and data acquisition.The proposed algorithm is validated through simulation and experimental data in Section 5.These results are further discussed and analyzed in Section 6, where the relationship between the acceleration information and residual phase error is explored.Conclusions are given in Section 7.

Theory
This section introduces the signal reconstruction model of airborne multi-beam CSAR and the motion error analysis of CSAR.

Azimuth Signal Reconstruction Model
The azimuth beam width is one of the factors that restrict the spotlight region of CSAR.Multi-beam CSAR essentially extends the concept of azimuth multi-channel SAR.This is achieved by configuring a beam directed toward different azimuthal directions, effectively expanding the observation area of CSAR.The observation geometry of multi-beam CSAR is depicted in Figure 1, whereas Figure 2 illustrates the observation geometry of CSAR.The SAR platform flies around a circle, and the radar points to a specific angle within each pulse repetition time (PRT), with N PRTs being scanned repeatedly in a group.The data collected at the same angle are considered as an azimuth channel, resulting in N channels of data.The pulse repetition frequency (PRF) of each channel is 1/N of the system's PRF.The multi-channel data essentially utilize the same set of transmitting and receiving antennas, with each channel pointing to different positions at different times.The amplitude and phase errors between the multi-channel data are small and can be considered negligible.This paper primarily focuses on the motion error analysis of CSAR.In this context, the flight direction of the carrier is referred to as the azimuth direction, whereas the radar line-of-sight direction corresponds to the range direction.The positive direction of the X-axis represents the relative starting point, and the azimuth angle is denoted as θ.The radar transmits a chirp signal, with the carrier frequency denoted as f c , pulse width denoted as T p , and frequency modulation denoted as K r .The expression for the chirp signal is as follows: Assuming channel n c as the reference channel and considering a side-looking antenna configuration, the echo signal from a single point target located at P can be expressed as follows: where A 0 is the backscattering coefficient of the target P, τ is the range fast time, η is the azimuth slow time, w r (•) is the range gate function, λ is the working wavelength of the radar, and R p (η) is the distance between the radar and the target P at the azimuth time η.
The signal after pulse compression can be expressed as )) After pulse compression, the data are subjected to Fourier transform in the azimuth direction to obtain S n c (τ, f a ), where f a ∈ [−pr f /2, pr f /2], pr f = PRF/N.The time interval between the adjacent channels of the multi-channel data is expressed as ∆δ = 1/(N • pr f ).
In the case of the azimuth multi-channel data, each channel has a different beam center pointing angle.When the pulse repetition frequency is small, there can be Doppler ambiguity.The multi-channel data are initially deblurred and then interpolated to match the original data length, resulting in the deblurred signal S n (τ, f ), Subsequently, the reconstructed pulse compression data can be obtained by compensating for the delay factor, with n c serving as the reference channel.

Motion Error Analysis of CSAR
This paper employs the BP algorithm for conducting two-dimensional imaging in CSAR.The BP algorithm, being a time-domain imaging algorithm, is not limited by the flight trajectory.It possesses qualities such as high imaging accuracy, broad applicability, and straightforward implementation, making it one of the primary imaging algorithms for CSAR [34].During CSAR imaging processing, it is crucial to select an appropriate reference plane as the focusing plane.The correct reference height enables target focusing, whereas an incorrect reference height results in the target defocusing into a circle [9,35].
The plane with height z re f is selected as the reference plane, and the reference plane is divided into uniform grids to obtain the imaging grid N x × N y .
⌢ s (τ, η), the timedomain signal after pulse compression, is obtained through an inverse Fourier transform of Equation (4).To determine the point P(x, y, z re f ) on the grid, we calculate the time delay ∆t p from the radar to the point at the azimuth time η.Subsequently, we employ interpolation within the signal ⌢ s (τ, η) to locate the echo data that correspond to the calculated time delay ∆t p .This interpolation operation enables us to retrieve the desired echo information for point P on the grid.After that, the interpolated data need to be compensated for the phase factor exp(j4πR zref (η) λ), where R zre f (η) is the distance from the radar at the azimuth time η to the reference plane.The phase factor exp(j4πR zref (η) λ) is related to the reference plane.Finally, the echo data obtained at all azimuth moments of the grid (x, y, z re f ) are coherently superimposed to obtain a two-dimensional image, as shown in Equation (5).
When the airborne platform is flying in a circular trajectory, it is susceptible to external factors that can introduce deviations between the radar measurement coordinates obtained from the navigation system and the actual coordinates.As shown in Figure 3, the radar position information obtained from the navigation system is A, and the actual radar position information is A ′ .For the point target P in the scene, its motion error can be expressed as the slant-range variation.At the azimuth time η, the slant distance from the actual radar coordinate to the target can be expressed as Here, R p (η) is the measured slant distance, and ∆R p (η) is the change in the slant distance caused by the motion error.Without considering the azimuth signal reconstruction, the echo after pulse compression can be expressed as Measurement errors between the actual trajectory and the measured trajectory are inevitable due to atmospheric disturbances and other factors.These errors introduce phase errors in the echo signal, preventing the coherent accumulation of energy and resulting in unfocused images.The phase error caused by residual motion errors in the cross-range resolution unit is the primary factor leading to image defocusing [24].
The multi-beam CSAR extends the spotlight region of CSAR through scanning.The BP algorithm is used to image the signal after azimuth reconstruction, providing robust imaging capabilities.However, the traditional BP algorithm is a time-domain imaging algorithm, lacking a clear Fourier transform relationship in the processed signal.As a result, the conventional autofocus methods cannot be directly applied to the BP imaging algorithm [36].Moreover, the flight trajectory of CSAR is complex, and the aircraft attitude changes frequently.The non-rigid connection between the IMU and GPS introduces significant errors in the POS's solution process.The inaccurate POS position information and complex beam pointing contribute to phase errors in the reconstructed signal, making compensation challenging.In light of these issues, this study focuses on motion compensation technology for high-resolution CSAR and leverages acceleration information from the POS to achieve high-resolution imaging in CSAR.

Method
This paper conducts two-dimensional imaging of multi-beam CSAR data and proposes a motion compensation algorithm for CSAR based on acceleration information from the airborne platform.The algorithm utilizes radial acceleration data from the POS to address phase errors stemming from deformation.Initially, azimuth reconstruction and coarse imaging of multi-beam CSAR data are undertaken.Following this, the ground calibrators within the image are chosen to extract deformation errors.Acceleration is employed for estimating and calibrating the system's deformation error parameters.Ultimately, leveraging the calibrated system parameters, the deformation error of the imaging scene can be computed based on acceleration, facilitating subsequent phase compensation.The proposed method in this paper is categorized into two main parts: data preprocessing and phase error compensation based on acceleration information.The primary processing flow is outlined as follows: 1.
The raw echo data undergo range-pulse compression, and the resulting data postpulse compression are extracted at 1/N intervals to obtain multi-channel data.

2.
The data are converted into the azimuth frequency domain, and the Doppler center frequency is estimated.Subsequently, the data undergo interpolation and de-ambiguity processing.

3.
Compensating for the delay factor in multi-channel data achieves time-domain alignment, followed by the reconstruction of the multi-channel signal with reference to the content in Section 2.1.

4.
The BP algorithm is employed for two-dimensional imaging.Calibrators are selected to mark the deformation error ε(θ) in the two-dimensional frequency domain.

5.
Acceleration information is calculated from the POS, the phase error ε(θ) is linearly estimated from the acceleration, and the system deformation error coefficient K is calibrated.6.
Using the system parameter K calibrated in the previous step, the phase error requiring compensation based on acceleration is calculated, and error compensation in the frequency domain is performed to obtain the focused image after compensation.
The overall flowchart is illustrated in Figure 4.

Raw echo Range pulse compression
Multi-channel data deinterleaving The complete algorithm for data preprocessing is given in Algorithm 1.

Deformation error calibration
Algorithm 1 Multi-beam CSAR data preprocessing.
Input: Radar raw echo Initialization: The raw echo is processed using pulse compression.Define the beam center pointing angle, radar system PRF, azimuth channel number N, and reference channel n c .
Step 1: The multi-beam CSAR data after pulse compression are extracted to obtain N sub-channel data; Step 2: The autocorrelation method is used to estimate and calculate the Doppler center frequency f dc of each channel; Step 3: The Doppler ambiguity m is calculated according to the oblique angle of the antenna beam center, and the ambiguity is solved to obtain Step 4: The multi-channel data are zero-padded in the frequency domain so that the data of each channel are equal to the length of the data before extraction; Step 5: Multiple sets of data are interpolated to compensate for the phase factor exp(j2π f • (n − n c )∆δ), where ∆δ = 1/PRF; Step 6: Multi-channel data are accumulated according to Formula (4) to obtain ⌢ S(τ, f ); An azimuth inverse Fourier transform is performed to obtain the reconstructed time-domain pulse-pressure data Step 7: Time-domain two-dimensional imaging using the BP algorithm.Output: CSAR two-dimensional images.
Assuming that the image with the phase error is represented as s ε (x, y), the spectrum obtained through the two-dimensional FFT is S ε (k x , k y ).
Assuming the target position is (x p , y p ), under ideal conditions, the two-dimensional spectral phase only includes the initial phase of the target and the linear phase related to the target position.By neglecting the initial phase and the envelope of the transmitted signal spectrum, we obtain Equation (11).
The phase difference between S ε (k x , k y ) and S 0 (k x , k y ) is the phase error.From Equation (7), it can be observed that the phase error at the central wavenumber primarily influences imaging.Therefore, the residual phase can be expressed as: The time-domain image obtained through the BP algorithm exhibits periodic blurring in the spectrum when the imaging grid is relatively large, i.e., when it does not satisfy the grid spacing required for full-aperture imaging.Therefore, before establishing the timefrequency relationship between the two-dimensional wavenumber domain (k x , k y ) and the angular wavenumber domain (θ, K), deblurring is necessary.The central frequencies of the two-dimensional spectrum of the radar signal can be expressed as , α is the look angle of the radar, and f c is the radar signal carrier frequency.Under full-aperture conditions, the plane resolution of CSAR can be expressed as where f max = f c + B/2, and B is the signal bandwidth.When the imaging grid spacing is ∆x = ∆y, the frequency range obtained in the f x and f y directions is [−c/(4∆x), c/(4∆x)], and the blur period can be represented as where M x is the blur number in the f x direction, and M y is the blur number in the f y direction.Finally, the residual phase at the central frequency is extracted from the frequencydomain image.However, there is no direct relationship between the azimuth time and phase error.We need to extract the phase error for each sub-aperture.Sometimes, there is no suitable isolated point target in the scene, and the phase error is difficult to obtain directly.Therefore, a more adaptive method is needed.This article proposes a CSAR motion compensation method based on acceleration information.Radial acceleration can, to some extent, reflect the degree of deformation of the lever arm between the IMU and GPS.Furthermore, radial acceleration can reflect motion errors.In this paper, a motion compensation method for CSAR based on acceleration information is proposed, which uses the acceleration to estimate the phase error of each sub-aperture and then compensates to obtain the focused image.The complete algorithm for CSAR motion compensation based on acceleration information is given in Algorithm 2.

Algorithm 2 CSAR motion compensation algorithm based on acceleration.
Input: Two-dimensional CSAR images.Initialization: Select the calibrator and apply windowing to reduce clutter and interference from nearby targets, obtaining data s ε (x, y).
Step 1: Perform a two-dimensional Fourier transform on s ε (x, y) to move to the frequency domain S ε (k x , k y ); Step 2: Deblur in the frequency domain.M x = round( f xc /(c/(2∆x))), M y = round( f yc /(c/(2∆y))); Step 3: Extract the phase at the central wavenumber k c and unwrap the phase; Step 4: Remove the linear phase to extract the residual phase errors ε(θ); obtain the residual motion errors ∆R(θ) based on Equation (12); Step 5: Calculate the radial acceleration ar(θ) using POS data; Step 6: The radial acceleration information ar(θ) and motion error ∆R(θ) are estimated using the least-square method; determine the system parameter K; Step 7: According to the system parameter K obtained in Step 6, the phase error is calculated and compensated through the acceleration, and the compensated spectrum S(k x , k y ) is obtained; Step 8: Perform a two-dimensional inverse Fourier transform on S(k x , k y ) to obtain the motion-compensated image s(x, y); Step 9: For the new scene, after transforming to the two-dimensional frequency domain, Steps 7 and 8 are repeated to obtain the compensated focusing image.Output: Two-dimensional CSAR image after motion compensation.

Airborne Experiments and Data Acquisition
In 2021, AIRCAS conducted airborne multi-beam CSAR flight experiments in Rizhao, Shandong, China.The experimental system is shown in Figure 5.The GPS antenna is located on top of the aircraft, and the phased-array radar closely integrated with the IMU is positioned on the belly of the aircraft.The system is configured with six beam positions, scanning from −10 deg to +10 deg.The aircraft flies counterclockwise at an altitude of 5500 m, with the radar operating in the left-looking mode and the working frequency band in the X-band.Table 1 provides the configuration parameters of the airborne multi-beam CSAR experimental system.Figure 6 illustrates the flight trajectory of the radar platform.Figure 7 shows an optical image and a SAR image of the imaging area.Before imaging, it is necessary to reconstruct the azimuth signal of the obtained multibeam CSAR.This involves deinterleaving the data after pulse compression to extract multibeam data, numbering channel 1 to channel 6, respectively.To estimate the center frequency of the signal in the baseband, the autocorrelation method is employed.Additionally, the Doppler ambiguity is resolved by considering the beam-pointing relationship.The multichannel data are then deblurred and interpolated in the frequency domain, ensuring that the data length matches that of the original data and compensating for the delay phase factor.Subsequently, the multi-channel data are combined into a single-channel signal, with attention given to the energy balance during spectrum splicing.As shown in Figure 8, four sub-apertures with varying angles are chosen for data preprocessing, leading to the generation of reconstructed signals and imaging results.The sub-aperture length is set to 6000 pulses.The reconstructed data can be imaged, but they still have the residual phase error.Therefore, it is necessary to study CSAR motion compensation technology based on acceleration information.

Experimental Results
The proposed method's effectiveness is initially validated through simulation experiments.The verification process involves introducing a known error into the simulation data and utilizing the proposed method for estimation and compensation.It is noteworthy that the introduced error should be associated with acceleration information.The method proposed in this paper primarily focuses on compensating for the system error arising from the deformation of the lever arm due to rapid changes in platform acceleration.In actual data processing, we chose a calibration field for the purpose of validity verification.The process involves coarse imaging of a sub-aperture, followed by deformation phase error extraction and system parameter calibration to establish the relationship between the acceleration and phase error.Subsequently, the acceleration is utilized to directly estimate and compensate for the motion errors of the remaining sub-apertures.An analysis of the results before and after compensation is conducted to validate the efficacy of the proposed method.

Simulation Experiments
This section employs real POS data and incorporates system parameter design to conduct semi-physical simulation experiments.The synthetic aperture size is set to 10 deg.Referring to the actual radial acceleration information, a slant-range error is added, as shown in Figure 9, with an error greater than λ/16.Table 2 provides the coordinates of the scene center and point targets.The imaging grid size is 400 × 400.To verify the effectiveness of the algorithm, the method proposed in this paper was compared with that in [25], as shown in Figure 10 and Table 3.It can be seen that both methods can achieve good focusing effects.However, the algorithm proposed in this article has lower time complexity and less dependence on point targets.The proposed method calculates the deformation error based on the acceleration in the wavenumber domain, with a time complexity of O(N 2 log N).Furthermore, in contrast to traditional methods, which involve coarse imaging and error calculation for each sub-aperture, the proposed approach requires the use of point targets only in the initial sub-aperture to establish the system parameters between the deformation error and acceleration.Subsequently, it can calculate and compensate for the phase errors of the remaining sub-apertures based on acceleration.The system parameters are equally applicable to different scenes, so deformation error compensation can be performed even if there are no suitable point targets in the new scene.This not only saves time in extracting each sub-aperture error but also reduces dependence on point targets.when determining system parameters. 1 The comparison method [25]. 2 The proposed method.
( It is noteworthy that, to some extent, acceleration reflects the degree of deformation between the IMU and GPS, so analyzing the relationship between the acceleration and residual phase error is meaningful.Due to the relatively complex flight trajectory of CSAR, the direction of the aircraft's velocity and acceleration constantly changes, and the non-rigid connection between the GPS and IMU introduces significant errors in the radar position solution of the POS.Therefore, referencing the radar's three-dimensional position coordinates for imaging introduces residual phase errors.This paper creatively introduces acceleration information and utilizes it for residual phase error estimation.Following parameter calibration, the phase error caused by deformation can be estimated and compensated through acceleration.

Real Data Experiments
This section validates the effectiveness of the proposed method using real data.The characteristics of the ground calibrators are analyzed.The ground calibrators used in the experiment are assembled using four triangular reflectors, as shown in Figure 11.Based on the opening directions of the ground calibrators, four sub-apertures are selected in this paper, each with an aperture length of 96,000 pulses.One sub-aperture is used for parameter estimation, and the other three sub-apertures are used to verify the effectiveness of the proposed method.The method in this article needs to determine the parameter relationship between the acceleration and phase error through sub-aperture sub1.It is worth noting that during the parameter determination stage, point targets need to be calibrated.Afterward, both the errors of the other sub-apertures and the errors of other new scenes can be calculated and compensated through acceleration, and point targets are no longer needed.Taking the north direction as an example, the aperture sub1 is used for two-dimensional imaging and parameter estimation.As depicted in Figure 12a, one group of ground calibrators is windowed to suppress the influence of surrounding clutter, and the image s ε (x, y) is obtained.A two-dimensional Fourier transform is performed to obtain the frequency-domain image S ε (k x , k y ).To facilitate subsequent phase extraction, the spectral phase variation should not be too dense.Here, we employ cyclic shifting to position the target point at a corner of the image.This approach effectively reduces errors in phase unwrapping.According to Section 3, we know that for the spectrum phase, we only need to calculate the residual phase at the center of the spectrum, that is, the phase at position ( f xc (θ), f yc (θ)).
As shown in Figure 13, the phase at the center frequency is extracted, and phase unwrapping is performed.Then, the linear phase is removed, as shown in Figure 13b, representing the phase errors that need compensation.According to Equation (12), ∆R(θ) can be calculated for motion compensation.
Figure 14a shows the speed of the carrier in the POS in three directions and derives it to obtain the acceleration in those directions.Figure 14b shows the radial acceleration after coordinate system transformation, and it can be seen that the change in the trend of the radial acceleration is consistent with the change in the trend of the phase error.The phase error is estimated using acceleration information through a least-squares method, and the estimated results are shown in Figure 15.Phase compensation is applied to the scene, and the results are shown in Figure 16. Figure 17 provides information on the calibrator, showing the results of residual phase compensation based on acceleration.The target is effectively compensated, and the experimental results verify the effectiveness of phase error estimation based on acceleration.
The parameter relationship between the acceleration information and motion errors has been established through least-squares estimation.Now, for different sub-apertures, radial acceleration can be calculated from the POS.This parameter is then used for phase error estimation, resulting in the estimation results for phase compensation.Figure 18 shows the extracted phase error and acceleration estimation results for the other three sub-apertures.These results validate the feasibility of the proposed method, and the system parameters estimated through sub-aperture sub1 can be effectively applied to sub-apertures sub2, sub3, and sub4.
Figure 19 shows the spectral phase results of the point target before and after three different sub-aperture motion error compensations.The target characteristics from Figures 17 and 19 are analyzed and the results are presented in Table 4.The effectiveness of the method has been verified with real data.the third row represents the target image after motion compensation using the method proposed in this paper; the fourth row represents the spectrum phase image of the target after motion compensation using the method proposed in this paper.

Discussion
The deformation error of the lever arm exhibits a strong correlation with acceleration.Consequently, this paper introduces a CSAR motion compensation method based on acceleration information.The method initially utilizes point targets to calibrate the system parameters, subsequently calculating the deformation error that requires compensation at each position based on the calibrated parameters and acceleration.Once the system parameters are calibrated, there is no longer a need to employ a point target to estimate the phase error for new scenes and new sub-apertures.This effectively eliminates the dependence on point targets.The method proposed in this paper has been validated through airborne flight experiments, and the experimental results are presented in Section 5.
Firstly, we reconstruct the azimuth signal of multi-beam CSAR based on Algorithm 1. Subsequently, motion compensation based on the acceleration information is executed, following the steps outlined in Algorithm 2. As depicted in Figure 12, the selected calibrator (or isolated strong point) is windowed in the image domain, and the synthetic aperture size is set to 15∼20 deg.Equation ( 14) is employed for deblurring, resulting in a frequency spectrum phase, as illustrated in Figure 12d.For a target, its frequency spectrum phase should only contain an initial phase and a linear phase, with the phase uniformly varying.We extract the phase at the center frequency, and after unwrapping the phase and removing the linear phase, we obtain a residual phase error.This phase error is challenging to express directly in an equation for compensation during the CSAR imaging process.There is a correlation between the phase error and radial acceleration.The parameter relationship between the acceleration information and phase error can be estimated from a sub-aperture, as shown in Figure 15.This allows us to calculate the phase error compensation required at various azimuth times based on the acceleration information.One source of error is the frequent changes in the flight attitude of CSAR.Coupled with the non-rigid connection between the IMU and GPS, the deformation error of the lever arm becomes significant.The variations in acceleration serve as an indicator of the degree of lever arm deformation.Consequently, utilizing acceleration to estimate residual motion errors is a rational and effective approach.Figure 18 illustrates the estimated results for sub-apertures in various directions.It is evident that the radial acceleration information aligns well with the residual motion error.Figures 17 and 19 showcase the motion compensation results for four directions.The compensated target's spectral phase becomes uniform, retaining only the linear phase associated with the target position, resulting in a focused target point.Further analysis of the target's azimuth characteristics reveals a significant improvement in the peak sidelobe ratio and integral sidelobe ratio, confirming the effectiveness of the proposed method.It is noteworthy that this paper introduces acceleration information creatively and verifies its application in high-resolution CSAR imaging.

Conclusions
Considering the complex trajectory of CSAR and the non-rigid connection between the IMU and GPS, the deformation error can cause defocusing in CSAR imaging.Based on the results of airborne flight experiments, this paper determines the correlation between the deformation error and acceleration information and proposes a high-resolution CSAR motion compensation method based on acceleration information.Initially, the parameters for the deformation error are calibrated using a sub-aperture.Subsequently, the deformation error of new scenes or new sub-apertures is computed through acceleration, followed by phase compensation.It is noteworthy that this method requires the use of point targets solely for parameter calibration in the initial sub-aperture.Afterward, both other sub-apertures and new scenes no longer require point targets, thereby reducing dependence on them.After compensation, both the peak sidelobe ratio and integral sidelobe ratio of the target are effectively improved, resulting in enhanced image quality.The effectiveness of the proposed method is validated using airborne measured data.Compared to traditional com-

Figure 1 .
Figure 1.The observation geometry of multi-beam circular SAR.During each pulse repetition time, the radar is directed to a specific angle, scanning repeatedly in three sets within one cycle while the carrier aircraft flies along a circular trajectory.

Figure 2 .
Figure 2. The observation geometry of circular SAR.

Figure 3 .
Figure 3. Circular SAR motion model, where A ′ is the actual position of the radar, and A is the measured radar position.

Figure 4 .
Figure 4. Motion compensation algorithm for high-resolution circular SAR.

Figure 5 .
Figure 5. Airborne multi-beam circular SAR experimental system.(a) System installation diagram.(b) When the acceleration changes, the lever arm undergoes deformation.

Figure 8 .
Figure 8. Multi-beam circular SAR data preprocessing.(a-d) The four sub-apertures with different angles, with an aperture length of 6000 pulses.The first row shows the Doppler center frequency estimation, the second row shows the azimuth signal reconstruction, and the third row shows the imaging results.

Figure 9 .Table 2 .
Figure 9. Motion error.Table 2. Coordinates of the scene center and point targets.Position Coordinate (m)

Figure 10 .
Figure 10.Comparison of CSAR motion compensation before and after.(a) Error image.(b) Phase error estimation.Method 1 is the comparison method [25], and Method 2 is the proposed method.(c) Compensation results of Method 1.(d) Compensation results of the proposed method.

Figure 11 .
Figure 11.Positions of ground calibrators and sub-apertures.(a) The set of ground calibrators.(b) The positions of the sub-apertures align with the opening directions of the ground calibrators.The sub-apertures are numbered from sub1 to sub4.

Figure 15 .
Figure 15.Phase error estimation based on acceleration information.

Figure 16 .
Figure 16.Result after motion error compensation.

Figure 17 .
Figure 17.The results of the calibrator after motion error compensation.(a) The target image before motion compensation.(b) The target image after motion compensation using the method proposed in this paper.(c) The spectrum phase image of the target before motion compensation.(d) The spectrum phase image of the target after motion compensation using the method proposed in this paper.

Figure 19 .
Figure 19.Comparison of the results before and after motion error compensation.(a-c) Sub-apertures sub2, sub3, and sub4, respectively.The first row represents the target image before motion compensation; the second row represents the spectrum phase image of the target before motion compensation;the third row represents the target image after motion compensation using the method proposed in this paper; the fourth row represents the spectrum phase image of the target after motion compensation using the method proposed in this paper.

Table 3 .
Simulation data, characteristic analysis of the target in the azimuth direction.

Table 4 .
Actual data, characteristic analysis of the target in the azimuth direction.