<div dir="ltr">G'Day Martin:<div><br></div><div>The GeoTools has been configured with GridShft files for NAD27 - this is not something I have personally tried before so I am not quite sure how it works.</div><div><br></div><div>I had a look with Betsey, the data being projected is in Texas (so slightly different from Cuba so the <span style="font-size:13px">TOWGS84 parameters sound legit</span>).</div><div><br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>--</div><div>Jody Garnett</div></div></div></div></div></div>
<br><div class="gmail_quote">On 24 May 2016 at 17:51, Martin Desruisseaux <span dir="ltr"><<a href="mailto:martin.desruisseaux@geomatys.com" target="_blank">martin.desruisseaux@geomatys.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello Betsy<br>
<br>
The difference that you get may depend on whether a datum shift has been<br>
applied or not. If both Proj.4 and GeoTools have applied a datum shift,<br>
then it may depend on the values in the TOWGS84 expression. There is<br>
many different possible TOWGS84 values for a transformations between<br>
NAD27 and WGS84, depending on the geographic area. If you are<br>
transforming coordinates in USA, you should have TOWGS84[-8, 160, 176]<br>
(for an error below 10 metres). If you are transforming coordinates in<br>
Canada, you should have TOWGS84[-10, 158, 187] (for an error below 20<br>
metres). The TOWGS84 parameters in the GeoTools CRS provided with your<br>
email are for Cuba (operation code EPSG:15978) except the sign of 0.526<br>
which should be -0.526.<br>
<br>
In order to compare with Proj.4, we would need to know which parameters<br>
Proj.4 has selected. It may also use the NADCON grids in USA, or NTv2<br>
grids in Canada. I do not know how to get this information from Proj.4.<br>
<br>
Alternatively you could also try Apache SIS 0.7 for coordinate<br>
transformations. It takes in account the geographic area of the points<br>
to transform, support GML and WKT version 1 and 2, tell you which<br>
parameters it selected, in which geographic area they are valid and what<br>
accuracy to expect, etc. Version 0.7 should be released this week.<br>
<br>
    Martin<br>
<br>
<br>
Le 24/05/16 à 18:11, Betsy Emmons a écrit :<br>
<div><div class="h5">> Hello,<br>
><br>
> I'm using Proj.4 v4.9.1 and GeoTools v13.x. When a WMS layer is<br>
> brought into QGIS as EPSG:4326 there is a noticeable offset compared<br>
> to the same local layer. The WMS layer source is EPSG:4267.<br>
><br>
> A WMS 4326 sample point has a coordinate of -97.2341220, 28.0611696<br>
> and the same local 4326 sample point has a coordinate of -97.2339531,<br>
> 28.0608752. This is roughly a 40m difference. Some local calculations<br>
> have shown the error margin to be <=111m.<br>
><br>
> Here is the source 4267 projection<br>
> definition: GEOGCS["GCS_North_American_1927",DATUM["D_North_American_1927",SPHEROID["Clarke_1866",6378206.4,294.9786982]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]].<br>
><br>
> Here's what geotools is using:<br>
> GEOGCS["NAD27",<br>
>   DATUM["North American Datum 1927",<br>
>     SPHEROID["Clarke 1866", 6378206.4, 294.9786982138982,<br>
> AUTHORITY["EPSG","7008"]],<br>
>     TOWGS84[2.478, 149.752, 197.726, 0.526, -0.498, 0.501, 0.685],<br>
>     AUTHORITY["EPSG","6267"]],<br>
>   PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],<br>
>   UNIT["degree", 0.017453292519943295],<br>
>   AXIS["Geodetic longitude", EAST],<br>
>   AXIS["Geodetic latitude", NORTH],<br>
>   AUTHORITY["EPSG","4267"]]<br>
><br>
> Proj.4 is really a great library! Thank you for maintaining and<br>
> contributing to it!<br>
><br>
> Thank you,<br>
> Betsy Emmons<br>
<br>
</div></div>_______________________________________________<br>
MetaCRS mailing list<br>
<a href="mailto:MetaCRS@lists.osgeo.org">MetaCRS@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/metacrs" rel="noreferrer" target="_blank">http://lists.osgeo.org/mailman/listinfo/metacrs</a></blockquote></div><br></div>