pro test_cov
Compile_opt idl2
file_path="C:\Users\86178\Desktop\地物类别\新建文本文档5.txt"
openr,lun,file_path,/get_lun
temp1=intarr(6,50);根据提取的数据个数更改
for i=0,49 do begin
temp=intarr(6)
readf,lun,temp
temp1[*,i]=temp
endfor
free_lun,lun
data=transpose(temp1)
dims=size(data,/dim)
data_mean=fltarr(dims[1])
for k=0,dims[1]-1 do begin
data_mean[k]=mean(data[*,k])
endfor
img_result=fltarr(dims[1],dims[1])
for i=0,dims[1]-1 do begin
for j=0,dims[1]-1 do begin
img_result[i,j]=total((data[*,i]-data_mean[i])*(data[*,j]-data_mean[j]))/(dims[0]*dims[0])
endfor
endfor
print,'平均值'
print,data_mean
print,'协方差'
print,img_result
end