<div dir="ltr"><div>There are some solutions to the problem of datum <Qatar_National_Datum_1995> not recognized by GRASS:</div><div><br></div><div>- add Qatar_National_Datum_1995 to the datums known by GRASS</div><div>- 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<br></div><div><br></div><div><div>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.</div><div><br></div><div>Markus M<br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Dec 3, 2019 at 8:52 PM Helmut Kudrnovsky <<a href="mailto:hellik@web.de">hellik@web.de</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Markus Neteler wrote<br>
> Hi,<br>
> <br>
> On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei <<br>
<br>
> zoltans@.co<br>
<br>
> > wrote:<br>
>><br>
>> Hi,<br>
>> I'm using EPSG 2932 and QGIS etc all are OK with it.<br>
>><br>
>> I have installed everything using OSGeo4W64, so how come Grass does not<br>
>> use the same projections database (as QGis etc)?<br>
>> (and please can I have a pointer as to how to introduce EPSG 2932 to<br>
>> Grass 7.8.1)<br>
>><br>
>> Thanks and regards,<br>
>> Zoltan<br>
>><br>
>> ******** Running against:<br>
>> D:\GDBroad_tiles\shp_road_poly_few\23103770_road_poly.shp<br>
>> WARNING: Datum <br>
> <Qatar_National_Datum_1995><br>
> not recognised by GRASS and<br>
>> no parameters found<br>
> <br>
> Is it possible that you have a PROJ software version mixup?<br>
> <br>
> I tried on my Linux box:<br>
> <br>
> grass78 -c epsg:2932 ~/grassdata/test_2932<br>
> Starting GRASS GIS...<br>
> Creating new GRASS GIS location <br>
> <test_2932><br>
> ...<br>
> Cleaning up temporary files...<br>
> <br>
> __________ ___ __________ _______________<br>
> / ____/ __ \/ | / ___/ ___/ / ____/ _/ ___/<br>
> / / __/ /_/ / /| | \__ \\_ \ / / __ / / \__ \<br>
> / /_/ / _, _/ ___ |___/ /__/ / / /_/ // / ___/ /<br>
> \____/_/ |_/_/ |_/____/____/ \____/___//____/<br>
> <br>
> Welcome to GRASS GIS 7.8.2dev (ad4836c73)<br>
> ...<br>
> <br>
> GRASS 7.8.2dev (test_2932):~ > g.proj -w<br>
> PROJCS["QND95 / Qatar National Grid",<br>
> GEOGCS["QND95",<br>
> DATUM["Qatar_National_Datum_1995",<br>
> SPHEROID["International 1924",6378388,297,<br>
> AUTHORITY["EPSG","7022"]],<br>
> <br>
> TOWGS84[-119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065],<br>
> AUTHORITY["EPSG","6614"]],<br>
> PRIMEM["Greenwich",0,<br>
> AUTHORITY["EPSG","8901"]],<br>
> UNIT["degree",0.0174532925199433,<br>
> AUTHORITY["EPSG","9122"]],<br>
> AUTHORITY["EPSG","4614"]],<br>
> PROJECTION["Transverse_Mercator"],<br>
> PARAMETER["latitude_of_origin",24.45],<br>
> PARAMETER["central_meridian",51.21666666666667],<br>
> PARAMETER["scale_factor",0.99999],<br>
> PARAMETER["false_easting",200000],<br>
> PARAMETER["false_northing",300000],<br>
> UNIT["metre",1,<br>
> AUTHORITY["EPSG","9001"]],<br>
> AXIS["Easting",EAST],<br>
> AXIS["Northing",NORTH],<br>
> AUTHORITY["EPSG","2932"]]<br>
<br>
when importing a shapefile with the same epsg projection, I get:<br>
<br>
v.in.ogr --verbose input=D:\wd\test_epsg2932.shp <br>
Using OGR driver 'ESRI Shapefile/ESRI Shapefile'<br>
WARNING: Datum <Qatar_National_Datum_1995> von GRASS nicht erkannt und keine<br>
Parameter gefunden.<br>
Die Projektionsinformationen des Eingabedatensatzes und der aktuellen<br>
Location scheinen übereinzustimmen.<br>
<br>
D:\wd>ogrinfo --version<br>
GDAL 3.0.2, released 2019/10/28<br>
<br>
D:\wd>ogrinfo test_epsg2932.shp -al -so<br>
INFO: Open of `test_epsg2932.shp'<br>
using driver `ESRI Shapefile' successful.<br>
<br>
Layer name: test_epsg2932<br>
Metadata:<br>
DBF_DATE_LAST_UPDATE=2019-12-03<br>
Geometry: Point<br>
Feature Count: 1<br>
Extent: (823205.438352, 25769.568359) - (823205.438352, 25769.568359)<br>
Layer SRS WKT:<br>
PROJCS["QND95 / Qatar National Grid",<br>
GEOGCS["QND95",<br>
DATUM["Qatar_National_Datum_1995",<br>
SPHEROID["International 1924",6378388,297,<br>
AUTHORITY["EPSG","7022"]],<br>
AUTHORITY["EPSG","6614"]],<br>
PRIMEM["Greenwich",0,<br>
AUTHORITY["EPSG","8901"]],<br>
UNIT["degree",0.0174532925199433,<br>
AUTHORITY["EPSG","9122"]],<br>
AUTHORITY["EPSG","4614"]],<br>
PROJECTION["Transverse_Mercator"],<br>
PARAMETER["latitude_of_origin",24.45],<br>
PARAMETER["central_meridian",51.2166666666667],<br>
PARAMETER["scale_factor",0.99999],<br>
PARAMETER["false_easting",200000],<br>
PARAMETER["false_northing",300000],<br>
UNIT["metre",1,<br>
AUTHORITY["EPSG","9001"]],<br>
AXIS["Easting",EAST],<br>
AXIS["Northing",NORTH],<br>
AUTHORITY["EPSG","2932"]]<br>
Data axis to CRS axis mapping: 1,2<br>
id: Integer64 (10.0)<br>
test: Integer (9.0)<br>
<br>
D:\wd>testepsg test_epsg2932.prj<br>
Validate Succeeds.<br>
WKT[test_epsg2932.prj] =<br>
PROJCS["QND95 / Qatar National Grid",<br>
GEOGCS["QND95",<br>
DATUM["Qatar_National_Datum_1995",<br>
SPHEROID["International 1924",6378388,297,<br>
AUTHORITY["EPSG","7022"]],<br>
AUTHORITY["EPSG","6614"]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["Degree",0.0174532925199433]],<br>
PROJECTION["Transverse_Mercator"],<br>
PARAMETER["latitude_of_origin",24.45],<br>
PARAMETER["central_meridian",51.2166666666667],<br>
PARAMETER["scale_factor",0.99999],<br>
PARAMETER["false_easting",200000],<br>
PARAMETER["false_northing",300000],<br>
UNIT["metre",1,<br>
AUTHORITY["EPSG","9001"]],<br>
AXIS["Easting",EAST],<br>
AXIS["Northing",NORTH]]<br>
<br>
Simplified WKT[test_epsg2932.prj] =<br>
PROJCS["QND95 / Qatar National Grid",<br>
GEOGCS["QND95",<br>
DATUM["Qatar_National_Datum_1995",<br>
SPHEROID["International 1924",6378388,297]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["Degree",0.0174532925199433]],<br>
PROJECTION["Transverse_Mercator"],<br>
PARAMETER["latitude_of_origin",24.45],<br>
PARAMETER["central_meridian",51.2166666666667],<br>
PARAMETER["scale_factor",0.99999],<br>
PARAMETER["false_easting",200000],<br>
PARAMETER["false_northing",300000],<br>
UNIT["metre",1]]<br>
<br>
Old Style WKT[test_epsg2932.prj] = PROJCS["QND95 / Qatar National<br>
Grid",GEOGCS["QND95",DATUM["Qatar_National_Datum_1995",SPHEROID["International<br>
1924",6378388,297]],PRIMEM["Greenwich",0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",24.45],PARAMETER["central_meridian",51.2166666666667],PARAMETER["scale_factor",0.99999],PARAMETER["false_easting",200000],PARAMETER["false_northing",300000],UNIT["metre",1]]<br>
ESRI'ified WKT[test_epsg2932.prj] =<br>
PROJCS["QND_1995_Qatar_National_Grid",<br>
GEOGCS["GCS_QND_1995",<br>
DATUM["D_QND_1995",<br>
SPHEROID["International_1924",6378388.0,297.0]],<br>
PRIMEM["Greenwich",0.0],<br>
UNIT["Degree",0.0174532925199433]],<br>
PROJECTION["Transverse_Mercator"],<br>
PARAMETER["False_Easting",200000.0],<br>
PARAMETER["False_Northing",300000.0],<br>
PARAMETER["Central_Meridian",51.2166666666667],<br>
PARAMETER["Scale_Factor",0.99999],<br>
PARAMETER["Latitude_Of_Origin",24.45],<br>
UNIT["Meter",1.0]]<br>
PROJ.4 rendering of [test_epsg2932.prj] = +proj=tmerc +lat_0=24.45<br>
+lon_0=51.2166666666667 +k=0.99999 +x_0=200000 +y_0=300000 +ellps=intl<br>
+units=m +no_defs<br>
<br>
maybe worth to open a ticket<br>
<br>
<br>
<br>
<br>
-----<br>
best regards<br>
Helmut<br>
--<br>
Sent from: <a href="http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html" rel="noreferrer" target="_blank">http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html</a><br>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a></blockquote></div>