gdalwarp creates correct ENVI hdr but incomplete raster

[piyushrpt]

I'm trying to warp some swath data using geolocation arrays. All input data is stored in flat binary files with appropriate VRTs. The methodology works for most cases but there were a couple of instances where it resulted in incomplete raster files (missing lines at the bottom).

Command used:

gdalwarp -t_srs EPSG:4326 -of ENVI -ot Float32  -srcnodata 0 -te -122.32 37.57-121.102 38.617 -geoloc -tr 0.0001 0.0001 swathRaster.dat.vrt  mapRaster.dat

The resulting mapRaster.dat.hdr is correct with the expected dimensions, but the generated mapRaster.dat has missing lines. One thing to note is that the missing lines should contain nodata (0) and lie outside the support of the geoloc arrays.

I can overcome the issue in two ways:

1. Create a Geotiff and then translate it to an ENVI file. This worked fine.
2. Increase working memory using "-wm 2000" and this worked as well.

Looks like there could be issues with nodata regions and memory settings when using the ENVI driver to warp raster files. I'm also attaching two log files generated with debug turned on - one for default call using gdalwarp and one using more working memory if it helps.

[piyushrpt]

Log for gdalwarp with default settings

[piyushrpt]

Log for gdalwarp with extra working memory "-wm 2000"

[Even Rouault]

In 35523:

ENVI: on closing, pad image file with trailing nul bytes if needed (fixes #6662)

[Even Rouault]

In 35524:

ENVI: on closing, pad image file with trailing nul bytes if needed (fixes #6662)