<br><br><div class="gmail_quote">On Fri, Oct 19, 2012 at 7:41 AM, jcurru <span dir="ltr"><<a href="mailto:jcurru@yahoo.es" target="_blank">jcurru@yahoo.es</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hello!<br>
<br>
I have an old problem with TOWGS84 and OGRSpatialReference. Most likely it's<br>
because I'm no expert in GDAL. It is important for my project, and limits<br>
many of its capabilities. The fact is that, for unknown reasons (to me),<br>
proj4 does not do well conversions between coordinate system if there are<br>
not a TOWGS84 nodes inside the OGRSpatialReference(s).<br>
<br>
First of all: is this so? Do I really need that node? If not properly added,<br>
clearly some conversions fail.<br></blockquote><div><br></div><div>JCurru,</div><div><br></div><div>PROJ.4 should still handle reprojection for coordinate systems without</div><div>a TOWGS84 or other datum relationship to WGS84 *but* it will ignore</div>
<div>datum shifting if the source and destination coordinate system do not</div><div>have well defined relationships to WGS84.  So the results will (often)</div><div>be somewhat inaccurate.  How big an issue this is will depend on your</div>
<div>application and the typical coordinate systems you work with.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I continue: Imagine getting an OGRSpatialReference a shapefile or a TIF. To<br>
assign the node TOWGS84, what I do is look into the file "EPSG", found in<br>
many sites related to GDAL, that contains a list of EPSGs with corresponding<br>
parameters. Knowing the EPSG, I can generate the TOWGS84 node smoothly.<br>
<br>
The problem comes when I do not know the OGRSpatialReference’s EPSG,<br>
something very common: AutoIdentifyEPSG () gets nothing most of the time, at<br>
least with the systems in which I work.<br>
<br>
So what to do then? I need the TOWGS84 node but, with my knowledge, I can<br>
only get it from the EPSG.<br></blockquote><div><br></div><div>If you want to handle a datum shift properly and you don't know how</div><div>the datum of your coordinate system relates to WGS84 then you are </div><div>
essentially out of luck.  </div><div><br></div><div>Best regards,</div></div>-- <br>---------------------------------------+--------------------------------------<br>I set the clouds in motion - turn up   | Frank Warmerdam, <a href="mailto:warmerdam@pobox.com" target="_blank">warmerdam@pobox.com</a><br>
light and sound - activate the windows | <a href="http://pobox.com/~warmerdam" target="_blank">http://pobox.com/~warmerdam</a><br>and watch the world go round - Rush    | Geospatial Software Developer<br><br>