<div dir="ltr"><div>Assuming that it does not happen in the main area of Ireland...I guess it is because it is outside of the area of usage of the datum transformation <a href="https://epsg.org/transformation_1641/TM65-to-WGS-84-2.html">https://epsg.org/transformation_1641/TM65-to-WGS-84-2.html</a></div><div>Probably (this is just a guess) gdalwarp is not taking it into consideration, making a ballpark transformation between TM65 and WGS84, while ogr2ogr doesn't.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, 8 Nov 2023 at 13:35, Sam Hall via gdal-dev <<a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</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">I have some data that covers the whole of the Republic of Ireland, including the Blasket Islands off the West Coast, situated around 24439, 95540 in EPSG:29902. These islands straddle the western boundary of the area of use of the coordinate system -10.56° W.<br>
<br>
When transforming a raster dataset from EPSG:29902 to EPSG:4326 using gdalwarp at gdal >= 3.5.2, there seems to be a roughly 50 m shift in the data at the area of use boundary, while if I use ogr2ogr on a vector copy of this data, I don't see this shift.<br>
<br>
I have a reproducible example of this at <a href="https://gitlab.com/Sam.Hall1/29902_gdalwarp_vs_ogr2ogr.git" rel="noreferrer" target="_blank">https://gitlab.com/Sam.Hall1/29902_gdalwarp_vs_ogr2ogr.git</a>. In this example I created a rectangular polygon with interpolated points and saved it to a shapefile, then used gdal_rasterize to create a GeoTiff. gdalwarp gives the desired result using a gdal 3.5.1 docker image. The full commands are in the docker-compose.yml in the repository but the parameter's I'm using are:<br>
<br>
gdalwarp -t_srs EPSG:4326 -tr 0.000045 0.000045 -tap -overwrite input.tif output.tif<br>
ogr2ogr -t_srs EPSG:4326 output.shp input.shp<br>
<br>
I've also seen similar results when transforming from EPSG:27700 to EPSG:4326 in Jersey and Guernsey in the English Channel but I'd have to use a much earlier version of gdal to reproduce the desired result.<br>
<br>
Is there a gridshift file or any additional parameters that I could use to get the output to look more continuous across the area of use line? <br>
<br>
Sam<br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</blockquote></div>