专利名称：基于量子遗传优化的流形降维医学图像检索方法
技术领域：
本发明涉及的是一种图像的检索方法，特别是一种属于医学图像的检索方 法。具体涉及一种基于量子遗传优化的流形降维医学图像检索方法。
(二)
背景技术：
随着医学影像设备的不断发展，临床上产生的多种模态的医学影像数据越来 越多，对影像内容的描述也变得越来越丰富，产生的特征数据集的维数也随之不 断增高。在基于内容的医学图像检索(Content-based Medical Image Retrieval, CBMIR)中，对医学影像内容的描述通常用特征向量的形式来表示，为了尽可 能准确、全面地描述图像信息，这些特征向量往往是从图像中抽取出来的几十维 甚至上百维的高维特征向量,这对数据库的组织、索引、基于距离函数的相似性 度量以及检索都是很大的挑战，即所谓的维数灾难。事实上，对于某一特定的高 维特征向量集合，以及它们所张成的高维空间，这些高维向量的重要程度是不同 的，其间仍然存在着很大的冗余。因此，通过降维技术去除次要信息的特征向量 以及相关性较强的特征向量，并保持其内在结构是CBMIR研究的关键问题之一。
流形学习是一种非线性的维数縮减方法，它的本质是将高维空间中的数据映 射到一个低维的非线性子空间(即低维流形)，从而实现降维。局部线性嵌入 (locally linear embedding, LLE)算法是S.T.Roweis和L.K.Saul.在2000年提出的 一种非线性流形降维方法，其用于医学图像特征降维的基本思想是利用局部线性 重建反映高维数据空间中的非线性结果，使降维后的数据保持原始空间的拓扑结 构。传统局部线性嵌入-LLE算法中最核心的问题是求取使重建误差最小的近邻 局部重建权值矩阵。然而，该算法是针对局部数据进行操作，而大多数的研究者 都采用与欧氏距离有关的变量来定义该权值矩阵，默认距离近的点相互影响大而 距离远的点间影响小，这使得该算法对样本中的噪音很敏感，此外该算法收敛速 度不够快。
量子遗传算法(Quantum Genetic Algorithm, QGA)是量子计算(Quantum Computing)与遗传算法(Genetic Algorithm, GA)相结合的产物。以量子计算的概念和理论为基础，用量子位编码表征染色体，通过量子门旋转更新来完成进 化搜索，能够表示出解的线性叠加态，获得更好的种群多样性、更快的收敛速度 和全局寻优的能力。但到目前为止，还没有人将该算法应用于局部线性嵌入降维 的优化方面。
发明内容
本发明的目的在于提供一种可以降低噪声点对样本数据的影响，提高运行速 度的基于量子遗传优化的流形降维医学图像检索方法。 本发明的目的是这样实现的
1、 提出査询请求首先由用户通过査询接口向医学图像数据库提出査询目 标图像的查询请求，所述的医学图像数据库是图像入库时，将各序列的关键图像 统一转换为512X512大小的JPG格式图像构成医学图像数据库；
2、 提取高维医学图像特征查询处理模块根据请求提取目标图像和医学图 像数据库中每幅图像的72维灰度特征和48维纹理特征，将这些高维特征以向量 的形式存入高维特征向量库；
3、 基于量子遗传优化的降维对描述图像内容的高维特征进行基于量子遗 传优化的降维，得到32维低维特征子集；
4、 相似性度量通过相似性度量将检索结果返回给用户。 本发明还可以包括
1、所述的基于量子遗传优化的降维的方法为
(1) 医学图像特征向量高维空间/^中的W个Z)维随机向量",X2,…;cJ， x^iT，通过降维把它们映射到低维空间^中，得到降维后的输出向量乂 ，
e [1, W] ， d《Z);
计算每一个样本点即特征向量jc,的;t个最近邻点；对高维空间中的样本点;c,, 计算该点与其他W-1个样本点之间的距离，将距离排序，选择前A:个与x,最近 的点作为其邻近点；
(2) 由每个样本点的近邻点计算出该样本点的局部重建权值矩阵『；用每 个特征向量的近邻点对该特征向量进行重建，求取使重建误差最小的近邻局部重 建权值矩阵^。
G)保持上步求取的权值『不变，求Xi在低维嵌入空间的输出特征向量y,;为了使输出数据在低维空间保持高维空间原有的拓扑结构需构造一个代价函数， 映射过程中必须使代价函数值达到最小，该函数形式为
'=1 7 = 1
其中，y,是x,的输出向量，乂, 01，2，…,"是》，,的A个近邻点； 为使低维重构误差最小化同样也要满足两个约束条件 1) Y,是一个平移不变量，即^.=0;
2 )为了避免产生退化解需使= / ，其中N是数据点的个数，/是"xd 的单位矩阵；
对e(r)进行求解，求取的最优解^可以使得代价函数sO0达到最小值；
构造一个WxW的稀疏矩阵『，W,(:1,2，…,7V)可以存储在其中，当、是;c,的
近邻点时，R;=W .，否则，^,,=0。 此时误差代价函数式可改写为
其中，M—M,.,)是一个WxW的对称矩阵，其表达式为
似=(/一『/(/一『)
利用Lagrange乘子法转化为在丄^X =/条件下的特征值和特征向量的问
题(M-八)77'=0，艮卩Ma"";
要最小化代价函数，则取M的最小d个特征值对应的向量为列向量组成矩
阵y，则r的列向量即为"维空间的降维向量表示。
2、所述的求取使重建误差最小的近邻局部重建权值矩阵的方法为
1) 初始化种群初始化医学图像特征高维向量中包含7V个个体的种群
其中《C/、l,2,…A0为种群中第f代的第J个个体，种群中
全部染色体的所有基因a、 e都被初始化为1/W，这意味着一个染色体所表达 的是其全部可能状态的等概率叠加；
2) 对种群中的个体进行测量对初始种群中的个体进行一次测量以获得一组确定的解尸。={",其中《(_/ = l,...，AO是第f代种群中的第j个
个体的测量值，其表现形式为量子为数目即染色体长度为附的二进制字符串； 其中每-"位为0或1是根据量子比特的概率| ,(|2或|/ /|2 (/ = 1,...，附)测量得到的；
具体实现方法为随机产生[o，i]上的一个数e，若^>|";|2，则"中相应的位取
值为1，否则取值为0;
3) 评价二进制染色体种群p(r)并保留最优解用适应度评价函数评价种群
尸(r)中每个染色体，并保留此代中的最优解；如算法获得满意解，算法停止； 否则，转入步骤4)继续进行；
4) 量子遗传操作使用量子旋转门f/②更新P(r);
5) 改变进化代数进化代数加1，如仍未满足最大进化代数rm^，算法转
至歩骤2)继续进行。
3、 所述的相似性度量的方法为直方图的交，
假设/和2是两个含有W个biri的颜色直方图，则它们之间的相交距离用下 式表示
|>11(/" 2》
产l
直方图的相交是指两个直方图在每个bin中共有的像素数量；有时，该值还 可以通过除以其中一个直方图中所有的像素数量来实现标准化，从而使其值处于的值域范围，如下式所示
s(/,e)《min(/^)/i:;x
4、 所述的相似性度量的方法为余弦距离法-
s(/，e)=/7'*2/(imi*iieii)
其中，I和Q分别表示查询图像和数据库中图像的特征向量，ih哮示向量范数，
计算得到的相似性度量值在[o， l]之间，该值越大，表示图像越相似。
5、 所述的相似性度量的方法，对于组合特征，相似性度量定义为各个特征 相似性度量的加权和，其公式为<formula>formula see original document page 10</formula>表示由m个特征组合而成，其中^表示第7'个 特征的权重系数，它表示第y个特征的重要性， 一般取各,相等，又S,(/，0表
示第/个特征的相似性度量函数值。
本发明将量子遗传算法应用于局部线性嵌入降维的优化。 局部线性嵌入(LLE)方法
LLE算法有一个前提假设采样数据所在的低维流形在局部是线性的，即每 个采样点可以用它的近邻点线性表示。
该算法的学习目标是在低维空间中保持每个邻域中的权值不变，即假设嵌 入映射在局部是线性的条件下，最小化重构误差。
设有高维空间i D中的TV个D维随机向量",X2，…xJ ，通过降维把它们映射 到低维空间W中，得到降维后的输出向量乂， ， = 1,2,…，W， "<<Z)。通常希望能 够用显示的映射F来表示，但当U》是非线性结构时，求出显示映射P相当困难。 基于流形学习的原理，可以将非线性结构的局部用线性结构来近似，即可以将非 线性结构的降维分解为局部的线性降维。这就是LLE算法的基本原理，该算法 的主要步骤如下。
(1) 计算每一个样本点(特征向量)x,的k个最近邻点。对高维空间中的样本 点x,,计算该点与其他TV- 1个样本点之间的距离，将距离排序，选择前A:个与x,
最近的点作为其邻近点。
(2) 由每个样本点的近邻点计算出该样本点的局部重建权值矩阵『。用每个特 征向量的近邻点对该特征向量进行重建，求取使重建误差最小的近邻局部重建权 值矩阵。
定义重建拟合误差函数如下
<formula>formula see original document page 10</formula>
其中，s(y",2,…,的为;c,的第7个近邻点，、是x,与x(,之间的权值在满足以下
两个约束条件时，通过最小化误差函数得到局部重建权值矩阵，即由样本点的近 邻点，构造出最优『矩阵使误差函数值达到最小。
a)每一个数据点《都只能由它的邻近点来表示，若x,不是近邻点，则w7 = 0 ;b)权值矩阵的每一行的和为1，即满足归一化约束1^=1。
片1
G)保持权值『不变，求;c,在低维嵌入空间的输出特征向量y,。为了使输出数
据在低维空间保持高维空间原有的拓扑结构需构造一个代价函数，映射过程中必 须使代价函数值达到最小。该函数形式为-
i>,& ii2
'=1 y=l
其中，y,是 的输出向量，A(X，2，…力是y,的A个近邻点。 为使低维重构误差最小化同样也要满足两个约束条件
a) y,是一个平移不变量，即！>,.=0;
b) 为了避免产生退化解需使^f凡少/ = / ，其中W是数据点的个数，/ 是dxd的单位矩阵。
此时，对s(J0进行求解，求取的最优解y,可以使得代价函数sOO达到最小值。
构造一个WxW的稀疏矩阵ff ， ；(, = 1,2，一^)可以存储在其中，当x,是;c,的近邻 点时，R, =w ，否则，K,,=0。 此时误差代价函数式可改写为
/=1 乂=1
其中，似=(似,:,)是一个7^^的对称矩阵，其表达式为
Af = (/ —『/■(/ —『)
利用Lagmnge乘子法可转化为在丄fx =/条件下的特征值和特征向量的
问题(M-A)yr=0，艮卩Mor"or。
要最小化代价函数，则取M的最小c/个特征值对应的向量为列向量组成矩 阵y，则F的列向量即为J维空间的降维向量表示。
量子遗传智能优化算法
传统局部线性嵌入-LLE算法中最核心的问题是求取使重建误差最小的近邻 局部重建权值矩阵，然而，该算法是针对局部进行操作，大多数的研究者都采用与欧氏距离有关的变量来定义该权值矩阵，默认距离近的点相互影响大而距离远 的点间影响小，这使得该算法对样本中的噪音很敏感，此外该算法收敛速度不够 快。量子遗传智能优化算法建立在量子态矢量表述基础上，将量子比特的几率幅 表示应用于染色体的编码，使得一条染色体可以表达多个态的叠加，并利用量子 旋转门实现染色体的更新操作，从而实现了目标的优化求解。所以，本发明利用 量子遗传算法从数据间的本征特性出发，寻找使医学图像特征向量降维重建误差 最小的权重值，从而达到符合使误差函数最小的目的。
本发明将LLE算法中样本点与其临近点的重建权值集合构造为量子遗传算 法的量子染色体，然后量子遗传中的个体在寻优的过程中找到最优位置，最终得 到近邻局部重建权值矩阵。 (1)量子比特编码
在量子遗传算法中，染色体用量子比特来表示，或者说是用随机概率方式表 示。量子比特是一个双态量子系统，与经典位的不同之处在于， 一个量子比特除 了处于量子位的两个基本态0态或1态外，还可以同时处在O, l两个本征态的迭加 状态， 一个量子比特状态可以表示为-
其中，"和"表示相应态出现的概率幅。对量子比特测量时得到O的概率为
1 |2，得到l的概率为l-l2。 二者满足归一化条件1"12+1/3|2=1。所以，在QGA中，/ 个量子比特的概率幅可表示为
《,=
卜"2、
〃2…-,
则w个/位基因构成的一个染色体可以表示为
仏'=
"12a21"22a2, m2、
/ 21A2A,
其中，z表示此染色体的编号，z表示此染色体当前进化的代数，w表示染色体 基因的个数，/表示基因里量子比特的个数。因此，如果存在一个具有w个量子比
特位的系统，则该系统能同时表达2"个状态。 (2)量子旋转门在量子遗传算法中，遗传操作主要是将构造的量子门c/("作用于量子叠加态
或纠缠的基态，使其相互干涉，相位发生改变，从而改变各基态的概率幅。在本 发明中，量子门主要采用量子旋转门，通过量子旋转门可以实现染色体的调整， 其调整操作如下
在
cos《 —sin《 sin《 cos《
式中(",，a)为第/个量子比特，a为旋转角，其幅度影响收敛速度，但是如果其
幅度太大，会导致早熟，因此，《的大小和方向本文根据通用的、与问题无关 的调整策略确定。
(3)适应度评价函数 适应度评价函数是衡量个体优劣的标志，其作用类似于度量自然界中生物适 应环境的能力。本发明根据医学图像特征向量作为个体模型的特殊性，所采用的 适应度评价函数为-
乂=1
式中，M^.表示种群中第i代的第y个个体，iV表示种群中个体的数量。
本发明对医学图像检索中高维特征向量的降维情况进行了研究，采用基于量 子遗传优化的流形降维方法实现了医学图像的有效检索。
量子遗传算法比传统遗传算法拥有更好的种群多样性和更少的种群规模，并 且，由于其采用量子比特编码也可以获得较好的收敛性。因而，本发明利用量子
遗传算法从数据间的本征特性出发，寻找在医学图像特征向量降维中使重建误差 最小的权重值，从而使得医学图像检索在保证精度的条件下达到较快的检索速
为了更好的发现医学图像检索中基于流形的非线性降维方法的本征维数，并 寻找使误差函数最小的权值矩阵，本发明提出了基于量子遗传优化的流形降维医 学图像检索方法。在本发明中，对LLE非线性降维方法进行了改进，引入量子 遗传优化局部重建权值矩阵，降低了噪声点对样本数据的影响，提高了算法收敛 速度，提高了医学图像检索速度和精度。(四)
图1.量子遗传优化算法流程图； 图2.医学图像检索系统原理图。
(五)
具体实施例方式
下面结合附图举例对本发明做更详细地描述
本发明针对不同解剖部位及各类不同疾病的多模态医学影像进行检索研究， 选取2种模态(CT与MR)， 3种解剖部位(脑、胸、腹)，6种疾病(血肿、垂 体瘤、乳腺癌、肺癌、胰腺癌、膀胱癌)共约800个病例的影像资信与正常或其 他非正常资料共约IOOO例来构成研究数据库。图像入库时，各序列的关键图像 统一转换为512X512大小的JPG格式图像。首先由用户通过査询接口提出査询 请求(目标图像)，查询处理模块根据请求提取目标图像和图像数据库中每幅图像 的72维灰度特征和48维纹理特征，这些高维特征以向量的形式存入高维特征向 量库；然后，对描述图像内容的高维特征进行基于量子遗传优化的降维，得到 32维低维特征子集；通过相似性度量将检索结果返回给用户。系统结构框图如 图2所示，具体实施方法如下-
1. 提取高维医学图像特征
对医学图像进行检索，首先要自动抽取图像特征，建立特征索引库。灰度特 征相对几何体特征而言，具有一定的稳定性，对大小，方向都不敏感，表现出相 当强的鲁棒性，因此是医学图像检索最常用的特征之一。由于人体各器官组织的 纹理结构具有一定的区别，并且，同一器官组织的纹结构相当明显，因此，对医
学图像进行检索也常选取纹理特征。
本发明提取医学图像的灰度和纹理特征，其中，提取的HSV颜色(灰度) 特征为72维，Gabor纹理特征为48维，因此检索的联合特征向量为72+48=120 维。
2. 基于量子遗传优化的特征向量的流形降维
(1)设医学图像特征向量高维空间/^中的iV个D维随机向量k,x,，…;cj，
;c,.ei^，通过降维把它们映射到低维空间i "中，得到降维后的输出向量y,， 计算每一个样本点(特征向量)x,的A个最近邻点。对高维空间中的样本点x,,计算该点与其他iV- 1个样本点之间的距离，将距离排序，选择前/t个与x,最 近的点作为其邻近点。
(2)由每个样本点的近邻点计算出该样本点的局部重建权值矩阵『。用每个特 征向量的近邻点对该特征向量进行重建，求取使重建误差最小的近邻局部重建权 值矩阵。
定义重建拟合误差函数如下
s,)=iik，i>,~.ii2
其中，x,01,2，…,/t)为x,的j个近邻点，^,是^与^之间的权值。
在满足以下两个约束条件时，通过最小化误差函数得到局部重建权值矩阵， 即由样本点的近邻点，构造出最优『矩阵使误差函数值达到最小。
1 )每一个数据点x,都只能由它的邻近点来表示，若X,不是近邻点，则w, = 0 ;
2)权值矩阵的每一行的和为1，即满足归一化约束gw,, = 1 。
在求取最优『矩阵的过程中本发明将LLE算法中样本点与其临近点的重建 权值向量集合构造为量子遗传算法的量子染色体，然后量子遗传中的个体在寻优 的过程中找到最优位置，最终得到近邻局部重建权值矩阵。
1) 初始化种群初始化医学图像特征高维向量中包含W个个体的种群 2(r) = W,其中《(j、l,2,…A0为种群中第,代的第y个个体。种群中 全部染色体的所有基因a、 e都被初始化为1/A，这意味着一个染色体所表达 的是其全部可能状态的等概率叠加。
2) 对种群中的个体进行测量对初始种群中的个体进行一次测量以获得一 组确定的解户。=沐,其中《0:1,…,7V)是第,代种群中的第y个
个体的测量值，其表现形式为量子为数目即染色体长度为W的二进制字符串。
其中每一位为o或i是根据量子比特的概率Kf或l/ /f ( / = l，...，w )观糧得到的。
具体实现方法为随机产生
上的一个数e，若6>卜,'|2，则"中相应的位取 值为1，否则取值为0;3) 评价二进制染色体种群p(r)并保留最优解用适应度评价函数评价种群 尸(r)中每个染色体，并保留此代中的最优解。如算法获得满意解，算法停止；
否则，转入歩骤4)继续进行；
4) 量子遗传操作使用量子旋转门W"更新尸(r);
5) 改变进化代数进化代数加1，如仍未满足最大进化代数rm^，算法转
至步骤2)继续进行。
按照以上步骤即可求取最优近邻局部重建权值矩阵『。
(3)保持上步求取的权值『不变，求X,在低维嵌入空间的输出特征向量y,。为
了使输出数据在低维空间保持高维空间原有的拓扑结构需构造一个代价函数，映 射过程中必须使代价函数值达到最小。该函数形式为<formula>formula see original document page 16</formula>其中，y,是、的输出向量，力01，2,…,/t)是x的/l个近邻点。 为使低维重构误差最小化同样也要满足两个约束条件
1) y,是一个平移不变量，即|;々=0;
2) 为了避免产生退化解需使<formula>formula see original document page 16</formula>，其中N是数据点的个数，/是dxd 的单位矩阵。
此时，对s(r)进行求解，求取的最优解;;,可以使得代价函数S(10达到最小值。 构造一个WxW的稀疏矩阵『，A(,、1,2,…,W)可以存储在其中，当x,是;c,的近邻 点时，^ ,=,,否则，Wj = 0 。
此时误差代价函数式可改写为<formula>formula see original document page 16</formula>
其中，M—M,》是一个WxW的对称矩阵，其表达式为<formula>formula see original document page 16</formula>利用Lagrange乘子法可转化为在1^X =/条件下的特征值和特征向量的
问题0W —A)"'=0，艮卩Mo^/la。
要最小化代价函数，则取M的最小d个特征值对应的向量为列向量组成矩 阵r，则r的列向量即为"维空间的降维向量表示。
3. 相似性度量
对医学图像数据库中的图像和待检索图像进行特征向量降维后就可以用该 低维特征向量来表征对应的图像，然后需要对待检索图像和数据库中图像进行相 似度测量，本发明采用如下相似测度方法
(1) 直方图的交
假设/和g是两个含有iV个bin的颜色直方图，则它们之间的相交距离用下 式表示
直方图的相交是指两个直方图在每个bin中共有的像素数量。有时，该值还 可以通过除以其中一个直方图中所有的像素数量来实现标准化，从而使其值处于
的值域范围，如下式所示
丰e)《m《e》/2:》
(2) 余弦距离法
5(/,2)=/、/(||/||*||冊
其中，I和Q分别表示査询图像和数据库中图像的特征向量，1HI表示向量范数。 计算得到的相似性度量值在[O， l]之间，该值越大，表示图像越相似。
对于组合特征，相似性度量定义为各个特征相似性度量的加权和。其公式为: S(/,2) = Z:%*^(/,0，表示由w个特征组合而成，其中w,表示第j'个
特征的权重系数，它表示第j'个特征的重要性， 一般取各w,相等。又^(/,2)表
示第_/个特征的相似性度量函数值。
4. 算法的性能评价本发明采用MPEG-7推荐的检索性能评价公式来评价本系统的检索性能，设 尸w为返回的前W个结果中的正确率，又设W为某一具有特定语义含义的图像集
合。现给出示例图像A, G' = l，2，…，M)， M为测试系统检索性能而预先设定的 示例图像个数，存在《e/ ，在对《的一次检索中系统返回的前7V个结果为 A,)二1，2,…,A^，那么正确率/ w说)定义如下式所示
前iV个结果的査全率~(仏)可以用下式表示
硫n
其中，IWII表示图像集R所含的图像数。
査准率和査全率越高，表明该检索系统的效果越好。 一般地，査准率和查全 率是一对矛盾体，当要求精度较高时，査全率较低，反之亦然。因此， 一般的检 索系统只要求在这两者之间达到一个最优的平衡点，就认为达到了较好的检索性
权利要求
1、一种基于量子遗传优化的流形降维医学图像检索方法，其特征是(1)提出查询请求首先由用户通过查询接口向医学图像数据库提出查询目标图像的查询请求，所述的医学图像数据库是图像入库时，将各序列的关键图像统一转换为512×512大小的JPG格式图像构成医学图像数据库；(2)提取高维医学图像特征查询处理模块根据请求提取目标图像和医学图像数据库中每幅图像的72维灰度特征和48维纹理特征，将这些高维特征以向量的形式存入高维特征向量库；(3)基于量子遗传优化的降维对描述图像内容的高维特征进行基于量子遗传优化的降维，得到32维低维特征子集；(4)相似性度量通过相似性度量将检索结果返回给用户。
2、 根据权利要求1所述的基于量子遗传优化的流形降维医学图像检索方法，其特征是所述的基于量子遗传优化的降维的方法为(1) 医学图像特征向量高维空间i D中的iV个D维随机向量^，X2,…JcJ，X,ei 通过降维把它们映射到低维空间/^中，得到降维后的输出向量乂 ，计算每一个样本点即特征向量X,的个最近邻点；对高维空间中的样本点X,， 计算该点与其他iV-l个样本点之间的距离，将距离排序，选择前A个与X,最近 的点作为其邻近点；(2) 由每个样本点的近邻点计算出该样本点的局部重建权值矩阵『；用每 个特征向量的近邻点对该特征向量进行重建，求取使重建误差最小的近邻局部重 建权值矩阵『；G)保持上步求取的权值『不变，求x,在低维嵌入空间的输出特征向量y,;为了使输出数据在低维空间保持高维空间原有的拓扑结构需构造一个代价函 数，映射过程中必须使代价函数值达到最小，该函数形式为<formula>formula see original document page 2</formula>其中，y,是;c,的输出向量，7,,(7 = 1,2,...,*)是^的*个近邻点；为使低维重构误差最小化同样也要满足两个约束条件 1) y,是一个平移不变量，即|>,.=0;2)为了避免产生退化解需使=/ ，其中N是数据点的个数，/是"的单位矩阵；对^(y)进行求解，求取的最优解y,可以使得代价函数f()O达到最小值； 构造一个WxW的稀疏矩阵『，0、1,2,…,W)可以存储在其中，当x,是;c,的近邻点时，R,;，否则，K,o; 此时误差代价函数式可改写为=结^力,="歴7.) 其中，似=(似,,)是一个^><^的对称矩阵，其表达式为M =(/-『)7'(/-『)利用Lagrange乘子法转化为在丄^>,乂 =/条件下的特征值和特征向量的 问题(M—A)f7'=0，艮卩Mq^Aq;;要最小化代价函数，则取M的最小t/个特征值对应的向量为列向量组成矩阵r，则r的列向量即为c/维空间的降维向量表示。
3、根据权利要求2所述的基于量子遗传优化的流形降维医学图像检索方法， 其特征是所述的求取使重建误差最小的近邻局部重建权值矩阵的方法为1) 初始化种群初始化医学图像特征高维向量中包含iV个个体的种群2(r) = WU(,...，A}，其中《C/、1,2,…A0为种群中第M戈的第J个个体，种群中全部染色体的所有基因a 、 e都被初始化为i/v^，这意味着一个染色体所表达的是其全部可能状态的等概率叠加；2) 对种群中的个体进行测量对初始种群中的个体进行一次测量以获得一组确定的解尸(r)—",其中《O一l,...,AO是第H戈种群中的第乂个 个体的测量值，其表现形式为量子为数目即染色体长度为W的二进制字符串；其中每一位为O或l是根据量子比特的概率hl或|《| (/ = l,...,m)测量得到的；具体实现方法为随机产生[O,1]上的一个数0，若6>>| ,'|2，则《中相应的 位取值为1，否则取值为O;3) 评价二进制染色体种群尸(:r)并保留最优解用适应度评价函数评价种群尸(r)中每个染色体，并保留此代中的最优解；如算法获得满意解，算法停止; 否则，转入歩骤4)继续进行；4) 量子遗传操作使用量子旋转门w"更新户(:r);5) 改变进化代数进化代数加l，如仍未满足最大进化代数rm^，算法转至步骤2)继续进行。
4、根据权利要求3所述的基于量子遗传优化的流形降维医学图像检索方 法，其特征是所述的相似性度量的方法为直方图的交，/和g是两个含有W个bin的颜色直方图，它们之间的相交距离用下式表直方图的相交是指两个直方图在每个bin中共有的像素数量；有时，该值还可以通过除以其中一个直方图中所有的像素数量来实现标准化，从而使其值处于[O， l]的值域范围，如下式所示
5、根据权利要求3所述的基于量子遗传优化的流形降维医学图像检索方法，其特征是所述的相似性度量的方法为余弦距离法其中，I和Q分别表示查询图像和数据库中图像的特征向量，IH哮示向量范数， 计算得到的相似性度量值在[O， l]之间，该值越大，表示图像越相似。
6、根据权利要求3所述的基于量子遗传优化的流形降维医学图像检索方法，其特征是所述的相似性度量的方法，对于组合特征，相似性度量定义为各个特征相似性度量的加权和，其公式为s(/,2)=/"e/(imi*.s(/,0 = Z:%*^(/,g)，表示由附个特征组合而成，其中^表示第y个特征的权重系数，它表示第J'个特征的重要性， 一般取各w,相等，又S,(/，g) 表示第/个特征的相似性度量函数值。
全文摘要
本发明提供的是一种基于量子遗传优化的流形降维医学图像检索方法。首先由用户通过查询接口向医学图像数据库提出查询目标图像的查询请求；查询处理模块根据请求提取目标图像和医学图像数据库中每幅图像的72维灰度特征和48维纹理特征，将这些高维特征以向量的形式存入高维特征向量库；对描述图像内容的高维特征进行基于量子遗传优化的降维，得到32维低维特征子集；通过相似性度量将检索结果返回给用户。本发明提出了基于量子遗传优化的流形降维医学图像检索方法。在本发明中，对LLE非线性降维方法进行了改进，引入量子遗传优化局部重建权值矩阵，降低了噪声点对样本数据的影响，提高了算法收敛速度，提高了医学图像检索速度和精度。
文档编号G06F17/30GK101546332SQ20091007196
公开日2009年9月30日 申请日期2009年5月7日 优先权日2009年5月7日
发明者望 丛, 冯耀宇, 吕美超, 金 李, 杨广达, 洪 梁, 汤连志, 磊 王, 胡文广 申请人:哈尔滨工程大学