초록 ▼

SAR imaging method that includes applying PRF decimation to range-compressed IQ data to generate PRF-decimated range-compressed IQ data for each image block of an image and applying motion compensation to the PRF-decimated range-compressed IQ data to generate motion-compensated data for each image b

SAR imaging method that includes applying PRF decimation to range-compressed IQ data to generate PRF-decimated range-compressed IQ data for each image block of an image and applying motion compensation to the PRF-decimated range-compressed IQ data to generate motion-compensated data for each image block. The method includes computing first stage image values at image grid point intersections of iso-range lines and vertical grid lines for each image bock based on the motion-compensated data for each image block. The method also includes computing second stage image values for the image grid point intersections by interpolation using the first stage image values at the image grid point intersections and correcting image phase of the second stage image values for the image grid point intersections to generate phase-corrected image values for each image block. The method includes generating a full-resolution SAR image by summing the phase-corrected image values for each image block.

대표청구항 ▼

1. A synthetic aperture radar (SAR) imaging method for processing moving air borne radar acquired images for multiple performance matrices, the method comprising: physically acquiring SAR in-plane and quadratic data from a scan of an image gridapplying range pulse compression to SAR in-phase and qua

1. A synthetic aperture radar (SAR) imaging method for processing moving air borne radar acquired images for multiple performance matrices, the method comprising: physically acquiring SAR in-plane and quadratic data from a scan of an image gridapplying range pulse compression to SAR in-phase and quadrature (IQ) data for each image block of a SAR image grid to generate range-compressed IQ data for each image block;applying pulse range frequency (PRF) decimation to the range-compressed IQ data to generate PRF-decimated range-compressed IQ data for each image block;applying motion compensation to the PRF-decimated range-compressed IQ data to generate motion-compensated data for each image block;computing first stage image values at image grid point intersections of iso-range lines and vertical grid lines for each image block based on the motion-compensated data for each image block;computing second stage image values for the image grid point intersections by interpolation using the first stage image values at the image grid point intersections;correcting image phase of the second stage image values for the image grid point intersections to generate phase-corrected image values for each image block; andgenerating a full-resolution SAR image by summing the phase-corrected image values for each image block. 2. The method of claim 1, wherein the step of applying motion compensation comprises performing a fast Fourier transform in the SAR image grid range dimension. 3. The method of claim 2, wherein the fast Fourier transform in the SAR image grid range dimension is performed in accordance with: S′⁡(n,i)=∑k=0Nfft-1⁢S⁡(k,i)·ⅇ-j2π·n⁢⁢kNfft⁢⁢for⁢⁢i∈[1,2,3,…⁢,Np] where S(k,i) is a two-dimensional range compressed SAR response with k being range bin index and i being pulse index, Np is number of pulses in a processing pulse group, and Nfft is the least integer power of 2 greater than the number of slant range bins covering the processed image block. 4. The method of claim 3, comprising calculating S1(n,i)=S′(n,i)·Ref(n) by multiplying S′(n,i) with a reference function of the form: Ref(n)==exp(−i2π·n1·ΔR(i)/Δrs), where n=[1,2,3 . . . , Nfft], n1=[0,1,2 . . . , Nfft/2−1,−Nfft/2,−Nfft/2+1, . . . , −1], ΔR(i)=R(i)−R(ic), where R(i) is distance between a reference scatter and the SAR at pulse i which is within a range of [1,2,3, . . . , Np], and ic=Np/2+1. 5. The method of claim 4, wherein the motion-compensated data for each image block is determined by performing a fast Fourier transform in the range dimension in accordance with: Smocomp⁡(k,i)=∑n=0Nfft-1⁢S1⁡(n,i)·ⅇ-j2π·k⁢⁢nNfft⁢⁢for⁢⁢i∈[1,2,3,…⁢,Np]. 6. The method of claim 1, wherein the first stage image values are computed in accordance with: Sstage_⁢1⁡(Rj,fd⁡(Rj,n))=∑i=1Np⁢Smocomp⁡(i,Rj)·exp⁡(j2π·f⁢⁢d⁡(Rj,n)·i′) where i′=[−Np/2, . . . , −1,0,2, . . . Np/2−1], fd(Rj,n) is Doppler frequency at (Rj,n), Rj=∥ac(ic)−(b+n·az+β·rg)∥, b is position of each image block center, ac(i) is position of the SAR aircraft at the i-th pulse in the pulse group, {circumflex over (p)}az is azimuth unit vector of the image grid, {circumflex over (p)}rg is range unit vector of the image grid, n is azimuth coordinate index, ic=Np/2+1, Np is number of pulses in a processing pulse group, where n is the along-track pixel index and β is cross-track pixel index associated with the pixel having a slant range of Rj. 7. The method of claim 6, wherein the second stage image values are computed in accordance with: Sstage_⁢2⁡(m,n)=∑k=1K⁢Sstage_⁢1⁡(Rj,n)·wk⁡(Rj,m) where R(m,n) is distance between grid point (m,n) and ac(ic), and wk(Rj,m) are the interpolation coffiicients determined based on the range of the point (m,n) and Rj. 8. The method of claim 7, wherein the full-resolution SAR image is generated in accordance with: SF_res⁡(m,n)=∑l=1L⁢S⁡(m,n,l)·exp⁢{-j4π⁢⁢Ric⁡(l),m,nλ}, where Ric(l),m,n=∥ac(ic(l))−m,n∥, ic(l) is pulse index of the center pulse in group l, ac(ic(l)) is the aircraft position at the center of each pulse group, and Ric(l),m,n is the distance between ac(ic(l)) and the image pixel (m,n). 9. A synthetic aperture radar (SAR) imaging system for processing moving images acquired by moving airborne radar apparatus to provide images in multiple performance metrices, the system comprising: (a) radar apparatus for acquiring from SAR in-plane and quadrature data from a scan of an image;(b) one or more computer hardware processors; and(c) a memory, the memory including executable code representing instructions that when executed cause the system to: (i) apply range pulse compression to SAR in-phase and quadrature (IQ) data for each image block of a SAR image grid to generate range-compressed IQ data for each image block;(ii) apply pulse range frequency (PRF) decimation to the range-compressed IQ data to generate PRF-decimated range-compressed IQ data for each image block;(iii) apply motion compensation to the PRF-decimated range-compressed IQ data to generate motion-compensated data for each image block;(iv) compute first stage image values at image grid point intersections of iso-range lines and vertical grid lines for each image bock based on the motion-compensated data for each image block;(v) compute second stage image values for the image grid point intersections by interpolation using the first stage image values at the image grid point intersections;(vi) correct image phase of the second stage image values for the image grid point intersections to generate phase-corrected image values for each image block; and(vii) generate a full-resolution SAR image by summing the phase-corrected image values for each image block. 10. The system of claim 9, wherein the memory includes executable code representing instructions that when executed cause the system to apply motion compensation by performing a fast Fourier transform in the SAR image grid range dimension. 11. The system of claim 10, wherein the fast Fourier transform in the SAR image grid range dimension is performed in accordance with: S′⁡(n,i)=∑k=0Nfft-1⁢S⁡(k,i)·ⅇ-j2π·n⁢⁢kNfft⁢⁢for⁢⁢i∈[1,2,3,…⁢,Np] where S(k,i) is a two-dimensional range compressed SAR response with k being range bin index and i being pulse index, Np is number of pulses in a processing pulse group, and Nfft is the least integer power of 2 greater than the number of slant range bins covering the processed image block. 12. The system of claim 11, comprising wherein the memory includes executable code representing instructions that when executed cause the system to calculate S1(n,i)=S′(n,i)·Ref(n) by multiplying S′(n,i) with a reference function of the form: Ref(n)==exp(−i2π·n1ΔR(i)/Δrs), where n=[1,2,3 . . . , Nfft], n1=[0,1,2 . . . , Nfft/2−1,−Nfft/2,−Nfft/2+1, . . . , −1], ΔR(i)=R(i)−R(ic), where R(i) is distance between a reference scatter and the SAR at pulse i which is within a range of [1,2,3, . . . , Np], and ic=Np/2+1. 13. The system of claim 12, wherein the memory includes executable code representing instructions that when executed cause the system to determine the motion-compensated data for each image block by performing a fast Fourier transform in the range dimension in accordance with Smocomp⁡(k,i)=∑n=0Nfft-1⁢S1⁡(n,i)·ⅇ-j2π·k⁢⁢nNfft⁢⁢for⁢⁢i∈[1,2,3,…⁢,Np]. 14. The system of claim 13, wherein the memory includes executable code representing instructions that when executed cause the system to compute the first stage image values in accordance with: Sstage_⁢1⁡(Rj,fd⁡(Rj,n))=∑i=1Np⁢Smocomp⁡(i,Rj)·exp⁡(j2π·f⁢⁢d⁡(Rj,n)·i′) where i′=[−Np/2, . . . , −1,0,2, . . . Np/2−1], fd(Rj,n) is Doppler frequency at (Rj,n), Rj=∥ac(ic)−(b+n·{circumflex over (p)}az+β·{circumflex over (p)}rg)∥, b is position of each image block center, ac(i) is position of the SAR aircraft at the i-th pulse in the pulse group, {circumflex over (p)}az is azimuth unit vector of the image grid, {circumflex over (p)}rg is range unit vector of the image grid, n is azimuth coordinate index, ic=Np/2+1, Np is number of pulses in a processing pulse group, where α is along-track pixel index n and β is cross-track pixel index associated with the pixel having a slant range of Rj. 15. The system of claim 14, wherein the memory includes executable code representing instructions that when executed cause the system to compute the second stage image values in accordance with: Sstage_⁢2⁡(m,n)=∑k=1K⁢Sstage_⁢1⁡(Rj,n)·wk⁡(Rj,m) where R(m,n) is distance between grid point (m,n) and ac(ic), and wk(Rj,m) are the interpolation coeffiicents determined based on the range of the point (m,n) and Rj. 16. The system of claim 15, wherein the full-resolution SAR image is generated in accordance with: SF_res⁡(m,n)=∑l=1L⁢S⁡(m,n,l)·exp⁢{-j4π⁢⁢Ric⁡(l),m,nλ}, where Ric(l),m,n=∥ac(ic(l))−m,n∥, ic(l) is pulse index of the center pulse in group l, ac(ic(l)) is the aircraft position at the center of each pulse group, and Ric(l),m,n is the distance between ac(ic(l)) and the image pixel (m,n).