<div dir="ltr"><div>Hi Markus, list</div><div><br></div><div>Maybe (un)related. </div><div><br></div><div>From a lat-long location, I'm trying to import TOPO-WMS from mundialis (je! I like the colors very much! :)) and so, I saved my region of interest (my_region) and run:</div><div><br></div><div>GRASS 7.1.svn (latlong_wgs84):~ > r.in.wms url=<a href="http://ows.mundialis.de/services/service" target="_blank">http://ows.mundialis.de/services/service</a>? output=mundialis_topo layers=TOPO-WMS styles=default srs=4326 region=my_region format=png<br></div><div><div>Downloading data from WMS server...</div><div>Importing raster map into GRASS...</div><div>ERROR: Projection of dataset does not appear to match current location.</div><div><br></div><div> Location PROJ_INFO is:</div><div> name: Lat/Lon</div><div> proj: ll</div><div> datum: wgs84</div><div> ellps: wgs84</div><div> no_defs: defined</div><div> towgs84: 0.000,0.000,0.000</div><div><br></div><div> Import dataset PROJ_INFO is:</div><div> cellhd.proj = 0 (unreferenced/unknown)</div><div><br></div><div> In case of no significant differences in the projection definitions,</div><div> use the -o flag to ignore them and use current location definition.</div><div> Consider generating a new location from the input dataset using the</div><div> 'location' parameter.</div><div>ERROR: r.in.gdal failed</div></div><div><br></div><div>but i get ERROR, even though EPSG:4326 is listed as available for the layer.</div><div><br></div><div>GRASS 7.1.svn (latlong_wgs84):~ > r.in.wms url=<a href="http://ows.mundialis.de/services/service" target="_blank">http://ows.mundialis.de/services/service</a>? -c<br></div><div>...<br></div><div> <Layer queryable="1"><br></div><div><div> <Name>TOPO-WMS</Name></div><div> <Title>Topographic WMS - by mundialis</Title></div><div> <LatLonBoundingBox minx="-180" miny="-89.999999" maxx="180" maxy="89.999999" /></div><div> <BoundingBox SRS="EPSG:4326" minx="-180" miny="-90" maxx="180" maxy="90" /><br></div><div>...</div></div><div><br></div><div>Also from the GUI I get the same error regarding projection (see screenshot). Maybe I'm doing it completely wrong, dunno... help is welcome :)<br></div><div><br></div><div>Testing examples from the manual page, I get the same error:</div><div><br></div><div><div>r.in.wms output=kn url=<a href="http://wms.cuzk.cz/wms.asp">http://wms.cuzk.cz/wms.asp</a> layers=prehledka_kraju-linie srs=4326 format=png</div><div>Downloading data from WMS server...</div><div>Importing raster map into GRASS...</div><div>ERROR: Projection of dataset does not appear to match current location.</div><div><br></div><div> Location PROJ_INFO is:</div><div> name: Lat/Lon</div><div> proj: ll</div><div> datum: wgs84</div><div> ellps: wgs84</div><div> no_defs: defined</div><div> towgs84: 0.000,0.000,0.000</div><div><br></div><div> Import dataset PROJ_INFO is:</div><div> cellhd.proj = 0 (unreferenced/unknown)</div><div><br></div><div> In case of no significant differences in the projection definitions,</div><div> use the -o flag to ignore them and use current location definition.</div><div> Consider generating a new location from the input dataset using the</div><div> 'location' parameter.</div><div>ERROR: r.in.gdal failed</div></div><div><br></div><div>My box is Fedora 23 and I'm using:<br></div><div>version=7.1.svn<br></div><div><div>revision=r68259<br></div><div>build_date=2016-04-14</div><div>build_platform=x86_64-pc-linux-gnu</div><div>proj4=4.9.1<br></div><div>gdal=2.0.2</div><div>geos=3.4.2</div><div>sqlite=3.11.0</div></div><div><br></div><div>Cheers, </div><div>Vero</div></div><div class="gmail_extra"><br><div class="gmail_quote">2016-04-14 11:54 GMT-03:00 Markus Neteler <span dir="ltr"><<a href="mailto:neteler@osgeo.org" target="_blank">neteler@osgeo.org</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
I have generated a new location from the EU DEM 25m which comes with<br>
these metadata:<br>
<br>
gdalinfo eudem_dem_3035_europe.tif<br>
Driver: GTiff/GeoTIFF<br>
Files: eudem_dem_3035_europe.tif<br>
eudem_dem_3035_europe.tif.ovr<br>
Size is 240000, 200000<br>
Coordinate System is:<br>
PROJCS["ETRS89 / LAEA Europe",<br>
GEOGCS["ETRS89",<br>
DATUM["European_Terrestrial_Reference_System_1989",<br>
SPHEROID["GRS 1980",6378137,298.2572221010002,<br>
AUTHORITY["EPSG","7019"]],<br>
AUTHORITY["EPSG","6258"]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4258"]],<br>
PROJECTION["Lambert_Azimuthal_Equal_Area"],<br>
PARAMETER["latitude_of_center",52],<br>
PARAMETER["longitude_of_center",10],<br>
PARAMETER["false_easting",4321000],<br>
PARAMETER["false_northing",3210000],<br>
UNIT["metre",1,<br>
AUTHORITY["EPSG","9001"]],<br>
AUTHORITY["EPSG","3035"]]<br>
Origin = (2000000.000000000000000,6000000.000000000000000)<br>
Pixel Size = (25.000000000000000,-25.000000000000000)<br>
[...]<br>
<br>
I have then imported the file, cut out a portion with r.mapcalc and<br>
packed up the result with r.pack since I have to deliver the data to<br>
another GRASS user.<br>
<br>
For a test, I created a new location with<br>
<br>
grass70 -c epsg:3035 ~/grassdata/eu_laea_extra<br>
<br>
and tried to import it again. Surprisingly, different spelling created an issue:<br>
<br>
r.unpack subregion_eudem_25m.pack<br>
WARNING: Difference between PROJ_INFO file of packed map and of current<br>
location:<br>
- name: ETRS89 / LAEA Europe<br>
? -------<br>
+ name: ETRS89 / ETRS-LAEA<br>
? +++++<br>
datum: etrs89<br>
ellps: grs80<br>
proj: laea<br>
lat_0: 52<br>
lon_0: 10<br>
x_0: 4321000<br>
y_0: 3210000<br>
- towgs84: 0,0,0,0,0,0,0<br>
no_defs: defined<br>
ERROR: Projection of dataset does not appear to match current location. In<br>
case of no significant differences in the projection definitions,<br>
use the -o flag to ignore them and use current location definition.<br>
<br>
Question:<br>
in lib/proj/convert.c we have a function papszDatumEquiv(). Do we also<br>
need a similar function for names?<br>
<br>
thanks,<br>
Markus<br>
<br>
PS: If it matters, here the GDAL version (a Linux Mint box):<br>
gdalinfo --version<br>
GDAL 1.10.1, released 2013/08/26<br>
_______________________________________________<br>
grass-dev mailing list<br>
<a href="mailto:grass-dev@lists.osgeo.org">grass-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-dev" rel="noreferrer" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-dev</a></blockquote></div><br></div>