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."]]
