[PROJ] [gdal-dev] OGRCreateCoordinateTransformation()
Even Rouault
even.rouault at spatialys.com
Thu Mar 24 14:18:53 PDT 2022
Actually my below statement is not completely true
>
> And the ABRIDGEDTRANSFORMATION[] syntax cannot be used for vertical
> transformations where the geographic CRS would not be the
> interpolation CRS.
>
https://github.com/OSGeo/PROJ/pull/3134 should improve this, by using
the approach used in the EPSG database using a PARAMETER["EPSG code for
Interpolation CRS",7886]
So something like below could be used to convey the interpolation CRS
BOUNDCRS[
SOURCECRS[
VERTCRS["unknown",
VDATUM["unknown using geoidgrids=@foo.gtx"],
CS[vertical,1],
AXIS["gravity-related height (H)",up,
LENGTHUNIT["metre",1,
ID["EPSG",9001]]]]],
TARGETCRS[
GEOGCRS["WGS 84",
DATUM["World Geodetic System 1984",
ELLIPSOID["WGS 84",6378137,298.257223563,
LENGTHUNIT["metre",1]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433]],
CS[ellipsoidal,3],
AXIS["latitude",north,
ORDER[1],
ANGLEUNIT["degree",0.0174532925199433]],
AXIS["longitude",east,
ORDER[2],
ANGLEUNIT["degree",0.0174532925199433]],
AXIS["ellipsoidal height",up,
ORDER[3],
LENGTHUNIT["metre",1]],
ID["EPSG",4979]]],
ABRIDGEDTRANSFORMATION["unknown to WGS84 ellipsoidal height",
METHOD["GravityRelatedHeight to Geographic3D"],
PARAMETERFILE["Geoid (height correction) model file","@foo.gtx",
ID["EPSG",8666]],
PARAMETER["EPSG code for Interpolation CRS",7886]]]
Even
>
> --
> http://www.spatialys.com
> My software is free, but my time generally not.
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev
--
http://www.spatialys.com
My software is free, but my time generally not.
More information about the PROJ
mailing list