多光谱医学图像关键波段选取方法、系统与可读存储介质技术方案

本发明专利技术公开了一种多光谱医学图像关键波段选取方法、系统与可读存储介质，具体包含以下步骤：收集多光谱医学图像数据，设定参数，迭代求解得到估计的系数矩阵，采用谱聚类方法得到选取的波段集合。与已有技术方案相比，具有以下优点：利用谱聚类方法，对于高维数据计算复杂度较低；使用高斯核函数计算相似性矩阵，可得到更好的聚类效果；基于低秩假设选取波段，可实现更准确的多光谱医学图像分类效果。可实现更准确的多光谱医学图像分类效果。可实现更准确的多光谱医学图像分类效果。

【技术实现步骤摘要】
多光谱医学图像关键波段选取方法、系统与可读存储介质


[0001]本专利技术涉及多光谱医学图像处理领域，具体涉及一种多光谱医学图像关键波段选取方法、系统与可读存储介质。

技术介绍

[0002]成像光谱技术起源于遥感领域，通过在卫星上安装可探测多种波段光的探测器，可采集海洋森林矿产地理位置、环境气象、地震火山活动等方面的数据，相对于单一波段成像技术，多光谱成像技术极大提高了目标探测的准确性和探测能力。近年来，多光谱成像技术被引入到医疗诊断中，利用一台成像设备即可获得病变部位多个波段光的图像，由于不同组织在不同波段图像中显著信息不同，从而增强了组织间的对比度，可实现疾病的高效诊断。然而由于波段较多，各个波段图像间存在冗余信息，如何选取关键波段用于后续的多光谱医学图像处理是多光谱医学图像在临床诊断中应用的前提。
[0003]在多光谱医学图像处理的相关研究中，李亚研究多光谱医学图像的组织分割和特征描述方法，其中采用Sheffield指数分析不同波段图像相关性，实现关键波段选取；通过图像类间方差确定分割阈值，实现多组织分割；并提出改进局部二进制模式的图像纹理描述方法进行血管特征提取。

技术实现思路

