如何在python中使用2d数组访问网格的单元格?

时间:2015-03-03 23:13:22

标签: python arrays numpy netcdf

我正在尝试使用python为一个NetCDF数据集将一个单元格变成一个值。我是python的新手,我无法理解我在这里做错了什么。我有一个2d数组,我通过迭代数据集dataArr给出了4个点。 我使用Lerp方法进行计算。不知何故,当我传递3个值时,它似乎认为它超过3。 我在这做错了什么?谢谢你。

这是我得到的错误>> self.Lerp((dataArr [i] [j]),(dataArr [i + 1] [j]),fraction)正好取3个参数(给定4个)

def compute(self,varval):      
        vars=self.data.variables
        for var in vars:
            if var==varval:
                ntimes, ny, nx=vars[var].shape #inherit the method above.
        print(ntimes, ny, nx)
        #create the old computational grid.
        computational_grid=np.zeros((ny,nx),dtype=int) 
        fraction=.5 
        newnx,newny =(nx*fraction,ny*fraction)
        new_computational_grid=np.zeros((newny,newnx),dtype=int)
        phy_value_arr=self.get_data(varval)
        t=10 #send this t value with coords
        dataArr=self.data.variables['tos'][t]     
        for i in range(0,(ny-1),1):
            for j in range(0,(nx-1),1):
               a=self.Lerp((dataArr[i][j]),(dataArr[i+1][j]),fraction)
               b=self.Lerp((dataArr[i][j]),(dataArr[i+1][j]),fraction)
               self.tempY.append(self.Lerp(a,b,fraction))
       tempY.reshape(newnx,newny)
       pcolormesh(self.tempY)
       colorbar()

def Lerp( _a, _b, _t) :
      return _a+(_b-_a)*_t

1 个答案:

答案 0 :(得分:1)

你打电话给Lerp就像是:

self.Lerp(a,b,fraction)

根据定义,它会将self作为第一个参数。因此,你将总共有四个。但是,您将Lerp定义为

def Lerp( _a, _b, _t) : #<- this takes only three arguments

我认为应该是:

def Lerp(self, _a, _b, _t) : #<- this takes four, and self is first one