[GRASS-user] Datum not recognized by Grass
Helmut Kudrnovsky
hellik at web.de
Tue Dec 3 14:30:59 PST 2019
> ignore the warning and use GRASS with PROJ6, granted that authority name
(e.g. EPSG) and authority code (e.g. 2932) are known for both CRS's in case
of reprojection
this issue is already by GRASS with PROJ6, see
####################
GRASS version: 7.8.1
Code revision: c865432c9
Build date: 2019-11-10
Build platform: x86_64-w64-mingw32
GDAL: 3.0.2
PROJ: 6.2.1 <=
GEOS: 3.8.0
SQLite: 3.29.0
Python: 3.7.0
wxPython: 4.0.7
Platform: Windows-10-10.0.18362-SP0 (OSGeo4W)
####################
and the output from the underlying PROJ: 6.2.1
####################
C:\>projinfo EPSG:2932 -o PROJ,WKT2_2018
PROJ.4 string:
+proj=tmerc +lat_0=24.45 +lon_0=51.2166666666667 +k=0.99999 +x_0=200000
+y_0=300000 +ellps=intl
+towgs84=-119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065
+units=m +no_defs +type=crs
WKT2_2018 string:
PROJCRS["QND95 / Qatar National Grid",
BASEGEOGCRS["QND95",
DATUM["Qatar National Datum 1995",
ELLIPSOID["International 1924",6378388,297,
LENGTHUNIT["metre",1]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433]],
ID["EPSG",4614]],
CONVERSION["Qatar National Grid",
METHOD["Transverse Mercator",
ID["EPSG",9807]],
PARAMETER["Latitude of natural origin",24.45,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8801]],
PARAMETER["Longitude of natural origin",51.2166666666667,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8802]],
PARAMETER["Scale factor at natural origin",0.99999,
SCALEUNIT["unity",1],
ID["EPSG",8805]],
PARAMETER["False easting",200000,
LENGTHUNIT["metre",1],
ID["EPSG",8806]],
PARAMETER["False northing",300000,
LENGTHUNIT["metre",1],
ID["EPSG",8807]]],
CS[Cartesian,2],
AXIS["(E)",east,
ORDER[1],
LENGTHUNIT["metre",1]],
AXIS["(N)",north,
ORDER[2],
LENGTHUNIT["metre",1]],
USAGE[
SCOPE["unknown"],
AREA["Qatar - onshore"],
BBOX[24.55,50.69,26.2,51.68]],
ID["EPSG",2932]]
####################
>The problem is that GRASS still assumes datum transformation from X to
WGS84, whereas with PROJ6,
>WGS84 is no longer required as pivot datum. Datum transformations from
datum X to datum Y can
>sometimes (often) be done without going through WGS84. The requirement is
to have EPSG or other
>authority codes and to use PROJ 6.
this warning is just from a simple v.in.ogr of a EPSG:2932-based shapefile
into a EPSG:2932-created GRASS location.
so, there may be some more adaptions needed for actions like v.in.ogr and
co?
-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html
More information about the grass-user
mailing list