擅长:python、mysql、java
<p>你忘了通知大地测量系统。下面这条路行得通。在</p>
<pre><code>import ogr, sys, os
import osgeo.osr as osr
line=ogr.Geometry(ogr.wkbLineString)
line.AddPoint(10,10)
line.AddPoint(20,20)
driver=ogr.GetDriverByName('ESRI Shapefile')
ds=driver.CreateDataSource('test.shp')
srs = osr.SpatialReference()
srs.ImportFromEPSG(4326)
layer=ds.CreateLayer('test', srs, ogr.wkbLineString)
fieldDefn_=ogr.FieldDefn('id', ogr.OFTInteger)
layer.CreateField(fieldDefn_)
featureDefn=layer.GetLayerDefn()
feature=ogr.Feature(featureDefn)
feature.SetGeometry(line)
feature.SetField('id',1)
layer.CreateFeature(feature)
</code></pre>