1.目的
虽然ArcGIS也提供了Arcpy这样的工具,但是只能运行在windows上,在Mac和Linux上是没法用的, 因此希望可以脱离ArcGIS平台通过一些开源python库完成一些简单的空间分析操作。
GIS中的开源库首推GDAL(Geographical Data Abstract Library),ogr是GDAL中用于矢量数据格式的模块。
2.主要内容
利用ogr包批量裁剪shapefile文件并保存为新的shapefile文件
3.代码
直接上代码,ogr的Layer可以通过SetAttributeFilter函数先根据属性进行筛选,表达式如"LAND_USE_T like '25%'"
import ogr
inputfilename = '/Users/yych97/data/2019-8-25-05/2019-8-25-05.shp'
clipfilename = '/Users/yych97/data/clip/clip.shp'
outputfilename = '/Users/yych97/data/2019-8-25-05/clipped.shp'
def clip(inputfilename, clipfilename, outputfilename):
'''
usage: clip a layer through an input Methodlayer like ArcMAP clip tool\n
parameters:
inputfilename type str,
clipfilename type str,
outputfilename type str
'''