[GRASS-user] Datum not recognized by Grass

GRASS GRASS user
,
1

Hi,
I'm using EPSG 2932 and QGIS etc all are OK with it.

I have installed everything using OSGeo4W64, so how come Grass does not use the same projections database (as QGis etc)?
(and please can I have a pointer as to how to introduce EPSG 2932 to Grass 7.8.1)

Thanks and regards,
Zoltan

******** Running against: D:\GDBroad_tiles\shp_road_poly_few\23103770_road_poly.shp
WARNING: Datum <Qatar_National_Datum_1995> not recognised by GRASS and no parameters found
Check if OGR layer <23103770_road_poly> contains polygons...
Creating attribute table for layer <23103770_road_poly>...
WARNING: Name <23103770_road_poly> is not SQL compliant. Must start with a letter.
Importing 28 features (OGR layer <23103770_road_poly>)...

--

=============================================
Zoltan Szecsei GPrGISc 0031
Geograph (Pty) Ltd.
GIS and Photogrammetric Services

P.O. Box 7, Muizenberg 7950, South Africa.

Mobile: +27-83-6004028 (WhatsApp only)
Qatar: +974 5083 2722 www.geograph.co.za

2

Hi,

On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei <zoltans@geograph.co.za> wrote:

Hi,
I'm using EPSG 2932 and QGIS etc all are OK with it.

I have installed everything using OSGeo4W64, so how come Grass does not
use the same projections database (as QGis etc)?
(and please can I have a pointer as to how to introduce EPSG 2932 to
Grass 7.8.1)

Thanks and regards,
Zoltan

******** Running against:
D:\GDBroad_tiles\shp_road_poly_few\23103770_road_poly.shp
WARNING: Datum <Qatar_National_Datum_1995> not recognised by GRASS and
no parameters found

Is it possible that you have a PROJ software version mixup?

I tried on my Linux box:

grass78 -c epsg:2932 ~/grassdata/test_2932
Starting GRASS GIS...
Creating new GRASS GIS location <test_2932>...
Cleaning up temporary files...

          __________ ___ __________ _______________
         / ____/ __ \/ | / ___/ ___/ / ____/ _/ ___/
        / / __/ /_/ / /| | \__ \\_ \ / / __ / / \__ \
       / /_/ / _, _/ ___ |___/ /__/ / / /_/ // / ___/ /
       \____/_/ |_/_/ |_/____/____/ \____/___//____/

Welcome to GRASS GIS 7.8.2dev (ad4836c73)
...

