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

我最近使用了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数组有问题,但是没有报错。完全不知道问题出在哪。请问各位有什么办法吗?已经困扰好几天了,在线等待解惑,这对我真的很重要!

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