<div dir="ltr">I don't pretend to be an expert on this subject although we certainly do have experts at Esri. <div><br></div><div>Anyways, it's been explained to me that the exact datum transformation depends on the geography of the source data, as well as the exact time. The datum conversion I offered was given to me as the most general transformation between NAD83 and WGS84.</div><div><br></div><div>One should expect the transformed data to be slightly off from the source. This is not a worry though because it is highly unlikely that the source data capture is rated to the precision level that matches the change.</div></div><br><div class="gmail_quote"><div dir="ltr">On Mon, Feb 8, 2016 at 11:17 AM Even Rouault <<a href="mailto:even.rouault@spatialys.com">even.rouault@spatialys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Le lundi 08 février 2016 15:54:11, Daniel Fenton a écrit :<br>
> Hi Maning,<br>
><br>
> For NAD83 you need to convert the Datum to WGS84. The WKT string you show<br>
> does not include that conversion. Proj.4 is missing these in a lot of<br>
> cases.<br>
><br>
> Try this string instead:<br>
><br>
> PROJCS["NAD83_California_zone_5_ftUS",GEOGCS["GCS_North_American_1983",DATU<br>
> M["D_North_American_1983",SPHEROID["GRS_1980",6378137,298.257222101]],TOWGS<br>
> 84[-0.9956,1.9013,0.5215,0.025915,0.009426,0.011599,-0.00062],PRIMEM["Green<br>
> wich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Lambert_Conformal<br>
> _Conic"],PARAMETER["standard_parallel_1",35.46666666666667],PARAMETER["stan<br>
> dard_parallel_2",34.03333333333333],PARAMETER["latitude_of_origin",33.5],PA<br>
> RAMETER["central_meridian",-118],PARAMETER["false_easting",6561666.667],PAR<br>
> AMETER["false_northing",1640416.667],UNIT["Foot_US",0.30480060960121924]]<br>
><br>
> Hopefully that works for you. Let me know if it doesn't.<br>
<br>
(CC'ing proj mailing list)<br>
<br>
That's interesting, as indeed, proj assumes that WGS84 == NAD83, if no extra parameters<br>
are specified. Comparing the TOWGS84 values you suggest with the EPSG database,<br>
I see you suggest using EPSG:6866 coordinate transformation "ITRF2000 to NAD83(CORS96)"<br>
(actually the opposite transformation, NAD83(CORS96) to ITRF2000).<br>
<br>
Wondering if it might be appropriate to change proj.4 defaults to the above ?<br>
<br>
Although it is not clear to me which realization of WGS84 is supposed to be the one used<br>
by proj.4 as its pivot when datum transformations are involved. I guess this must depends<br>
on the semantics of the grid or towgs84 parameters used...<br>
<br>
There's also the EPSG:6864 (ITRF96 to NAD83(CORS96)) values that could be used if we rely rather<br>
on the original WGS84 realization : -0.9910, 1.9072,0.5129,0.02579,0.00965,0.01166,0<br>
<br>
<br>
Using the HTDP grids ( <a href="https://trac.osgeo.org/proj/wiki/HTDPGrids" rel="noreferrer" target="_blank">https://trac.osgeo.org/proj/wiki/HTDPGrids</a> ) might also be an option.<br>
<br>
<br>
Indeed Daniel's suggestion seems to do the trick, or at least better match<br>
other data from the site <a href="http://geohub.lacity.org/datasets/813fcefde1f64b209103107b26a8909f_0" rel="noreferrer" target="_blank">http://geohub.lacity.org/datasets/813fcefde1f64b209103107b26a8909f_0</a><br>
. I see it has also a geojson output :<br>
<a href="http://geohub.lacity.org/datasets/813fcefde1f64b209103107b26a8909f_0.geojson" rel="noreferrer" target="_blank">http://geohub.lacity.org/datasets/813fcefde1f64b209103107b26a8909f_0.geojson</a><br>
<br>
And the first feature returned is OBJECT_ID=4001 whose first vertex is :<br>
-118.276946349107,33.7920287419671<br>
<br>
If I do,<br>
ogr2ogr -f kml /vsistdout/ /vsizip/Building_Footprints.zip  -where "OBJECTID=4001"<br>
<br>
I get :<br>
-118.276934177326,33.7920239166625<br>
<br>
which is off by a bit more than 1 meter from the position of the geojson output.<br>
<br>
But If I do :<br>
<br>
ogr2ogr -f kml /vsistdout/ /vsizip/Building_Footprints.zip  -where "OBJECTID=4001" \<br>
 -s_srs "+proj=lcc +lat_1=35.46666666666667 +lat_2=34.03333333333333 \<br>
+lat_0=33.5 +lon_0=-118 +x_0=2000000.000101601 +y_0=500000.0001016002 \<br>
+ellps=GRS80 +towgs84=-0.9956,1.9013,0.5215,0.025915,0.009426,0.011599,-0.00062 +units=us-ft +no_defs" \<br>
-t_srs EPSG:4326<br>
<br>
I get:<br>
<br>
-118.276946349109,33.7920287419681<br>
<br>
which is identical to the geojson output (difference is in the 10 micrometer precision range)<br>
<br>
><br>
> Daniel Fenton<br>
> Software Engineer<br>
> Esri R&D<br>
><br>
> On Mon, Feb 8, 2016 at 6:56 AM maning sambale <<a href="mailto:emmanuel.sambale@gmail.com" target="_blank">emmanuel.sambale@gmail.com</a>><br>
><br>
> wrote:<br>
> > Hi,<br>
> ><br>
> > I'm trying to convert LA City building data from State Plane to<br>
> > EPSG:4326. I converted the data using ogr2ogr.  Visually inspecting<br>
> > using Bing imagery, I'm seeing an offset.<br>
> ><br>
> > Details of SRS info and sample images here:<br>
> > <a href="https://gist.github.com/maning/093a866fd465ea880032" rel="noreferrer" target="_blank">https://gist.github.com/maning/093a866fd465ea880032</a><br>
> ><br>
> > It is possible that Bing has the offset, but the source of Bing comes<br>
> > from LARIAC (Los Angeles Region Imagery Acquisition Consortium)<br>
> > imagery, the same source as the building footprints were traced.<br>
> ><br>
> > Any hints on the "proper" conversion from State Plane to EPSG:4326.<br>
> ><br>
> > --<br>
> > cheers,<br>
> > maning<br>
> > ------------------------------------------------------<br>
> > "Freedom is still the most radical idea of all" -N.Branden<br>
> > <a href="https://epsg4253.wordpress.com/" rel="noreferrer" target="_blank">https://epsg4253.wordpress.com/</a><br>
> > <a href="http://twitter.com/maningsambale" rel="noreferrer" target="_blank">http://twitter.com/maningsambale</a><br>
> > ------------------------------------------------------<br>
> > _______________________________________________<br>
> > gdal-dev mailing list<br>
> > <a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
> > <a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
<br>
--<br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
</blockquote></div>