Yes, it works if I fix the longitude sign, and make it first. Thank you
On Wed, Dec 11, 2019 at 2:56 AM Martin Desruisseaux < [email protected]> wrote: > Hello Thierry > > Indeed from the WKT below, longitude is first. I do not know yet why we > got different string (I should look again at the test code). Have you tried > with longitude first again, but fixing only the longitude sign? > > Martin > > Le 11/12/2019 à 00:50, Thierry Danard a écrit : > > Thank you > It does work using your method. > > I was wondering, I got this operation string, how do you know from this > operation string whether longitude or latitude is first? I would have > assumed longitude is first from the string. > > CoordinateOperation["AGD66 to GDA94 (11)", > SourceCRS[GeodeticCRS["AGD66", > Datum["Australian Geodetic Datum 1966", > Ellipsoid["Australian National Spheroid", 6378160.0, 298.25]], > CS[ellipsoidal, 2], > Axis["Geodetic longitude (Lon)", east], > Axis["Geodetic latitude (Lat)", north], > Unit["degree", 0.017453292519943295]]], > TargetCRS[GeodeticCRS["GDA94", > Datum["Geocentric Datum of Australia 1994", > Ellipsoid["GRS 1980", 6378137.0, 298.257222101]], > CS[ellipsoidal, 2], > Axis["Geodetic longitude (Lon)", east], > Axis["Geodetic latitude (Lat)", north], > Unit["degree", 0.017453292519943295]]], > Method["NTv2"], > Parameter["Latitude and longitude difference file", "A66 National > (13.09.01).gsb"], > OperationAccuracy[0.1], > Scope["0.1m accuracy."], > Area["Australia - Australian Capital Territory; New South Wales; > Northern Territory; Queensland; South Australia; Tasmania; Western > Australia; Victoria."], > BBox[-43.70, 112.85, -9.86, 153.69], > Remark["Replaces AGD66 to GDA94 variants 6, 7 and 10 (codes 1506 1507 > 1596). Input expects longitudes to be positive west; EPSG GeogCRS AGD66 > (code 4202) and GDA94 (code 4283) both have longitudes positive east. May > be used as tfm to WGS 84 - see code 15786."]] > > >
