24小时热门版块排行榜    

查看: 740  |  回复: 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.55.1.O.54,科目齐全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 18:47 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:39 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 17:35 by CgyNCDVNhGVg
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 17:23 by CgyNCDVNhGVg
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-15 15:35 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-15 15:27 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:50 by VMDqgnlKimtZ
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 14:38 by VMDqgnlKimtZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:34 by VMDqgnlKimtZ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:38 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 05:22 by 5BDX0d0WFp7t
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见