提取满足另一个矩阵中阈值的数据

Question

我在Matlab遇到了一些麻烦，这可能更像是一个方法论/逻辑问题。

所以我有三个netCDF文件，我将其加载到三个单独的矩阵中。这些矩阵中的两个是两个不同实验的测试运行的值（所有三个文件都是网格化的），另一个是两个实验之间的数据的p值。

我想创建一个新的netCDF文件，如果该特定网格的pvalue＆lt; = 0.05，则该文件具有两个实验之间差异的值。

我认为这很容易。这是我实施的内容：

%find the difference between experiments, averaged over the 3rd dimension
yr_diff = mean(mat_2012-mat_1980,3);
yr_diff(pvalue >= 0.05) = NaN;
final_grid = yr_diff;

然后我使用了一个相当简单的netCDF导出方案。据我所知，mat_2012和mat_1980矩阵的组织方式与pvalue矩阵相同。

这产生了一种奇怪的网格，并没有真正显示任何东西。下面的第一张图显示了差异的填充轮廓，其中区域非常重要的黑色轮廓线。第二个图显示了使用此逻辑创建的netCDF文件的填充轮廓，具有相同的重要轮廓线。

所以我的问题：逻辑是否错误，我应该采取新的方向吗？这与我的出口计划有关吗？或者这完全是另一回事吗？