Fast Gpu-Based Interpolation for SAR Backprojection

We introduce and discuss a parallel SAR backprojection algorithm using a Non-Uniform FFT (NUFFT) routine implemented on a GPU in CUDA language. The details of a convenient GPU implementation of the NUFFT-based SAR backprojection algorithm, amenable to further generalizations to a multi-GPU architecture, are also given. The performance of the approach is analyzed in terms of accuracy and computational speed by comparisons to a \standard", parallel version of the backprojection algorithm exploiting FFT+interpolation instead of the NUFFT. Difierent interpolators have been considered for the latter processing scheme. The NUFFT-based backprojection has proven signiflcantly more accurate than all the compared approach, with a computing time of the same order. An analysis of the computational burden of all the difierent steps involved in both the considered approaches (i.e., standard and NUFFT backprojections) has been also reported. Experimental results against the Air Force Research Laboratory (AFRL) airborne data delivered under the \challenge problem for SAR-based Ground Moving Target Identiflcation (GMTI) in urban environments" and collected over circular ∞ight paths are also shown.

[1]  Amedeo Capozzoli,et al.  GPU-BASED Ï-K TOMOGRAPHIC PROCESSING BY 1D NON-UNIFORM FFTS , 2012 .

[2]  A. F. Yegulalp.  Fast backprojection algorithm for synthetic aperture radar , 1999, Proceedings of the 1999 IEEE Radar Conference. Radar into the Next Millennium (Cat. No.99CH36249).

[3]  Maurizio Migliaccio,et al.  Knab Sampling Window for InSAR Data Interpolation , 2007, IEEE Geoscience and Remote Sensing Letters.

[4]  Achim Hein.  Processing of SAR Data: Fundamentals, Signal Processing, Interferometry , 2003 .

[5]  LeRoy A. Gorham,et al.  SAR image formation toolbox for MATLAB , 2010, Defense + Commercial Sensing.

[6]  Martin Blom,et al.  VHF SAR image formation implemented on a GPU , 2005, Proceedings. 2005 IEEE International Geoscience and Remote Sensing Symposium, 2005. IGARSS '05..

[7]  R.L. Moses,et al.  Polar format algorithm for bistatic SAR , 2004, IEEE Transactions on Aerospace and Electronic Systems.

[8]  R. Keys.  Cubic convolution interpolation for digital image processing , 1981 .

[9]  A. R. Brenner,et al.  PAMIR - a wideband phased array SAR/MTI system , 2003 .

[10]  Karsten Fourmont.  Non-Equispaced Fast Fourier Transforms with Applications to Tomography , 2003 .

[11]  Satoshi Matsuoka,et al.  Auto-tuning 3-D FFT library for CUDA GPUs , 2009, Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis.

[12]  Paul Suetens,et al.  Efficient GPU-Based Texture Interpolation using Uniform B-Splines , 2008, J. Graph. Tools.

[13]  V. Pascazio,et al.  WASAR: a wide-angle SAR processor , 1992 .

[14]  Richard Bamler,et al.  Evaluation of interpolation kernels for SAR interferometry , 1999, IEEE Trans. Geosci. Remote. Sens..

[15]  Fuk K. Li,et al.  Synthetic aperture radar interferometry , 2000, Proceedings of the IEEE.

[16]  LeRoy A. Gorham,et al.  A challenge problem for SAR-based GMTI in urban environments , 2009, Defense + Commercial Sensing.

[17]  Andreas Reigber,et al.  SAR tomography with optimized constellation and its application to forest scenes , 2009 .

[18]  A. Capozzoli,et al.  Fast CPU / GPU Pattern Evaluation of Irregular Arrays , 2010 .

[19]  Jørgen Dall,et al.  SAR focusing of P-band ice sounding data using back-projection , 2010, 2010 IEEE International Geoscience and Remote Sensing Symposium.

[20]  Maurizio Migliaccio,et al.  A new interpolation kernel for SAR interferometric registration , 2003, IEEE Trans. Geosci. Remote. Sens..

[21]  G. Schirinzi,et al.  'Traveling' sampling of scattered fields , 1985, 1985 Antennas and Propagation Society International Symposium.

[22]  Richard Bamler,et al.  A comparison of range-Doppler and wavenumber domain SAR focusing algorithms , 1992, IEEE Trans. Geosci. Remote. Sens..

