Python 我想知道，如何绘制异常海平面温度的T检验结果？（SST）数据文件=SST.mean.anom.nc

我想知道，如何映射异常海表面温度（SST）数据文件=SST.mean.anom.nc的T检验结果。 我写了两个代码，但我不知道这两个代码中哪一个效果更好

代码1

代码2

meanSumaX = sumaX / len(X)
meanSumaHx = sumaHx / len(Hx)
difMeanSumaHxX = meanSumaHx - meanSumaX

ddof_MeanSumaHx = meanSumaHx.var(ddof=1)

ddof_MeanSumaX = meanSumaX.var(ddof=1)

s = np.sqrt((ddof_MeanSumaX + ddof_MeanSumaHx)/2)

t = difMeanSumaHxX / (s * np.sqrt(2/len(sst.loc['1948-03-01':'1998-03-01':12])))

# Using geophysical units. `robust` disregards outliers for colour map creation.
fig = plt.figure(5, figsize=(15., 12.))
ax = plt.axes(projection=ccrs.PlateCarree(central_longitude=0.0))
ax.coastlines()
ax.add_feature(cf.LAND)  # utilizar con internet
t.plot()
plt.title("Mean SST anomali (C) 1948-03 - 1998-03 T - test  ")
ax.gridlines(draw_labels=True)
plt.show()

meanSumaX = sumaX / len(X)
meanSumaHx = sumaHx / len(Hx)
difMeanSumaHxX = meanSumaHx - meanSumaX
std1, std2 = std(meanSumaHx, ddof=1), std(meanSumaX, ddof=1)

se1, se2 = std1 / np.sqrt(len(meanSumaHx)), std2 / np.sqrt(len(meanSumaX))

sed = np.sqrt(se1 ** 2 + se2 ** 2)

t = difMeanSumaHxX / sed

# Using geophysical units. `robust` disregards outliers for colour map creation.
fig = plt.figure(5, figsize=(15., 12.))
ax = plt.axes(projection=ccrs.PlateCarree(central_longitude=0.0))
ax.coastlines()
ax.add_feature(cf.LAND)  # utilizar con internet
t.plot()
plt.title("Mean SST anomali (C) 1948-03 - 1998-03 T - test  ")
ax.gridlines(draw_labels=True)
plt.show()