获取arcgis克里格插值后的每个点的插值结果

      最近在arcgis克里格插值后需要得到插值后的矩阵。找了很久才找到工具箱中的calculateZvalue,但是一次只能计算一个点,因此需要用脚本去提取。


import arcpy
arcpy.CheckOutExtension("GeoStats")


total=0
points=[]
for width in range(37578564,37636464,50):
    for row in range(3279256,3354160,50):
        point="%d"%width+"  "+"%d"%row#整数转化为字符串
        outCZV = arcpy.GACalculateZValue_ga("G:\科研立项\kriginghigher.lyr", point)
       
        addprint=point+"   "+str(outCZV)+"\n"
        points.append(addprint)
           
    myfile = open("C:\\Users\\Administrator\\Desktop\\resultofkrige.txt", "a")
    myfile.writelines( points)
    myfile.close()
    total+=1
    print "第"+str(total)
    points=[]
            
    


    

      

转载自:https://blog.csdn.net/sisterhuang/article/details/16995109

You may also like...