[发明专利]结合空间信息的高光谱分类结果优化方法

权利要求书

1.结合空间信息的高光谱分类结果优化方法，其特征在于该方法包括以下步骤：

1）高光谱谱域分类；

首先对数据做归一化处理，并依据先验知识，在每个地物类别区域内，随机的选取一定比例的训练样本构成训练样本集；然后进行分类器的训练，采用的分类器为支持向量机；最后用训练好的分类器进行数据的测试，即能得到谱域分类的结果；

2）自适应阈值边缘提取；

自适应阈值边缘提取分为提取多波段梯度图和用一个自适应阈值二值化处理两部分；

2.1）提取多波段梯度图，具体是：

第一步：选择n个目标边缘清晰的波段组成用于目标边缘提取的数据集；

并计算各波段每个像素点在横坐标和纵坐标下的偏导数：

un=∂f(n)∂x,n=1,2,3......m---(1)]]>

vn=∂f(n)∂y,n=1,2,3....m---(2)]]>

其中f(n)表示每波段数据；

第二步：计算各波段偏导数点乘的和：

gxx=Σn=1m(un·un)=Σn=1m(unTun)=∂f(1)∂x·∂f(1)∂x+∂f(2)∂x·∂f(2)∂x·+......+∂f(n)∂x·∂f(n)∂x---(3)]]>

gyy=Σn=1m(vn·vn)=Σn=1m(vnTvn)=∂f(1)∂x·∂f(1)∂x+∂f(2)∂x·∂f(2)∂x·+......+∂f(n)∂x·∂f(n)∂x---(4)]]>

gxy=Σn=1m(un·vn)=Σn=1m(unTvn)=∂f(1)∂x·∂f(1)∂y+∂f(2)∂x·∂f(2)∂y·+......+∂f(n)∂x·∂f(n)∂y---(5)]]>

即可得到3个矩阵；

第三步：计算最大变化率方向θ和该方向上的变化率F(θ)：

θ=12arctan[2gxy(gxy-gyy)]---(6)]]>

F(θ)={12[(gxx+gyy)+(gxx-gyy)cos2θ+2gxysin2θ]}12---(7)]]>

即可得到多波段数据边缘梯度图F(θ)；

2.2）自适应阈值二值化处理，具体是：

首先，根据同一靶标的像素点应该是成区域分布这一特性，把光谱域分类结果点自动分为多块区域，并计算每一块区域的中心点；

然后，以每个中心点为圆心，画一个圆，半径根据每块区域目标点聚类的个数求得，选择圆的面积是分块区域像素点的3倍，这样就能基本覆盖每个目标，提取圆内的梯度图，从而能得到N幅目标梯度图；

同时计算每个圆内像素梯度的均值K_N,根据公式

T_N＝aK_N+b    （8）

计算得到每个圆内梯度图的阈值T_N，其中a,b是两个变量，需根据实际数据大小调节；

最后，用每个阈值T_N对每个圆内梯度图进行二值化处理和边缘细化处理，即能得到N个目标的边缘；求并集，即是所有目标的边缘；

3）内部膨胀,具体是：

第一步：提取目标现有的像素点，并判断每一个像素点（i,j）是否在每个轮廓区域内,判断的步骤为：（1）选取与该像素点同行的轮廓点坐标(X_1n,Y_1n)，n=1,2,...和与该像素点同列的轮廓点坐标(X_2m,Y_2m),m=1,2,...；（2）如果存在且则可判断该像素点在目标区域内部；

第二步：如果该像素点在目标内部，则对该点进行一次膨胀,新得到的点归为目标点；其中，膨胀操作采用的模板为：

p=010111010]]>

第三步：重复步骤第一步和第二步，直到没有点增加时结束；

第四步：对结果做并集就能得到目标和轮廓的叠加图；减掉所有轮廓点即能得到纯目标点，对图像再做一次膨胀；即可得到最终的结果图。