So, this means that I'm not able to use g.proje -jf output to properly define a coordinate system right?<br><br><div class="gmail_quote">2010/6/29 Markus Neteler <span dir="ltr"><<a href="mailto:neteler@osgeo.org">neteler@osgeo.org</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hi Kim,<br>
<div class="im"><br>
On Tue, Jun 29, 2010 at 6:37 PM, Kim Besson <<a href="mailto:kimbesson1981@gmail.com">kimbesson1981@gmail.com</a>> wrote:<br>
> Hello Markus<br>
> Thanks for the reply<br>
> To do what you suggest, I would have to retrieve LOCATION's CRS EPSG code<br>
> and I'm not able to get it automatically. I only can get that PROJ4 string<br>
> by using g.proj. Is there a way to automatically get a LOCATIONS's EPSG<br>
> code?<br>
<br>
</div>That not possible (yet?) - also GDAL/OGR does not have a reverse<br>
lookup - too complicated or too many exceptions.<br>
<br>
The string<br>
<div class="im">+proj=utm +south +no_defs +zone=24 +a=6378137 +rf=298.257223563<br>
+towgs84=0.000,0.000,0.000 +to_meter=1<br>
<br>
</div>looks like<br>
<a href="http://spatialreference.org/ref/epsg/32724/" target="_blank">http://spatialreference.org/ref/epsg/32724/</a><br>
to me.<br>
<br>
cs2cs -le | grep 'WGS 84'<br>
WGS84 a=6378137.0 rf=298.257223563 WGS 84<br>
<font color="#888888"><br>
Markus<br>
</font></blockquote></div><br>