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

最新推荐文章于 2023-06-06 13:04:04 发布

m0_62393891

最新推荐文章于 2023-06-06 13:04:04 发布

阅读量737

点赞数 1

分类专栏： IDL 文章标签：大数据

本文链接：https://blog.csdn.net/m0_62393891/article/details/125643751

版权

本文介绍了如何在IDL环境中使用极大似然分类器进行遥感图像分类。通过两种不同的方法，计算每个像元落在不同类别上的概率，并根据最大概率确定像元的分类。最终，两种方法得到的分类结果相同。

摘要由CSDN通过智能技术生成

暂且写出了两种方法

第一种;

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