[Proj] Dutch correction grid

[Jan Hartmann]

(About creating a grid shift file for the Dutch RD projection):
> On 04/08/11 18:29, Frank Warmerdam wrote:
>> As input we would require an input raster with two "bands"
>> representing the
>> offset in longitude and latitude.   I'm not sure off hand what the units
>> are - possibly arcseconds?   The input file would need to be
>> georeferenced
>> in latitude and longitude.
>>
>> If your input points are sparse and irregular you could likely use
>> gdal_grid to interpolate a grid from the point data.
>>
>> My point was mainly that you don't necessarily have to write an NTv2
>> "writer" - it is sufficient to get the datum shift data into any GDAL
>> support grid format.
>>
>> Best regards,
>>
> I made a 400*400 gdal raster called dxy.tif with two bands, the first
> with offsets for longitude, the second for latitude, all in arcseconds.
>
OK, now I get the following error (gdal svn 1.9dev, released 2011/01/18)

$ gdal_translate -of NTv2 -ot Float32  dxy.tif dxy.gsb
Input file size is 400, 400
0ERROR 1: Input and output dataset sizes or band counts do not
match in GDALDatasetCopyWholeRaster()

I had a quick look at the source (gdal/frmts/raw/ntv2dataset.cpp), but 
perhaps someone could point out the error, or where I should look.

Jan

> Driver: GTiff/GeoTIFF
> Files: dxy.tif
> Size is 400, 400
> Coordinate System is:
> GEOGCS["RT38",
>       DATUM["Stockholm_1938",
>           SPHEROID["Bessel 1841",6377397.155,299.1528128000008,
>               AUTHORITY["EPSG","7004"]],
>           AUTHORITY["EPSG","6308"]],
>       PRIMEM["Greenwich",0],
>       UNIT["degree",0.0174532925199433],
>       AUTHORITY["EPSG","4308"]]
> Origin = (3.386235000000000,50.762069400000001)
> Pixel Size = (0.009627452750000,0.007066402750000)
> Metadata:
>     AREA_OR_POINT=Area
> Image Structure Metadata:
>     INTERLEAVE=PIXEL
> Corner Coordinates:
> Upper Left  (   3.3862350,  50.7620694) (  3d23'10.45"E, 50d45'43.45"N)
> Lower Left  (   3.3862350,  53.5886305) (  3d23'10.45"E, 53d35'19.07"N)
> Upper Right (   7.2372161,  50.7620694) (  7d14'13.98"E, 50d45'43.45"N)
> Lower Right (   7.2372161,  53.5886305) (  7d14'13.98"E, 53d35'19.07"N)
> Center      (   5.3117256,  52.1753499) (  5d18'42.21"E, 52d10'31.26"N)
> Band 1 Block=400x1 Type=Float64, ColorInterp=Gray
> Band 2 Block=400x1 Type=Float64, ColorInterp=Undefined
>