<div dir="ltr"><div><div><div>Steve,<br><br></div>Essentially (leaving out the double inversions), your &quot;case 2&quot; does this:<br><br>&lt;case2&gt;<br>proj=pipeline<br>step proj=tmerc lat_0=49 lon_0=-2 k=0.9996012717<br>                          x_0=400000 y_0=-100000 ellps=airy<br>step proj=vgridshift grids=OSGM15_GTX_OSGBtoGRS80.gtx<br>step proj=hgridshift grids=OSTN15_NTv2_OSGBtoETRS.gsb inv<br><br></div>But your grid is given in geographic coordinates, not transverse mercator<br>(cf. gdalinfo output [1] below), so starting your pipeline going from geo<br></div>to transverse mercator takes you way outside the grid coverage.<br><div><br></div><div>So start doing the gridshifts, then do the projection. When doing that, you<br></div><div>get this:<br><br>$ echo -3 55 100 | cct -t0 +init=OSGBtrans:case2<br></div><div>  335948.6065    567743.2288       47.5270        0.0000<br><br></div><div>/Thomas<br></div><div><br></div><div><br><br><br>[1]  gdalinfo OSGM15_GTX_OSGBtoGRS80.gtx <br><br>Driver: GTX/NOAA Vertical Datum .GTX<br>Files: OSGM15_GTX_OSGBtoGRS80.gtx<br>Size is 661, 1441<br>Coordinate System is:<br>GEOGCS[&quot;WGS 84&quot;,<br>    DATUM[&quot;WGS_1984&quot;,<br>        SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,<br>            AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],<br>        AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],<br>    PRIMEM[&quot;Greenwich&quot;,0,<br>        AUTHORITY[&quot;EPSG&quot;,&quot;8901&quot;]],<br>    UNIT[&quot;degree&quot;,0.0174532925199433,<br>        AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],<br>    AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]]<br>Origin = (-9.008333333333329,61.004166666666627)<br>Pixel Size = (0.016666666666667,-0.008333333333333)<br>Corner Coordinates:<br>Upper Left  (  -9.0083333,  61.0041667) (  9d 0&#39;30.00&quot;W, 61d 0&#39;15.00&quot;N)<br>Lower Left  (  -9.0083333,  48.9958333) (  9d 0&#39;30.00&quot;W, 48d59&#39;45.00&quot;N)<br>Upper Right (   2.0083333,  61.0041667) (  2d 0&#39;30.00&quot;E, 61d 0&#39;15.00&quot;N)<br>Lower Right (   2.0083333,  48.9958333) (  2d 0&#39;30.00&quot;E, 48d59&#39;45.00&quot;N)<br>Center      (  -3.5000000,  55.0000000) (  3d30&#39; 0.00&quot;W, 55d 0&#39; 0.00&quot;N)<br>Band 1 Block=661x1 Type=Float32, ColorInterp=Undefined<br>  NoData Value=-88.8888000000000034<br><br><br><br><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">2018-02-23 8:46 GMT+01:00 stevejking <span dir="ltr">&lt;<a href="mailto:stevejking@talktalk.net" target="_blank">stevejking@talktalk.net</a>&gt;</span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Kristian,<br>
Thanks for the reply, I&#39;ll take a look at the link.<br>
<br>
For the moment, I&#39;ve added the missing +step to case 3, and it now works :)<br>
I also removed the command line -t0&#39;s which I had in error put in the init<br>
file.<br>
<br>
Unfortunately case 2 is still giving the error. I even tried removing the<br>
breakline \&#39;s.<br>
<br>
Will look more when time allows.<br>
<br>
...<br>
<div class="HOEnZb"><div class="h5"><br>
Steve<br>
<br>
<br>
<br>
--<br>
Sent from: <a href="http://osgeo-org.1560.x6.nabble.com/PROJ-4-f3840930.html" rel="noreferrer" target="_blank">http://osgeo-org.1560.x6.<wbr>nabble.com/PROJ-4-f3840930.<wbr>html</a><br>
______________________________<wbr>_________________<br>
Proj mailing list<br>
<a href="mailto:Proj@lists.maptools.org">Proj@lists.maptools.org</a><br>
<a href="http://lists.maptools.org/mailman/listinfo/proj" rel="noreferrer" target="_blank">http://lists.maptools.org/<wbr>mailman/listinfo/proj</a><br>
</div></div></blockquote></div><br></div>