资源描述：
《Python与开源GIS：使用OGR创建数据集的几何形状》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、Python与开源GIS:使用OGR创建数据集的几何形状2012-05-2514:52作者:bukunu来源:本站浏览648views我要评论字号:大中小摘要：下面我们来看一下使用OGR创建数据集的更常用的情形。这一部分重点说明一下如何使用OGR创建数据集的几何形状，对于字段、属性的处理，放到后面。对于创建Geometry来讲，wkt是最直观的。wkt是最简单的字符串格式，Python提供了丰富的函数来对字符串进行…F面我们来看一下使用OGR创建数据集的更常用的情形。这一SB分重点说明一下如彳可使用OGR创建数据集的几何形
2、状,对于字段、属性的处理，放到后面。对于创建Geometry来讲，wkt是最直观的。wkt是最简单的字符串格式，Python提供了丰富的函数来对字符串进行处理。使用wkt创建矢量数据集与使用拷贝的方法创建数据集的基本步骤是一样的。首先是创建Driver,然后创建dataSource,再创建Layer#接下来就创建Feature了。这是分别创建点、线、与多边形，使用extent变量来存储多边形的四个角点，点、线的坐标根据extent变量来生成。下面先看一下生成的数据的图，可能会更加有助于理解。创建点状数据集下面是生成点状要素
3、的Python代码。一般情况不会生成f单独的点，需要使用For循环来对数组参数进行遍历,得到想要的结果。dI1fromosgeoimportogr2importos,math3driver=ogr.GetDriverByName(HESRIShapefile11)44extfile='point_demo・shp‘5point_coors=[[300,450],[750,700],[1200,450],[750,200],[750,450]]6printpoint_coors7driver=ogr.GetDriverByN
4、ame(HESRIShapefile11)8ifos.access(extfile,os.F_OK):1()driver.DeleteDataSource(extfile)11newds=driver.CreateDataSource(extfile)12laycmcw=ncwds.CrcatcLaycrC^oint^Nonc^ogr.wkbPoint)13foraainpoint_coors:14wkt='POINT(*+str(aa[0])+fromosgeoimportogrimportos,math*+str(aa
5、[1])+y15geom=ogr.CreateGeometryFromWkt(wkt)16feat=ogr.Feature(layernew.GetLayerDefn())17fcat.SctGcomctry(gcom)18layernew.CreateFeature(feat)1919newds.Destroy()创建线状数据集下面解释一下如何单独创建一个线状数据集。3driver=ogr.GetDriverByName(HESRIShapefile11)44extfile=line_demo.shp‘5point_co
6、ors=[300,450,750,700,1200,450,750,200]6printpoint_coors3driver=ogr.GetDriverByName(HESRIShapefile11)4ifos.access(extfile,os.F_OK):1()driver.DeleteDataSource(extfile)11newds=driver.CreateDataSource(extfile)12laycmcw=ncwds.CrcatcLaycr(pointNonc^ogr.wkbLincString)13
7、13wkt=LINESTRING(%f%f,%f%f,%f%f,%f%f,%f%f)'%(point_coorslO],point_coorsl1J,14point_coors[2],point_coors[3],point_coors[4],point_coors[5],15point_coors[6],point_coors[7],poinCcoors[0]5point_coors[1])16print(wkt)17geom=ogr.CreateGeometryFromWkt(wkt)18feat=ogr.Featur
8、e(layernew.GetLayerDefn())19feat.SetGeometry(geom)20layemew.CreateFeature(feat)21newds.DcstroyO注意坐标的格式：不同点的坐标之间用T分割，坐标的不同维度用空格分割。创建多边形数据集下面解释一下如何单独创建一个矢
