<div dir="ltr">When I load both of them into QGIS they align pretty well, with the output being slightly smaller than the input file.<div><br></div><div><img src="cid:ii_jhgf4nst0" alt="Screenshot from 2018-05-21 16-40-32.png" width="562" height="298"></div><div><br></div><div>I suspect this is due to  the pixel size changing, but I'm not sure why it is so large as QGis is using GDAL internally to reproject the input to 3857 here.</div><div><br></div><div>Ian<br><br></div></div><br><div class="gmail_quote"><div dir="ltr">On Mon, 21 May 2018 at 16:17, Ian Hannah <<a href="mailto:ihannah@meniscus.co.uk">ihannah@meniscus.co.uk</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I am using the C# GDAL interfact to create a GeoTiff file ( test.tiff<br>
<<a href="http://osgeo-org.1560.x6.nabble.com/file/t383179/test.tiff" rel="noreferrer" target="_blank">http://osgeo-org.1560.x6.nabble.com/file/t383179/test.tiff</a>>  ).<br>
<br>
This file has the EPSG:27700 projection set within it and a GeoTransform as<br>
follows:<br>
<br>
                    double[] transform = new double[6];<br>
                    transform[0] = -24000;<br>
                    transform[1] = 1000;<br>
                    transform[2] = 0;<br>
                    transform[3] = 996000;<br>
                    transform[4] = 0;<br>
                    transform[5] = -1000;<br>
                    dataset.SetGeoTransform(transform);<br>
<br>
So the GeoTiff image goes from -24000 to 679000 in the Eastings and 996000<br>
to -14000 in the Northings.<br>
<br>
I then use the GDAL warp application using this command:<br>
<br>
gdalwarp -s_srs EPSG:27700 -t_srs EPSG:3857 -overwrite test.tiff output.tiff<br>
<br>
This creates me the output  output.tiff<br>
<<a href="http://osgeo-org.1560.x6.nabble.com/file/t383179/output.tiff" rel="noreferrer" target="_blank">http://osgeo-org.1560.x6.nabble.com/file/t383179/output.tiff</a>>  . The large<br>
square in the middle should be over Leeds when placed onto a map but it is<br>
not. The right hand side does not extend out far enough and I have no idea<br>
why.<br>
<br>
I just want to know if anyone has had any experience in doing this<br>
transformation or can see anything that I am not doing correctly.<br>
<br>
Thanks<br>
Ian<br>
<br>
<br>
<br>
<br>
--<br>
Sent from: <a href="http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html" rel="noreferrer" target="_blank">http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html</a><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></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature">Ian Turton<br></div>