代码任务
已知有一个TXT文件,每行格式:country_name:X1 Y1,X2 Y2,…Xn,Yn,其中,Xn,Yn指折点坐标,把它导出为shp文件,TXT如下:
代码思路
1、导入相关库包,切换到当前文件夹
2、注册驱动,创建矢量文件,获取图层
3、创建输出文件,并获取图层(没有属性定义)
4、定义输出图层字段属性:
(1)TXT文件中只能得到一个字段,命名为name,设置属性
(2)打开TXT,读取遍历每一行:
创建ring类型geom,分割每行内容为两部分(name和x,y值),将这行所有x,y值添加为ring坐标点
创建polygon,将含有坐标的ring添加进去构成一个polygon
创建这一行对应的一个要素,并写入name字段值
要素写入图层
清除此循环中ring、polygon及feature缓存,便于下次循环
(3)清除DataSource并关闭TXT文件
关键点
每行创建一个要素
polygon创建思路:ring-addpoint to ring -polygon–addring to polygon
split分割字段和坐标点
代码
#导入相关库包
import osgeo,gdal,os,ogr,osr
#创建输出shp
driver = ogr.GetDriverByName('ESRI Shapefile')
os.chdir(r'F:\PQ\python_exercise\7weeks数据\ospy_data2\ospy_data2')
if os.path.exists<