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

问题描述

我想知道，如何绘制异常海面温度 (sst) 数据文件 = sst.mean.anom.nc 的 T 检验结果。 我写了两个代码，但我不知道这两个代码哪个效果更好。

代码1

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

代码 2

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

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）