Distributed ISAR Subimage Fusion of Nonuniform Rotating Target Based on Matching Fourier Transform

In real applications, the image quality of the conventional monostatic Inverse Synthetic Aperture Radar (ISAR) for the maneuvering target is subject to the strong fluctuation of Radar Cross Section (RCS), as the target aspect varies enormously. Meanwhile, the maneuvering target introduces nonuniform rotation after translation motion compensation which degrades the imaging performance of the conventional Fourier Transform (FT)-based method in the cross-range dimension. In this paper, a method which combines the distributed ISAR technique and the Matching Fourier Transform (MFT) is proposed to overcome these problems. Firstly, according to the characteristics of the distributed ISAR, the multiple channel echoes of the nonuniform rotation target from different observation angles can be acquired. Then, by applying the MFT to the echo of each channel, the defocused problem of nonuniform rotation target which is inevitable by using the FT-based imaging method can be avoided. Finally, after preprocessing, scaling and rotation of all subimages, the noncoherent fusion image containing all the RCS information in all channels can be obtained. The accumulation coefficients of all subimages are calculated adaptively according to the their image qualities. Simulation and experimental data are used to validate the effectiveness of the proposed approach, and fusion image with improved recognizability can be obtained. Therefore, by using the distributed ISAR technique and MFT, subimages of high-maneuvering target from different observation angles can be obtained. Meanwhile, by employing the adaptive subimage fusion method, the RCS fluctuation can be alleviated and more recognizable final image can be obtained.


