请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 1302|回复: 2

[经验总结] 月平均求季节平均权重

[复制链接]

新浪微博达人勋

发表于 2023-7-25 00:10:56 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 一大碗年糕 于 2023-7-25 00:18 编辑

试了下geocat的month_to_season(原ncl month_to_season)函数,与自己计算的两种重考虑权重的月平均均有差异。说明该函数确实没考虑不同月天数的权重,与考虑权重时的误差在E-2,如果不在意这点误差可以考虑使用
  1. ### gc ###
  2. gc_ds = gc.climatologies.month_to_season(ds.sel(time=slice("1940","2022")),"JJA")


  3. ### weighted ###
  4. ds_season = ds.groupby('time.season')['JJA']
  5. month_length = ds_season.time.dt.days_in_month
  6. month_length
  7. # 按“ time.season”分组来计算权重。
  8. weights = month_length.groupby('time.season')['JJA'].groupby('time.year') / month_length.groupby('time.season')['JJA'].groupby('time.year').sum()
  9. weights
  10. # 计算加权平均值
  11. ds_season_weighted = (ds_season * weights).groupby('time.year').sum(dim='time')
复制代码


密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2023-7-25 12:53:18 | 显示全部楼层
好的,谢谢楼主
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2023-7-25 23:27:13 | 显示全部楼层

密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表