[GRASS-user] Importing ESRI .hdr Labelled raster via r.in.gdal
displays strange aliasing pattern
Patton, Eric
epatton at nrcan.gc.ca
Fri May 11 13:35:28 EDT 2007
I'm importing an Arc Floating-point raster using the 'EHdr/ESRI .hdr Labelled' driver in gdal:
$ gdalinfo 45_04_N_66_44_W_BATHY.flt
Driver: EHdr/ESRI .hdr Labelled
Size is 2641, 3721
Coordinate System is `'
Origin = (2481609.500000000000000,7340725.500000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Corner Coordinates:
Upper Left ( 2481609.500, 7340725.500)
Lower Left ( 2481609.500, 7337004.500)
Upper Right ( 2484250.500, 7340725.500)
Lower Right ( 2484250.500, 7337004.500)
Center ( 2482930.000, 7338865.000)
Band 1 Block=2641x1 Type=Byte, ColorInterp=Undefined
NoData Value=-99999.999
The projection is New Brunswick Double Stereographic (EPSG #2953). I have set up a Location with these projection parameters.
Th import seems to progress fine with no errors, but the imported raster doesn't display any data, just a jagged, greyscale grid. (A screenshot of the imported raster was too big to attach, but I can provide it in a separate email attachment if anyone is interested.)
I suspect, but am not sure, that I need to use gdal_translate instead and pass some CREATOPT parameters?
Also noteworthy is the fact that I had to do a tiny hack to the hdr file as per Frank Warmerdam's suggestion through the GDAL bugtracker:
(Background: I was received a malloc error when trying to import the ESRI Floating-point raster. See http://trac.osgeo.org/gdal/ticket/1573)
>From Frank Warmerdam:
"The problem is related to initialization when the NBITS keyword is missing
from the .hdr file. In 1.4.0 a bunch of new code was introduced for NBITS
< 8, but it also added a dependency on having the keyword present.
The workaround is to add "NBITS 8" in the .hdr file.
Fix is to default to 8 in the code. The fix has been committed in trunk,
and 1.4 branch and it should be released in 1.4.2."
I have added the above fix to the hdr file prior to using r.in.gdal to import the raster.
Any help appreciated!
~ Eric.
More information about the grass-user
mailing list