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

紫藤萝

西红柿毛毛雨 发表于 2022-9-27 16:34
average()函数中axis所表示在元组中指定的所有轴上执行平均，而不是像以前那样在单个轴或所有轴上执行， ...

还是没懂

紫藤萝

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

如果这里的axis=1的话（sst_mean=np.average(area,axis=1,weights=wgts, returned=False)）会报错：ValueError: Length of weights not compatible with specified axis.

紫藤萝

紫藤萝 发表于 2023-1-1 00:53
如果这里的axis=1的话（sst_mean=np.average(area,axis=1,weights=wgts, returned=False)）会报错：Value ...

前提是数据的格式是lat x lon

西红柿毛毛雨

紫藤萝 发表于 2023-1-1 01:01
前提是数据的格式是lat x lon

我用的是NCEP的NC数据，这样计算应该没问题

puck66

是否有比较过直接mean(lat,lon)的结果差异有多大？

Pure

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

就是经向才是lat啊，怎么会对lon那一维进行不同权重的平均呢？

西红柿毛毛雨

Pure 发表于 2023-3-28 21:21
就是经向才是lat啊，怎么会对lon那一维进行不同权重的平均呢？

我又重新看了一下我原代码的数据，我的数据是三维的，所以axis=1，就是对lat求平均

西红柿毛毛雨

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

最近又有人问我这个问题，我看了一下我的数据，是三维的，包括时间，纬度，经度~~~所以第一次求平均时候axis=1，就是对lat求加权平均，然后第二次求平均，axis=1是对经度求平均，最后剩下时间维，得到一串时间序列。是这样的

西红柿毛毛雨

puck66 发表于 2023-2-7 23:10
是否有比较过直接mean(lat,lon)的结果差异有多大？

中高纬度地区差异还是比较明显的

wangzhe1993

area是什么呢