Ошибка индекса: индекс 206893 выходит за пределы оси 0 с размером 206893, проблема с данными сетки ⇐ Python

[Anonymous]

Последние 4 дня у меня возникла проблема с расшифровкой ошибки Python:

Код: Выделить всё

IndexError: index 206893 is out of bounds for axis 0 with size 206893

при применении данных сетки и «ближайшего» метода интерполяции используются следующие строки:

Код: Выделить всё

# create a matrix where I will store the first interpolated file
tempnew    = np.ones((np.asarray(w1[0,0,:,:]).shape))*np.nan
# The lon, lat coordinate points of the original grid
lonl,latl  = np.meshgrid(lon,lat)
points     = np.vstack((np.array(lonl).flatten(),np.array(latl).flatten())).transpose()

# The values of the original file
values     = np.array([np.asarray(temp[0,0,:,:])]).flatten()

# The dimensions of the grid that I want to interpolate to
lons       = np.array(nav_lon)
lats       = np.array(nav_lat)
X,Y        = np.meshgrid(lons,lats)

# Interpolation
tempnew    = griddata(points,values, (X,Y), method = "nearest",fill_value=-3)

Размер каждой из переменных, которые я использовал выше:

Код: Выделить всё

#tempnew.shape: (728, 312) #(Dimensions of tempnew is (lats,lons))
#lat.shape: (661,) #(original latitude)
#lon.shape: (313,) #(original longitude)
#points.shape: (206893, 2)
#values.shape: (206893,)
#X.shape: (728, 312)
#Y.shape: (728, 312)

Здесь я хотел бы отметить, что исходная сетка файла представляет собой обычные данные сетки (типа A), тогда как сетка, в которую я хочу интерполировать, не является регулярной (данные сетки C)
Ошибка выглядит следующим образом:

Код: Выделить всё

In [36]: tempnew    = sp.interpolate.griddata(points,values, (X,Y), method = "nearest
...: ",fill_value=-3)
---------------------------------------------------------------------------
IndexError                                Traceback (most recent call last)
 in 
----> 1 tempnew    = sp.interpolate.griddata(points,values, (X,Y), method =
"nearest",fill_value=-3)

~/software/anaconda3/envs/mhw/lib/python3.7/site-packages/scipy/interpolate/ndgriddata.py in
griddata(points, values, xi, method, fill_value, rescale)
217     elif method == 'nearest':
218         ip = NearestNDInterpolator(points, values, rescale=rescale)
--> 219         return ip(xi)
220     elif method == 'linear':
221         ip = LinearNDInterpolator(points, values, fill_value=fill_value,

~/software/anaconda3/envs/mhw/lib/python3.7/site-packages/scipy/interpolate/ndgriddata.py in
__call__(self, *args)
79         xi = self._scale_x(xi)
80         dist, i = self.tree.query(xi)
---> 81         return self.values[i]
82
83

IndexError: index 206893 is out of bounds for axis 0 with size 206893

Заранее спасибо,
Софи

Подробнее здесь: https://stackoverflow.com/questions/609 ... 3-griddata