爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 20407|回复: 12

[其他] cloudsat轨道数据2B-GEOPROF怎么网格化成1°×1°的格点数据

[复制链接]

新浪微博达人勋

发表于 2020-5-28 13:55:07 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 liweier123 于 2020-5-28 13:56 编辑

我有2006-2017年cloudsat的2B-GEOPROF轨道数据,想处理成月平均的1°×1°网格点数据,请求大神帮忙





Snipaste_2020-05-28_13-55-56.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2020-6-3 11:10:15 | 显示全部楼层
请问用ncl能不能实现呢,哪位大神帮帮忙,急~~!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-6-4 20:01:58 | 显示全部楼层
我也在想这个问题
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-6-5 16:35:42 | 显示全部楼层
甜酒果 发表于 2020-6-4 20:01
我也在想这个问题

关于这个网格化可以参考这个例子https://www.ncl.ucar.edu/Applications/Scripts/binning_1.ncl
但怎么批量处理按Julian day 存储的文件毫无头绪
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-7-2 22:07:54 | 显示全部楼层
最近也在做这个,楼主实现了吗?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-7-2 23:36:06 | 显示全部楼层
轨道数据网格化可以参考这个例子 :http://www.ncl.ucar.edu/Applications/Scripts/hdf5eos_3a.ncl

;---Set up options for regridding
    Opt                 = True
    Opt@SrcGridLat      = lat               ;
  • unstructured
        Opt@SrcGridLon      = lon               ;

  •     Opt@InterpMethod    = "bilinear"

    ;---If you don't set these two, the regridding will be VERY slow
        if (DstGridType.eq."1x1") then
            Opt@DstURCorner     =  (/ 89.5d,  179.5d/)    ; for 2x2
        end if
        if (DstGridType.eq."2x2") then
            Opt@DstURCorner     =  (/ 89d,  179d/)    ; for 2x2
        end if
        if (DstGridType.eq."5x5") then
            Opt@DstURCorner     = (/ 87.5d,  177.5d/)
        end if
        Opt@DstLLCorner     =  -Opt@DstURCorner

        Opt@DstGridType     = DstGridType ; destination grid
        Opt@SrcRegional     = False
        Opt@DstRegional     = False

        Opt@RemoveSrcFile   = True        ; remove SCRIP grid destination files
        Opt@RemoveDstFile   = True        ; remove SCRIP grid destination files

        Opt@ForceOverwrite  = True
        Opt@Debug           = True

    ;---Extract desired level for regridding

        xp                  = x(:,{plvl})
        stat_xp = stat_dispersion(xp, opt )    ; ALWAYS look at satellite data
                                               ; Here, just 'xp'

        Opt@SrcGridMask     = where(ismissing(lat) .or. ismissing(xp) .or. \
                                    abs(lat).gt.90 .or. abs(lon).gt.180, 0,1)
        Opt@WgtFileName     = "MLS-Aura_L2GP-"+app+"."+Opt@DstGridType+".nc"
          
        xp_regrid           = ESMF_regrid(xp,Opt)

    把每天的轨道数据转为格点数据,再求格点数据的月平均。
  • 密码修改失败请联系微信:mofangbao
    回复 支持 反对

    使用道具 举报

    新浪微博达人勋

     楼主| 发表于 2020-7-3 21:20:23 | 显示全部楼层
    Slouvast111 发表于 2020-7-2 23:36
    轨道数据网格化可以参考这个例子 :http://www.ncl.ucar.edu/Applications/Scripts/hdf5eos_3a.ncl

    ;-- ...

    好的,谢谢!
    密码修改失败请联系微信:mofangbao
    回复 支持 反对

    使用道具 举报

    新浪微博达人勋

    发表于 2020-7-26 15:38:46 | 显示全部楼层
    请问一下cloudsat网站的数据集要如何下载?
    密码修改失败请联系微信:mofangbao
    回复 支持 反对

    使用道具 举报

    新浪微博达人勋

     楼主| 发表于 2020-7-29 12:14:10 | 显示全部楼层
    zhangjiaqing 发表于 2020-7-26 15:38
    请问一下cloudsat网站的数据集要如何下载?

    cloudDPC上注册,登录ftp://ftp.cloudsat.cira.colostate.edu下载
    密码修改失败请联系微信:mofangbao
    回复 支持 反对

    使用道具 举报

    新浪微博达人勋

    发表于 2020-7-30 15:39:38 | 显示全部楼层
    liweier123 发表于 2020-7-29 12:14
    cloudDPC上注册,登录ftp://ftp.cloudsat.cira.colostate.edu下载

    多谢,已经注册下载
    密码修改失败请联系微信:mofangbao
    回复 支持 反对

    使用道具 举报

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

    本版积分规则

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

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

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