Introduction
Conventional monostatic ISAR images are obtained under the condition that the change of the radar observation angle is small (<5 o ) [1]. For modern high maneuvering aircraft with high speed, stealth and other characteristics, even a slight change of the observation angle can cause a fluctuation of 10 to 15 dB on the RCS [2]. The strong fluctuation of the RCS will cause a deterioration in the ISAR image quality. Therefore, the imaging quality of monostatic ISAR for a high-maneuvering target is easily influenced by the observation angle. Meanwhile, after translation motion compensation, the nonuniform rotation caused by the maneuvering motion will degrade the imaging performance of the conventional FT-based method.
Unlike the conventional monostatic ISAR, the distributed ISAR technique can utilize the data acquired from multiple observation angle to improve the image quality [3][4][5][6]. Each radar sensor is characterized by either transmitting capability or receiving capability. Moreover, the receiving sensors can receive and separate all the transmitting signals. Therefore, any transmitting sensor and any receiving sensor can form a transmitting/receiving channel (or can be considered to form an equivalent self-transmitting and self-receiving sensor). With an appropriate formation of the radar sensors, the target can be observed from multiple observation angles, which provides the ways to overcome the RCS fluctuation.
Recently, the research on the distributed ISAR is making rapid growth. Pastina et al. [3][4][5][6] analyzed the potential of the distributed ISAR to increase the cross-range resolution by exploiting multiple equivalent sensors to increase the global variation of the view angle. However, its effectiveness relies on the assumption that the change of the observation angle is small so that the range-compressed echoes of different channels have stable phases. However, in order to get the target image that overcomes the RCS fluctuation, the observation angles must be very different and the phase stability cannot be guaranteed. Thus, the proposed methods in [3][4][5] are not easy to apply in this scenario. Furthermore, for the high-maneuvering target which is inevitable in real scenarios, there is a nonuniform rotation after translation motion compensation. As a result, the change of observation angle of each equivalent sensor is no longer a linear function with respect to the slow time, which means that the combined view angle is not continuous and cannot provide better resolution. Thus, the methods proposed in [3][4][5] have some limitation for radar formations involving large variation among the individual perspectives and high-maneuvering targets.
Image fusion using subimages obtained from the distributed ISAR system is a solution to overcome the RCS flucation problem of the high maneuvering target. The distributed ISAR system enables the target to be observed from multiple perspectives, and the image fusion instead of echo fusion can reduce the limitation of the change of observation angle. References [7,8] acquire the fusion image after estimating the rotation rate and the bistatic angle from two subimages. Without more observation channels and more RCS information, the improvement of the image quality of the bistatic ISAR is limited. Reference [9] proposed a image fusion method for the uniform rotating target via distributed ISAR, but without consideration for the nonuniform rotating target.
In this work, the authors are interested in the potential of the distributed ISAR to acquire the image of the high-maneuvering target. Thus, the two aforementioned key problems must be considered, namely RCS fluctuation in different observation angles and target nonuniform rotation after translation compensation. For the first problem, the distributed ISAR technique is applied to obtain subimages from different observation angles with different RCS values. Then, the fusion image which contains all the RCS information and with improved quality can be acquired by the non-coherent fusion method.
For the second problem, it is proved that the radar echos can be modelled as a polynomial phase signal (PPS), and the Fourier transform is inappropriate for the azimuth focusing. To solve this problem, the Range-Instantaneous-Doppler (RID) algorithm has been proposed to improve the ISAR image quality, where the Fourier transform is substituted by the time-frequency representations (TFRs). The first approach of RID algorithm is based on the (TFRs) with high concentration and reduced cross-terms, such as the imaging methods based on Short Time Fourier Transform (STFT) [10], the Wigner-Ville distribution (WVD) [11] and so on. Though the RID algorithms performs well in terms of the computational efficiency, they still suffer from the tradeoff between the time frequency concentration and the cross-terms. Also, the obtained RID images are not very stable which will cause diffculty to the subimage fusion step. The second approach is based on the parameters estimation technique, the high quality ISAR images can be obtained by the parameters estimation of the coefficients of the PPS. The parameterized imaging methods (such as [12,13]) are effective for the enhancement of ISAR image. However, for the distributed ISAR system, it is very computational expensive to estimate the signal parameters of all range bins in echoes from all observation channels. However, the rotational parameters of the target can be regarded as invariable in the whole observation time, and the ratio of the rotational acceleration to rotational speed is also fixed. Therefore, estimating the ratio of rotational parameters in one or several range bins, rather than estimating signal parameters in all range bins, may greatly reduce the amount of computation. Wu [14] describes the rotational nonuniformity-relative angular acceleration (RAA) and relative angular jerk (RAJ), and with the estimated RAA and RAJ, rotational nonuniformity compensation is carried out. This method is effective to get the high quality ISAR image, but for the distributed ISAR, it is still cumbersome to construct compensation matrix for each equivalent sensor imaging. Thus, the MFT imaging method is used here. By estimating the rotation parameters, the ISAR images of all observation channels can be obtained directly through MFT. This paper is organized as follows. After presenting the signal model of distributed ISAR in Section 2, the imaging method of the nonuniform rotation target based on MFT is introduced in Section 3. Then, both simulation and experimental results are presented to validate the effectiveness of the proposed method in Section 4. Finally, we conclude this paper in Section 5.