GRASS 7.8.2dev (test_2932):~ > g.proj -w
PROJCS["QND95 / Qatar National Grid",
    GEOGCS["QND95",
        DATUM["Qatar_National_Datum_1995",
            SPHEROID["International 1924",6378388,297,
                AUTHORITY["EPSG","7022"]],
            TOWGS84[-119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065],
            AUTHORITY["EPSG","6614"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4614"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",24.45],
    PARAMETER["central_meridian",51.21666666666667],
    PARAMETER["scale_factor",0.99999],
    PARAMETER["false_easting",200000],
    PARAMETER["false_northing",300000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH],
    AUTHORITY["EPSG","2932"]]

# The versions here:

GRASS 7.8.2dev (test_2932):~ > g.version -rge
version=7.8.2dev
date=2019
revision=ad4836c73
build_date=2019-11-28
build_platform=x86_64-pc-linux-gnu
build_off_t_size=8
libgis_revision=00000
libgis_date="?"
proj=5.2.0
gdal=2.3.2
geos=3.7.1
sqlite=3.30.0

Looks all fine.
Please check if an old PROJ 4 is installed which might confuse things.

Best
Markus

3

On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei <zoltans@geograph.co.za> wrote:

Hi,
I'm using EPSG 2932 and QGIS etc all are OK with it.

I have installed everything using OSGeo4W64, so how come Grass does not
use the same projections database (as QGis etc)?

Please post the output of

g.version -rge

Best
Markus

4

Markus Neteler wrote

Hi,

On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei &lt;

zoltans@.co

&gt; wrote:

Hi,
I'm using EPSG 2932 and QGIS etc all are OK with it.

I have installed everything using OSGeo4W64, so how come Grass does not
use the same projections database (as QGis etc)?
(and please can I have a pointer as to how to introduce EPSG 2932 to
Grass 7.8.1)

Thanks and regards,
Zoltan

******** Running against:
D:\GDBroad_tiles\shp_road_poly_few\23103770_road_poly.shp
WARNING: Datum

<Qatar_National_Datum_1995>
not recognised by GRASS and

no parameters found

Is it possible that you have a PROJ software version mixup?

I tried on my Linux box:

grass78 -c epsg:2932 ~/grassdata/test_2932
Starting GRASS GIS...
Creating new GRASS GIS location
<test_2932>
...
Cleaning up temporary files...

          __________ ___ __________ _______________
         / ____/ __ \/ | / ___/ ___/ / ____/ _/ ___/
        / / __/ /_/ / /| | \__ \\_ \ / / __ / / \__ \
       / /_/ / _, _/ ___ |___/ /__/ / / /_/ // / ___/ /
       \____/_/ |_/_/ |_/____/____/ \____/___//____/

Welcome to GRASS GIS 7.8.2dev (ad4836c73)
...

GRASS 7.8.2dev (test_2932):~ > g.proj -w
PROJCS["QND95 / Qatar National Grid",
    GEOGCS["QND95",
        DATUM["Qatar_National_Datum_1995",
            SPHEROID["International 1924",6378388,297,
                AUTHORITY["EPSG","7022"]],
           
TOWGS84[-119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065],
            AUTHORITY["EPSG","6614"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4614"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",24.45],
    PARAMETER["central_meridian",51.21666666666667],
    PARAMETER["scale_factor",0.99999],
    PARAMETER["false_easting",200000],
    PARAMETER["false_northing",300000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH],
    AUTHORITY["EPSG","2932"]]

# The versions here:

GRASS 7.8.2dev (test_2932):~ > g.version -rge
version=7.8.2dev
date=2019
revision=ad4836c73
build_date=2019-11-28
build_platform=x86_64-pc-linux-gnu
build_off_t_size=8
libgis_revision=00000
libgis_date="?"
proj=5.2.0
gdal=2.3.2
geos=3.7.1
sqlite=3.30.0

Looks all fine.
Please check if an old PROJ 4 is installed which might confuse things.

with

C:\>g.version -rge
version=7.8.1
date=2019
revision=c865432c9
build_date=2019-11-10
build_platform=x86_64-w64-mingw32
build_off_t_size=8
libgis_revision=00000
libgis_date="?"
proj=6.2.1
gdal=3.0.2
geos=3.8.0
sqlite=3.29.0

I get with the location wizzard:

Starting GRASS GIS...
Traceback (most recent call last):
  File
"C:\OSGEO4~1\apps\grass\grass78\gui\wxpython\location_wizard\wizard.py",
line 1607, in OnPageChanging
    self.parent.parent, transforms=ret)
  File
"C:\OSGEO4~1\apps\grass\grass78\gui\wxpython\location_wizard\dialogs.py",
line 617, in __init__
    transforms = '---\n\n0\nDo not apply any datum transformations\n\n' +
transforms
TypeError: can only concatenate str (not "bytes") to str

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

5

Markus Neteler wrote

Hi,

On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei &lt;

zoltans@.co

&gt; wrote:

Hi,
I'm using EPSG 2932 and QGIS etc all are OK with it.

I have installed everything using OSGeo4W64, so how come Grass does not
use the same projections database (as QGis etc)?
(and please can I have a pointer as to how to introduce EPSG 2932 to
Grass 7.8.1)

Thanks and regards,
Zoltan

******** Running against:
D:\GDBroad_tiles\shp_road_poly_few\23103770_road_poly.shp
WARNING: Datum

<Qatar_National_Datum_1995>
not recognised by GRASS and

no parameters found

Is it possible that you have a PROJ software version mixup?

I tried on my Linux box:

grass78 -c epsg:2932 ~/grassdata/test_2932
Starting GRASS GIS...
Creating new GRASS GIS location
<test_2932>
...
Cleaning up temporary files...

          __________ ___ __________ _______________
         / ____/ __ \/ | / ___/ ___/ / ____/ _/ ___/
        / / __/ /_/ / /| | \__ \\_ \ / / __ / / \__ \
       / /_/ / _, _/ ___ |___/ /__/ / / /_/ // / ___/ /
       \____/_/ |_/_/ |_/____/____/ \____/___//____/

Welcome to GRASS GIS 7.8.2dev (ad4836c73)
...

GRASS 7.8.2dev (test_2932):~ > g.proj -w
PROJCS["QND95 / Qatar National Grid",
    GEOGCS["QND95",
        DATUM["Qatar_National_Datum_1995",
            SPHEROID["International 1924",6378388,297,
                AUTHORITY["EPSG","7022"]],
           
TOWGS84[-119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065],
            AUTHORITY["EPSG","6614"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4614"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",24.45],
    PARAMETER["central_meridian",51.21666666666667],
    PARAMETER["scale_factor",0.99999],
    PARAMETER["false_easting",200000],
    PARAMETER["false_northing",300000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH],
    AUTHORITY["EPSG","2932"]]

C:\>g.proj -w
PROJCS["QND95 / Qatar National Grid",
    GEOGCS["QND95",
        DATUM["Qatar_National_Datum_1995",
            SPHEROID["International 1924",6378388,297,
                AUTHORITY["EPSG","7022"]],
           
TOWGS84[-119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065],
            AUTHORITY["EPSG","6614"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4614"]],
    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,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH],
    AUTHORITY["EPSG","2932"]]

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

6

Markus Neteler wrote

Hi,

On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei &lt;

zoltans@.co

&gt; wrote:

Hi,
I'm using EPSG 2932 and QGIS etc all are OK with it.

I have installed everything using OSGeo4W64, so how come Grass does not
use the same projections database (as QGis etc)?
(and please can I have a pointer as to how to introduce EPSG 2932 to
Grass 7.8.1)

Thanks and regards,
Zoltan

******** Running against:
D:\GDBroad_tiles\shp_road_poly_few\23103770_road_poly.shp
WARNING: Datum

<Qatar_National_Datum_1995>
not recognised by GRASS and

no parameters found

Is it possible that you have a PROJ software version mixup?

I tried on my Linux box:

grass78 -c epsg:2932 ~/grassdata/test_2932
Starting GRASS GIS...
Creating new GRASS GIS location
<test_2932>
...
Cleaning up temporary files...

          __________ ___ __________ _______________
         / ____/ __ \/ | / ___/ ___/ / ____/ _/ ___/
        / / __/ /_/ / /| | \__ \\_ \ / / __ / / \__ \
       / /_/ / _, _/ ___ |___/ /__/ / / /_/ // / ___/ /
       \____/_/ |_/_/ |_/____/____/ \____/___//____/

Welcome to GRASS GIS 7.8.2dev (ad4836c73)
...

GRASS 7.8.2dev (test_2932):~ > g.proj -w
PROJCS["QND95 / Qatar National Grid",
    GEOGCS["QND95",
        DATUM["Qatar_National_Datum_1995",
            SPHEROID["International 1924",6378388,297,
                AUTHORITY["EPSG","7022"]],
           
TOWGS84[-119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065],
            AUTHORITY["EPSG","6614"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4614"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",24.45],
    PARAMETER["central_meridian",51.21666666666667],
    PARAMETER["scale_factor",0.99999],
    PARAMETER["false_easting",200000],
    PARAMETER["false_northing",300000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH],
    AUTHORITY["EPSG","2932"]]

when importing a shapefile with the same epsg projection, I get:

v.in.ogr --verbose input=D:\wd\test_epsg2932.shp
Using OGR driver 'ESRI Shapefile/ESRI Shapefile'
WARNING: Datum <Qatar_National_Datum_1995> von GRASS nicht erkannt und keine
Parameter gefunden.
Die Projektionsinformationen des Eingabedatensatzes und der aktuellen
Location scheinen übereinzustimmen.

D:\wd>ogrinfo --version
GDAL 3.0.2, released 2019/10/28

D:\wd>ogrinfo test_epsg2932.shp -al -so
INFO: Open of `test_epsg2932.shp'
      using driver `ESRI Shapefile' successful.

Layer name: test_epsg2932
Metadata:
  DBF_DATE_LAST_UPDATE=2019-12-03
Geometry: Point
Feature Count: 1
Extent: (823205.438352, 25769.568359) - (823205.438352, 25769.568359)
Layer SRS WKT:
PROJCS["QND95 / Qatar National Grid",
    GEOGCS["QND95",
        DATUM["Qatar_National_Datum_1995",
            SPHEROID["International 1924",6378388,297,
                AUTHORITY["EPSG","7022"]],
            AUTHORITY["EPSG","6614"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4614"]],
    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,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH],
    AUTHORITY["EPSG","2932"]]
Data axis to CRS axis mapping: 1,2
id: Integer64 (10.0)
test: Integer (9.0)

D:\wd>testepsg test_epsg2932.prj
Validate Succeeds.
WKT[test_epsg2932.prj] =
PROJCS["QND95 / Qatar National Grid",
    GEOGCS["QND95",
        DATUM["Qatar_National_Datum_1995",
            SPHEROID["International 1924",6378388,297,
                AUTHORITY["EPSG","7022"]],
            AUTHORITY["EPSG","6614"]],
        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,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH]]

Simplified WKT[test_epsg2932.prj] =
PROJCS["QND95 / Qatar National Grid",
    GEOGCS["QND95",
        DATUM["Qatar_National_Datum_1995",
            SPHEROID["International 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]]

Old Style WKT[test_epsg2932.prj] = PROJCS["QND95 / Qatar National
Grid",GEOGCS["QND95",DATUM["Qatar_National_Datum_1995",SPHEROID["International
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]]
ESRI'ified WKT[test_epsg2932.prj] =
PROJCS["QND_1995_Qatar_National_Grid",
    GEOGCS["GCS_QND_1995",
        DATUM["D_QND_1995",
            SPHEROID["International_1924",6378388.0,297.0]],
        PRIMEM["Greenwich",0.0],
        UNIT["Degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["False_Easting",200000.0],
    PARAMETER["False_Northing",300000.0],
    PARAMETER["Central_Meridian",51.2166666666667],
    PARAMETER["Scale_Factor",0.99999],
    PARAMETER["Latitude_Of_Origin",24.45],
    UNIT["Meter",1.0]]
PROJ.4 rendering of [test_epsg2932.prj] = +proj=tmerc +lat_0=24.45
+lon_0=51.2166666666667 +k=0.99999 +x_0=200000 +y_0=300000 +ellps=intl
+units=m +no_defs

maybe worth to open a ticket

-----
best regards
Helmut
--
Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

7

There are some solutions to the problem of datum <Qatar_National_Datum_1995> not recognized by GRASS:

  • add Qatar_National_Datum_1995 to the datums known by GRASS
  • 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

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.

Markus M

On Tue, Dec 3, 2019 at 8:52 PM Helmut Kudrnovsky <hellik@web.de> wrote:

Markus Neteler wrote

Hi,

On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei <

zoltans@.co

> wrote:

Hi,
I’m using EPSG 2932 and QGIS etc all are OK with it.

I have installed everything using OSGeo4W64, so how come Grass does not
use the same projections database (as QGis etc)?
(and please can I have a pointer as to how to introduce EPSG 2932 to
Grass 7.8.1)

Thanks and regards,
Zoltan

******** Running against:
D:\GDBroad_tiles\shp_road_poly_few\23103770_road_poly.shp
WARNING: Datum
<Qatar_National_Datum_1995>
not recognised by GRASS and
no parameters found

Is it possible that you have a PROJ software version mixup?

I tried on my Linux box:

grass78 -c epsg:2932 ~/grassdata/test_2932
Starting GRASS GIS…
Creating new GRASS GIS location
<test_2932>
…
Cleaning up temporary files…

/ / __ / | / / / / __/ / /
/ / __/ // / /| | _ \ \ / / __ / / _
/ // / , / ___ |/ // / / // // / / /
_// |// |/// _///_/

Welcome to GRASS GIS 7.8.2dev (ad4836c73)
…

GRASS 7.8.2dev (test_2932):~ > g.proj -w
PROJCS[“QND95 / Qatar National Grid”,
GEOGCS[“QND95”,
DATUM[“Qatar_National_Datum_1995”,
SPHEROID[“International 1924”,6378388,297,
AUTHORITY[“EPSG”,“7022”]],

TOWGS84[-119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065],
AUTHORITY[“EPSG”,“6614”]],
PRIMEM[“Greenwich”,0,
AUTHORITY[“EPSG”,“8901”]],
UNIT[“degree”,0.0174532925199433,
AUTHORITY[“EPSG”,“9122”]],
AUTHORITY[“EPSG”,“4614”]],
PROJECTION[“Transverse_Mercator”],
PARAMETER[“latitude_of_origin”,24.45],
PARAMETER[“central_meridian”,51.21666666666667],
PARAMETER[“scale_factor”,0.99999],
PARAMETER[“false_easting”,200000],
PARAMETER[“false_northing”,300000],
UNIT[“metre”,1,
AUTHORITY[“EPSG”,“9001”]],
AXIS[“Easting”,EAST],
AXIS[“Northing”,NORTH],
AUTHORITY[“EPSG”,“2932”]]

when importing a shapefile with the same epsg projection, I get:

v.in.ogr --verbose input=D:\wd\test_epsg2932.shp
Using OGR driver ‘ESRI Shapefile/ESRI Shapefile’
WARNING: Datum <Qatar_National_Datum_1995> von GRASS nicht erkannt und keine
Parameter gefunden.
Die Projektionsinformationen des Eingabedatensatzes und der aktuellen
Location scheinen übereinzustimmen.

D:\wd>ogrinfo --version
GDAL 3.0.2, released 2019/10/28

D:\wd>ogrinfo test_epsg2932.shp -al -so
INFO: Open of test_epsg2932.shp' using driver ESRI Shapefile’ successful.

Layer name: test_epsg2932
Metadata:
DBF_DATE_LAST_UPDATE=2019-12-03
Geometry: Point
Feature Count: 1
Extent: (823205.438352, 25769.568359) - (823205.438352, 25769.568359)
Layer SRS WKT:
PROJCS[“QND95 / Qatar National Grid”,
GEOGCS[“QND95”,
DATUM[“Qatar_National_Datum_1995”,
SPHEROID[“International 1924”,6378388,297,
AUTHORITY[“EPSG”,“7022”]],
AUTHORITY[“EPSG”,“6614”]],
PRIMEM[“Greenwich”,0,
AUTHORITY[“EPSG”,“8901”]],
UNIT[“degree”,0.0174532925199433,
AUTHORITY[“EPSG”,“9122”]],
AUTHORITY[“EPSG”,“4614”]],
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,
AUTHORITY[“EPSG”,“9001”]],
AXIS[“Easting”,EAST],
AXIS[“Northing”,NORTH],
AUTHORITY[“EPSG”,“2932”]]
Data axis to CRS axis mapping: 1,2
id: Integer64 (10.0)
test: Integer (9.0)

D:\wd>testepsg test_epsg2932.prj
Validate Succeeds.
WKT[test_epsg2932.prj] =
PROJCS[“QND95 / Qatar National Grid”,
GEOGCS[“QND95”,
DATUM[“Qatar_National_Datum_1995”,
SPHEROID[“International 1924”,6378388,297,
AUTHORITY[“EPSG”,“7022”]],
AUTHORITY[“EPSG”,“6614”]],
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,
AUTHORITY[“EPSG”,“9001”]],
AXIS[“Easting”,EAST],
AXIS[“Northing”,NORTH]]

Simplified WKT[test_epsg2932.prj] =
PROJCS[“QND95 / Qatar National Grid”,
GEOGCS[“QND95”,
DATUM[“Qatar_National_Datum_1995”,
SPHEROID[“International 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]]

Old Style WKT[test_epsg2932.prj] = PROJCS[“QND95 / Qatar National
Grid”,GEOGCS[“QND95”,DATUM[“Qatar_National_Datum_1995”,SPHEROID[“International
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]]
ESRI’ified WKT[test_epsg2932.prj] =
PROJCS[“QND_1995_Qatar_National_Grid”,
GEOGCS[“GCS_QND_1995”,
DATUM[“D_QND_1995”,
SPHEROID[“International_1924”,6378388.0,297.0]],
PRIMEM[“Greenwich”,0.0],
UNIT[“Degree”,0.0174532925199433]],
PROJECTION[“Transverse_Mercator”],
PARAMETER[“False_Easting”,200000.0],
PARAMETER[“False_Northing”,300000.0],
PARAMETER[“Central_Meridian”,51.2166666666667],
PARAMETER[“Scale_Factor”,0.99999],
PARAMETER[“Latitude_Of_Origin”,24.45],
UNIT[“Meter”,1.0]]
PROJ.4 rendering of [test_epsg2932.prj] = +proj=tmerc +lat_0=24.45
+lon_0=51.2166666666667 +k=0.99999 +x_0=200000 +y_0=300000 +ellps=intl
+units=m +no_defs

maybe worth to open a ticket

best regards
Helmut

Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html

grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

8

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

9

Hi Markus,
Thanks for the interest in my problem.

Herewith g.version output:

Kind regards,
Zoltan

(Wed Dec 4 10:23:18 2019)
D:\sandbox\\tests\PERMANENT\.tmp/unknown\1608.0.py
version=7.8.1
date=2019
revision=c865432c9
build_date=2019-11-10
build_platform=x86_64-w64-mingw32
build_off_t_size=8
libgis_revision=00000
libgis_date="?"
proj=6.2.1
gdal=3.0.2
geos=3.8.0
sqlite=3.29.0
(Wed Dec 4 10:23:18 2019) Command finished (0 sec)

On 2019/12/03 20:06, Markus Neteler wrote:

On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei <zoltans@geograph.co.za> wrote:

Hi,
I'm using EPSG 2932 and QGIS etc all are OK with it.

I have installed everything using OSGeo4W64, so how come Grass does not
use the same projections database (as QGis etc)?

Please post the output of

g.version -rge

Best
Markus

--

=============================================
Zoltan Szecsei GPrGISc 0031
Geograph (Pty) Ltd.
GIS and Photogrammetric Services

P.O. Box 7, Muizenberg 7950, South Africa.

Mobile: +27-83-6004028 (WhatsApp only)
Qatar: +974 5083 2722 www.geograph.co.za

10

Hi Helmut,
Thanks for your comments.

I installed everything with OSGeo4W64, and QGIS get the EPSG:2932 but Grass not.

Perhaps I have a PATH or some other setting problem?
Perhaps let me know what search paths Grass uses for proj4, and what proj files and locations I should scan for.

Regards,
Zoltan

On 2019/12/04 01:30, Helmut Kudrnovsky wrote:

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
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

--

=============================================
Zoltan Szecsei GPrGISc 0031
Geograph (Pty) Ltd.
GIS and Photogrammetric Services

P.O. Box 7, Muizenberg 7950, South Africa.

Mobile: +27-83-6004028 (WhatsApp only)
Qatar: +974 5083 2722 www.geograph.co.za

11

On Wed, Dec 4, 2019 at 12:23 PM Zoltan Szecsei <zoltans@geograph.co.za> wrote:

Hi Helmut,
Thanks for your comments.

I installed everything with OSGeo4W64, and QGIS get the EPSG:2932 but
Grass not.

Perhaps I have a PATH or some other setting problem?
Perhaps let me know what search paths Grass uses for proj4, and what
proj files and locations I should scan for.

In this case where the EPSG code is known, there is no need to do anything but to ignore the warning. GRASS will use the EPSG code if available and passes it to PROJ when it comes to reprojection.

Markus M

Regards,
Zoltan

On 2019/12/04 01:30, Helmut Kudrnovsky wrote:

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

grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

–

=============================================
Zoltan Szecsei GPrGISc 0031
Geograph (Pty) Ltd.
GIS and Photogrammetric Services

P.O. Box 7, Muizenberg 7950, South Africa.

Mobile: +27-83-6004028 (WhatsApp only)
Qatar: +974 5083 2722 www.geograph.co.za

grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

12

Within GRASS:

g.proj epsg=2932 -p
-PROJ_INFO-------------------------------------------------
name : QND95 / Qatar National Grid
ellps : international
proj : tmerc
lat_0 : 24.45
lon_0 : 51.2166666666667
k : 0.99999
x_0 : 200000
y_0 : 300000
towgs84 : -119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065
no_defs : defined
-PROJ_EPSG-------------------------------------------------
epsg : 2932
-PROJ_UNITS------------------------------------------------
unit : meter
units : meters
meters : 1

looks fine.

Outside GRASS, using projinfo from PROJ 6.2.1:
projinfo -o PROJ -s EPSG:2932 -t EPSG:4326
Candidate operations found: 1

Operation n°1:

unknown id, Inverse of Qatar National Grid + QND95 to WGS 84 (1), 0 m, Qatar - onshore

PROJ string:
+proj=pipeline
+step +inv +proj=tmerc +lat_0=24.45 +lon_0=51.2166666666667 +k=0.99999 +x_0=200000 +y_0=300000 +ellps=intl
+step +proj=push +v_3 +step +proj=cart +ellps=intl
+step +proj=helmert +x=-119.4248 +y=-303.65872 +z=-11.00061 +rx=1.164298 +ry=0.174458 +rz=1.096259 +s=3.657065 +convention=position_vector
+step +inv +proj=cart +ellps=WGS84 +step +proj=pop +v_3
+step +proj=unitconvert +xy_in=rad +xy_out=deg +step
+proj=axisswap +order=2,1

Datum transformation happens with +proj=helmert which is the same as the towgs84 parameters in GRASS. This does not matter, GRASS does not use these parameters for reprojection, it relies on PROJ to select an appropriate datum transformation depending on the source and target CRS.

Markus M

On Mon, Dec 9, 2019 at 10:15 PM Markus Metz <markus.metz.giswork@gmail.com> wrote:

On Wed, Dec 4, 2019 at 12:23 PM Zoltan Szecsei <zoltans@geograph.co.za> wrote:

Hi Helmut,
Thanks for your comments.

I installed everything with OSGeo4W64, and QGIS get the EPSG:2932 but
Grass not.

Perhaps I have a PATH or some other setting problem?
Perhaps let me know what search paths Grass uses for proj4, and what
proj files and locations I should scan for.

In this case where the EPSG code is known, there is no need to do anything but to ignore the warning. GRASS will use the EPSG code if available and passes it to PROJ when it comes to reprojection.

Markus M

Regards,
Zoltan

On 2019/12/04 01:30, Helmut Kudrnovsky wrote:

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

grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

–

=============================================
Zoltan Szecsei GPrGISc 0031
Geograph (Pty) Ltd.
GIS and Photogrammetric Services

P.O. Box 7, Muizenberg 7950, South Africa.

Mobile: +27-83-6004028 (WhatsApp only)
Qatar: +974 5083 2722 www.geograph.co.za

grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user