我有两个降水地理数据样本,每个样本都是一个数组,维度为时间、纬度、经度坐标。我对这两个样本进行了统计测试,我想得到的结果是一个包含lat-lon网格上的测试统计数据(有两个统计数据(的系列。然后,我想在lat-lon网格上绘制一个测试统计值。我所做的是定义为一个名为test的字典,test={},然后我应用了如下测试:
for ii in range(96):
for jj in range(144):
test[ii,jj]=stats.ks_2samp(year_avgprec1[:,ii,jj],year_avgprec2[:,ii,jj])
测试的结果是形式的字典
dict_keys([(0, 0), (0, 1), (0, 2), (0, 3), ... up to (96,144)
和
dict_values([Ks_2sampResult(statistic=0.25, pvalue=0.571336004933722), Ks_2sampResult(statistic=0.25, pvalue=0.571336004933722), Ks_2sampResult(statistic=0.25, pvalue=0.571336004933722),...
如何获取我的p值并关联作为字典关键字的维度lat-lon(96x144,以便在lat-lon网格上绘制(?
谢谢,
在等待建议时,我能想出的最好的办法是:
matrix = []
for item in test:
print(item,round(test[item][0],5))
matrix.append(test[item][0])
x = np.reshape(matrix, (96, 144))
print(x.shape)
print(type(x))
这产生了一个96x144(lat-lon(矩阵,然后我可以绘制它。如果还有其他关于蟒蛇的建议,我很乐意阅读。