关于GIS中的Python使用的方法其实可以直接参考GIS的帮助文章,利用里边的python语言直接捷星编程,这个过程其实十分简单,就是根据提供的代码进行改写,这里举一个例子,方法的目的是多值提取至点(Extract MultiValues To Points),具体的代码如下:
# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
#输入你的工作目录,注意路径的表达方式,为了简便,这个路径也是所有数据存放的路径,数据会在这个文件夹中进行运算。
env.workspace ="G:\\soil_erosion_calculate\\small_watershed"
inPointFeatures = "dem_forest.shp" #点文件的名称
inRasterList = [["forest_04.tif","aprf"],["forest_05.tif","mayf"],["forest_07.tif","julyf"],["forest_08.tif","augf"],["forest_09.tif","sepf"]] #栅格数据的名称以及栅格数据提取后在shapefile中数据属性的显示名称。这里需要注意的是名称中必须都是字符,不能出现数字,否则会报错。当然如果非要以数字命名,会有相对应的方法,这里就不再讲解了。
arcpy.CheckOutExtension("Spatial")
ExtractMultiValuesToPoints(inPointFeatures, inRasterList, "NONE") #这个里边的NONE可以用“BILINEAR”代替,具体看需要。
其实GIS中的python工具是和现有的工具箱中的设置是对应的关系,如果有哪些代码不清楚的可以对比这现有的工具进行理解。
网友评论