error in CRS-list??

classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|

error in CRS-list??

marc.vandereycken
Dear Sir,


I just started using Qgis- GRASS and might have found an error in the CRS-data.
See the pdf mentioned above.


Kind regards,

Marc Vandereycken.

_______________________________________________
grass-web mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-web

Error in GRASS.pdf (1M) Download Attachment
Reply | Threaded
Open this post in threaded view
|

Re: error in CRS-list??

Markus Neteler
Dear Mr Vandereycken,

On Sat, Oct 6, 2018 at 6:39 PM <[hidden email]> wrote:
>
> Dear Sir,
>
> I just started using Qgis- GRASS and might have found an error in the CRS-data.
> See the pdf mentioned above.

I have checked locally:
EPSG 31300 shown in WKT notation:

grass76 -c epsg:31370 --tmp-location --exec g.proj -w
Cleaning up temporary files...
Starting GRASS GIS...
Creating new GRASS GIS location <tmploc>...
Executing <g.proj -w> ...
PROJCS["Belge 1972 / Belgian Lambert 72",
    GEOGCS["Belge 1972",
        DATUM["Reseau_National_Belge_1972",
            SPHEROID["International 1924",6378388,297,
                AUTHORITY["EPSG","7022"]],
            TOWGS84[-106.8686,52.2978,-103.7239,0.3366,-0.457,1.8422,-1.2747],
            AUTHORITY["EPSG","6313"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4313"]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",51.16666723333333],
    PARAMETER["standard_parallel_2",49.8333339],
    PARAMETER["latitude_of_origin",90],
    PARAMETER["central_meridian",4.367486666666666],
    PARAMETER["false_easting",150000.013],
    PARAMETER["false_northing",5400088.438],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["X",EAST],
    AXIS["Y",NORTH],
    AUTHORITY["EPSG","31370"]]
Execution of <g.proj -w> finished.
Cleaning up temporary files...

and EPSG 31300 shown in PROJ notation:

grass76 -c epsg:31370 --tmp-location --exec g.proj -p
Cleaning up temporary files...
Starting GRASS GIS...
Creating new GRASS GIS location <tmploc>...
Executing <g.proj -p> ...
-PROJ_INFO-------------------------------------------------
name       : Belge 1972 / Belgian Lambert 72
datum      : bel72
ellps      : international
proj       : lcc
lat_1      : 51.16666723333333
lat_2      : 49.8333339
lat_0      : 90
lon_0      : 4.367486666666666
x_0        : 150000.013
y_0        : 5400088.438
no_defs    : defined
-PROJ_EPSG-------------------------------------------------
epsg       : 31370
-PROJ_UNITS------------------------------------------------
unit       : meter
units      : meters
meters     : 1
Execution of <g.proj -p> finished.
Cleaning up temporary files...

The code representation looks identical to the official register at
https://www.epsg-registry.org/
--> Retrieve by code: 31370

https://www.epsg-registry.org/export.htm?wkt=urn:ogc:def:crs:EPSG::31370

###################

Now the second reported code:
EPSG 31300 shown in WKT notation:

grass76 -c epsg:31300 --tmp-location --exec g.proj -w
Cleaning up temporary files...
Starting GRASS GIS...
Creating new GRASS GIS location <tmploc>...
Executing <g.proj -w> ...
PROJCS["Belge 1972 / Belge Lambert 72",
    GEOGCS["Belge 1972",
        DATUM["Reseau_National_Belge_1972",
            SPHEROID["International 1924",6378388,297,
                AUTHORITY["EPSG","7022"]],
            TOWGS84[-106.8686,52.2978,-103.7239,0.3366,-0.457,1.8422,-1.2747],
            AUTHORITY["EPSG","6313"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4313"]],
    PROJECTION["Lambert_Conformal_Conic_2SP_Belgium"],
    PARAMETER["standard_parallel_1",49.83333333333334],
    PARAMETER["standard_parallel_2",51.16666666666666],
    PARAMETER["latitude_of_origin",90],
    PARAMETER["central_meridian",4.356939722222222],
    PARAMETER["false_easting",150000.01256],
    PARAMETER["false_northing",5400088.4378],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["X",EAST],
    AXIS["Y",NORTH],
    AUTHORITY["EPSG","31300"]]
Execution of <g.proj -w> finished.
Cleaning up temporary files...

and EPSG 31300 shown in PROJ notation:

grass76 -c epsg:31300 --tmp-location --exec g.proj -p
Cleaning up temporary files...
Starting GRASS GIS...
Creating new GRASS GIS location <tmploc>...
Executing <g.proj -p> ...
-PROJ_INFO-------------------------------------------------
name       : Belge 1972 / Belge Lambert 72
datum      : bel72
ellps      : international
proj       : lcc
lat_1      : 49.83333333333334
lat_2      : 51.16666666666666
lat_0      : 90
lon_0      : 4.356939722222222
x_0        : 150000.01256
y_0        : 5400088.4378
no_defs    : defined
-PROJ_EPSG-------------------------------------------------
epsg       : 31300
-PROJ_UNITS------------------------------------------------
unit       : meter
units      : meters
meters     : 1
Execution of <g.proj -p> finished.
Cleaning up temporary files...

Also here I see the same representation as in the online epsg registry, see:
https://www.epsg-registry.org/export.htm?wkt=urn:ogc:def:crs:EPSG::31300

"https://www.epsg-registry.org/
-> code 31300

"Remarks:
Replaces 21500 (Belge 1950 / Belge Lambert 50). An alternative - Belge
1972 / Belgian Lambert 72 (code 31370) - was introduced in 2000 to use
the standard Lambert Conic Conformal (2SP) projection method (code
9802). EPSG recommends this alternative.
"

I cannot spot any differences to the EPSG registry.
Which source of information did you use?

Best regards,
Markus Neteler

--
Markus Neteler, PhD
http://www.mundialis.de - free data with free software
http://grass.osgeo.org
http://courses.neteler.org/blog
_______________________________________________
grass-web mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-web