[Gdal-dev] gdal_translate ArcGrid to ENVI error
Shoaib Burq (VPAC)
sab at vpac.org
Tue Apr 5 18:03:06 EDT 2005
ok the data is now sitting at
http://www.vpac.org/~sab/data/arc2envi_error.tar.gz
The other thing I forgot to mention is that I am able to convert this
ArcGrid without problems in ENVI (Win32) using the Save As functionality.
regards
shoaib
On Wed, 6 Apr 2005, Shoaib Burq (VPAC) wrote:
>
> 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
>
> _______________________________________________
> Gdal-dev mailing list
> Gdal-dev at xserve.flids.com
> http://xserve.flids.com/mailman/listinfo/gdal-dev
>
--
Shoaib Burq
--
VPAC - Geospatial Applications Developer
Building 91, 110 Victoria Street,
Carlton South, Vic 3053, Australia
_______________________________________________________________
w: www.vpac.org | e: sab_AT_vpac_DOT_org | mob: +61.431-850039
More information about the Gdal-dev
mailing list