python +气象+ T显著性检验时的困扰

Question

我最近使用了MERRA-2的气溶胶再分析数据，需要做T显著性检验的过程中遇到了麻烦。

t,p = ttest_ind(aod_volc,aod_clim,equal_var=False)

我用ttest_ind输入的两个场的维度分别是aod_volc（41，361，576）和aod_clim（1，361，576）。出来的P值打印一看，发现都是mask为True的值，画图的时候一个点也没有。

masked_array(
  data=[[--, --, --, ..., --, --, --],
        [--, --, --, ..., --, --, --],
        [--, --, --, ..., --, --, --],
        ...,
        [--, --, --, ..., --, --, --],
        [--, --, --, ..., --, --, --],
        [--, --, --, ..., --, --, --]],
  mask=[[ True,  True,  True, ...,  True,  True,  True],
        [ True,  True,  True, ...,  True,  True,  True],
        [ True,  True,  True, ...,  True,  True,  True],
        ...,
        [ True,  True,  True, ...,  True,  True,  True],
        [ True,  True,  True, ...,  True,  True,  True],
        [ True,  True,  True, ...,  True,  True,  True]],
  fill_value=1e+20,
  dtype=float64)

请问这种情况是是一个点也没有通过显著性检验吗？
【正常情况下，mask里应该有False；但是我的都是True，不知道是不是因为这个原因导致画不出图】

之后，我又发现不管是

hatches=['...', None]

还是

hatches=['...', '...']

都是没有点的，我认为应该是P数组有问题，但是没有报错。完全不知道问题出在哪。请问各位有什么办法吗？已经困扰好几天了，在线等待解惑，这对我真的很重要！

Answer 1

补充：当两个场的维度是（41，361，576）和（1，361，576），得到的p是无效数组
但是当我调整了一下变量，即增加了一年，维度变为（41，361，576）和（2，361，576），得到的p是有效的，可以进行画图的。但是这样的调整并不是我想要的，我想要的输入场维度仍然是（1，361，576）。
我猜想：（1，361，576）三维数组是不是等同于（361，576）二维数组，因此使用ttest_ind方法失败，但是我想不到任何替代方法了。
想请教更好的间接和方法，谢谢大家了。