[Gdal-dev] gdal_translate ArcGrid to ENVI error

[Shoaib Burq (VPAC)]

Hi

I'm using gdal-dev 1.2.5 and have a set of data that I batch convert from
ArcGrid to ENVI format. Most of the data works out OK.

However a couple of grids die near the end of the conversion with the 
error:

processing: ./data/geol_epp/hdr.adf ....
Input file size is 2544, 3419
0...10...20...30...40...50...60...70...80...90...ERROR 1: IReadBlock 
failed at X offset 6, Y offset 854
ERROR 1: GetBlockRef failed at X block offset 6, Y block offset 854
Open GDAL Datasets:
  1 N ENVI   2544x3419x1 ./ENVI_format/geol_epp.img

processing: ./data/new_lsys/hdr.adf ....
Input file size is 2544, 3419
0...10...20...30...40...50...60...70...80...90...ERROR 1: IReadBlock 
failed at X offset 7, Y offset 854
ERROR 1: GetBlockRef failed at X block offset 7, Y block offset 854
Open GDAL Datasets:
  1 N ENVI   2544x3419x1 ./ENVI_format/new_lsys.img


Had a look at the pix in question but don't appear to have anything 
unusual.

here is the output of gdalinfo on both the source images:

[sab at uga arc2envi]$ gdalinfo ./data/geol_epp/
Driver: AIG/Arc/Info Binary Grid
Size is 2544, 3419
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            TOWGS84[0,0,0,0,0,0,0],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9108"]],
        AXIS["Lat",NORTH],
        AXIS["Long",EAST],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",147],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",10000000],
    UNIT["METERS",1]]
Origin = (255994.767802,5919365.411465)
Pixel Size = (20.00000000,-20.00000000)
Corner Coordinates:
Upper Left  (  255994.768, 5919365.411) (  0d 0'0.00"E, 90d 0'0.00"S)
Lower Left  (  255994.768, 5850985.411) (  0d 0'0.00"E, 90d 0'0.00"S)
Upper Right (  306874.768, 5919365.411) (  0d 0'0.00"E, 90d 0'0.00"S)
Lower Right (  306874.768, 5850985.411) (  0d 0'0.00"E, 90d 0'0.00"S)
Center      (  281434.768, 5885175.411) (  0d 0'0.00"E, 90d 0'0.00"S)
Band 1 Block=256x4 Type=Int32, ColorInterp=Undefined
  Min=2004.000 Max=99999.000
  NoData Value=-2.14748e+09

[sab at uga arc2envi]$ gdalinfo ./data/new_lsys/
Driver: AIG/Arc/Info Binary Grid
Size is 2544, 3419
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            TOWGS84[0,0,0,0,0,0,0],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9108"]],
        AXIS["Lat",NORTH],
        AXIS["Long",EAST],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",147],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",10000000],
    UNIT["METERS",1]]
Origin = (255994.767802,5919365.411465)
Pixel Size = (20.00000000,-20.00000000)
Corner Coordinates:
Upper Left  (  255994.768, 5919365.411) (  0d 0'0.00"E, 90d 0'0.00"S)
Lower Left  (  255994.768, 5850985.411) (  0d 0'0.00"E, 90d 0'0.00"S)
Upper Right (  306874.768, 5919365.411) (  0d 0'0.00"E, 90d 0'0.00"S)
Lower Right (  306874.768, 5850985.411) (  0d 0'0.00"E, 90d 0'0.00"S)
Center      (  281434.768, 5885175.411) (  0d 0'0.00"E, 90d 0'0.00"S)
Band 1 Block=256x4 Type=Byte, ColorInterp=Undefined
  Min=1.000 Max=60.000
  NoData Value=255

Again nothing too strange, Right?.... When the tremendously-able
gdal_translate finishes it does create an output image although
incomplete. here is the output of gdalinfo on the ENVI images. Note the
difference in the Center coordinate :

[sab at uga arc2envi]$ gdalinfo ./ENVI_format/geol_epp.img
Driver: ENVI/ENVI .hdr Labelled
Size is 2544, 3419
Coordinate System is:
LOCAL_CS["arbitrary",
    UNIT["Meter",1]]
Origin = (255994.767802,5919365.411470)
Pixel Size = (20.00000000,20.00000000)
Corner Coordinates:
Upper Left  (  255994.768, 5919365.411)
Lower Left  (  255994.768, 5987745.411)
Upper Right (  306874.768, 5919365.411)
Lower Right (  306874.768, 5987745.411)
Center      (  281434.768, 5953555.411)
Band 1 Block=2544x1 Type=Int32, ColorInterp=Undefined

[sab at uga arc2envi]$ gdalinfo ./ENVI_format/new_lsys.img
Driver: ENVI/ENVI .hdr Labelled
Size is 2544, 3419
Coordinate System is:
LOCAL_CS["arbitrary",
    UNIT["Meter",1]]
Origin = (255994.767802,5919365.411470)
Pixel Size = (20.00000000,20.00000000)
Corner Coordinates:
Upper Left  (  255994.768, 5919365.411)
Lower Left  (  255994.768, 5987745.411)
Upper Right (  306874.768, 5919365.411)
Lower Right (  306874.768, 5987745.411)
Center      (  281434.768, 5953555.411)
Band 1 Block=2544x1 Type=Byte, ColorInterp=Undefined

Any guidence is appreciated. Also quite happy to post the rasters to look
at.

Thanks heap
regards
shoaib