24小时热门版块排行榜    

查看: 715  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天务委会开完了,明天出结果吗 +7 angus9576 2026-08-25 7/350 2026-08-25 15:50 by Weed114114
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 18/900 2026-08-25 14:50 by Weed114114
[基金申请] 什么时候开奖? +12 CrisMessi 2026-08-18 13/650 2026-08-25 14:40 by ZJTJZ
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 2026国自然函评费到账 +19 羊腰板 2026-08-21 22/1100 2026-08-25 10:38 by popular289
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
信息提示
请填处理意见