<div dir="ltr"><div>It doesn't surprise me that you have difficulties getting a better transformation between your WGS84-based custom CRS and the GDA94/GDA2020 CRS's. WGS84 is rather ill-defined, especially when using deprecated constructs such as "+datum=WGS84". If your data came from a RTK GNSS receiver I would expect it to already be in either GDA94 or GDA2020.</div><div><br></div><div>It is likely that you need to construct a custom transformation pipeline to achieve your goal. You can use projinfo to get more insight to the various transformations you are getting from cs2cs. <br></div><div><br></div><div>/Kristian<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, 11 Oct 2024 at 12:27, Nigel Stewart <<a href="mailto:nigel.stewart@emesent.io">nigel.stewart@emesent.io</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">> First of all I think you are being way too unrealistic about what you can use ChatGPT for. <div>>If you'd tried the command it's suggested to you, you would know that it doesn't work.</div><div>> It has given you an ill-formed CRS, as well as suggesting you use deprecated syntax.<br><br></div><div>It didn't smell good, I must admit.<br> <br>> It would be helpful if you can be more explicit about what you are trying to do. </div><div>> What EPSG-codes are you trying to get a transformation between?<br><br>In a nutshell:<br><br># to GDA94 / MGA zone 56<br><br></div><div>$ echo "0 0 0" | cs2cs +proj=tmerc +lat_0=-34.0 +lon_0=150.0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs +type=crs +to 'EPSG:28356'<br>222908.70       6233785.28 0.00<br><br># to GDA2020 / MGA zone 56<br><br>$ echo "0 0 0" | cs2cs +proj=tmerc +lat_0=-34.0 +lon_0=150.0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs +type=crs +to 'EPSG:7856'<br>222908.70       6233785.28 0.00<br><br># to GDA2020 / MGA zone 56<br><br>$ echo "222908.70       6233785.28 0.00" | cs2cs 'EPSG:28356' +to 'EPSG:7856'<br>222909.17       6233786.71 0.00<br><br>Our expectation is that GDA2020 is within a meter or two of GDA94 due to the tectonic drift of the Australian plate towards the north east.<br>Our concern is that from these custom transverse mercator projections, this difference isn't reflected.<br>We can see the Helmert transformation being applied in the third case, but not in the first two.<br>ChatGPT doesn't seem completely off track pointing to EPSG:8048 for GDA94 to GDA2020 Helmert.<br><br>Our application is RTK-based georeferencing, possibly WGS84 is not appropriate for our application.<br><br>$ cs2cs<br>Rel. 6.3.1, February 10th, 2020<br><br>(But we're seeing similar for our stack using PROJ 9.4.1)<br><br>Thanks in advance,<br><br>- Nigel Stewart</div></div>
</blockquote></div>