[Proj] Dutch correction grid
Jan Hartmann
j.l.h.hartmann at uva.nl
Fri Apr 29 15:10:50 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
> 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.
If I use this file in the +nadgrids clause of cs2cs, I don't see any
difference in the output. If I change the filename into something
non-existing, I don 't see any differences either and do not get an
error message. Any idea how I can test this further?
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
More information about the Proj
mailing list