24小时热门版块排行榜    

查看: 773  |  回复: 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一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-14 13:02 by KXypaGgBBlw7
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 FChKX8meu1ds 2026-11-09 10/500 2026-11-13 11:08 by Equinoxhua
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 C6CRn66RYYXS 2026-11-10 6/300 2026-11-13 11:04 by Equinoxhua
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 cdqGD6ftHmbl 2026-11-11 5/250 2026-11-12 12:28 by ECLTXnJnRSXf
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 cdqGD6ftHmbl 2026-11-11 4/200 2026-11-12 12:26 by ECLTXnJnRSXf
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:33 by HajghKOMynmz
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 6/300 2026-11-12 11:22 by HajghKOMynmz
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 5/250 2026-11-12 11:14 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:30 by HajghKOMynmz
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 8/400 2026-11-12 10:27 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 11/550 2026-11-12 09:57 by HajghKOMynmz
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-12 09:21 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 YHiw1gou1Eda 2026-11-10 7/350 2026-11-12 05:41 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 05:30 by 9fOTqj074tIg
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 05:30 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 10/500 2026-11-12 05:19 by 9fOTqj074tIg
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 cqv1zKnzd7hr 2026-11-09 8/400 2026-11-12 05:08 by 9fOTqj074tIg
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 02:23 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-12 00:44 by 9fOTqj074tIg
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
信息提示
请填处理意见