Distributed ISAR Echo Model
Consider a 3D coordinate system (O, X, Y, Z) with the origin in the target's fulcrum, and the target's motion can be decomposed into a translation of the fulcrum and a rotation of the target body. We assume here that any relative translation motion between the distributed sensors and the target fulcrum have been already compensated. Therefore, we can focus on the target rotation. To simply the processing algorithm, the dominant rotation around the vertical axis is considered only. Also, we model the target as a rigid body consisting of Q scatterers. The radar formation and the target rotation are shown in Figure 1. The distributed ISAR system consists of M transmitting sensors and N receiving sensors as shown in Figure 1. m is used to represent the mth transmitting sensor and n is used to represent the nth receiving sensor. Each of them can placed in the ground and carry an antenna appropriately steered toward the moving target in the air. Also, another possible application is that each radar sensor is carried by an aircraft to observe the target on the ground. The detailed placement of the sensors will be introduced later.
Assume the M transmitting signals s m (t) are orthogonal, and each receiving sensor has the ability to receive and separate the signals from different transmitting sensors. Thus, I = MN transmitting/receiving channels can be formed. Meanwhile, we assume that all sensors have achieved time synchronization to ensure accurate matching of the transmitting signal and receiving signal. After demodulation and range compression, the received backscattered signal in the m-nth observation channel is denoted as wheret is the fast time, t p is the slow time, c is the wave velocity, λ is the carrier wavelength, σ mn,q is the scattering coefficient of the qth scatterer in the m-nth observation channel, and R mn,q (t p ) is the propagation distance of the signal in the m-nth channel of the qth scatterer. p m (t) denotes the point spread function of s m (t). ∀m ∈ {1, 2, · · · , M}, p m (t) ≈ p(t) where p(t) is the sinc function [15]. The position vector of the qth scatterer can be written as where r q is the distance of the scatterer from the fulcrum O, θ 0 q is the initial azimuth angle, κ q0 is the elevation angle of the qth scatterer above the XOY plane, and φ(t p ) is the rotation angle at slow time t p measured in clockwise. So, θ 0 q + φ(t p ) is the real azimuth angle at t p . The position unit vectors of the mth transmitting sensor and the nth receiving sensor are denoted as R 0m and R 0n : where R 0m and R 0n are the distance between the fulcurm and the transmitting or receiving sensor, ζ n (ζ m ) is the angle between the projection of the line connecting the receiving sensor (the transmitting sensor) and the target fulcrum on the XOY plane and the positive direction of Y-axis, which is drawn in Figure 1. ψ m and ψ n are grazing angle of the transmitting and receiving sensors. For the sake of simplicity, the grazing angles of all sensors are assumed to be the same as ψ 0 , which is reasonable when the sensors are not very far from each other and the target locates far away. Therefore, under the far-field assumption, the propagation distance R mn,q (t p ) can be expressed approximately as = 2(R mn +r q cos κ q0 cos ψ 0 cos(θ 0 q +φ(t p )−α mn ) cos β mn + sin ψ 0 sin κ q0 ), where R mn = (R 0m + R 0n )/2, α mn = (ζ m + ζ n )/2, and β mn = (ζ m − ζ n )/2 are the mean distance, mean angle, and the half difference angle of the m-nth transmitting/receiving channel pair, respectively [4]. The transmitting/receiving pair (m, n) can be regarded as the ith equivalent sensor. By setting α i = α mn and β i = β mn and neglecting the constant distance R mn under the assumption that the translation motion has been compensated, R mn,q (t p ) can be written as Therefore, the received signal of the ith equivalent sensor can be expressed as Due to the spatial separation of each equivalent sensor, the imaging projection plane (IPP) may not be the same. Since each equivalent sensor can be regarded as working independently, and the analysis of them are similar, we only use the analysis of the ith equivalent sensor to illustrate. As shown in Figure 2 , since the target rotates around the Z-axis, the rotation vector ω is very simple as ω = |ω| [0, 0, −1]. R i represents the range unit vector of the ith IPP and is pointed from the fulcrum O to the ith equivalent sensor.
In the ISAR imaging, the azimuth of the IPP is the cross-product of the range unit vector and the effective rotation vector ω ei , so the red line in Figure 2 represents the direction of azimuth. The ith IPP is the plane containing R i and When the IPPs are not the same, neither are the subimages of different observation angles. In order to make all the subimages in the same plane, it needs to project them into the unified IPP [6]. Therefore, we assume that ψ 0 = 0, according to the rotation vector ω, the IPP is the XOY plane in Figure 1. However, we must declare that in the real scene, the radar sensors should be placed reasonably so that the IPPs can be approximated as the same plane. Once the radar formation does not meet the requirement, the method in this paper need some modification.
Based on the above assumption, the imaging of the target is the projection in that plane. In this paper, the IPP is the XOY plane in Figure 1. Although any value of α i and β i can be set at this case, in reality, we still have to rationalize the radar and limit the value of ψ 0 , α i and β i . At the same time, assuming that the target has been projected onto the XOY plane, then let κ q0 = 0, so r iq can be re-expressed as Equations (7) and (8) are the new echo expression, which will be analyzed to get the target's image. For the distributed ISAR system, the observation time is short and φ(t p ) is small. By using approximations cos φ(t p ) ≈ 1 and sin φ(t p ) ≈ 0 in the range dimension, which are reasonable as the range error caused by these approximations is negigible compared with the range resolution, the compressed range of the qth scatterer can be expressed as While in the cross-range dimension, we use more accurate approximations sin(φ(t p )) ≈ φ(t p ), cos(φ(t p )) ≈ 1. Based on these approximations, S i (t, t p ) can be rewritten as In Equation (10), the first exponent term is a constant related to the scattering point and transmitting/receiving channel, which has no effect on the imaging result. The cross-range imaging information is contained in the second exponential phase. For nonuniform rotation target, φ(t p ) is a polynomial function and S i (t, t p ) is a polynomial phase signal (PPS). If the cross-range compression is achieved by using the FT, the second order terms (corresponding to the rotational acceleration) or even higher order terms (corresponding to the high order rotational motion) of φ(t p ) will make the image blurred in the cross-range dimension. To avoid this effect, instead of using FT-based imaging method, the MFT is used here for cross-range compression. The MFT is a generalization of the FT and can effectively deal with the PPS.

