Two-dimensional signal processing from hexagonal rasters

Two-dimensional signals are normally sampled and processed as rectangular arrays. For signals which are bandlimited with a circular region of support in the Fourier plane, however, it has been known for some time that a savings in the number of samples required for an exact reconstruction can be realized by sampling the signal on a hexagonal sampling raster. In this presentation we show that a substantial savings in computation can result as well. Included are methods for signal representation, linear system implementation, Fourier transform computation and FIR filter design. Some comparisons between systems defined over rectangular and hexagonal sampling rasters will also be given.