[gdal-dev] Which Proj.4 transforms are available in GDAL?
Even Rouault
even.rouault at spatialys.com
Tue Sep 5 11:50:08 PDT 2017
> Except that they have a scale factor with LCC 2SP, but it gets silently
> dropped by GDAL/PROJ.4.
Which is OK. The LCC 2SP formulation just needs the 2 standard parallels and the latitude of
origin. Basically when GDAL imports a ESRI WKT and sees that there's 2 standard parallel, it
considers it is a LCC_2SP. The scale factor is redundant there. proj.4 doesn't need it when the
2 standard parallels are provided; You need either the 2 standard parallels, or 1 standard
parallel + the scale factor, and with some maths you can go from one form to the equivalent
other one (the importFromESRI() method could probably be enhanced to check the
consistency of the 4 parameters when present and emit a warning in case they are not
consistent. Exercice left to contributors)
Actually looking at proj.4 code, you must *not* include a +k_0 (other than 1) when using the
+lat_1 +lat_2 formulation, otherwise the scale factor would be applied twice.
> ... and that looks nothing similar to the Ocotepeque datum shift. They
> are about 400 meters apart.
Well that's expected if NAD27 and Octopeque are indeed 2 different datums: there must be
a shift between. The EPSG database has also shifts from Octopeque to NAD27. Apparently
they also offer the datasets in WGS 84, so you can perhaps check which hypothesis is right
from that.
--
Spatialys - Geospatial professional services
http://www.spatialys.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20170905/3f77a473/attachment.html>
More information about the gdal-dev
mailing list