[GRASS-user] projection issue

Vincent Bain bain at toraval.fr
Mon Mar 24 05:10:03 PDT 2014


Comment on my previous message :

if we have a closer look at parameters, we notice a difference between
the current (grass location) projection parameters and epsg:3857
parameters :

Current :
        a (semi-major axis) = 6378137
        rf (flattening factor inverse) =298.257223563
am I wrong or as a consequence 
        b (semi-minor axis) = 6356752 ?

epsg:3857 :
        a (semi-major axis) = 6378137
        b (semi-minor axis) = 6378137 (a sphere ?)
        
Now if I run cs2cs with these distinct params :

         echo "6.46311951 45.78037144" | cs2cs +init=epsg:4326 -f "%.8f"
        +to +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0
        +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext
        +no_defs
returns
        719471.17278927	5745223.17292760 0.00000000

and
        echo "6.46311951 45.78037144" | cs2cs +init=epsg:4326 -f "%.8f"
        +to +proj=merc +a=6378137 +b=6356752 +lat_ts=0.0 +lon_0=0.0
        +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext
        +no_defs
returns
        719471.17278927	5714587.36406099 0.00000000
that is what m.proj initially said.

It looks as if grass was taking into account the flatten-shape of the
ellipsoid, where epsg registry considers it is a sphere.

Sorry for the noise, but I would actually enjoy understanding where I am
wrong ;-(

Bye,
Vincent



More information about the grass-user mailing list