Skip to content

Conversation

@himsve
Copy link
Contributor

@himsve himsve commented Oct 2, 2025

Cause of all the merge messes I had a new PR for #148. #148 is closed.

@rouault
Copy link
Member

rouault commented Oct 2, 2025

@himsve As this grid is about a vertical shift between 2 vertical CRS, the use of TYPE=VERTICAL_OFFSET_GEOGRAPHIC_TO_VERTICAL and geoid_undulation for the band description is not appropriate. You may recreate it with "vertoffset_grid_to_gtiff.py --type VERTICAL_TO_VERTICAL --source-crs EPSG:5776 --target-crs EPSG:5941 --interpolation-crs EPSG:4937" if you use vertoffset_grid_to_gtiff.py
If you create it manually, here's an example of another similar grid to take inspiration from:

$ gdalinfo /home/even/proj/PROJ-data/nz_linz/nz_linz_stisht1977-nzvd2016.tif
Driver: GTiff/GeoTIFF
Files: /home/even/proj/PROJ-data/nz_linz/nz_linz_stisht1977-nzvd2016.tif
Size is 49, 31
Coordinate System is:
GEOGCRS["NZGD2000",
    DATUM["New Zealand Geodetic Datum 2000",
        ELLIPSOID["GRS 1980",6378137,298.257222101004,
            LENGTHUNIT["metre",1]]],
    PRIMEM["Greenwich",0,
        ANGLEUNIT["degree",0.0174532925199433]],
    CS[ellipsoidal,2],
        AXIS["geodetic latitude (Lat)",north,
            ORDER[1],
            ANGLEUNIT["degree",0.0174532925199433]],
        AXIS["geodetic longitude (Lon)",east,
            ORDER[2],
            ANGLEUNIT["degree",0.0174532925199433]],
    ID["EPSG",4167]]
Data axis to CRS axis mapping: 2,1
Origin = (167.183333333333337,-46.483333333333334)
Pixel Size = (0.033333333333333,-0.033333333333333)
Metadata:
  TIFFTAG_IMAGEDESCRIPTION=NZVD2016 height (EPSG:7839) to Stewart Island 1977 height (EPSG:5772). Converted from stisht1977-nzvd2016.gtx (last modified at 2020/01/12)
  TIFFTAG_DATETIME=2020:02:02 00:00:00
  TIFFTAG_COPYRIGHT=Derived from work by Land Information New Zealand. Creative Commons Attribution 4.0 International licence https://creativecommons/licenses/by/4.0/
  area_of_use=New Zealand - Stewart Island
  source_crs_epsg_code=7839
  target_crs_epsg_code=5772
  TYPE=VERTICAL_OFFSET_VERTICAL_TO_VERTICAL
  AREA_OR_POINT=Point
Image Structure Metadata:
  COMPRESSION=DEFLATE
  INTERLEAVE=BAND
  PREDICTOR=3
Corner Coordinates:
Upper Left  ( 167.1833333, -46.4833333) (167d11' 0.00"E, 46d29' 0.00"S)
Lower Left  ( 167.1833333, -47.5166667) (167d11' 0.00"E, 47d31' 0.00"S)
Upper Right ( 168.8166667, -46.4833333) (168d49' 0.00"E, 46d29' 0.00"S)
Lower Right ( 168.8166667, -47.5166667) (168d49' 0.00"E, 47d31' 0.00"S)
Center      ( 168.0000000, -47.0000000) (168d 0' 0.00"E, 47d 0' 0.00"S)
Band 1 Block=49x31 Type=Float32, ColorInterp=Gray
  Description = vertical_offset
  Unit Type: metre

@rouault
Copy link
Member

rouault commented Oct 20, 2025

@himsve little ping w.r.t. my above comments

@rouault rouault merged commit 05009f0 into OSGeo:master Oct 22, 2025
1 check passed
@rouault rouault added this to the 1.24 milestone Oct 22, 2025
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants