FastSAR: High-Throughput, High-Fidelity Open-Source SAR Image Formation on CPUs, GPUs and TPUs
Abstract: Exact backprojection of a spaceborne SAR collection needs one interpolated read per pixel for every pulse; for the 12,207 by 8,808 pixel Umbra image of the Panama Canal this balloons to 1.6 trillion reads. FastSAR is an open-source library for SAR image formation on CPUs, NVIDIA GPUs and Google TPUs, with factorized backprojection, exact backprojection and polar format behind one Python interface and readers for the vendors' phase-history files. Its factorized kernels form the Panama image in 11.2 s on a 16-vCPU instance, 3.0 s on an NVIDIA L4, 0.6 s on an H100 and 2.9 s on a TPU v6e, with errors of -56 to -63 dB relative to a float64 exact backprojection in three regions over six cloud devices; the L4, a modest GPU, gives the lowest cost of the six at 0.59 US dollars per thousand images. At exact backprojection FastSAR is as fast as the MatX backprojection example, by far the fastest open-source implementation found (10.9 against 11.6 s at the same accuracy, with MatX's pixel coordinates in float64); its factorized algorithm is 3.8 times faster than MatX, and its polar format forms the image in 1.4 s on the L4 at -32 dB. On five Capella collections in three modes the worst-region error is below -59 dB.