The +nadgrid command expects a NTv2 file, with extension gsb, not a tiff file. You must use gdal_translate to convert it to the NTv2 format. I think that&#39;s what Frank Warmerdam referred in a previous email.<br><br>Regards<br>
<br>Jose&#39; A. Gonçalves<br><br><br><br><br><div class="gmail_quote">2011/4/29 Jan Hartmann <span dir="ltr">&lt;<a href="mailto:j.l.h.hartmann@uva.nl">j.l.h.hartmann@uva.nl</a>&gt;</span><br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
(About creating a grid shift file for the Dutch RD projection):<br>
<div class="im"><br>
On 04/08/11 18:29, Frank Warmerdam wrote:<br>
&gt; As input we would require an input raster with two &quot;bands&quot;<br>
&gt; representing the<br>
&gt; offset in longitude and latitude.   I&#39;m not sure off hand what the units<br>
&gt; are - possibly arcseconds?   The input file would need to be<br>
&gt; georeferenced<br>
&gt; in latitude and longitude.<br>
&gt;<br>
&gt; If your input points are sparse and irregular you could likely use<br>
&gt; gdal_grid to interpolate a grid from the point data.<br>
&gt;<br>
&gt; My point was mainly that you don&#39;t necessarily have to write an NTv2<br>
&gt; &quot;writer&quot; - it is sufficient to get the datum shift data into any GDAL<br>
&gt; support grid format.<br>
&gt;<br>
&gt; Best regards,<br>
&gt;<br>
</div>I made a 400*400 gdal raster called dxy.tif with two bands, the first<br>
with offsets for longitude, the second for latitude, all in arcseconds.<br>
If I use this file in the +nadgrids clause of cs2cs, I don&#39;t see any<br>
difference in the output. If I change the filename into something<br>
non-existing, I don &#39;t see any differences either and do not get an<br>
error message. Any idea how I can test this further?<br>
<br>
Jan<br>
<br>
Driver: GTiff/GeoTIFF<br>
Files: dxy.tif<br>
Size is 400, 400<br>
Coordinate System is:<br>
GEOGCS[&quot;RT38&quot;,<br>
     DATUM[&quot;Stockholm_1938&quot;,<br>
         SPHEROID[&quot;Bessel 1841&quot;,6377397.155,299.1528128000008,<br>
             AUTHORITY[&quot;EPSG&quot;,&quot;7004&quot;]],<br>
         AUTHORITY[&quot;EPSG&quot;,&quot;6308&quot;]],<br>
     PRIMEM[&quot;Greenwich&quot;,0],<br>
     UNIT[&quot;degree&quot;,0.0174532925199433],<br>
     AUTHORITY[&quot;EPSG&quot;,&quot;4308&quot;]]<br>
Origin = (3.386235000000000,50.762069400000001)<br>
Pixel Size = (0.009627452750000,0.007066402750000)<br>
Metadata:<br>
   AREA_OR_POINT=Area<br>
Image Structure Metadata:<br>
   INTERLEAVE=PIXEL<br>
Corner Coordinates:<br>
Upper Left  (   3.3862350,  50.7620694) (  3d23&#39;10.45&quot;E, 50d45&#39;43.45&quot;N)<br>
Lower Left  (   3.3862350,  53.5886305) (  3d23&#39;10.45&quot;E, 53d35&#39;19.07&quot;N)<br>
Upper Right (   7.2372161,  50.7620694) (  7d14&#39;13.98&quot;E, 50d45&#39;43.45&quot;N)<br>
Lower Right (   7.2372161,  53.5886305) (  7d14&#39;13.98&quot;E, 53d35&#39;19.07&quot;N)<br>
Center      (   5.3117256,  52.1753499) (  5d18&#39;42.21&quot;E, 52d10&#39;31.26&quot;N)<br>
Band 1 Block=400x1 Type=Float64, ColorInterp=Gray<br>
Band 2 Block=400x1 Type=Float64, ColorInterp=Undefined<br>
<div><div></div><div class="h5">_______________________________________________<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" target="_blank">http://lists.maptools.org/mailman/listinfo/proj</a><br>
</div></div></blockquote></div><br>