On Thu, Jan 22, 2004 at 08:15:19PM -0500, SWlab wrote:
Hello list,
I have some problems with "r.in.gdal":
I run grass 5.0.3, compiled with the "--with-gdal" option, on a SuSE8.2 box.
The gdal version is 1.1.9, compiled with "--without-jasper" and
"--without-grass" options. [Both compiled OK, from 'clean' sources (viz,
after a make distclean eachtime)]
I try to import this ASCII DEM:
http://new-gis.mannlib.cornell.edu/CUGIR_Data/data1/x42elu.dem.gz
When I run
r.in.gdal i=x42elu.dem o=dem lo=x42
I get a segmentation fault.
Did you try 'gdalinfo'?
I have downloaded the file, and it reports:
gdalinfo x42elu.dem
Driver: USGSDEM/USGS Optional ASCII DEM
Size is 1041, 1395
Coordinate System is:
PROJCS["UTM Zone 18, Northern Hemisphere",
GEOGCS["NAD27",
DATUM["North_American_Datum_1927",
SPHEROID["Clarke 1866",6378206.4,294.978698213898,
AUTHORITY["EPSG","7008"]],
TOWGS84[-3,142,183,0,0,0,0],
AUTHORITY["EPSG","6267"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9108"]],
AXIS["Lat",NORTH],
AXIS["Long",EAST],
AUTHORITY["EPSG","4267"]],
PROJECTION["Transverse_Mercator"],
[...]
But, in fact I also get a segmentation fault (r.in.gdal from 5.3-CVS and
GDAL-CVS).
It crashes in wkt_to_grass() in main.c of r.in.gdal (crash
in line 1054: free( pszUnitsName ); ). Actually I don't know
why, maybe it does not fetch the projection info properly.
Other formats seem to be imported without trouble (GeoTIFF, for example).
FYI, the same file is imported without trouble on a SuSE7.1 box, with
grass5.0.2 and gdal1.1.9.
It seems to be related to USGS DEM only.
gdal_translate -of GTiFF x42elu.dem x42elu.tif
Input file size is 1041, 1395
0...10...20...30...40...50...60...70...80...90...100 - done.
r.in.gdal in=x42elu.tif out=dem loc=demdem
100%
CREATING SUPPORT FILES FOR dem
This works fine. So something runs odd when importing USGS DEM with
r.in.gdal.
Any idea of troubleshooting more than welcome, it's rather frustrating...
Thanks in advance
Unless it is fixed, you may use gdal_translate to convert the map to
GeoTIFF.
Regards
Markus