I quickly ran this in a python shell:<div><br></div><div>>>> from osgeo import osr</div><div><br></div><div>>>> wkt = 'GEOGCS["GCS_ITRF_2000",DATUM["D_ITRF_2000",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]'</div>
<div><br></div><div>>>> srs = osr.SpatialReference()</div><div><br></div><div>>>> srs.ImportFromWkt(wkt)</div><div><br></div><div>>>> srs.ExportToProj4()</div><div>'+proj=longlat +ellps=GRS80 +no_defs '</div>
<div><br></div><div>- Jamie<br><br><div class="gmail_quote">On Thu, May 13, 2010 at 1:32 PM, Stephen Woodbridge <span dir="ltr"><<a href="mailto:woodbri@swoodbridge.com">woodbri@swoodbridge.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hi all,<br>
<br>
I have a shapefile with the following in the prj file:<br>
<br>
GEOGCS["GCS_ITRF_2000",DATUM["D_ITRF_2000",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]]<br>
<br>
Anyone know how to convert this to a proj4 definition?<br>
<br>
Thanks,<br>
-Steve<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" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</blockquote></div><br></div>