Skip to content
This repository has been archived by the owner on May 18, 2022. It is now read-only.

Add USA geoid models 1999, 2003, 2006, 2009 and 2018 #60

Merged
merged 5 commits into from
Nov 4, 2019

Conversation

rouault
Copy link
Member

@rouault rouault commented Nov 2, 2019

All of them are referenced, at least partially, by EPSG and/or LIDAR WKT1 CompoundCRS

rouault added a commit to rouault/PROJ that referenced this pull request Nov 2, 2019
…D99, 03, 06, 09 and 18. Related to OSGeo/proj-datumgrid#60 that add those grids in proj-datumgrid-north-america
@rouault
Copy link
Member Author

rouault commented Nov 4, 2019

@kbevers @hobu Does that look good to you ? After doing that, I realized that https://download.osgeo.org/proj/vdatum/ had .zip with previous grids. The names betwen those .zip and this pull requests are not always the same:

  • I've tried to align on the name of the grids referenced by the EPSG dataset as much as possible. If they refer to geoid06_ak.bin, then I create a geoid06_ak.gtx (potentially merging the 4 NOAA .bin files together)
  • For other files not referenced by EPSG, I keep the name of the upstream NOAA file and replace its .bin extension by .gtx

@kbevers kbevers added this to the datumgrid-north-america-1.3 milestone Nov 4, 2019
@kbevers
Copy link
Member

kbevers commented Nov 4, 2019

It looks good to me. I haven't tested any of the grids but I trust that you have converted them correctly.

@rouault
Copy link
Member Author

rouault commented Nov 4, 2019

I haven't tested any of the grids but I trust that you have converted them correctly.

I trust the GDAL NGSGEOID driver that has been used recently to do the job
Also cross-checking on a random point among the various GEOID realizations and egm96_15.gtx shows consistency

For CONUS

$ gdallocationinfo -wgs84 ../egm96_15.gtx -77 40
Report:
  Location: (412P,200L)
  Band 1:
    Value: -34.5047569274902

$ gdallocationinfo -wgs84 g1999u03.gtx 283 40
Report:
  Location: (1140P,1080L)
  Band 1:
    Value: -33.820499420166

$ gdallocationinfo -wgs84 geoid03_conus.gtx 283 40
Report:
  Location: (3180P,1080L)
  Band 1:
    Value: -33.8349685668945

$ gdallocationinfo -wgs84 geoid09_conus.gtx 283 40
Report:
  Location: (3180P,1080L)
  Band 1:
    Value: -33.8376998901367

$ gdallocationinfo -wgs84 g2012bu0.gtx 283 40
Report:
  Location: (3180P,1080L)
  Band 1:
    Value: -33.8595008850098

$ gdallocationinfo -wgs84 g2018u0.gtx 283 40
Report:
  Location: (3180P,1080L)
  Band 1:
    Value: -33.880199432373

Alaska

$ gdallocationinfo -wgs84 ../egm96_15.gtx 175 65
Report:
  Location: (1420P,100L)
  Band 1:
    Value: 6.03143215179443

$ gdallocationinfo -wgs84 g1999a01.gtx  175 65
Report:
  Location: (180P,420L)
  Band 1:
    Value: 6.23759984970093

$ gdallocationinfo -wgs84 g2003a01.gtx  175 65
Report:
  Location: (180P,420L)
  Band 1:
    Value: 5.28455686569214

$ gdallocationinfo -wgs84 geoid06_ak.gtx  175 65
Report:
  Location: (180P,420L)
  Band 1:
    Value: 3.61653804779053

$ gdallocationinfo -wgs84 geoid09_ak.gtx  175 65
Report:
  Location: (180P,420L)
  Band 1:
    Value: 4.35279989242554

$ gdallocationinfo -wgs84 g2012ba0.gtx  175 65
Report:
  Location: (180P,420L)
  Band 1:
    Value: 4.43069982528687

@rouault rouault merged commit 6ecbb6e into OSGeo:master Nov 4, 2019
Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants