24小时热门版块排行榜    

查看: 654  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 材料工程,326分,求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 0856化工原理 +6 z2839474511 2026-03-10 6/300 2026-03-13 10:41 by houyaoxu
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
信息提示
请填处理意见