在IDL中实现极大似然分类器

暂且写出了两种方法

第一种;

For i=0,ns-1 do begin
      For j=0,nl-1 do begin

        tm_ij = Transpose(tm_data[i,j,*])   
    ;根据极大似然分类器求个像元落在各类的“概率”
    ;alog 自然对数
    ;determ 行列数求值
    ;transpose 矩阵转置
    ;invert 求逆矩阵
    ;A#B  A的列乘以B的行 相当于B乘A
    P_urban=-(alog(determ(cov_urban))+transpose(tm_ij-mean_urban)#invert(cov_urban)#(tm_ij-mean_urban))/2
    P_forest=-(alog(determ(cov_forest))+transpose(tm_ij-mean_forest)#invert(cov_forest)#(tm_ij-mean_forest))/2
    P_water=-(alog(determ(cov_water))+transpose(tm_ij-mean_water)#invert(cov_water)#(tm_ij-mean_water))/2
    P_crop=-(alog(determ(cov_crop))+transpose(tm_ij-mean_crop)#invert(cov_crop)#(tm_ij-mean_crop))/2
    P_barren=-(alog(determ(cov_barren))+transpose(tm_ij-mean_barren)#invert(cov_barren)#(tm_ij-mean_barren))/2
    pro_mtx=[P_urban,P_forest,P_water,P_crop,P_barren]
    max_p