Thanks Helli and Moritz for help :) [...]
> I can confirm the issue when using v.in.ogr + v.proj, when using the file > to define the location: > > grass73 -c /data/GEODATA/OSM/coastlines-generalized-3857/coastlines_z8.shp > /data/GRASS/DATA7/test_3857/ > v.in.ogr /data/GEODATA/OSM/coastlines-generalized-3857/coastlines_z8.shp > > and then in the EPSG 3035 location: > > v.proj location=test_3857 mapset=PERMANENT input=coastlines_z8 > Reprojecting primitives ... > ATTENTION: pj_transform() a échoué: failed to load datum shift file > ERREUR :Unable to re-project vector map <coastlines_z8@PERMANENT> from > <test_3857> > > > However, when I create the 3857 location using the EPSG code, it works: > > grass73 -c epsg:3857 /data/GRASS/DATA7/test_3857_epsg/ > [...] > v.proj location=test_3857_epsg mapset=PERMANENT input=coastlines_z8 > > success. > yeah, I just tested here with the same result... so, creating the 3857 location with the coastlines .shp file then produces the error when reprojecting to 3035 because of this extra "nadgrids: @null" in the .prj... good to know :) There is a difference in the projection definition between the two: > > diff -u test_3857/PERMANENT/PROJ_INFO test_3857_epsg/PERMANENT/PROJ_INFO > --- test_3857/PERMANENT/PROJ_INFO 2016-08-02 18:24:25.965809631 +0200 > +++ test_3857_epsg/PERMANENT/PROJ_INFO 2016-08-02 18:28:18.427085934 +0200 > @@ -7,7 +7,5 @@ > x_0: 0.0 > y_0: 0 > k: 1.0 > -nadgrids: @null > wktext: defined > no_defs: defined > -over: defined > > Erasing the 'nadgrids: @null' from the first allows to reproject. > > The .prj file of the original shapefile contains: > > EXTENSION["PROJ4","+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 +over"] > > I guess the issue comes from how this is handled... > It seems so... Thank you very much to both for your time!!! Vero
_______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
