多特征整合的宫颈细胞图像特征提取与识别方法及宫颈细胞特征识别装置的制造方法

多特征整合的宫颈细胞图像特征提取与识别方法及宫颈细胞特征识别装置的制造方法
【技术领域】
[0001] 本发明属于医学细胞图像处理领域，具体设及一种多特征整合的宫颈细胞图像特 征提取与识别方法及宫颈细胞特征识别装置。
【背景技术】
[0002] 宫颈癌是女性中最常见的恶性肿瘤之一。宫颈癌从早期到晚期一般需要8-10年， 而且没有显著的症状。因此，宫颈癌的早期诊断对于及时发现宫颈癌早期病变，将宫颈癌消 灭在萌芽状态从而挽救女性患者的生命具有重要的意义。在宫颈癌的筛查方法中，宫颈细 胞学筛查方法被认为是预防宫颈癌最有效的手段之一。宫颈细胞诊断技术已经成为女性宫 颈癌临床诊断、病理分析W及治疗的重要辅助技术。临床已经证实每年检查一次宫颈癌恶 化致死率可W降低92%。人工宫颈细胞学检查是一项单调乏味、耗时、易出错的工作。随着 计算机和医疗技术的不断发展，对宫颈细胞检测提出了新的要求一一 '陕速化、高精度、智能 化。因此，开发一种计算机辅助自动宫颈细胞筛选与诊断系统十分重要，有助于帮助女性患 者预防宫颈癌的发生。
[0003] 目前的宫颈细胞图像特征提取方法是仅在细胞图像的空间域中提取其特征，该类 方法造成了部分宫颈细胞图像信息的丢失W及宫颈细胞图像特征提取的不充分，导致识别 精度低，且不利于实现对宫颈细胞图像的自动识别分类。

【发明内容】

