24小时热门版块排行榜    

查看: 709  |  回复: 1

cczzw

捐助贵宾 (初入文坛)

[求助] 求教一个IDL程序的问题

FUNCTION read_modis_cloud_mask, filename


    sd_names = 'Cloud_Mask'
    sd_names_lat = 'Latitude'
    sd_names_lon = 'Longitude'
   
    sd_id = HDF_SD_START(filename, /READ)
   
    sds_index_cm = HDF_SD_NAMETOINDEX(sd_id, sd_names)
    sds_index_lat = HDF_SD_NAMETOINDEX(sd_id, sd_names_lat)
    sds_index_lon = HDF_SD_NAMETOINDEX(sd_id, sd_names_lon)
   
    sds_id_cm = HDF_SD_SELECT(sd_id, sds_index_cm)
    sds_id_lat = HDF_SD_SELECT(sd_id, sds_index_lat)
    sds_id_lon = HDF_SD_SELECT(sd_id, sds_index_lon)
   
    ;Get information of hdffile.
    HDF_SD_GETINFO, sds_id_cm, NAME = name, DIMS = dims, TYPE = data_type
   
    ;Read byte 1.
    count = [dims[0], dims[1], 1]
    HDF_SD_GETDATA, sds_id_cm, data_cm, COUNT = count
   
    data_cm_size = SIZE(data_cm, /DIMENSIONS)
   
    ;Read latitude
    HDF_SD_GETDATA, sds_id_lat, data_lat
    data_lat = CONGRID(data_lat, data_cm_size[0], data_cm_size[1], /INTERP)
   
    ;Read longitude
    HDF_SD_GETDATA, sds_id_lon, data_lon
    data_lon = CONGRID(data_lon, data_cm_size[0], data_cm_size[1], /INTERP)
   
    HDF_SD_ENDACCESS, sds_id_cm
    HDF_SD_ENDACCESS, sds_id_lat
    HDF_SD_ENDACCESS, sds_id_lon
   
    HDF_SD_END, sd_id
   
    ;BIT 1,2 - unobstructed FOV quality flag
    ; 0 - Cloud
    ; 1 - 66% probality of clear
    ; 2 - 95% probality of clear
    ; 3 - 99% probality of clear
     
    bit1 = bit_get(data_cm, 1)
    bit2 = bit_get(data_cm, 2)
   
    clear99_prob_idx = WHERE(bit2 AND bit1)
    clear95_prob_idx = WHERE(bit2 AND (~ bit1))
    clear66_prob_idx = WHERE((~ bit2) AND bit1)
    cloud_idx = WHERE((~ bit2) AND (~ bit1))
   
    mask = BYTARR(dims[0], dims[1])
    mask[clear99_prob_idx] = 3
    mask[clear95_prob_idx] = 2
    mask[clear66_prob_idx] = 1
    mask[cloud_idx] = 0
   
    data = {mask : mask, lat : data_lat, lon : data_lon}
   
    RETURN, data
END
程序运行不了,不懂为什么,有没有懂IDL编程的朋友帮忙解答,万分感谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

金金金音欠

铁虫 (初入文坛)

楼主您好,请问您的问题解决了吗?想向您请教一下MOD35云掩膜提取的方法
这个人很懒,没有个性签名
2楼2020-09-02 12:38:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 学员wurF9B 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 8/400 2026-08-02 15:01 by jHcrHMhbKCFV
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 14/700 2026-08-02 14:58 by jHcrHMhbKCFV
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 18/900 2026-08-02 14:52 by jHcrHMhbKCFV
[基金申请] 面上提前没消息,有中的吗 +7 archvillain 2026-08-02 9/450 2026-08-02 14:36 by biomichael
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 12/600 2026-08-02 13:43 by jHcrHMhbKCFV
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +17 低垂的野花 2026-07-31 23/1150 2026-08-02 13:41 by RRRKKK
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cBGhfJFuNMca 2026-08-02 16/800 2026-08-02 13:34 by jHcrHMhbKCFV
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cBGhfJFuNMca 2026-08-02 14/700 2026-08-02 13:34 by jHcrHMhbKCFV
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 21/1050 2026-08-02 10:56 by iaeyuan
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 10/500 2026-08-02 10:00 by 锦衣卫寒战
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 7/350 2026-08-02 07:59 by 88817753
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
信息提示
请填处理意见