From: Even Rouault
Sent: Thursday, March 24, 2022 10:50 AM
To: Steve Riddell ; proj
Cc: gdal-dev@lists.osgeo.org
Subject: Re: OGRCreateCoordinateTransformation()
Steve,
(adding proj mailing list since it's actually purely a PROJ topic
Does the fact that there isn’t a consensus mean
does it fix these problems?
Thanks,
steve
From: Even Rouault
Sent: Wednesday, March 23, 2022 6:06 PM
To: Steve Riddell ; gdal-dev@lists.osgeo.org
Subject: Re: OGRCreateCoordinateTransformation()
Steve,
you may try PR https://github.com/OSGeo/PROJ/pull/3126 or
https://github.com/OSGeo/PROJ/
p +inv +proj=vgridshift +grids=HT2_0.gtx +multiplier=1\
+step +proj=tmerc +lat_0=49.351346659616 +lon_0=-123.20266499149 +k=1\
+x_0=15307.188 +y_0=6540.975 +ellps=GRS80";
The call to Transform() gives: (15307.1880, 6540.9748, 308.042) -- XY matches
above, but now I get the e
Thanks very much, Even. For reporting the problem, I “packaged up” a case with
CRS’s with EPGS codes. For the actual problem I’m working, the target CRS is a
custom projection –no EPSG.
From: Even Rouault
Sent: Friday, March 18, 2022 5:05 PM
To: Steve Riddell ; gdal-dev@lists.osgeo.org
1,ID["EPSG",8611,BOUNDCRS[SOURCECRS[VERTCRS["CGVD28 height -
HT2_0",VDATUM["Canadian Geodetic Vertical Datum of
1928"],CS[vertical,1],AXIS["gravity-related
height",up,LENGTHUNIT["metre",1]],ID["EPSG",5713]]],TARGETCRS[GEOGCRS[&q
When trying to transform a projected coordinate (UTM 14) from WGS84 to
NAD83, using the OGRCoordinateTransformation class, my transformed output
point exactly matches my input point, even though I have specified valid,
non-zero +towgs84 parameters for the NAD83 coordinate system.
If I do the s