[0004] 本发明旨在解决现有技术中存在的技术问题之一。
[0005] 1998年化ang根据单频分量信号本身的特征提出了一种提取单频分量的方 法--经验模态分解方法（EmpiricalModeDecomposition,EMD),经验模态分解能够从时 变的信号中分解出保持信号自身性质不变的单频信号分量。
[0006] 经验模态分解是一种完全自适应的分解过程。EMD方法中的信号要满足一定的条 件：
[0007] (1)数据信号中至少包含一个极大值和一个极小值两个极值点；
[000引 （2)特征尺度用极值点之间的时间间隔定义；
[0009] 做当整个数据序列没有极值点，而只有拐点时，它能够在进行一阶或多阶的求导 运算后重建极值点，最后可W将各模态积分得到各成分。
[0010] EMD方法的目的就是将数据中的不同特征尺度加W分解区分开来，把复杂 的信号分解为从高频到低频的若干阶简单的单分量固有模态函数（IntrinsicMode 化nction，IM巧组合。各IMF需要满足如下两个条件；
[ocm] (1)在整个信号长度上，一个IMF的极值点和过零点数目必须相等或至多只相差 一个；
[0012] (2)在任意时刻，由极大值点定义的上包络线和由极小值点定义的下包络线的平 均值为零。
[0013] EMD方法可W从一维推广到二维，法国J.Nunes提出了对二维信号直接分解的二 维经验模态分解方法炬idimensionalEmpiricalModeDecomposition，BEMD)。二维经验 模态分解是对EMD的进一步发展，在图像处理领域中发挥了其优良的特性，得到了广泛的 应用。本发明将二维经验模态分解方法引入到宫颈细胞图像处理领域中，为宫颈细胞图像 特征提取与识别提供了新方法。
[0014] 为此，本发明一方面提供一种多特征整合的宫颈细胞图像特征提取与识别方法， 可W充分提取宫颈细胞图像特征，提高识别精度，而且便于自动识别分类。
[0015] 本发明提供的一种多特征整合的宫颈细胞图像特征提取与识别方法，包括如下步 骤：
[0016] S100 ;对宫颈细胞图像进行二维经验模态分解，得到1个残差分量和频率由高到 低依次递减的前K个二维本征模态函数分量IMF ;
[0017]S200 ;对步骤S100中获得到的所述K个二维本征模态函数分量IMF和所述1个残 差分量进行二维傅里叶变换，开展幅频和相频特性分析W得到幅频特性和相频特性；
[001引S300 ;对步骤S200中获得到的所述幅频特性、所述相频特性W及步骤S100中获得 的所述K个二维本征模态函数分量及所述1个残差值分量作奇异值分解，获得奇异值分解 后的特征值；
[0019] S400 ;将获得的所述奇异值分解后的特征值构成宫颈细胞图像特征向量；
[0020] S500 ;依据最小欧氏距离检索方法对宫颈细胞测试样本图像的特征向量与细胞特 征数据库中的基准训练细胞图像特征进行相似性比较W进行测试细胞图像分类识别。
[0021] 在步骤100中所采用的二维经验模态分解是一种在时频域都具有良好表征信号 局部特征能力的时频分析方法，采用二维经验模态分解提取宫颈细胞图像的特征既能表征 宫颈细胞的频率信息又能兼顾空间位置信息，该分解方法是从信号本身的尺度特征出发分 解信号，是完全由数据驱动分解的方法，无需像小波变换需要选取基函数；二维经验模态分 解即具有多尺度多分辨率的优势，又具有自适应分解信号的能力。在步骤200中进行二维 傅里叶变换又可W准确提取细胞图像的幅频与相频特性，因此获得的特征多。在步骤S300 中，奇异值分解达到了对细胞图像特征向量的有效降维。奇异值分解（SingularValue Decomposition,SVD)是一种有效的代数特征提取方法，矩阵的奇异值是唯一的，同时奇异 值具有稳定性，比例不变性，旋转和平移不变性等性质。因此，二维经验模态分解、二维傅里 叶变换和奇异值分解的多特征整合识别技术对宫颈细胞图像特征的充分、精确提取，提高 了宫颈细胞识别的速度与准确率。
[0022] 进一步的，步骤S100中的二维经验模态分解过程包括如下步骤：
[002引 S110 ;利用像素点的8领域值比较法确定输入图像信号hk，i (m, n)第k级第1次筛 选后被所有的极大值max(m，n)与极小值min(m，n);
[0024] S120 ;将所有的所述极大值与所述极小值分别构造DenaulayS角形网格，由双S 次样条插值求取所述极大值构成的极大值点包络面(m，n)和所述极小值构成的极小值 点包络面emh，i(m，n);
[0025] S130;通过算术平均值计算所述极大值点包络面与所述极小值点包络面的平均包 络面；
[0026]S140 ;按公式；hk，w(m,n) =hk，i(m,n)-emean，i(m,n)计算第k个IMF分量的第 1 次 筛选值；
[0027] S150 ;计算筛选终止条件S，并判断S《T是否成立，若不成立则返回步骤S110进 行第1 + 1次筛选，直至S《T是成立，其中
T为阔值；
[002引 S160 ;获取第k个二维本征模态函数分量IMFIMFk(m，n) =hk，i(m,n)，并获取第k 个残差分量fk(m,n)=町_1(m,n)-IMFk(m,n);
[0029] S170 ;判断是否达到预定分解级数k，如果为否，则返回步骤SllO进行第k+1级第 1次筛选，直至达到预定分解级数k，完成宫颈细胞图像的二维经验模态分解，获得k个二维 本征模态函数IMF分量和1个残差分量。
[0030] 进一步的，步骤S300中的奇异值分解过程为
[0031] 首先设矩阵a为奇异值分解的矩阵，计算获得矩阵a的大小，行列数分别为 sizeal，sizea2,初始化最大循环次数loopmax= 100*max(sizeal,sizea2)，初始化计数 器loopcount= 1，初始化误差值化r=realmax，初始化矩阵S=a'，初始化矩阵V= 巧e(sizea2)，初始化矩阵U=巧e(sizeal);
[0032] 然后按如下步骤进行分解：
[003引 S310 ;对矩阵s'执行QR分解；虹S] =qr(s'），获取矩阵U=u*q，对矩阵s'再 次执行QR分解；虹S] =qr(s'），获取矩阵V=v*q;
[0034] S320 ;利用步骤S310中获取的矩阵S，提取矩阵S的上S角矩阵，计算公式如下：
[0035] e=triu(s, 1)，对所述上S角矩阵e，求取其二范数E;
[0036] S330 ;计算步骤S310中获取的矩阵S主对角线的二范数值F，判断F是否为0,判 断结果为是，F= 1，执行步骤S340 ;判断结果为否，执行步骤S340 ;
[0037] S340 ;利用步骤S320、S330中获取的二范数值E、F计算误差值化r=E/F，计数器 loopcount方口 1
[003引 S350 ;判断所述误差值Err是否小于预设误差tol或计数器计数大于预设最大循 环次数，判断结果均为是，执行步骤S360 ;否则，返回执行步骤S310[0039] S360 ;提取所述矩阵S的主对角线，生成sizeal行、sizea2列的零矩阵SS，初始化 n= 1 ；
[0040] S370 ;提取步骤S360中获取的矩阵SS的第n个值ssn，赋值S(n，n)为ssn的绝 对值，判断ssn是否为0,判断结果为是，u(:，n) = -u(:，n)，n=n+1，并执行步骤S380 ;判 断结果为否，n=n+1,并执钉步骤S380 ;
[004US380 ;判断n是否小于矩阵SS的长度，判断结果为是，输出矩阵s即为奇异值分解 的特征值；判断结果为否，返回执行步骤S370。
[0042] 进一步的，所述步骤S500中的细胞特征数据库中的基准训练细胞图像特征是由 各种不同类别状态的宫颈细胞图像按权利要求1中的步骤S100、S200、S300、S400获得的。
[0043]本发明另一方面提供宫颈细胞特征识别装置，能够快速准确的对宫颈细胞进行筛 选判断，所述宫颈细胞特征识别装置包括细胞图像输入接收模块、细胞图像分析模块、细胞 图像比较判断模块，其中所述细胞图像输入接收模块接收细胞图像的输入并输出至所述细 胞图像分析模块，所述细胞图像分析模块采用上述多特征整合的宫颈细胞图像特征提取与 识别方法中的步骤S100、S200、S300、S400对细胞图像进行分析W获得宫颈细胞图像特征 向量，所述细胞图像分析模块与所述细胞图像比较判断模块连接并将所获得宫颈细胞图像 特征向量输送至所述细胞图像比较判断模块，所述细胞图像比较判断模块包括存储有基准 训练细胞图像特征的细胞特征数据库，并按照上述多特征整合的宫颈细胞图像特征提取与 识别方法中的步骤S500对输入的所述宫颈细胞图像进行分类判断。
[0044] 进一步的，供宫颈细胞特征识别装置还包括与细胞图像比较判断模块连接用于提 示细胞图像分类判断结果的判断结果提示装置，所述判断结果提示装置包括语音或/和图 像提不。
【附图说明】
[0045] 本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变 得明显和容易理解，其中：
[0046] 图1为本发明的多特征整合的宫颈细胞图像特征提取与识别方法原理示意图；
[0047] 图2为本发明的多特征整合的宫颈细胞图像特征提取与识别方法中二维经验模 态分解原理示意图；
[0048] 图3为本发明的多特征整合的宫颈细胞图像特征提取与识别方法中奇异值分解 原理示意图；
[0049] 图4为7类宫颈细胞图像；其中图4(a)为鱗状细胞癌细胞图像（carcinoma ins;Uu)，图4(b)为轻度鱗状上皮内病变细胞图像（li曲tdysplastic)，图4(c)为 中度鱗状上皮内病变细胞图像（moderatedysplastic);图4(d)为正常柱状细胞图像 (normalcolumnar),4(e)为正常中层细胞图像（normalintermediate),图 4(f)为正常 表层细胞图像（normalsuperficiel)，图4(g)为重度鱗状上皮内病变细胞图像（severe dysplastic)；
[0化0] 图5为7类宫颈细胞图像二维经验模态分解的=级分解结果；其中，每行从左到右 分别为图4中对应类型宫颈细胞的灰度图，第1级IMF分量，第2级IMF分量，第3级IMF分量W及残差分量；
[0051] 图6为7类宫颈细胞图像灰度图与二维经验模态分解的=级分解分量的幅频特性 图。其中，每行从左到右分别为图4中对应类型宫颈细胞的灰度图幅频特性图，第1级IMF 分量幅频特性图，第2级IMF分量幅频特性图，第3级IMF分量幅频特性图W及残差分量幅 频特性图；
[0052]图7为7类宫颈细胞图像灰度图与二维经验模态分解的S级分解分量的相频特性 图。其中，每行从左到右分别为图4中对应类型宫颈细胞的灰度图相频特性图，第1级IMF 分量相频特性图，第2级IMF分量相频特性图，第3级IMF分量相频特性图W及残差分量相 频特性图；
[005引图8为7类宫颈细胞图像二维经验模态分解IMF1分量的奇异值分解提取的特征 值对比曲线图；
[0054]图9为本发明宫颈细胞识别实验曲线图；
[0化5]图10为本发明的宫颈细胞特征识别装置组成示意图。
【具体实施方式】
[0056] 为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和具体实 施方式对本发明进行进一步的详细描述。需要说明的是，在不冲突的情况下，本申请的实施 例及实施例中的特征可W相互组合。
[0057] 在下面的描述中阐述了很多具体细节W便于充分理解本发明，但是，本发明还可 W采用其他不同于在此描述的方式来实施，因此，本发明的保护范围并不受下面公开的具 体实施例的限制。
[0化引下面参照图1-3对本发明实施例的多特征整合的宫颈细胞图像特征提取与识别 方法做进一步的描述。
[0化9]如图1所示，本发明的多特征整合的宫颈细胞图像特征提取与识别方法的原理示 意图，包括如下步骤：
[0060] S100 ;对宫颈细胞图像进行二维经验模态分解，得到1个残差分量和频率由高到 低依次递减的前K个二维本征模态函数分量IMF;
[0061] S200 ;对步骤S100中获得到的所述K个二维本征模态函数分量IMF和所述1个残 差分量进行二维傅里叶变换，开展幅频和相频特性分析W得到幅频特性和相频特性；
[0062]S300;对步骤S200中获得到的所述幅频特性、所述相频特性W及步骤S100中获得 的所述K个二维本征模态函数分量及所述1个残差值分量作奇异值分解，获得奇异值分解 后的特征值；
[0063] S400 ;将获得的所述奇异值分解后的特征值构成宫颈细胞图像特征向量；
[0064] S500 ;依据最小欧氏距离检索方法对宫颈细胞测试样本图像的特征向量与细胞特 征数据库中的基准训练细胞图像特征进行相似性比较W进行测试细胞图像分类识别。
[00化]在步骤100中所采用的二维经验模态分解是一种在时频域都具有良好表征信号 局部特征能力的时频分析方法，采用二维经验模态分解提取宫颈细胞图像的特征既能表征 宫颈细胞的频率信息又能兼顾空间位置信息，该分解方法是从信号本身的尺度特征出发分 解信号，是完全由数据驱动分解的方法，无需像小波变换需要选取基函数；二维经验模态分 解即具有多尺度多分辨率的优势，又具有自适应分解信号的能力。在步骤200中进行二维 傅里叶变换又可W准确提取细胞图像的幅频与相频特性，因此获得的特征多。在步骤S300 中，奇异值分解达到了对细胞图像特征向量的有效降维。因此，二维经验模态分解、二维傅 里叶变换和奇异值分解的多特征整合识别技术对宫颈细胞图像特征的充分、精确提取，提 高了宫颈细胞识别的速度与准确率。
[0066] 步骤一中获得1个残差分量和频率由高到低依次递减的前K个二维本征模态函数 分量IMF的过程是对宫颈细胞图像进行二维经验模态分解的过程，具体的分解过程为：
[0067] 设输入宫颈细胞图像为I(m，n)，其中m，n分别为宫颈细胞图像像素点的横坐标与 纵坐标，初始化残差项r。(m，n) =I(m，n)，第k级分解第1次筛选值为hk,i(m，n)，其中k为 二维经验模态分解级数，k= 1，2,. . .，K，1为筛选次数，初始化k= 1，1 = 1，然后按W下 步骤进行分解
[0068] S110 ;利用像素点的8领域值比较法确定输入图像信号hk,i(m，n)第k级第1次筛 选后被所有的极大值max(m,n)与极小值min(m,n);
[0069]S120 ;将所有的所述极大值与所述极小值分别构造DenaulayS角形网格，由双S 次样条插值求取所述极大值构成的极大值点包络面(m，n)和所述极小值构成的极小值 点包络面emh，i(m，n);
[0070] S130;通过算术平均值计算所述极大值点包络面与所述极小值点包络面的平均包 络面，即
[0071] 514〇;按公式屯，14(111，]1)=11、1(111，]1)-6。6。。，1加，11)计算第1^个加。分量的第1次 筛选值；
[007引 S150;计算筛选终止条件S，并判断S《T是否成立，若不成立则返回步骤S110进 行第1+1次筛选，直至S《T是成立，其中
T为阔值；
[007引 S160 ;获取第k个二维本征模态函数分量IMFIMFk(m，n) =hk，i(m,n)，并获取第k 个残差分量fk(m,n)=町_1(m,n)-IMFk(m,n);
[0074] S170 ;判断是否达到预定分解级数k，如果为否，则返回步骤SllO进行第k+1级第 1次 筛选，直至达到预定分解级数k，完成宫颈细胞图像的二维经验模态分解，获得k个二维 本征模态函数IMF分量和1个残差分量。
[0075] 经过二维经验模态分解，宫颈细胞图像I(m，n)被分解成下式：
[0076]
[0077] 本实施例中，图4所示，7类宫颈细胞图像作为训练细胞图像，其中a为鱗状细胞癌 细胞图像（carcinomainsitu),b为轻度鱗状上皮内病变细胞图像（li曲tdysplastic)，c为中度鱗状上皮内病变细胞图像（moderatedysplastic) ;d为正常柱状细胞图像 (normalcolumnar)，e为正常中层细胞图像（normalintermediate)，f为正常表层细胞图 像（normalsuperficiel)，g为重度鱗状上皮内病变细胞图像（severedysplastic);具体 对7类宫颈细胞图像二维经验模态分解的=级分解，其分解结果如图5所示，其中从左至右 数，第1列为宫颈细胞的灰度图、第2列为宫颈细胞的第1级IMF分量、第3列为宫颈细胞的 第2级IMF分量、第4列为宫颈细胞的第3级IMF分量、第5列为宫颈细胞的残差分量，即 al为鱗状细胞宫颈癌细胞的灰度图，a2为鱗状细胞宫颈癌细胞的第1级IMF分量，a3为鱗 状细胞宫颈癌细胞的第2级IMF分量、a4为鱗状细胞宫颈癌细胞的第3级IMF分量、a5为 鱗状细胞宫颈癌细胞的残差分量，依次类推其他细胞图像的二维经验模态分解结果bl、b2、 b3......
[007引在步骤S200中，对宫颈细胞图像的前K个二维本征模态函数分量IMF和1个残差 分量进行幅频（Ampl;Uude化equency，AF)和相频特性（Phase化equency，P巧分析。本实 施方式中采用二维离散傅里叶变换对二维分量进行时频分析，各二维分量信号的离散傅里 叶变换形式如下式：
[0079]
[0080]
[00 川 式中，u= 0, 1，2，...，M-1，V= 0, 1，2，...，N-1，i= 1，2，...，K，M、N分别为图像 的宽与高，IMFi(m，n)、rK(m，n)分别为步骤一中获得的宫颈细胞图像的二维本征模态函数分 量IMF和残差分量，F,V，)、f；Mw，v.')分别为步骤一中获得的宫颈细胞图像的二维本征模 态函数分量IMF和残差分量二维傅里叶变换的结果。
[0082] 傅里叶变换结果为复数形式，因此厂("，v)、f;.'、("，，..')可W写成：
[0085] 本实施例中，二维本征模态函数分量IMF和残差分量的幅频特性为：
[008引二维本征模态函数分量IMF和残差分量的相频特性为：
[0091] 7类宫颈细胞图像二维经验模态分解=级分解分量的幅频特性与相频特性分别如 图6、图7所示，其中图6从左至右分别为宫颈细胞灰度图幅频特性、第1级IMF分量幅频特 性、第2级IMF分量幅频特性、第3级IMF分量幅频特性、残差分量幅频特性；图7从左至右 分别为宫颈细胞灰度图相频特性、第1级IMF分量相频特性、第2级IMF分量相频特性、第 3级IMF分量相频特性、残差分量相频特性。
[0092] 本实施例中，步骤S300中的奇异值分解过程为：
[0093] 首先设矩阵a为奇异值分解的矩阵，计算获得矩阵a的大小，行列数分别为 sizeal，sizea2,初始化最大循环次数loopmax= 100*max(sizeal,sizea2)，初始化计数 器loopcount= 1，初始化误差值化r=realmax，初始化矩阵S=a'，初始化矩阵V= eye(sizea2)，初始化矩阵U=eye(sizeal);
[0094] 然后按如下步骤进行分解：
[00巧]S310 ;对矩阵s'执行QR分解；虹S] =qr(s'），获取矩阵U=u*q，对矩阵s'再 次执行QR分解；虹S] =qr(s'），获取矩阵V=v*q;
[0096] S320;利用步骤S310中获取的矩阵S，提取矩阵S的上S角矩阵，计算公式如下：
[0097]e=triu(s，1)，对所述上S角矩阵e，求取其二范数E;
[0098] S330 ;计算步骤S310中获取的矩阵S主对角线的二范数值F，判断F是否为0,判 断结果为是，F= 1，执行步骤S340 ;判断结果为否，执行步骤S340 ;
[0099]S340 ;利用步骤S320、S330中获取的二范数值E、F计算误差值化r=E/F，计数器 loopcount方日 1
[0100] S350 ;判断所述误差值Err是否小于预设误差tol或计数器计数大于预设最大循 环次数，判断结果均为是，执行步骤S360 ;否则，返回执行步骤S310[010US360;提取所述矩阵S的主对角线，生成sizeal行、sizea2列的零矩阵SS，初始化 n= 1 ；
[0102] S370 ;提取步骤S360中获取的矩阵SS的第n个值ssn，赋值s(n，n)为ssn的绝 对值，判断ssn是否为0,判断结果为是，u(:,n) = -u(:,n)，n=n+1，并执行步骤S380 ;判 断结果为否，n=n+1，并执行步骤S380 ;
[0103]S380;判断n是否小于矩阵SS的长度，判断结果为是，输出矩阵S即为奇异值分解 的特征值；判断结果为否，返回执行步骤S370。
[0104] 奇异值分解（SingularValueDecomposition,SVD)是一种有效的代数特征提取 方法，矩阵的奇异值是唯一的，同时奇异值具有稳定性，比例不变性，旋转和平移不变性等 性质。因此，本发明结合奇异值分解方法对宫颈细胞图像的二维经验模态分解分量及其幅 频、相频特征进行特征值提取。
[0105] 本发明采用奇异值分解对宫颈细胞图像的特征提取，图8为7类宫颈细胞图像的 二维经验模态分解IMF1分量特征曲线，由图8可W看出本实施例提取的特征曲线具有较大 的类间距离和明显的区分度。
[0106] 步骤S400中，执行二维经验模态分解分量幅频、相频特征及二维经验模态分解分 量奇异值分解提取的特征值构成的宫颈细胞图像的特征向量。【具体实施方式】为：
[0107] 2；'二レ/V巧；，?化V站化,,，?尸尸,J，其中C为宫颈细胞类别，i为宫颈细胞样本编 号，IMFi_3为IMF1-3分量特征向量，Resi化e4为残差分量特征向量，AF1_4为幅频特性特征向 量，PFi_4相频特性特征向量。
[0108] 如图1至图3所示，本实施例将7类宫颈细胞图像作为基准训练细胞图像按照上 述步骤S100、S200、S300、S400获得的图像特征向量形成细胞特征数据库，将测试宫颈细胞 的细胞图像照上述步骤S100、S200、S300、S400获得的图像特征向量，
[0109] 在步骤S500中，欧氏距离巧uclideanDistance)是一个通常采用的距离定义，指 在m维空间中两个点之间的距离，本实施中采用最小欧氏距离的方法对宫颈细胞图像的检 索识别，两宫颈细胞图像的特征向量的欧氏距离越小，说明两宫颈细胞图像的相似性越大。 令检索宫颈细胞图像特征向量为乂，k为宫颈细胞类别，i为宫颈细胞图像特征维数，则宫 颈细胞图像最小欧氏距离检索方式为：
[0110]
[01川式中，X为测试宫颈细胞图像的特征向量，cell_k为最小欧氏距离即相似性最大 检索宫颈细胞类别。
[0112] 本发明实验数据来于Herlev宫颈细胞图像数据集(hUp://I油s.fme.aegean. gr/decision/downloads)，Herlev宫颈细胞图像数据集由丹麦技术大学（Technical UniversityofDenmark)和Herlev大学医院化erlevUniversityHospital)联合开发， 图像的分辨率为0. 21微米/像素，共有917个宫颈单细胞图像，数据集中包含7类宫颈单 细胞，分别为；正常柱状细胞、正常中层细胞、正常表层细胞、轻度鱗状上皮内病变细胞、中 度鱗状上皮内病变细胞、重度鱗状上皮内病变细胞、鱗状细胞癌细胞，7类宫颈单细胞样本 图像如图4所示。本发明共设计五种实验，实验一；BEMD+AF+PF+SVD宫颈细胞识别；实验 二；BEMD+AF+SVD宫颈细胞识别；实验S;BEMD+PF+SVD宫颈细胞识别；实验四；BEMD+SVD宫 颈细胞识别；实验五；AF+PF+SVD宫颈细胞识别，其中实验四、五识别方法已在图像纹理识 别中应用，实验一、二、=为本实施例宫颈细胞图像识别所设计实验。每种实验分别将10%、 20%、30%、40%、50%、60%、70%、80%、90%宫颈细胞图像作为训练样本，其余作为测试 样本进行测试 实验。实验测试结果如图9所示。实验一宫颈细胞平均识别率为87.31%，实 验二宫颈细胞平均识别率为72. 31 %，实验=宫颈细胞平均识别率为72. 30%，实验四宫颈 细胞平均识别率为66. 94%，实验五宫颈细胞平均识别率为75. 52%。由实验数据可W得出 实验一方法相比实验二、S、四、五宫颈细胞平均识别率提高了 11. 79% -20. 37%。因此，本 发明最终采用实验一识别方法即多特征整合的宫颈细胞识别方法，达到了女性宫颈癌早期 诊断，及时发现宫颈癌的早期病变的目的。
[0113] 本发明另一方面提供宫颈细胞特征识别装置，如图10所示，所述宫颈细胞特征识 别装置包括细胞图像输入接收模块10、细胞图像分析模块20、细胞图像比较判断模块30， 其中所述细胞图像输入接收模块10接收细胞图像的输入并输出至所述细胞图像分析模块 20,所述细胞图像分析模块20采用上述多特征整合的宫颈细胞图像特征提取与识别方法 中的步骤S100、S200、S300、S400对细胞图像进行分析W获得宫颈细胞图像特征向量，所述 细胞图像分析模块20与所述细胞图像比较判断模块30连接并将所获得宫颈细胞图像特征 向量输送至所述细胞图像比较判断模块30,所述细胞图像比较判断模块30包括存储有基 准训练细胞图像特征的细胞特征数据库，并按照上述多特征整合的宫颈细胞图像特征提取 与识别方法中的步骤S500对输入的所述宫颈细胞图像进行分类判断。
[0114] 另外，如图10所示，所述供宫颈细胞特征识别装置还包括与细胞图像比较判断模 块30连接用于提示细胞图像分类判断结果的判断结果提示模块40,所述判断结果提示模 块40包括语音或/和图像提示。
[0115] W上实施例仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的 技术人员来说，本发明可W有各种更改和变化。凡在本发明的创造性精神和原则之内，所作 的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种多特征整合的宫颈细胞图像特征提取与识别方法，包括如下步骤： SlOO :对宫颈细胞图像进行二维经验模态分解，得到1个残差分量和频率由高到低依 次递减的前K个二维本征模态函数分量IMF ; S200 :对步骤SlOO中获得到的所述K个二维本征模态函数分量MF和所述1个残差分 量进行二维傅里叶变换，开展幅频和相频特性分析以得到幅频特性和相频特性； S300 :对步骤S200中获得到的所述幅频特性、所述相频特性以及步骤SlOO中获得的所 述K个二维本征模态函数分量及所述1个残差值分量作奇异值分解，获得奇异值分解后的 特征值； S400 :将获得的所述奇异值分解后的特征值构成宫颈细胞图像特征向量； S500 :依据最小欧氏距离检索方法对宫颈细胞测试样本图像的特征向量与细胞特征数 据库中的基准训练细胞图像特征进行相似性比较以进行测试细胞图像分类识别。2. 如权利要求1所述的多特征整合的宫颈细胞图像特征提取与识别方法，其特征在于 步骤SlOO中的二维经验模态分解过程包括如下步骤： SllO :利用像素点的8领域值比较法确定输入图像信号htl (m，η)第k级第1次筛选后 被所有的极大值max (m, η)与极小值min (m, η); S120 :将所有的所述极大值与所述极小值分别构造Denaulay三角形网格，由双三次样 条插值求取所述极大值构成的极大值点包络面ema!U (m，η)和所述极小值构成的极小值点包 络面 emin，！ (m, η); S130:通过算术平均值计算所述极大值点包络面与所述极小值点包络面的平均包络 面； S140 :按公式：hk，1+1(m,n) = huO^rO-e^dm，]!)计算第k个IMF分量的第1次筛选 值；S150:计算筛选终止条件S，并判断S彡τ是否成立，若不成立则返回步骤SllO进行 第1+1次筛选，直至S < τ是成立，其中」 τ为阈值； S160 :获取第k个二维本征模态函数分星IMtiMtkUii, n; = nk l uii, η;，并获取第k个残 差分量 rk (m，n) = Iv1 (m，n) -IMFk (m，η); S170 :判断是否达到预定分解级数k，如果为否，则返回步骤SllO进行第k+1级第1次 筛选，直至达到预定分解级数k，完成宫颈细胞图像的二维经验模态分解，获得k个二维本 征模态函数IMF分量和1个残差分量。3. 如权利要求1所述的多特征整合的宫颈细胞图像特征提取与识别方法，其特征在于 步骤S300中的奇异值分解过程 首先设矩阵a为奇异值分解的矩阵，计算获得矩阵a的大小，行列数分别为sizeal， sizea2,初始化最大循环次数loopmax = 100*max (sizeal, sizea2)，初始化计数器 Ioopcount = 1，初始化误差值Err = realmax，初始化矩阵s = a'，初始化矩阵V = eye(sizea2)，初始化矩阵 u = eye(sizeal); 然后按如下步骤进行分解： S310 :对矩阵s'执行QR分解：[q, s] = qr (s'），获取矩阵u = u*q，对矩阵s'再次执 行 QR 分解：[q, s] = qr (s'），获取矩阵 V = v*q ; S320 :利用步骤S310中获取的矩阵s，提取矩阵s的上三角矩阵，计算公式如下： e = triu (s, 1)，对所述上三角矩阵e，求取其二范数E ; S330 :计算步骤S310中获取的矩阵s主对角线的二范数值F，判断F是否为0,判断结 果为是，F = 1，执行步骤S340 ;判断结果为否，执行步骤S340 ; S340 :利用步骤S320、S330中获取的二范数值E、F计算误差值Err = E/F，计数器 Ioopcount 加 1 S350 :判断所述误差值Err是否小于预设误差tol或计数器计数大于预设最大循环次 数，判断结果均为是，执行步骤S360 ;否则，返回执行步骤S310 S360 :提取所述矩阵s的主对角线，生成sizeal行、sizea2列的零矩阵ss，初始化η = 1 ; S370 :提取步骤S360中获取的矩阵ss的第η个值ssn，赋值s (η，η)为ssn的绝对值， 判断ssn是否为0,判断结果为是，u(:，n) = -u(:，η)，η = n+1，并执行步骤S380 ;判断结 果为否，η = n+1，并执行步骤S380 ; S380 :判断η是否小于矩阵ss的长度，判断结果为是，输出矩阵s即为奇异值分解的特 征值；判断结果为否，返回执行步骤S370。4. 如权利要求1所述的多特征整合的宫颈细胞图像特征提取与识别方法，其特征在于 步骤S500中的细胞特征数据库中的基准训练细胞图像特征是由各种不同类别状态的宫颈 细胞图像按权利要求1中的步骤S100、S200、S300、S400获得的。5. -种宫颈细胞特征识别装置，包括细胞图像输入接收模块、细胞图像分析模块、细 胞图像比较判断模块，其中所述细胞图像输入接收模块接收细胞图像的输入并输出至所述 细胞图像分析模块，所述细胞图像分析模块采用如权利要求1中的步骤S100、S200、S300、 S400对细胞图像进行分析以获得宫颈细胞图像特征向量，所述细胞图像分析模块与所述 细胞图像比较判断模块连接并将所获得宫颈细胞图像特征向量输送至所述细胞图像比较 判断模块，所述细胞图像比较判断模块包括存储有基准训练细胞图像特征的细胞特征数据 库，并按照权利要求1中的步骤S500对输入的所述宫颈细胞图像进行分类判断。6. 如权利要求5所述的宫颈细胞特征识别装置，其特征在于还包括与细胞图像比较判 断模块连接用于提示细胞图像分类判断结果的判断结果提示装置，所述判断结果提示装置 包括语音或/和图像提示。
【专利摘要】本发明提供一种多特征整合的宫颈细胞图像特征提取与识别方法及宫颈细胞特征识别装置,多特征整合的宫颈细胞图像特征提取与识别方法，包括如下步骤：对宫颈细胞图像进行二维经验模态分解；进行二维傅里叶变换，开展幅频和相频特性分析以得到幅频特性和相频特性；作奇异值分解，获得奇异值分解后的特征值；将获得的所述奇异值分解后的特征值构成宫颈细胞图像特征向量；依据最小欧氏距离检索方法对宫颈细胞测试样本图像的特征向量与细胞特征数据库中的基准训练细胞图像特征进行相似性比较以进行测试细胞图像分类识别。本发明可以充分提取宫颈细胞图像特征，提高识别精度，而且便于自动识别分类。
【IPC分类】G06K9/00, G06K9/62
【公开号】CN104881631
【申请号】CN201510181052
【发明人】陈锦, 罗晓曙, 刘艳红
【申请人】广西师范大学
【公开日】2015年9月2日
【申请日】2015年4月16日