t.pro

在监测土壤湿度时

PRO T
i='0'XL
j='0'XL
;将ENVI图像的数据读入数组变量a1,a2,a3,并将数据保留小数点后两位
 envi_open_file,'I:\data\MODISdata\jiangsu\03\9-6to9-21\cut\evi', r_fid=fid
     envi_file_query,fid, ns=ns, nl=nl, nb=nb
      dims = [-1, 0, ns-1, 0, nl-1]
      pos  = lindgen(nb)
     a1 = ENVI_GET_DATA(fid=fid,dims=dims,pos=pos)
     b1=round(a1*100)
     array1=float(b1)/100
 envi_open_file,'I:\data\MODISdata\jiangsu\03\9-6to9-21\cut\day', r_fid=fid
     envi_file_query,fid, ns=ns, nl=nl, nb=nb
      dims = [-1, 0, ns-1, 0, nl-1]
      pos  = lindgen(nb)
     a2 = ENVI_GET_DATA(fid=fid,dims=dims,pos=pos)
     b2=round(a2*100)
     array2=float(b2)/100
 envi_open_file,'I:\data\MODISdata\jiangsu\03\9-6to9-21\cut\day', r_fid=fid
     envi_file_query,fid, ns=ns, nl=nl, nb=nb
      dims = [-1, 0, ns-1, 0, nl-1]
      pos  = lindgen(nb)
      a3 = ENVI_GET_DATA(fid=fid,dims=dims,pos=pos)
      b3=round(a3*100)
      array3=float(b3)/100
;将三辐影像数据读入一个3列矩阵
myarray=fltarr(3,ns*nl)
myarray[0,0:*]=array1[0:ns-1,0:nl-1]
myarray[1,0:*]=array2[0:ns-1,0:nl-1]
myarray[2,0:*]=array3[0:ns-1,0:nl-1]

;计算不为零的数据总数
x=where(array1 ne 0,n1)
y=where(array2 ne 0,n2)
z=where(array3 ne 0,n3)
n=n1>n2>n3

;创建一个新的数组，不含零值
newarray=fltarr(3,n-1)
while i lt ns*nl do begin

  if myarray[0,i] ne 0 and myarray[1,i] ne 0 and myarray[2,i] ne 0 then begin
     newarray[0,j]=myarray[0,i] & newarray[1,j]=myarray[1,i] & newarray[2,j]=myarray[2,i]
     j=j+1
  endif
 i=i+1
endwhile

;将数据读入txt文档
openw,1,'I:\data\MODISdata\jiangsu\03\9-6to9-21\statistic\9-6to9-21evi.txt'
printf,1,newarray
close,1
END