from osgeo import ogr # Create test polygon ring = ogr.Geometry(ogr.wkbLinearRing) ring.AddPoint(1179091.1646903288, 712782.8838459781) ring.AddPoint(1161053.0218226474, 667456.2684348812) ring.AddPoint(1214704.933941905, 641092.8288590391) ring.AddPoint(1228580.428455506, 682719.3123998424) ring.AddPoint(1218405.0658121984, 721108.1805541387) ring.AddPoint(1179091.1646903288, 712782.8838459781) poly = ogr.Geometry(ogr.wkbPolygon) poly.AddGeometry(ring) # Create the output Driver outDriver = ogr.GetDriverByName('GeoJSON') # Create the output GeoJSON outDataSource = outDriver.CreateDataSource('test.geojson') outLayer = outDataSource.CreateLayer('test.geojson', geom_type=ogr.wkbPolygon ) # Get the output Layer's Feature Definition featureDefn = outLayer.GetLayerDefn() # create a new feature outFeature = ogr.Feature(featureDefn) # Set new geometry outFeature.SetGeometry(poly) # Add new feature to output Layer outLayer.CreateFeature(outFeature) # destroy the feature outFeature.Destroy # Close DataSources outDataSource.Destroy() from osgeo import ogr # Create test polygon ring = ogr.Geometry(ogr.wkbLinearRing) ring.AddPoint(1179091.1646903288, 712782.8838459781) ring.AddPoint(1161053.0218226474, 667456.2684348812) ring.AddPoint(1214704.933941905, 641092.8288590391) ring.AddPoint(1228580.428455506, 682719.3123998424) ring.AddPoint(1218405.0658121984, 721108.1805541387) ring.AddPoint(1179091.1646903288, 712782.8838459781) poly = ogr.Geometry(ogr.wkbPolygon) poly.AddGeometry(ring) geojson = poly.ExportToJson() print geojson
有两个选项创建从一个几何GeoJSON。你可以创建一个新文件或简单的出口GeoJSON几何JSON和打印。这两个选项都在下面解释。
Copyright © Since 2014.
开源地理空间基金会中文分会
吉ICP备05002032号
Powered by TorCMS