[Proj] Dutch correction grid
j.l.h.hartmann at uva.nl
Tue May 3 08:34:48 EST 2011
(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
>> 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.
> Driver: GTiff/GeoTIFF
> Files: dxy.tif
> Size is 400, 400
> Coordinate System is:
> SPHEROID["Bessel 1841",6377397.155,299.1528128000008,
> Origin = (3.386235000000000,50.762069400000001)
> Pixel Size = (0.009627452750000,0.007066402750000)
> Image Structure Metadata:
> 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
More information about the Proj