SAR images are improved by a method for acquiring a synthetic aperture image from a sequence of periodic pulse returns where the sequence of periodic pulse returns is interspersed with interrupts, i.e. missing pulses. The interrupts mark the start and end of one or more segments, where the segments contain the periodic pulse returns form the SAR image. The method comprises the steps of: converting said pulse returns into a digital stream; performing an azimuth deskew on said digital stream to obtain a deskewed digital stream; forming a forward-backward data matrix from the deskewed digital stream for one or more segments; forming an average segment covariance from the forward-backward data matrix; computing a model order for the average segment covariance; computing one or more linear prediction coefficients using data contained in the forward backward data matrix, and model order; using the linear prediction coefficients to compute missing pulse returns belonging within the interrupts. The computation for extrapolating the missing pulse returns is introduced after the Stolt interpolator in RMA processing. In computing the model order, eigenvalues are found and compared to a threshold. Roots of a linear prediction polynomial are computed, then stabilized to obtain stabilized roots. Linear prediction coefficients are reconstituted using the stabilized roots. Sub-bands are used to decrease computing time for the missing pulse returns.
A moving radar (405) generates a synthetic aperture image from an incomplete sequence of periodic pulse returns. The incomplete sequence of periodic pulse returns has one or more missing pulses. The radar converts the incomplete sequence of pulse returns into a digital stream. A computer (403) processes the digital stream by computing an along track Fourier transform (402), a range compression (408), an azimuth deskew (410) and an image restoration and auto focus (412). The image restoration and autofocus (412) utilizes a low order autofocus (501), a gap interpolation using a Burg algorithm (503), and a high order autofocus (505) for generating an interpolated sequence. The interpolated sequence contains a complete sequence of periodic pulse returns with uniform spacing for generating the synthetic aperture image. The image restoration and autofocus (412) computes a linear prediction coefficients estimate using the Burg Algorithm (606). The linear prediction coefficients estimate (606) is used to compute a weighted forward-backward interpolation to generate the complete sequence of periodic pulse returns (608).
A moving radar generates a search mode synthetic aperture image of a patch from a sequence of periodic pulse returns having one or more missing pulses. An azimuth and range interpolation generates an interpolated sequence having samples oriented in range and azimuth frequency with uniform spacing. Range compression is performed using an IFFT. Azimuth deskew, an autofocus and pulse restore generates a focused and restored sequence. Azimuth reskew, and gain phase equalization generates an equalized sequence. A first linear phase is summed to the equalized sequence for applying a fractional sample shift in range frequency. A range FFT and Along Track IFFT is further applied to obtain a domain changed sequence. A second linear phase is summed to the domain changed sequence. A CT FFT of the result generates an image of the patch. The azimuth interpolation and range interpolation also include a Stolt interpolation after a matched filter function.
