[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