统计质量改善后的遥感影像的结果——第一部分，云层检测二值化

李百柏

已于 2024-06-09 21:15:06 修改

阅读量376

点赞数 4

分类专栏：基于IDL与ENVI API的Landsat数据简单处理文章标签：数据库

于 2024-06-06 20:33:42 首次发布

本文链接：https://blog.csdn.net/aoaoooo/article/details/139509487

版权

基于IDL与ENVI API的Landsat数据简单处理专栏收录该内容

11 篇文章 0 订阅

订阅专栏

pro Part_6
compile_opt idl2
e=envi(/headless)
for k=1, 69, 1 DO BEGIN
kk=strtrim(string(k),1)
file_mkdir,"D:\Lansat_Data_Cloud_Final\"+kk+"\"
endfor
for i=1,69,1 do begin
ii=strtrim(string(i),1)
folder = "D:\Lansat_Data_Composition\"
input = folder + ii + ".dat"
out_dir=strmid("D:\Lansat_Data_Cloud_Final\",0)+ii+"\"
out_name=out_dir+ii+".dat"
ENVI_OPEN_FILE,input,r_fid=fid
IF (fid[0] EQ -1) THEN CONTINUE
ENVI_FILE_QUERY, fid,dims=dims
t_fid = [fid,fid,fid]
pos = [0,1,2]
exp = "(((long(b1)+b2+b3) gt 0)*((long(b1)+b2+b3) lt 50000)*1 + ((long(b1)+b2+b3) ge 50000)*((long(b1)+b2+b3) lt 60000)*((float(b1)/b3) ge 0.9)*((float(b1)/b3) lt 1.1)*5 + ((long(b1)+b2+b3) ge 50000)*((long(b1)+b2+b3) lt 60000)*((float(b1)/b3) gt 0)*((float(b1)/b3) lt 0.9)*1 + ((long(b1)+b2+b3) ge 50000)*((long(b1)+b2+b3) lt 60000)*((float(b1)/b3) gt 1.1)*1 + ((long(b1)+b2+b3) ge 60000)*((long(b1)+b2+b3) lt 70000)*((float(b1)/b3) ge 0.82)*((float(b1)/b3) lt 1.75)*5 + ((long(b1)+b2+b3) ge 60000)*((long(b1)+b2+b3) lt 70000)*((float(b1)/b3) gt 0)*((float(b1)/b3) lt 0.82)*1 + ((long(b1)+b2+b3) ge 60000)*((long(b1)+b2+b3) lt 70000)*((float(b1)/b3) gt 1.75)*1 + ((long(b1)+b2+b3) ge 70000)*((long(b1)+b2+b3) lt 80000)*((float(b1)/b3) ge 0.9)*((float(b1)/b3) lt 1.9)*5 + ((long(b1)+b2+b3) ge 70000)*((long(b1)+b2+b3) lt 80000)*((float(b1)/b3) gt 0)*((float(b1)/b3) lt 0.9)*1 + ((long(b1)+b2+b3) ge 70000)*((long(b1)+b2+b3) lt 80000)*((float(b1)/b3) gt 1.9)*1 + ((long(b1)+b2+b3) ge 80000)*((long(b1)+b2+b3) lt 90000)*((float(b1)/b3) ge 0.9)*((float(b1)/b3) lt 2.5)*5 + ((long(b1)+b2+b3) ge 80000)*((long(b1)+b2+b3) lt 90000)*((float(b1)/b3) gt 0)*((float(b1)/b3) lt 0.9)*1 + ((long(b1)+b2+b3) ge 80000)*((long(b1)+b2+b3) lt 90000)*((float(b1)/b3) gt 2.5)*1 + ((long(b1)+b2+b3) ge 90000)*((long(b1)+b2+b3) lt 100000)*((float(b1)/b3) ge 0.95)*((float(b1)/b3) lt 3.5)*5 + ((long(b1)+b2+b3) ge 90000)*((long(b1)+b2+b3) lt 100000)*((float(b1)/b3) gt 0)*((float(b1)/b3) lt 0.95)*1 + ((long(b1)+b2+b3) ge 90000)*((long(b1)+b2+b3) lt 100000)*((float(b1)/b3) gt 3.5)*1 + ((long(b1)+b2+b3) ge 100000)*((long(b1)+b2+b3) lt 150000)*((float(b1)/b3) ge 0.95)*((float(b1)/b3) lt 4.0)*5 + ((long(b1)+b2+b3) ge 100000)*((long(b1)+b2+b3) lt 150000)*((float(b1)/b3) gt 0)*((float(b1)/b3) lt 0.95)*1 + ((long(b1)+b2+b3) ge 100000)*((long(b1)+b2+b3) lt 150000)*((float(b1)/b3) gt 4.0)*1)"
ENVI_DOIT,'math_doit',fid=t_fid,pos=pos,dims=dims,exp=exp,out_name=out_name,r_fid=r_fid
print,input,out_name
print,"-------------------------------------",i
endfor
end