<div dir="ltr"><div>Hi</div><div><br></div><div>I have a GeoTIFF in EPSG:27700, "OSGB36 / British National Grid"</div><div>For consistency reasons, I want to run gdal2tiles with an equivalent transformation than the one I will use later to transform points using a PROJ pipeline. For consistency reasons I will use that pipeline along the time to transform the points, so an update in EPSG/PROJ does not change the results (I generate the gdal2tiles at project creation time, but run the transformation on the marker points along the time).<br></div><div><br></div><div>Getting the PROJ pipeline is easy with something like</div><div><span style="font-family:monospace">projinfo -s EPSG:27700 -t EPSG:3857 --spatial-test intersects -o PROJ --single-line -q </span></div><div>(the option <span style="font-family:monospace">--spatial-test intersects</span> is very important. Otherwise you get a ballpark. I am not using grids. I can use the area of use of the GeoTIFF as well) It produces</div><div><div><span style="font-family:monospace">+proj=pipeline +step +inv 
+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000
 +ellps=airy +step +proj=push +v_3 +step +proj=cart +ellps=airy +step 
+proj=helmert +x=446.448 +y=-125.157 +z=542.06 +rx=0.15 +ry=0.247 
+rz=0.842 +s=-20.489 +convention=position_vector +step +inv +proj=cart 
+ellps=WGS84 +step +proj=pop +v_3 +step +proj=webmerc +lat_0=0 +lon_0=0 
+x_0=0 +y_0=0 +ellps=WGS84</span></div><div><span style="font-family:monospace"><br></span></div><div></div></div><div>However this pipeline does not work in gdal2tiles with the "-s" option, complaining with these errors, and placing the data in the wrong place (not very far from Null Island)<br></div><div></div><div><span style="font-family:monospace">ERROR 1: PROJ: proj_crs_get_coordinate_system: Object is not a SingleCRS<br>ERROR 1: PROJ: proj_as_wkt: PROJBasedOperation can only be exported to WKT2<br>ERROR 1: PROJ: proj_crs_get_coordinate_system: Object is not a SingleCRS<br>ERROR 1: PROJ: proj_as_wkt: PROJBasedOperation can only be exported to WKT2<br></span></div><br><div>Running projinfo on EPSG:27700 gives me the ballpark one<br></div><div><span style="font-family:monospace">projinfo EPSG:27700 -o PROJ --single-line<br>PROJ.4 string:<br>+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +units=m +no_defs +type=crs</span></div><div><br></div><div>The difference using this pipeline (without Helmert parameters) is about 100 m near London. Perfectly noticeable ;)<br></div><div><br></div><div>(I think it does not include the Helmert transformation, because the area of use of that transformation is a bit smaller than the geographic crs, EPSG:4277)<br></div><div><br></div><div>Of course I could create the command manually, but I would like to have a generic solution for multiple projects.</div><div><br></div><div>What can I do?</div><div><br></div><div>I'm using GDAL 3.0.4 (Ubuntu 20.04)<br></div><div><br></div><div>Thanks<br></div><div><div><div dir="ltr" data-smartmail="gmail_signature">.___ ._ ..._ .. . ._.  .___ .. __ . _. . __..  ... .... ._ .__<br>Entre dos pensamientos racionales <br>hay infinitos pensamientos irracionales.<br><br></div></div></div></div>