24小时热门版块排行榜

返回列表

kainan001

铁杆木虫 (正式写手)

EPI: 1
应助: 48 (小学生)
金币: 6418.6
散金: 652
红花: 51
帖子: 869
在线: 259.8小时
虫号: 3088580
注册: 2014-03-25
性别: GG
专业: 环境化工

[求助] 关于python for Arcgis的循环问题

大致说明：我要对某个栅格文件.tif，开一个n*n的窗口，从栅格第一个像元开始算它周围（包括本身）的均值。以下是我写的代码(这里n=1，即开一个3*3的窗口)：
>>> import arcpy
... import math
... import numpy as np
... alocation='a'
... aRaster=arcpy.Raster(alocation)
... a_array=arcpy.RasterToNumPyArray(aRaster,'#','#','#',9999)
... n=1          #'用户输入窗口大小,size',这里Num=2n+1，即窗口为3*3
... rows=a_array.shape[0]
... cols=a_array.shape[1]
... amean_array=np.zeros((rows,cols))
... for i in range(1,rows-1): #从第1行到倒数第2行
...    for j in range(1,cols-1):
...       for x in range(i-n,i+n):    #3*3的窗口，则左右扩1，5*5则扩2，依次类推
...          for y in range(j-n,j+n):
...                a_array[x,y]+=pow(2*n+1,2) #这里是个步数的迭代，例如3*3的窗口那么a_array[x,y]要迭代累加9次
...                mean=a_array[x,y]/pow(2*n+1,2) #算出来的累加值并除以窗口个数，即均值
...       amean_array[i,j]=mean
... mean_raster=arcpy.NumPyArrayToRaster(amean_array)
... mean_raster.save('D:/amean.tif')

这个可以运行，但是出来的结果不对，我的tif文件a_array列表几乎每个value都在30左右，但是最后算出来的3*3窗口的均值为4左右，这个结果肯定有问题。
请大家帮我看看代码哪有问题呢？

回复此楼