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