Блог

У меня есть два образца географических данных об осадках, каждый образец представляет собой массив с размерами как время, широта, долгота. координаты. Я применил статистический тест к этим двум выборкам, и в качестве результата я хочу получить серию, содержащую статистику тестирования (есть две статистики) в широтной сетке. Затем я хочу отобразить одно из значений статистики тестирования в широтной сетке. То, что я сделал, определяется как словарь с именем 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, которые являются ключами словаря (96×144, чтобы отобразить его в широтной сетке)?

Спасибо,

Ожидая предложений, лучшее, что я мог придумать, было:

 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))
  

Это создает матрицу 96×144 (широта-долгота), которую я затем могу построить.
Если есть какие-либо другие предложения pythoninc, я с удовольствием их прочитаю.