footprint function

Compute footprint of a raster