gdal.warp clipping the extent of a raster

[mluck]

When replacing rasterio with gdal, gdal.warp() seems to be producing a raster that is clipped on one of the edges. I've tried a number of parameters including outputBounds with no success.

For example, comparing the previous BLE benchmark inundation raster for 12040101 (left) with the new raw BLE depth raster (middle) and the GDAL transformed (using gdal.warp() overlaid on the first image) shows that the bottom of the grayscale image is missing ((right).

This issue might be ignored. It will mean that create_ble_benchmark.py will import both rasterio and gdal, but there do not appear to be any conflicts and this is a seldom-used script that doesn't impact FIM pipeline at all.