8000 geodetic CRS not found with cs2cs · Issue #1982 · OSGeo/PROJ · GitHub
[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Skip to content

geodetic CRS not found with cs2cs #1982

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Closed
snowman2 opened this issue Feb 25, 2020 · 4 comments
Closed

geodetic CRS not found with cs2cs #1982

snowman2 opened this issue Feb 25, 2020 · 4 comments
Labels

Comments

@snowman2
Copy link
Contributor

From: pyproj4/pyproj#549

Example of problem

cs2cs 'PROJCS["PCS_ITRF2000_TM",GEOGCS["ITRF2000",DATUM["International_Terrestrial_Reference_Frame_2000",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6656"]],PRIMEM["Greenwich",0],UNIT["Degree",0.0174532925199433],AUTHORITY["EPSG","4656"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",38],PARAMETER["central_meridian",127.5],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",1000000],PARAMETER["false_northing",2000000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH]]' 'EPSG:4326'
proj_create_operations: geodeticCRS not found
Rel. 6.3.1, February 10th, 2020
<cs2cs>: 
cannot initialize transformation
cause: generic error of unknown origin
program abnormally terminated

Environment Information

  • PROJ version 6.3.1
  • Ubuntu 18.04

Installation method

  • conda
@jorisvandenbossche
Copy link
Contributor

And as some additional context: the above WKT is what is produced by GDAL from reading a file with a slightly different WKT:

PROJCS["PCS_ITRF2000_TM",GEOGCS["GCS_ITRF_2000",DATUM["D_ITRF_2000",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",1000000.0],PARAMETER["False_Northing",2000000.0],PARAMETER["Central_Meridian",127.5],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",38.0],UNIT["Meter",1.0]]

With the above one, transformations work. But with the one from GDAL (through fiona) not.

@snowman2
Copy link
Contributor Author
snowman2 commented Feb 25, 2020

Looks like it worked in PROJ 6.1.1:

$ echo '197628.390 452882.814' | cs2cs 'PROJCS["PCS_ITRF2000_TM",GEOGCS["ITRF2000",DATUM["International_Terrestrial_Reference_Frame_2000",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6656"]],PRIMEM["Greenwich",0],UNIT["Degree",0.0174532925199433],AUTHORITY["EPSG","4656"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",38],PARAMETER["central_meridian",127.5],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",1000000],PARAMETER["false_northing",2000000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH]]' 'EPSG:4326'
23d50'18.824"N	119d38'16.144"E 0.000
$ cs2cs
Rel. 6.1.1, July 1st, 2019
usage: cs2cs [-dDeEfIlrstvwW [args]] [+opt[=arg] ...]
                   [+to +opt[=arg] ...] [file ...]

@snowman2
Copy link
Contributor Author

And in PROJ 6.2.1:

$ echo '197628.390 452882.814' | cs2cs 'PROJCS["PCS_ITRF2000_TM",GEOGCS["ITRF2000",DATUM["International_Terrestrial_Reference_Frame_2000",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6656"]],PRIMEM["Greenwich",0],UNIT["Degree",0.0174532925199433],AUTHORITY["EPSG","4656"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",38],PARAMETER["central_meridian",127.5],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",1000000],PARAMETER["false_northing",2000000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH]]' 'EPSG:4326'
23d50'18.824"N	119d38'16.144"E 0.000
$ cs2cs
Rel. 6.2.1, November 1st, 2019
usage: cs2cs [-dDeEfIlrstvwW [args]] [+opt[=arg] ...]
                   [+to +opt[=arg] ...] [file ...]

rouault added a commit to rouault/PROJ that referenced this issue Feb 25, 2020
rouault added a commit that referenced this issue Feb 25, 2020
createOperations(): be robust to a GeographicCRS having a wrong ID attached to it (fixes #1982)
@snowman2
Copy link
Contributor Author

Thanks @rouault!

rouault added a commit that referenced this issue Feb 25, 2020
rouault added a commit that referenced this issue Feb 25, 2020
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

2 participants
0