Single ISAR Imaging by Matching Fourier Transform
The MFT can focus signals with nonlinear phase changes [16,17]. Consider a continuous signal where ω is the MFT frequency. As aforementioned, for a nonuniform rotation target, φ(t p ) is a polynomial function. We use ω 0 , ω 1 , · · · to represent the different order components of the rotational angular velocity, ω 0 represents the uniform rotational component, ω 1 represents the first-order rotational acceleration and so on. Therefore, the rotation angle φ(t p ) can be expanded as where e d−1 = w d−1 /w 0 (e 0 = 1). ϑ(t p ) is determined by the target rotation characteristics and is the same for all scatterers. Substituting (12) into (10) and neglecting the first constant term in (10), the cross-range signal in one range bin can be rewritten as and is the equivalent cross-range position. It is clear that (13) is the sum of plurality of signals with the same frequency modulation function ϑ(t p ) and ϑ(t p ) is the monotonic bounded function with ϑ(0) = 0, which meets the definition of f (t). Thus, with the estimation of e d−1 , the MFT can be applied to S i (t p ). T is defined as the observation time, so the MFT expression is: Using the linear properties of the MFT and letting f d = ω/2π, we can get S i ( f d ) is a set of narrow sinc pulses after the MFT, and the equivalent cross-range position of the scattering point can be calculated from the peak position f dq of the pulse with the formula x iq = f dq λ/2ω 0 . Thus, the scatters are all focused both in range and cross range dimension. The width of the sinc pulse is 1/ϑ(T), and the equivalent cross-range resolution is ρ e = λ/2ω 0 ϑ(T) = λ/2φ(T). Consider two scatterers in the same range bin, i.e., their coordinate difference is (∆x, 0). By expanding x iq as (x q0 cos α i − y q0 sin α i ) cos β i ( where x q0 = r q sin θ 0 q , y q0 = r q cos θ 0 q ), it is easy to get that the distance of these two scatterers in the image is ∆x cos α i cos β i . Thus, in order to resolve these two scatterers, the following condition should be satisfied Obviously, after obtaining the estimation of the ratio e d−1 of the target rotation parameters, MFT is very succinct and can directly obtain the azimuthal focused results for all range bins. This article assumes that e d−1 has been estimated by other estimation methods. When the target's maneuvering is not severe, the PPS model degrades into an LFM signal. Accoding to previous studies, fractional Fourier transform (FrFT) [18], Radon-Wigner transform [19], adaptive Chirplet decomposition [20] and centroid frequency-chirp rate distribution(CFCRD) [21] are effective approaches for parameter estimation of LFM signals. When the target's maneuveirng is severe, several algorithms for cubic coefficient estimation can be retrieved, such as higher-order ambiguity function-integrated cubic phase function [22], scaled Fourier transform [23], keystone time-chirp rate distribution (KTCRD) [24]. Meanwhile, the rotational parameters can also be estimated according to the echoes from multi-channels [1,25,26].