[0004]为解决多光谱医学图像波段较多引起的存在冗余信息的技术问题，本专利技术提供一种多光谱医学图像关键波段选取方法、系统与可读存储介质，采用如下技术方案：
[0005]一种多光谱医学图像关键波段选取方法，包括以下步骤：
[0006]步骤一：收集多光谱医学图像数据
[0007]用矩阵表示一幅多光谱医学图像，其中x
i
,i＝1,2,
…
,N表示矩阵的第i列，M和N分别为矩阵X的行数和列数，M和N分别表示多光谱医学图像的像素点个数和波段数量；
[0008]步骤二：设定参数，包括所选取波段子集包含的波段数量k为正整数，误差阈值τ为正数，最大迭代次数t
max
为正整数，参数β的最大值β
max
为正数；设置带宽参数σ为正数；设置谱特征维度K为正整数，K>k；
[0009]步骤三：迭代求解得到估计的系数矩阵具体为：
[0010]3.1:初始化变量J
(0)
＝Z
(0)
＝0，E
(0)
＝X，系数系数参数β
(0)
＝10-6
；其中变量的上标(0)表示迭代开始前的初始值；其中Z表示系数矩阵，J为中间变量,E表示噪声矩阵；
[0011]3.2:令当前迭代次数t＝0；
[0012]3.3:计算J
(t+1)
，其中表示奇异值软阈值算子；
[0013]3.4:计算Z
(t+1)
，，其中E表示噪声矩阵，上标(t)表示第t次迭代的变量值；
[0014]3.5:计算E
(t+1)
，其中||
·
||
2,1
表示L21范数,||
·
||
F
表示F范数；
[0015]3.6:计算3.6:计算β
(t+1)
＝min(1.1β
(t)
,β
max
)；
[0016]3.7:如果||X-XZ
(t+1)-E
(t+1)
||
∞
>τ，且||Z
(t+1)-J
(t+1)
||
∞
>τ，或者当前迭代次数t<t
max
，则令t自增1，并转至3.3；否则转至3.8；
[0017]3.8:估计的系数矩阵即为
[0018]步骤四：采用谱聚类方法得到选取的波段集合，具体为：
[0019]4.1：根据步骤三得到的计算相似性矩阵W，其中W的第i行第j列元素计算方法为：
[0020][0021]其中i＝1,2,
…
,N，j＝1,2,
…
,N；σ为带宽参数，表示估计的系数矩阵的第i行第l列元素，l＝1,2,
…
,N；
[0022]4.2：根据相似性矩阵W，计算标准化拉普拉斯矩阵其中D为一个对角矩阵，其第i个对角元素
[0023]4.3：对L进行特征值分解，L的特征值从小到大排列，将前K个特征值所各自对应的特征向量合并写作矩阵其中u
m
表示第m个特征值对应的特征向量，m＝1,2,
…
,K，将矩阵U
K
按行标准化，得到特征矩阵F，其中
[0024]4.4：将U
K
的每个行向量作为一个样本，将N个样本聚为k个簇，每个簇的中心样本的行号对应的波段即构成所选取的波段子集。
[0025]一种多光谱医学图像关键波段选取系统，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述的多光谱医学图像关键波段选取方法。
[0026]一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述程序被处理器执行时实现上述的多光谱医学图像关键波段选取方法。
[0027]与现有技术相比，本专利技术的有益技术效果是：
[0028]利用谱聚类方法，对于高维数据计算复杂度较低；使用高斯核函数计算相似性矩阵，可得到更好的聚类效果；基于低秩假设选取波段，可实现更准确的多光谱医学图像分类效果。
附图说明
[0029]图1为本专利技术的流程图；
具体实施方式
[0030]下面结合附图对本专利技术的一种优选实施方式作详细的说明。
[0031]如图1所示的一种多光谱医学图像关键波段选取方法，包括以下步骤：
[0032]步骤一：收集多光谱医学图像数据
[0033]用矩阵表示一幅多光谱医学图像，其中x
i
,i＝1,2,
…
,N表示矩阵的第i列，M和N分别为矩阵X的行数和列数，M和N分别表示多光谱医学图像的像素点个数和波段数量；
[0034]步骤二：设定参数，包括所选取波段子集包含的波段数量k为正整数，误差阈值τ为正数，最大迭代次数t
max
为正整数，参数β的最大值β
max
为正数；设置带宽参数σ为正数；设置谱特征维度K为正整数，K>k；
[0035]步骤三：迭代求解得到估计的系数矩阵具体为：
[0036]3.1:初始化变量J
(0)
＝Z
(0)
＝0，E
(0)
＝X，系数系数参数β
(0)
＝10-6
；其中变量的上标(0)表示迭代开始前的初始值；其中Z表示系数矩阵，J为中间变量,E表示噪声矩阵；
[0037]3.2:令当前迭代次数t＝0；
[0038]3.3:计算本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种多光谱医学图像关键波段选取方法，其特征在于，包括以下步骤：步骤一：收集多光谱医学图像数据用矩阵表示一幅多光谱医学图像，其中x
i
，i＝1，2，
…
，N表示矩阵的第i列，M和N分别为矩阵X的行数和列数，M和N分别表示多光谱医学图像的像素点个数和波段数量；步骤二：设定参数，包括所选取波段子集包含的波段数量k为正整数，误差阈值τ为正数，最大迭代次数t
max
为正整数，参数β的最大值β
max
为正数；设置带宽参数σ为正数；设置谱特征维度K为正整数，K＞k；步骤三：迭代求解得到估计的系数矩阵具体为：3.1：初始化变量J
(0)
＝Z
(0)
＝0，E
(0)
＝X，系数系数参数β
(0)
＝10-6
；其中变量的上标(0)表示迭代开始前的初始值；其中Z表示系数矩阵，J为中间变量，E表示噪声矩阵；3.2：令当前迭代次数t＝0；3.3：计算J
(t+1)
，其中表示奇异值软阈值算子；3.4：计算Z
(t+1)
，，其中E表示噪声矩阵，上标(t)表示第t次迭代的变量值；3.5：计算E
(t+1)
，其中||
·
||
2，1
表示L21范数，||
·
||
F
表示F范数；3.6：计算β
(t+1)
，，β
(t+1)
＝min(1.1β
(t)
，β
max
)；3.7：如果||X-...

【专利技术属性】
技术研发人员：梅明亮，
申请(专利权)人：安徽威奥曼机器人有限公司，
类型：发明
国别省市：