遥感图像除云程序——转载自科学网

hillside

以下源码转载自科学网，本人对于内容并不熟悉，仅供有意者参考
http://blog.sciencenet.cn/blog-507873-619280.html
[转载] 除云程序源码

已有 346 次阅读 2012-10-4 23:16 |系统分类:[color=#999999 !important]科普集锦|关键词:程序

；===============================================
；this programe write by gxf1221@126.com，
；=============================================
  FUNCTION  cloud_only,b1,b2,b3
  temp1=(float(b1)-b2)/b1
  temp2=(float(b3)-b2)/b3
        result = (temp1 gt 0.2 and temp2 gt 0.2)*b2 $
        + (temp1 le 0.2 or temp2 le 0.2)*0
   return ,result
END
；===============================================

FUNCTION  cloud_done,b1,b2,b3
  temp1=(float(b1)-b2)/b1
  temp2=(float(b3)-b2)/b3
        result = (temp1 gt 0.2 and temp2 gt 0.2)*(b1+b3)/2 +(temp1le0.2ortemp2le0.2)∗b2  return,resultENDprocloudeffect1out,event    envi,/restorebasesavefiles    ENVISELECT,fid=fid,pos=pos,
  /no_dims, title='Please Input File'
  if (fid[0] eq -1) then return
   ENVI_FILE_QUERY, fid, ns=ns, nl=nl, nb=nb, dims=dims,bnames=bnames,datatype=datatype,
  interleave=interleave

  FOR i=0,nb DO BEGIN
      t_fid=[fid,fid,fid]
      j=i+1
      k=j+1

      tpos  = [i,j,k]
       snames=bnames[j]
      tnames=bnames+' '+bnames[j]+' '+bnames[k]

    exp='cloud_done(b1,b2,b3)'

      base = widget_auto_base(title='File Selection test')
       ws = WIDGET_SLABEL(base, /FRAME, PROMPT=tnames, XSIZE=40 , YSIZE=3 )
   wo = widget_outf(base, uvalue='outf', /auto)
   result = auto_wid_mng(base)
      if (result.accept eq 0) then return
     out_name =result.outf
      out_bname=out_name

    envi_doit, 'math_doit', fid=tfid,pos=tpos,dims=dims,
    exp=exp, out_name=out_name, out_bname=out_bname, $
    r_fid=r_fid
   
  ENDFOR

  envi_batch_exit
end

kongfeng0824

高端大气上档次啊

angelinacy

这个很实用啊