Запись определенного интеграла по времени трехмерной матрицы - PullRequest
Новая
       44

Запись определенного интеграла по времени трехмерной матрицы

0 голосов
/ 08 января 2020

У меня есть 3D-матрица, состоящая из времени, у и х. Я хочу интегрировать по определенному временному диапазону для каждой отдельной точки сетки, как бы я go сделал это?

Существует значение для каждого временного и пространственного местоположения (x и y), но я просто хотите интегрироваться с течением времени, но сохранить пространственные расположения одинаковыми.

Это то, что у меня пока есть, но есть ошибки: 

#PARSEVAL: integrate the square of the time series from 0 to the nyquist frequency

sla_squared = np.square(sla_standard) #where sla_squared is the 3D matrix with time, longitude, and latitude

t = np.linspace(0,177,710)
PD = np.zeros(sla_squared.shape)

for x in range(lon.size):
    for y in range (lat.size):
       PD = (np.trapz(sla_squared(t,x,y),t))
       print (PD)

Error: 'numpy.ndarray' object is not callable
...