IDL求图像的协方差

程序代码：
pro test_cov
   file_dir='C:\Users\86178\Desktop\遥感原理与应用及IDL二次开发\Default\ch03_IDL图像统计程序'
   tm_files=file_search(file_dir,'jz512-*.TIF',count=filenumb)；对文件里的图像进行检索
   
      tm_b1=read_tiff(tm_files[0])
      dims=size(tm_b1,/dimensions)
      tm_data=intarr(dims[0],dims[1],filenumb)
      tm_mean=fltarr(filenumb)
      tm_mean[0]=mean(tm_b1)
     
     
   for k=0,filenumb-1 do begin
       tm_data[*,*,k]=read_tiff(tm_files[k])
       tm_mean[k]=mean(tm_data[*,*,k])
   endfor
   
   img_result=fltarr(filenumb,filenumb)
   
   for i=0,filenumb-1,1 do begin
     for j=0,filenumb-1,1 do begin
       img_result[i,j]=total((tm_data[*,*,i]-tm_mean[i])*(tm_data[*,*,j]-tm_mean[j]))/(dims[0]*dims[1])
     endfor
   endfor
   
   print,img_result
   
end