[Gdal-dev] converting from NAD 83 to NAD 83 HARN
jhayes at lizardtech.com
Thu Jul 7 13:49:57 EDT 2005
Thanks for the ideas. I found some HARN grid shift files the other day so I will give this a try in the next few days.
From: Frank Warmerdam [mailto:fwarmerdam at gmail.com]
Sent: Wednesday, July 06, 2005 5:09 PM
To: John Hayes
Cc: gdal-dev at xserve.flids.com
Subject: Re: [Gdal-dev] converting from NAD 83 to NAD 83 HARN
On 7/6/05, John Hayes <jhayes at lizardtech.com> wrote:
> I'm trying to reproject an image from "Arizona Central State Plane NAD 83" to "Arizona Central State Plane NAD 83 HARN". This should shift the image by (X=-0.100, Y=0.182) meters according to http://www.ngs.noaa.gov/cgi-bin/nadcon2.prl. PROJ4 does not seem to be applying the shift.
> Does PROJ4 support HARN corrections?
> I have the "normal" NADCOM grid shift files installed and a 2005-04-25 snapshot of GDAL.
PROJ/GDAL is not distributed with any HARN related grid shift files
as far as I know, and in any event, the OGRSpatialReference needs
to have special cases in the WKT<->PROJ.4 translation to use grid
shift files, which it does not have for HARN.
How are you specifying the HARN coordinate systems? By EPSG
I think the solution to this issue would be for GDAL to support a WKT
extension for specifying grid shift files, as opposed to the TOWGS84
method of defining a datum. I think ESRI WKT already has this feature,
so I could do it similarly. However, this isn't really high on my
priority list at this time.
If you just wanted to hack something in, you would need to look in the
gdal/ogr/ogr_srs_proj.cpp code circa line 1670 (search for
+datum=nzgd49). You would also need to read up on use of
grid shift files on the PROJ.4 web site (the +nadgrids option) and
dig up suitable HARN grid shift files.
All in all quite alot of work!
I set the clouds in motion - turn up | Frank Warmerdam, warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush | Geospatial Programmer for Rent
More information about the Gdal-dev