【Python】利用ogr包进行空间分析（一）：shapefile裁剪操作

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
    '''