| 查看: 951 | 回复: 2 | ||
莲子清如水8木虫 (职业作家)
|
[求助]
哪位大神能详细解释一下这段idl程序，根据QA进行modis云检测及云处理
|
|
FUNCTION cloudmask, b1 R0 = float(b1) mod 2 R1 = float(b1/2) mod 2 R2 = float(b1/2^2) mod 2 R3 = float(b1/2^3) mod 2 R4 = float(b1/2^4) mod 2 R5 = float(b1/2^5) mod 2 R6 = float(b1/2^6) mod 2 R7 = float(b1/2^7) mod 2 R8 = float(b1/2^8) mod 2 R9 = float(b1/2^9) mod 2 R10 = float(b1/2^10) mod 2 R11 = float(b1/2^11) mod 2 R12 = float(b1/2^12) mod 2 R13 = float(b1/2^13) mod 2 R14 = float(b1/2^14) mod 2 R15 = float(b1/2^15) mod 2 ; land = (R3 eq 0) and (R4 eq 0) and (R5 eq 1) cloud = ((R0 eq 0) and (R1 eq 1)) or ((R0 eq 1) and (R1 eq 0)) shadow = (R2 eq 1) snow = (R12 eq 1) cloudmask = (((cloud eq 1) or (snow eq 1)))*1+ ((cloud ne 1) and (snow ne 1))*0 RETURN, cloudmask END pro qacloud For j=2011,2011,1 do begin time=strtrim(String(j),2) ; year=strtrim(string(j),2)+'\' For k=169,169,1 do begin day=strtrim(String(k, FORMAT='(I03)'),2) ; path='I:\'+'MOD09GA_'+time+'\' auxiliary='I:\'+'auxiliary data'+'\' prefix='clip'+'\' suffix1='NE_qacloud' suffix1$='cloud'+'\' result$='clip'+'\' file=path+result$+'NE'+'_GAstate'+time+day resultfile=path+suffix1$+suffix1+time+day shape=auxiliary+'NErang_Project.shp' envi_open_File,file,r_fid=fid envi_file_query,fid,dims=dims,nb=nb,nl=nl,ns=ns map_info = envi_get_map_info(fid=fid) b1=uindgen(nl,ns) b1=envi_get_data(fid=fid,dims=dims,pos=0)*0 result=cloudmask(b1) openw,lun,resultfile,/get_lun writeu,lun, result ENVI_SETUP_HEAD, fname=resultfile,$ ns=ns, nl=nl, nb=1, $ interleave=0, data_type=2, $ MAP_INFO = map_info , $ offset=0, /write free_lun,lun ; endfor endfor end 哪位大神能详细解释一下这段代码，看不懂；我要用mod09A1的话，这段程序该如何修改？ |
回复此楼
windwqf
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1053.3
- 散金: 100
- 红花: 2
- 帖子: 262
- 在线: 165.8小时
- 虫号: 1336685
- 注册: 2011-07-03
- 性别: GG
- 专业: 遥感机理与方法
2楼2013-10-15 11:41:04
绿叶丹
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1020.7
- 散金: 10
- 红花: 1
- 帖子: 191
- 在线: 38.9小时
- 虫号: 3622056
- 注册: 2014-12-31
- 专业: 全球变化生态学
3楼2015-03-26 10:29:46
