python怎么实现计算区域平均时加上面积权重

西红柿毛毛雨

求助各路大神指导，计算区域平均时，说是还要算面积加权平均，请问python有直接的函数或方法可以实现吗？？？

西红柿毛毛雨

louke 发表于 2021-10-20 21:37
同问，我也想了解这些

#计算面积权重
lat=f['latitude'][60:90]#纬度区域
lat = np.array(lat)
coslat = np.cos(np.deg2rad(lat)) #np.deg2rad将角度转化为弧度，np.cos求某个函数的cos值
wgts =np.array(coslat)
sst_mean=np.average(area,axis=1,weights=wgts, returned=False)#“对纬度进行加权平均，我的原始数据是三维的，axis=1就是对lat进行加权平均”
sst_mean=np.nanmean(sst_mean,axis=1)  ##第二次平均axis=1，是对lon进行平均，因为axis=0是我的时间维，最后剩下时间维，得到一串时间序列

这是我能想到的笨办法~~~

louke

同问，我也想了解这些

cnzkun

同问，我也想知道

紫藤萝

西红柿毛毛雨 发表于 2021-10-21 19:07
#计算面积权重
lat=f['latitude'][60:90]#纬度区域
lat = np.array(lat)

为什么要对axis=1进行加权呢，不应该对纬向加权吗，就是axis=0.
另外，为什么要计算两次平均，np.average和np.mean。谢谢了

紫藤萝

np.average可以处理缺测值吗

西红柿毛毛雨

紫藤萝 发表于 2022-9-8 11:33
为什么要对axis=1进行加权呢，不应该对纬向加权吗，就是axis=0.
另外，为什么要计算两次平均，np.averag ...

首先面积权重，主要是纬度的差异，也就是经向上的差别，所以axis=1，来算经向上的面积权重。
计算两次平均，是因为np.average（）函数只能对一个轴做加权平均，意思就是先对纬度经向加权平均后，再做一次平均，得到最后你想要的区域平均值。~~这个办法有点笨，但也不知道其他办法

西红柿毛毛雨

紫藤萝 发表于 2022-9-8 11:36
np.average可以处理缺测值吗

应该不行吧，不太清楚，我在带入函数时，已经处理过数据了

紫藤萝

西红柿毛毛雨 发表于 2022-9-13 13:41
首先面积权重，主要是纬度的差异，也就是经向上的差别，所以axis=1，来算经向上的面积权重。
计算两次平 ...

第二次平均应该是axis=0吧

西红柿毛毛雨

紫藤萝 发表于 2022-9-22 14:01
第二次平均应该是axis=0吧

average()函数中axis所表示在元组中指定的所有轴上执行平均，而不是像以前那样在单个轴或所有轴上执行，所以第一次加权平均后得到的是一列一维数据，所以第二次用mean（）计算平均时还是axis=1