你不是登录.您的编辑将被放置在队列中,直到完成同行评议

我们欢迎编辑,使文章更容易理解和更有价值的读者。因为社区成员会审查编辑,所以请尽量让帖子比你发现它的时候更好,例如,通过修改语法或添加额外的资源和超链接。

掩码阵列相关计算:

我有一个海洋表面温度时间序列,我想计算热带地区(30.0N到30S)节点之间的(Pearson)相关系数。在时间序列中,土地信息被掩盖。我不知道在相关性计算中如何处理蒙面数据。请帮助。我在这里使用的数据在这个链接中https://drive.google.com/file/d/1SVKQ4uBDEZOuN7_ftd5tqpF9_NGKY3pZ/view?usp=sharing

我尝试了下面的代码,但没有工作:

Temp5 = 'sst.day.mean.1983。fh5 =数据集(temp5, mode = 'r') sst5 = fh5。变量['sst'][:365] time = fh5。变量['time'][:] lat = fh5。变量['lat'][210:510][::35] #热带纬度lon = fh5。变量['lon'][::45] mar_05=[] for I in range(len(lat)): for j in range(len(lon)): for m in range(len(lat)): for n in range(len(lon)): mar_05.append(np。corrcoef (sst5 [59:90, i, j], sst5 [59:90, m, n] [0,1])) df = pd。DataFrame(data = mar_05)

回答

取消
1
  • 1
    \ begingroup美元 谢谢你的方便回复,它把我从嵌套循环中拯救了出来。之前,我计算了地球仪的相关性,没有实现lat = [..和lon =[]在嵌套循环条件中没有实现目的。谢谢你让我意识到我犯的错误。 \ endgroup美元
    - - - - - -Ruby
    3月9日10:44

Baidu
map