最近教师要求实现大量的二值图转换向量文件,发现用python调用arcgis的arcpy库非常方便,就记录下来。
我的数据都是0,1分布的二值图。 试着实际使用看看需求来修正吧。 提取值为1的区域,因此代码中的网格代码值为1。
#--coding:utf-8----#arcgis中导入的pythonIDE是ArcGIS附带的python.exeimportarcpyfromarcpyimportenvimportenvimportortosimporton Raster_to_polygonnew2_veg'#二值图路径output _ path=r ' e 3360 _ jw _ work _ GIS _ raster GIS raster _ to result/'#gridcode=1的shp file list=OS.list dir (env.workspace ) TIF _ list=[ ] forfileinfilelist 3360 ifos.path='.shp ' shp _ path=output _ pathshpfileresultpath=result _ pathshpfilearcpy.rastertopolygon _ conversion (' value #这里可以向你修改(round(time.time )-stime,2 ) )据说可以直接安装Arcpy库,我在这里直接安装arcgis时附带的python解释器那里附带了Arcpy。