Luke, you are so full of awesome. That page answers my question and has other good tips. Thank you so much.<br><br>Sent from my Verizon Wireless Phone<br><br>----- Reply message -----<br>From: &quot;Pinner, Luke&quot; &lt;Luke.Pinner@environment.gov.au&gt;<br>Date: Tue, Nov 23, 2010 7:04 pm<br>Subject: [gdal-dev] Unable to use geom.ExportToWkt() in Python ? [SEC=UNCLASSIFIED]<br>To: &quot;Gregor at HostGIS&quot; &lt;gregor@hostgis.com&gt;<br>Cc: &lt;gdal-dev@lists.osgeo.org&gt;<br><br><br>You are running into a python &quot;gotcha&quot; -<br><a href="http://trac.osgeo.org/gdal/wiki/PythonGotchas">http://trac.osgeo.org/gdal/wiki/PythonGotchas</a><br><br>Your reference to the feature is going out of scope when you use the<br>&quot;firstgeom &nbsp;= layer.GetFeature(0).GetGeometryRef()&quot; syntax. The<br>following works syntax fine:<br>feat &nbsp; &nbsp; &nbsp; = layer.GetFeature(0)<br>firstgeom &nbsp;= feat.GetGeometryRef()<br>print firstgeom.ExportToWkt()<br><br>Luke<br><br>________________________________<br><br>From: gdal-dev-bounces@lists.osgeo.org<br>[mailto:gdal-dev-bounces@lists.osgeo.org] On Behalf Of Gregor at HostGIS<br>Sent: Wednesday, 24 November 2010 1:50 PM<br>To: Frank Warmerdam<br>Cc: gdal-dev@lists.osgeo.org<br>Subject: Re: [gdal-dev] Unable to use geom.ExportToWkt() in Python ?<br><br><br>That&#39;s strange, the docs show ExportAsWkt as being a method of a<br>Geometry, as does the source code for the Python wrapper. Well, you<br>would know best...<br><br>If I do spatialref.ExportAsWkt() it gives the WKT of the SRS, right?<br>What if I want a feature&#39;s geometry? Would I need to generate it myself<br>from the geometry type and iterating over the vertices?<br><br>Sent from my Verizon Wireless Phone<br><br>----- Reply message -----<br>From: &quot;Frank Warmerdam&quot; &lt;warmerdam@pobox.com&gt;<br>Date: Tue, Nov 23, 2010 6:31 pm<br>Subject: [gdal-dev] Unable to use geom.ExportToWkt() in Python ?<br>To: &quot;Gregor at HostGIS&quot; &lt;gregor@hostgis.com&gt;<br>Cc: &lt;gdal-dev@lists.osgeo.org&gt;<br><br><br>Gregor at HostGIS wrote:<br>&gt; Hey all. I am using the OGR binding for python, to open a SHP and<br>print <br>&gt; a geometry&#39;s WKT. This works until I try to call geom.ExportToWkt() At<br><br>&gt; that time I simply get &quot;Premature end of script headers&quot; Running it<br>from <br>&gt; the command line I get a segmentation fault.<br>&gt; <br>&gt; The code snippet:<br>&gt; <br>&gt; &nbsp; &nbsp; driver &nbsp; &nbsp; = ogr.GetDriverByName(&#39;ESRI Shapefile&#39;)<br>&gt; &nbsp; &nbsp; datasource = driver.Open(shapefile,0)<br>&gt; &nbsp; &nbsp; layer &nbsp; &nbsp; &nbsp;= datasource.GetLayer()<br>&gt; &nbsp; &nbsp; spatialRef = layer.GetSpatialRef()<br>&gt; &nbsp; &nbsp; firstgeom &nbsp;= layer.GetFeature(0).GetGeometryRef()<br>&gt; &nbsp; &nbsp; wkt &nbsp; &nbsp; &nbsp; &nbsp;= firstgeom.ExportToWkt()<br><br>Gregor,<br><br>Geometries do not have the ExportToWkt() method. &nbsp;This is a method on<br>the ogr.SpatialReference class.<br><br>Perhaps try:<br> &nbsp;wkt = spatialRef.ExportToWkt()<br><br>Best regards,<br>-- <br>---------------------------------------+--------------------------------<br>------<br>I set the clouds in motion - turn up &nbsp; | Frank Warmerdam,<br>warmerdam@pobox.com<br>light and sound - activate the windows | <a href="http://pobox.com/~warmerdam">http://pobox.com/~warmerdam</a><br>and watch the world go round - Rush &nbsp; &nbsp;| Geospatial Programmer for Rent<br><br><br><br><br><br>If you have received this transmission in error please notify us immediately by return e-mail and delete all copies. If this e-mail or any attachments have been sent to you in error, that error does not constitute waiver of any confidentiality, privilege or copyright in respect of information in the e-mail or attachments.<br><br>Please consider the environment before printing this email.<br><br><br>