Subimage Fusion
By applying the range compression and the MFT to the echoes of all transmitting/receiving channels, all subimages can be obtained. Due to the different values of RCS in different transmitting/receiving channels, the quality of these subimages has large fluctuations. In some subimages, the scatterers can be clearly distinguished. While in other subimages, some scatterers are submerged in noise. To obtain a stable imaging result, accumulation of all subimages is performed in the image domain.

Subimage Resampling
For subsequent image processing, the first operation is to equalize the sampling grid in the range and cross-range dimension. Denote the sampling grid of the original image as (∆l s , ∆l a ), where ∆l s and ∆l a are the range sampling interval and the cross-range sampling interval: where f s is the sampling frequency in range dimension, PRF is the pulse repetition frequency, and N f is the number of points of the MFT. Denote the sampling grid of the image after resampling as (∆l s , ∆l a ). According to the relation between ∆l s and ∆l a , in order to make ∆l s = ∆l a , the original image can be up-sampled or down sampled in the range dimension or cross-range dimension.

Subimage Scaling and Rotation
Use y iq to reformulate the range compression position r iq of the qth scatterer in the ith equivalent sensor, according to the expression of r iq , y iq can be expressed as After the MFT, the cross-range compression position can be expressed as Therefore, the imaging position of a scatterer (x q0 , y q0 ) in the ith subimage is (x iq , y iq ). (x iq , y iq ) is related to α i and β i , which are different from a subimage to another subimage. To combine all subimages, the influence of α i , β i must be removed so that the imaging positions of the same scatterer in different subimages are aligned. By examing (20) and (21), the relationship between (x q0 , y q0 ) and (x iq , y iq ) can be represented by the rotation and scaling of the coordinate system as shown in Figure 3, which can be expressed as Thus, the alignment of all subimages can be achieved by applying an inverse scaling transform and an inverse rotation transform. The image transform is carried out in the homogeneous coordinate system. Assume the subimage I i has G pix pixels where G pix = L × K (L, K means the sampling length in range and cross-range dimension). The homogeneous coordinate of the (l, k) pixel is represented as i ilk = [l, k, 1] T . By stacking the column vectors {i ilk }, l = 1, 2, · · · , L, k = 1, 2, · · · , K in the row dimension, the coordinates of I i is converted into a 3 × G pix dimensional matrix. Therefore, the corrected subimage matrix I i can be expressed as I i = T r T s I i [27], where the scaling matrix T s and rotation matrix T r are expressed as After scaling and rotation, the coordinates of some pixel points of the subimage may not be integers, and the values are not defined for such coordinates which should be estimated from their neighbors. To achieve this, the bilinear interpolation algorithm for image is adopted here.

Subimage Fusion
After the above processing, each scattering point is located at the same pixel in all subimages. Therefore, the final image can be obtained by the accumulation of all the subimages. The weight coefficients of subimages are calculated adaptively according to the image qualities of different subimages. Herein, the entropy is used as a quality metric.
The entropy of subimage I i [28][29][30] is defined as where is the value of the pixel (l, k). The image entropy reflects the sharpness of the image, and the image with a small entropy value is clearer. Therefore, the image with a smaller entropy value is given more weight coefficients. The final fusion image can be accumulated by Figure 4 shows the complete processing chain. Assume the translation motion of the target has been compensated, and after range compression and neglecting the range migration, the MFT is used for each equivalent channel's echo to obtain the subimage. After the preprocessing, scaling and rotating of all subimages, the fusion image can be acquired by accumulating of all the subimages with coefficients determined by the subimages entropies.

