Python+OGR
文章平均质量分 58
Learner_P
山东科技大学遥感科学与技术本科毕业,现工作于西安
展开
-
Python+OGR库学习(一):读取点矢量文件属性值和坐标,并保存为TXT(一行一个要素值)
代码思路:1、导入相关库包,切换到当前文件夹2、注册驱动,打开点矢量文件,获取图层3、打开待写入TXT文件4、遍历要素:(1)获取当前要素‘ID’和‘cover’字段属性(2)获取当前点要素对应几何对象和其坐标值X,Y(3)将ID、cover、X、Y写入TXT文件(4)清除当前feature缓存并获取下一要素5、清除DataSource并关闭TXT文件关键点:遍历要素,对每一个...原创 2018-12-26 22:32:12 · 6128 阅读 · 6 评论 -
Python+OGR库学习(二):读取点矢量文件,复制特定属性值点并另存为shp文件
代码思路:1、导入相关库包,切换到当前文件夹2、注册驱动,打开点矢量文件,获取图层3、创建输出文件,并获取图层(没有属性定义)4、定义输出图层字段属性:假设已知文件所有属性字段定义(即ID和cover)(1)读取输入文件中某一要素(2)获取ID、cover字段定义(3)添加定义到输出图层(如果不知道属性表具体字段,可以参考前一篇文章获取:https://blog.csdn.net/...原创 2018-12-28 08:25:35 · 4378 阅读 · 2 评论 -
Python+OGR库学习(三):将含有字段值的TXT文件写入shpfile
代码任务已知有一个TXT文件,每行格式:country_name:X1 Y1,X2 Y2,…Xn,Yn,其中,Xn,Yn指折点坐标,把它导出为shp文件,TXT如下:代码思路1、导入相关库包,切换到当前文件夹2、注册驱动,创建矢量文件,获取图层3、创建输出文件,并获取图层(没有属性定义)4、定义输出图层字段属性:(1)TXT文件中只能得到一个字段,命名为name,设置属性(2)打...原创 2018-12-29 08:33:27 · 1915 阅读 · 1 评论 -
Python+OGR库学习(七):生成渔网
代码目的输入四角范围,指定单个网格大小,自动生成渔网代码关键点1、整个渔网的生成和数组组成一样,只不过用一网格代替一个数,还是用两个嵌套循环完成2、在1的基础上,从左上角开始,用四角范围中xmin和网格宽度控制列变化,ymax和网格高度控制行变化进行遍历代码思路()1、导入相关包2、输入渔网范围,格网宽度和高度3、计算行列数4、初始化左上角格网四角范围,创建输出文件5、遍历每一...原创 2019-01-05 18:36:20 · 3785 阅读 · 0 评论 -
Python+OGR库学习(八):关于面矢量文件的一些小操作
对于OGR库对shp文件的基本读写没有问题,整理一些对geometry简单的小操作代码功能1、面矢量转线:主要理解面由ring构成,提取ring直接写入线矢量文件就OK(或者构建多线几何体,将所有ring都添加进去,完成后可以在arcmap中用拆分多部件要素拆分)2、读取面矢量文件,输出四角多边形:获取图层四角范围,直接写入ring3、输入面矢量,输出它的凸多边形:关键在于创建多面几何体,...原创 2019-01-05 19:31:57 · 5481 阅读 · 19 评论 -
Python+OGR库碎片整理——ExecuteSQL筛选要素
整理一下SQL语句筛选shp要素,关于SQL语句详细使用参考:[OGR:SQL]https://www.gdal.org/ogr_sql.html1、获取cover是grass的要素并id降序返回,然后打印idimport ogr,osos.chdir(r'F:\Python+gdal\7weeks数据\7weeks数据\ospy_data3\ospy3_data')driver = og...原创 2019-01-01 22:03:16 · 1169 阅读 · 0 评论 -
Python+OGR库学习(六):批量投影转换(调用自定义函数和模块)
关键思路1、批量转换:定义投影转换函数,参数设置输入输出文件名,输入输出坐标系2、内部功能:针对每一个输入文件,字段属性批量定义(操作对象是layer)、字段值批量写入(操作对象是feature)代码思路1、定义.py文件,含3个函数:字段属性批量定义字段值批量写入投影转换:其中需要调用上面2个函数3、导入以上自定义模块,遍历文件夹,批量投影转换代码——投影转换模块#! usr...原创 2019-01-01 17:41:49 · 2064 阅读 · 7 评论