基于主成分分析(PCA)的EPFs(PCA-EPFs)方法在边缘保留特征在高光谱图像分类中的应用研究（Matlab代码实现）

 💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

 ⛳️赠与读者

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、数据、文章

---

 ⛳️赠与读者

👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。当哲学课上老师问你什么是科学，什么是电的时候，不要觉得这些问题搞笑。哲学是科学之母，哲学就是追究终极问题，寻找那些不言自明只有小孩子会问的但是你却回答不出来的问题。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能让人胸中升起一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它居然给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。

     或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎

💥1 概述

摘要— 通过将边缘保持滤波器应用于高光谱图像（HSI）获得的边缘保持特征（EPFs）在表征场景中物体的重要光谱和空间结构方面非常有效。然而，直接使用EPFs可能不足以提供对空间信息的完整描述，特别是在考虑到不同尺度的物体存在于图像中时。此外，边缘保持平滑操作不可避免地会降低不同类别物体之间的光谱差异，这可能会影响后续的分类。为了解决这些问题，本文提出了一种基于主成分分析（PCA）的EPFs（PCA-EPFs）方法用于HSI分类，包括以下步骤。首先，通过将不同参数设置的边缘保持滤波器应用于考虑的图像，构建标准EPFs，并将结果EPFs堆叠在一起。接下来，使用PCA减少堆叠EPFs的光谱维度，这不仅可以在均方意义上表示EPFs，还可以突出EPFs中像素的可分离性。最后，使用支持向量机（SVM）分类器对所得到的PCA-EPFs进行分类。在几组真实高光谱数据集上进行的实验表明，所提出的PCA-EPFs方法的有效性，大大提高了SVM分类器的准确性，相对于基于标准边缘保持滤波的特征提取方法以及其他广泛使用的光谱-空间分类器。
关键词— 边缘保持滤波、高光谱图像（HSI）、图像分类、主成分分析（PCA）、支持向量机（SVM）。
 

与旨在理解场景主要内容的一般土地利用分类方法[8]–[10]不同，HSI分类的目标是为HSI中的每个像素分配一个唯一的类标签，这是一个更具挑战性的任务。为实现这一目标，传统方法，如贝叶斯估计方法[11]、支持向量机（SVM）[12]和基于稀疏表示的技术[13]，[14]已成功应用于HSI分类。然而，当只有很少标记样本存在时，由于维度灾难，大多数先前提到的分类器无法获得令人满意的分类性能。此外，相邻的无噪声高光谱波段通常高度相关，并且高光谱维度也意味着分类过程的计算负担将相应增加。

为解决这些问题，特征提取已被证明是降低数据维度的有效方法，同时保留或增加不同对象的类别可分离性[15]–[18]。一些经典的信号分析工具，如主成分分析（PCA）[15]、独立分量分析（ICA）[16]、奇异谱分析（SSA）[21]和流形学习[17]方法，已成功应用于HSI的特征提取。然而，大多数这些方法在特征提取过程中仅利用光谱信息，因此通常表现出较低的性能。

为了在高光谱图像中结合光谱和空间信息，近年来对光谱-空间特征提取方法进行了广泛研究[20]，[22]，[23]。例如，Benediktsson等人[24]引入了扩展形态剖面（EMPs）来模拟HSI中对象的几何和纹理特性，在遥感领域引起了极大关注。EMP是通过依次应用开运算和闭运算构建的，对于表示图像中结构的多尺度变异性非常有效。然而，EMP严重依赖结构元素的形状，因此不适合模拟所有几何特征。

📚2 运行结果

部分代码：

addpath ('.\functions')
addpath (genpath('.\libsvm-3.22'))
%% load original image
path='.\Datasets\';
inputs = 'IndiaP';%145*145*200/10249/16 
location = [path,inputs];
load (location);

%%% size of image 
[no_lines, no_rows, no_bands] = size(img);
GroundT=GroundT';
load (['.\training_indexes\in_1.mat'])
%% Spectral dimension Reduction
 img2=average_fusion(img,15);
for i=1:10
    indexes=XX(:,i);
%% Normalization
no_bands=size(img2,3);
fimg=reshape(img2,[no_lines*no_rows no_bands]);
[fimg] = scale_new(fimg);
fimg=reshape(fimg,[no_lines no_rows no_bands]);
%% Feature extraction
 fimg1=spatial_feature(fimg,115,0.6);

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

🌈4 Matlab代码、数据、文章