Simulation and Experimental Results
In this section, simulations are conducted to demonstrate the effectiveness of the proposed method firstly. The distributed ISAR system are composed of four transmitting sensors and four receiving sensors (i.e., M = 4, N = 4). Therefore, sixteen equivalent sensors can be obtained.

Target Model and Echo Analysis
The transmitting signal is a set of orthogonal signals with the same center frequency (10 GHz) and the same bandwidth (300 MHz) which can achieve a range resolution of 0.5 m. Assume the target rotates nonuniformly with an angular speed of 0.01 rad/s and an acceleration rate of 0.04 rad/s 2 , so we have e 1 is 4. The accumulation time is 0.5 s to achieve a cross-range resolution of 1.5 m. The target is composed of eleven scatterers (see Figure 5a) which are isotropic and independent with each other and obey the Swerling I distribution model. Under this model, each scatterer has the same RCS value in the same observation channel, but it is distributed identically across all channels. Setting the noise power to ς 2 n , for a given signal to noise ratio (SNR), the signal power ς 2 s can be calculated by ς 2 s = ς 2 n · 10 SNR/10 . I random numbers σ i that satisfy the Rayleigh distribution with ς s as the standard deviation parameter are generated. Then the RCS value of the ith observation channel is set to σ i [9]. Figure 5b shows the range compressed echoes of all sixteen equivalent senors. Each equivalent sensor has 101 slow time sampling points, and there are 16 pieces of echo data. According to (20), the range compression position of each scatterer is related to α i and β i . Since the α i and β i are different from each other, the cross-range signal of each scatterer point spreads over different range bins and there are range jump between different equivalent sensors. Therefore, the observation angle in one range bin cannot be increased by directly splicing the cross-range signal from different equivalent sensors.
As the RCS fluctuation of high-maneuvering target and the nonuniform rotation are the main problem that affects the image quality, we mainly consider overcoming these problems by imaging the nonuniform rotating target from multi-channels.

Comparision of FT and MFT
Assume we have already estimated the ratio of rotation parameters by parametric estimation methods [18][19][20][21][22][23][24], or we can estimate the rotation parameters w 0 , w 1 [1,25,26] to calculate e 1 . Figure 6a,b show the subimages (without noise) of the twelfth equivalent channel with mean angle α i = 0 o and half difference angle β i = 18 o , obtained from the FT and the MFT, respectively. Figure 6c,d are the cross-range profiles of the range bin of (a) and (b) at range 0 m. There are 5 scatterers in the range bin. It is evident that the MFT has better performance than the FT when dealing with nonuniform rotation target. Without a prior knowledge, it is impossible to recognize the 5 scatterers in Figure 6a. Compared with Figure 6b, the scatterers in Figure 6a have wider mainlobe and the power of their sidelobes are comparable with that of the mainlobes. Thus, the resolution of FT in this case is limited and it is easy subject to the sidelobe interference. Conversely, the cross-range profile obtained by the MFT is easily recognized and has lower sidelobe (about −25 dB).
What needs to be pointed out here is that when the estimation of ϑ(t p ) in (13) is inaccurate, there will be errors in the MFT results. However, as long as the error of ϑ(t p ) is not too large, acceptable results can be still obtained by applying MFT.

Subimage Fuison
After applying MFT to range compressed echoes of all observation channels, the target images from different observation angles can be obtained. Figure 7 are the subimages from all the observation angles (SNR = 8 dB, ς 2 n = 4), and each subimage has a different degree of scaling and rotation. Due to the RCS fluctuation, in some subimages (like (h),(i),(j)), the scattering points can be easily identified, while this is not true for other subimages (like (b),(k),(l)). The image quality of the corrected subimages can be judged directly from each subimage, and they are not drawn here. Figure 8 is the fusion image, in which the scattering points can be more clearly recognized. Table 1 shows the entropies of all subimages and the final fusion image. We can see that the final fusion image has the smallest entropy and the best image quality.
To further analyze the performance of image fusion, the cross-range profiles of the range bin at range 0 m are drawn. Three corrected subimages are selected.  Figure 9d, it can be seen that the noise power has been reduced significantly and the positions of scatterers are clearly distinguishable. Although the non-coherent fusion of the subimages does not improve the image resolution, the image quality after the fusion is indeed improved.