[23]  P. C. Gao,et al.  Fast RCS prediction using multiresolution shooting and bouncing ray method on the GPU , 2010 .

[24]  Paolo Valerio Testa,et al.  NUFFT-based SAR backprojection on multiple GPUs , 2012, 2012 Tyrrhenian Workshop on Advances in Radar and Remote Sensing (TyWRRS).

[25]  Amedeo Capozzoli,et al.  NUFFT- & GPU-Based Fast Imaging of Vegetation , 2011, IEICE Trans. Commun..

[26]  Kan Xu,et al.  GPU Accelerated Unconditionally Stable Crank-Nicolson FDTD Method for the Analysis of Three-Dimensional Microwave Circuits , 2010 .

[27]  Hyeokho Choi,et al.  Direct-Fourier reconstruction in tomography and synthetic aperture radar , 1998, Int. J. Imaging Syst. Technol..

[28]  Yoram Bresler,et al.  O(N2log2N) filtered backprojection reconstruction algorithm for tomography , 2000, IEEE Trans. Image Process..

[29]  Riccardo Lanari,et al.  Synthetic Aperture Radar Processing with GPGPU , 2010, IEEE Signal Processing Magazine.

[30]  Jeffrey A. Fessler,et al.  Nonuniform fast Fourier transforms using min-max interpolation , 2003, IEEE Trans. Signal Process..

[31]  Lars M. H. Ulander,et al.  Synthetic-aperture radar processing using fast factorized back-projection , 2003 .

[32]  Juan M. Lopez-Sanchez,et al.  Efficient Interpolation of SAR Images for Coregistration in SAR Interferometry , 2007, IEEE Geoscience and Remote Sensing Letters.

[33]  Giorgio Franceschetti,et al.  A SAR processor based on two-dimensional FFT codes , 1990 .

[34]  A. Macovski,et al.  Selection of a convolution function for Fourier inversion using gridding [computerised tomography application]. , 1991, IEEE transactions on medical imaging.

[35]  Aining Li.  Algoritthms For The Implementation Of STOLT Interpolation In SAR Processing , 1992, [Proceedings] IGARSS '92 International Geoscience and Remote Sensing Symposium.

[36]  Veysel Demir,et al.  Graphics processor unit (GPU) acceleration of finite-difference frequency-domain (FDFD) method , 2012 .

[37]  N. Stevens.  Processing of sar data: fundamentals, signal processing, interferometry , 2004 .

[38]  Ahmed Fasih,et al.  GPU-accelerated synthetic aperture radar backprojection in CUDA , 2010, 2010 IEEE Radar Conference.

[39]  Emre Ertin,et al.  Sparse multipass 3D SAR imaging: applications to the GOTCHA data set , 2009, Defense + Commercial Sensing.

[40]  Marc Rodriguez-Cassola,et al.  Processing of Circular SAR trajectories with Fast Factorized Back-Projection , 2011, 2011 IEEE International Geoscience and Remote Sensing Symposium.

[41]  Alberto Moreira,et al.  A novel optimization approach to forest height reconstruction from multi-baseline data , 2007, 2007 IEEE International Geoscience and Remote Sensing Symposium.

[42]  F. Rocca,et al.  SAR data focusing using seismic migration techniques , 1991 .

[43]  Vito Pascazio,et al.  An efficient SAR parallel processor , 1991 .

[44]  Michal Mrozowski,et al.  A Memory Efficient and Fast Sparse Matrix Vector Product on a GPU , 2011 .

[45]  John J. Knab,et al.  Interpolation of band-limited functions using the approximate prolate series (Corresp.) , 1979, IEEE Trans. Inf. Theory.

[46]  H. Bao,et al.  FROM CPU TO GPU : GPU-BASED ELECTROMAGNETIC COMPUTING ( GPUECO ) , 2008 .

[47]  Michal Mrozowski,et al.  FINITE ELEMENT MATRIX GENERATION ON A GPU , 2012 .

[48]  John J. Knab,et al.  The sampling window , 1983, IEEE Trans. Inf. Theory.

[49]  W. Kenneth Jenkins,et al.  Convolution backprojection image reconstruction for spotlight mode synthetic aperture radar , 1992, IEEE Trans. Image Process..