envi提取纹理信息并进行地物可分离性指数分析,选取最合适窗口的纹理特征
2021/12/26 23:11:35
本文主要是介绍envi提取纹理信息并进行地物可分离性指数分析,选取最合适窗口的纹理特征,对大家解决编程问题具有一定的参考价值,需要的程序猿们随着小编来一起学习吧!
分别表示提取纹理的窗口大小和提取纹理的参数距离和角度。
glcm0=double(imread('GLCM\GLCM3_0_TF.tif'));%读取3*3窗口的四个角度的纹理数据
glcm1=double(imread('\GLCM\GLCM3_1_TF.tif'));
glcm2=double(imread('\GLCM\GLCM3_2_TF.tif'));
glcm3=double(imread('\GLCM\GLCM3_3_TF.tif'));
glcm00=reshape(glcm0,[],size(glcm0,3));
glcm11=reshape(glcm0,[],size(glcm1,3));
glcm22=reshape(glcm0,[],size(glcm2,3));
glcm33=reshape(glcm0,[],size(glcm3,3));
glcm_zh=cat(3,glcm00,glcm11,glcm22,glcm33);
glcm=mean(glcm_zh,3);%四个角度的纹理求平均
roi=imread('samples3.tif');%读取样本数据
SIab=zeros(8,8,48);
mea_si=zeros(48,1);
for i=1:48
for j=1:7
ind1=find(roi==j);
ind2=find(roi==j+1);
a=glcm(ind1,i);
b=glcm(ind2,i);
SIab(j,j+1,i)=abs(mean(a)-mean(b))/(std(a)+std(b));%计算SI,可分离性指数
end
ind3=find(SIab(:,:,i)~=0);
siab=SIab(ind3,i);
mea_si(i,1)=mean(siab);%计算可分离性指数的平均值
disp(i);
disp(mea_si(i,1));
end
这篇关于envi提取纹理信息并进行地物可分离性指数分析,选取最合适窗口的纹理特征的文章就介绍到这儿,希望我们推荐的文章对大家有所帮助,也希望大家多多支持为之网!
- 2024-05-15鸿蒙生态设备数量超8亿台
- 2024-05-13TiDB + ES:转转业财系统亿级数据存储优化实践
- 2024-05-09“2024鸿蒙零基础快速实战-仿抖音App开发(ArkTS版)”实战课程已上线
- 2024-05-09聊聊如何通过arthas-tunnel-server来远程管理所有需要arthas监控的应用
- 2024-05-09log4j2这么配就对了
- 2024-05-09nginx修改Content-Type
- 2024-05-09Redis多数据源,看这篇就够了
- 2024-05-09Google Chrome驱动程序 124.0.6367.62(正式版本)去哪下载?
- 2024-05-09有没有大佬知道这种数据应该怎么抓取呀?
- 2024-05-09这种运行结果里的10.100000001,怎么能最快改成10.1?