Application to Live Data
The aforementioned technique has been applied to the measured data of the conical target in the microwave anechoic chamber. The target is placed in the center of the turntable and remain motionless. The radar rotates around the target with the turntable and emits 8 GHz to 10 GHz stepped frequency modulated (SFM) signal with step of 20 MHz. The size and the shape of the target are shown in Figure 10 and the initial radar LOS is indicated by the dotted line. The PRF is 400 Hz. We simulate multiple equivalent sensors simultaneously observe the target by selecting radar echoes from different observation angles. In addition, we non-uniformly select the slow time sampling points to simulate nonuniform rotation of the target. Meanwhile, we add different levels of noise to the echoes of different observation channels to simulate the fluctuation of the background noise.   Figure 11a-e show the target with different rotation angles (since the half-difference angle of each observation angle is zero in this experiment, there is no scaling). Meanwhile, it also shows that in some observation angles, due to the mutual obstructions between scatterers, not all scatterers can be identified. Figure 11f is the non-coherent fused image, according to Table 2, the fused image has the smallest entropy and makes the shape of the target clearly visible from noise.

Conclusions
The distributed ISAR technique can utilize the data acquired from multiple observation angles. In this paper, a method which combines the distributed ISAR technique and the MFT is proposed to obtain the image of a high-maneuvering target. Two main problems, i.e., the RCS fluctuation and nonuniform rotation of high-maneuvering target, are solved.
In this paper, we assume that all IPPs are the same, which simplifies the processing chain. It should be pointed out that in order to satisfy this requirement, the radar sensors should not be placed far apart from each other in real scene. Based on the assumption, the multiple channel echoes of nonuniform rotation target are acquired from different observation angles firstly. Secondly, using the MFT for all channel echoes can avoid the azimuthal defocusing problem caused by the FT and get well-focused subimages. At the same time, by estimating the rotation parameters once, the MFT can apply directly to all the range bins of all echoes from different observation channel, which is computational effective. Thirdly, after the processing and accumulating of all subimages, the final fusion image can be acquired. To reduce the influence of RCS fluctuations, the accumulation coefficients are determined adaptively according to the subimage entropies. The simulations and experiment show the effectiveness of the proposed method to overcome the RCS fluctuation and can result in a final image with improved quality.
Thus, the innovativeness and contribution of this article are: • Based on the characteristics of the high-maneuvering target, the distributed ISAR technique is used to observe the target from multi-channels with different RCS values.

•
The MFT is applied to the echo of each channel to acquire well-focused subimages, which is computational efficient compared with other imaging methods.

•
Subimages fusion with adaptive coefficients calculated according to subimage entropies can effectively overcome the RCS fluctuations.
In this paper, the rotation around the Z-axis is considered only. However, in real scene, the target always has three-dimensional rotation which is more complicated. In this case, obtaining stable and recognized ISAR image is challenge. Meanwhile, the IPPs are assumed to be the same in this paper.
If the radar formation is not strictly limited, this assumption does not hold. Therefore, the subimages from different observation channels cannot be accumulated directly. Further research is necessary to increase the applicability of the method in this paper. Another research area of interest is how to use the raw data of each channel to get a ISAR image with higher resolution when the target rotates nonuniformly.
Author Contributions: Y.L. and W.Z. conceived and designed the experiments. Y.L. wrote the paper. Y.F. provided valuable discussions, substantially help as well as giving important comments for improving the presentation.