<div dir="ltr">Hello list, <div><br></div><div>I downloaded generalized coastlines from: <a href="http://openstreetmapdata.com/data/generalized-coastlines">http://openstreetmapdata.com/data/generalized-coastlines</a> and as they are in mercator projection (epsg:3857) and I'm working in eu_laea, i wanted to use v.import to reproject (and subset to my region) on the fly... but I can't get v.import nor v.proj to work...</div><div><br></div><div>This is the command I'm using (copied from the GUI) and the error:</div><div><br></div><div><div>v.import input=coastlines_z8.shp epsg=3857 datum_trans=0 --> the same with whatever value > 0 that I put here</div><div><br></div><div>WARNING: All available OGR layers will be imported into vector map <coastlines_z8></div><div>Projection information updated</div><div>WARNING: All available OGR layers will be imported into vector map <coastlines_z8></div><div>Importing </home/veroandreo/gisdata/coastlines-generalized-3857/coastlines_z8.shp> ...</div><div>WARNING: All available OGR layers will be imported into vector map <coastlines_z8></div><div>Check if OGR layer <coastlines_z8> contains polygons...</div><div>Default driver / database set to:</div><div>driver: sqlite</div><div>database: $GISDBASE/$LOCATION_NAME/$MAPSET/sqlite/sqlite.db</div><div>Importing 30179 features (OGR layer <coastlines_z8>)...</div><div>-----------------------------------------------------</div><div>Building topology for vector map <coastlines_z8@PERMANENT>...</div><div>Registering primitives...</div><div>30179 primitives registered</div><div>2796082 vertices registered</div><div>Building areas...</div><div>[...]</div><div>Reprojecting <coastlines_z8>...</div><div>Reprojecting primitives ...</div><div>WARNING: pj_transform() failed: failed to load datum shift file</div><div>ERROR: Unable to re-project vector map <coastlines_z8@PERMANENT> from <temp_import_location_12909></div><div>ERROR: Unable to to reproject vector <coastlines_z8></div></div><div><br></div><div>Moreover, if I use datum_trans=-1 which is supposed to throw the list of parameters, I get the following error:<br></div><div><br></div><div><div>v.import input=coastlines_z8.shp epsg=3857 <b>datum_trans=-1</b></div><div><b><br></b></div><div>ERROR: No output format specified, define one of flags -p, -g, -j, or -w</div><div>Traceback (most recent call last):</div><div>  File "/home/veroandreo/software/grass-7.3.svn/dist.x86_64</div><div>-pc-linux-gnu/scripts/v.import", line 290, in <module></div><div>    sys.exit(main())</div><div>  File "/home/veroandreo/software/grass-7.3.svn/dist.x86_64</div><div>-pc-linux-gnu/scripts/v.import", line 158, in main</div><div>    datum_trans=options['datum_trans'])</div><div>  File "/home/veroandreo/software/grass-7.3.svn/dist.x86_64</div><div>-pc-linux-gnu/etc/python/grass/script/core.py", line 410, in</div><div>run_command</div><div>    return handle_errors(returncode, returncode, args,</div><div>kwargs)</div><div>  File "/home/veroandreo/software/grass-7.3.svn/dist.x86_64</div><div>-pc-linux-gnu/etc/python/grass/script/core.py", line 329, in</div><div>handle_errors</div><div>    returncode=returncode)</div><div>grass.exceptions.CalledModuleError: Module run None</div><div>['g.proj', 'epsg=3857', 'datum_trans=-1'] ended with error</div><div>Process ended with non-zero return code 1. See errors in the</div><div>(error) output.</div></div><div><br></div><div>Then, I created a new location from the .prj and imported the file. So, I went to my target location to reproject with the good old v.proj but I get the same warning message at the end and the file is not projected:</div><div><br></div><div>v.proj location=mercator mapset=PERMANENT input=coastlines_z8<br></div><div><br></div><div><div>WARNING: pj_transform() failed: failed to load datum shift file</div><div>ERROR: Unable to re-project vector map <coastlines_z8@PERMANENT> from <mercator></div></div><div><br></div><div>In conclusion:</div><div><br></div><div>1. can someone illustrate on what datum_trans means in v.import, what number goes there and/or provide a working example, please?</div><div><br></div><div>2. what is the datum shift file and where is it supposed to be?<br></div><div><br></div><div>Thanks a million in advance for any help/hint!</div><div><br></div><div>Best, <br>Vero</div><div><br></div><div>ps1: I'm using grass 73 svn r69040M under Fedora 24</div><div>ps2: my way-around to get the file imported was: opening the mercator file in qgis, saving as eu_laea and importing with v.in.ogr into my eu_laea location (suboptimal, but works)</div></div>