Skip to content
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

GEOIDMODEL support #1841

Closed
snowman2 opened this issue Jan 11, 2020 · 1 comment
Closed

GEOIDMODEL support #1841

snowman2 opened this issue Jan 11, 2020 · 1 comment

Comments

@snowman2
Copy link
Contributor

The GEOIDMODEL tag is in the WKT here:

" VERTCRS[\"NAVD88 height\",\n"
" VDATUM[\"North American Vertical Datum 1988\"],\n"
" CS[vertical,1],\n"
" AXIS[\"gravity-related height (H)\",up,\n"
" LENGTHUNIT[\"metre\",1]],\n"
" GEOIDMODEL[\"GEOID12B\"]]]";

And it is in the PROJ database:

---- Geoid models -----
INSERT INTO "geoid_model" SELECT 'GEOID99', auth_name, code FROM grid_transformation WHERE auth_name = 'EPSG' AND grid_name LIKE 'g1999%' AND deprecated = 0;
INSERT INTO "geoid_model" SELECT 'GEOID03', auth_name, code FROM grid_transformation WHERE auth_name = 'EPSG' AND grid_name LIKE 'geoid03%' AND deprecated = 0;
INSERT INTO "geoid_model" SELECT 'GEOID06', auth_name, code FROM grid_transformation WHERE auth_name = 'EPSG' AND grid_name LIKE 'geoid06%' AND deprecated = 0;
INSERT INTO "geoid_model" SELECT 'GEOID09', auth_name, code FROM grid_transformation WHERE auth_name = 'EPSG' AND grid_name LIKE 'geoid09%' AND deprecated = 0;
-- Geoid12A and Geoid12B are identical
INSERT INTO "geoid_model" SELECT 'GEOID12A', auth_name, code FROM grid_transformation WHERE auth_name = 'EPSG' AND grid_name LIKE 'g2012b%' AND deprecated = 0;
INSERT INTO "geoid_model" SELECT 'GEOID12B', auth_name, code FROM grid_transformation WHERE auth_name = 'EPSG' AND grid_name LIKE 'g2012b%' AND deprecated = 0;
INSERT INTO "geoid_model" SELECT 'GEOID18', auth_name, code FROM grid_transformation WHERE auth_name = 'EPSG' AND grid_name LIKE 'g2018%' AND deprecated = 0;

Thoughts on adding:

  1. A C-API method to query for GEOIDMODEL in PROJ DB.
  2. A C-API method to get the GEOIDMODEL from the CRS if it exists

Number 2 can currently be done via PROJ JSON:

{"$schema":"https://proj.org/schemas/v0.2/projjson.schema.json","type":"VerticalCRS","name":"NAVD88 height","datum":{"type":"VerticalReferenceFrame","name":"North American Vertical Datum 1988"},"coordinate_system":{"subtype":"vertical","axis":[{"name":"Gravity-related height","abbreviation":"H","direction":"up","unit":"metre"}]},"geoid_model":{"name":"GEOID12B"}}

Not really a big need. Currently can be done with PROJ JSON. Just mostly raising for discussion/thoughts.

@snowman2
Copy link
Contributor Author

Going to close this as it isn't really needed :)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

1 participant