[GRASS-user] r.in.gdal global raster incorrect extents
Hamish
hamish_b at yahoo.com
Tue Oct 14 01:07:11 EDT 2008
> > gdalinfo W180N10.tif
> ....
> > Origin = (-180.000416666999996,10.000416666700000)
> > Pixel Size = (0.000833333333333,-0.000833333333333)
Hamish:
> note in this case the amount it exceeds 180 by is 1/2 the
> cell resolution.
....
> So in this case, the SRTM tile is not broken, it is what it
> is and we just have to process it correctly to compensate.
oops, -180.000416666999996 is broken after all, "180.000416666 999996"
it seems that 180.000416667 is asked for more precision than it has.
(but it works fine )
I see that gdalinfo 1.5.2 uses %.15f:
printf( "Origin = (%.15f,%.15f)\n",
adfGeoTransform[0], adfGeoTransform[3] );
printf( "Pixel Size = (%.15f,%.15f)\n",
adfGeoTransform[1], adfGeoTransform[5] );
but why then does it work fine for me?
I downloaded a 180W tile, that looks fine too:
ftp://e0srp01u.ecs.nasa.gov/srtm/version2/SRTM3/Australia/S20W180.hgt.zip
....
Origin = (-180.000416666666666,-18.999583333333334)
Pixel Size = (0.000833333333333,-0.000833333333333)
where do your .tif files come from?
(see also this discussion:
http://thread.gmane.org/gmane.comp.gis.grass.devel/29622/ )
Hamish
More information about the grass-user
mailing list