CN103366183A - 一种局灶性病灶的非参数自动检测方法 - Google Patents

一种局灶性病灶的非参数自动检测方法 Download PDF

Info

Publication number
CN103366183A
CN103366183A CN2013103021240A CN201310302124A CN103366183A CN 103366183 A CN103366183 A CN 103366183A CN 2013103021240 A CN2013103021240 A CN 2013103021240A CN 201310302124 A CN201310302124 A CN 201310302124A CN 103366183 A CN103366183 A CN 103366183A
Authority
CN
China
Prior art keywords
pixel
zone
image
focus
local description
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2013103021240A
Other languages
English (en)
Other versions
CN103366183B (zh
Inventor
阳维
黄美燕
冯前进
佘广南
卢振泰
陈武凡
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southern Medical University
Original Assignee
Southern Medical University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southern Medical University filed Critical Southern Medical University
Priority to CN201310302124.0A priority Critical patent/CN103366183B/zh
Publication of CN103366183A publication Critical patent/CN103366183A/zh
Application granted granted Critical
Publication of CN103366183B publication Critical patent/CN103366183B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

一种局灶性病灶的非参数自动检测方法，预先以已经勾画好病灶的医学病灶图像作为样本构建样本集A1数据库，分别提取样本集A1中每一幅图像的目标区域和背景区域的局部描述子，得到样本集A1中所有目标区域局部描述子组成的目标特征库F1和样本集A1中所有背景区域局部描述子组成的背景特征库F2；处理时具体包括下列步骤：（1）将待处理的医学病灶图像I分为多个子区域；（2）用NBNN分类器对每个子区域进行分类；（3）计算目标函数，得到病灶检测结果。本发明不需预先建立参数化模型，且不需要对图像局部描述子进行量化，应用灵活，图像局部描述子鉴别能力强，能精确检测局灶性病灶。

Description

一种局灶性病灶的非参数自动检测方法
技术领域
本发明涉及医学图像中的病灶自动检测方法，具体来说涉及一种医学图像中局灶性病灶的非参数自动检测方法。
背景技术
在临床诊断中，医学图像中的病灶检测和分割十分重要，因为所提取的病灶区域可以为医生提供病灶的解剖结构信息，并为后续的治疗和病人跟踪提供依据。然而，人工检测和分割病灶区域非常费时，由于一些病灶组织与其周围的正常组织对比度不明显，不同的临床专家所分割的病灶区域可能存在差异，容易增加后续诊断和跟踪的难度。
针对人工检测和分割病灶区域的不足，研究人员提出了很多自动检测和分割病灶的方法。其中大多数自动检测病灶的方法都是基于建立参数化模型，通过训练样本来学习和调整模型中的参数，然后使用模型对图像中的像素点或区域进行分类判别，判断其是否属于病灶边界或病灶区域。这种基于参数化模型的病灶检测和分割方法有几点不足：首先，针对不同的训练样本（如：不同的病灶类型图像），该方法都需要重新训练模型中的参数，不够简便；其次，由于训练样本一般比较少，从而参数模型会产生过拟合的现象，即对于训练样本能很好地进行分类判断，但是对于样本集以外的数据分类性能较差。
相对于参数模型，非参数模型不需要经过训练阶段，因此，非参数模型使用简便，而且可以避免参数模型中的过拟合问题。现有技术中，在计算机视觉领域，基于词袋（Bag of Visual Words, BoW）模型的方法是进行图像分类的常用方法。词袋模型需要对图像局部描述子进行量化，而量化会降低图像局部描述子的鉴别能力。
因此，针对现有技术不足，提供一种能精确检测局灶性病灶的非参数自动检测方法甚为必要。
发明内容
本发明的目的在于针对现有技术不足，提供一种医学图像中局灶性病灶的非参数自动检测方法，该方法不需预先建立参数化模型，且不需要对图像局部描述子进行量化，应用灵活，图像局部描述子鉴别能力强，能精确检测局灶性病灶。
本发明的上述目的通过如下技术方案实现。
一种局灶性病灶的非参数自动检测方法，预先以已经勾画好病灶的医学病灶图像作为样本构建样本集A1数据库,分别提取样本集A1中每一幅图像的目标区域和背景区域的局部描述子，得到样本集A1中所有目标区域局部描述子组成的目标特征库F1和样本集A1中所有背景区域局部描述子组成的背景特征库F2；
处理时具体包括下列步骤：
（1）将待处理的医学病灶图像I分为多个子区域；
（2）用NBNN分类器对每个子区域进行分类；
（3）计算目标函数，得到病灶检测结果。
步骤（1）具体是将待处理的医学病灶图像I进行过分割，第i个小区域记为R i 。
优选的，步骤（1）具体是采用超像素算法将待处理的医学病灶图像I进行过分割。
步骤（2）用NBNN分类器对每个小区域进行分类，具体包括：
（2.1）对区域R i 内的每个像素点提取其局部描述子，记为f j ，用区域R i 内的所有局部描述子组成R i 的特征集                                                
Figure 2013103021240100002DEST_PATH_IMAGE001
， ，其中，n为区域R i 内局部描述子的个数；
（2.2）用NBNN分类器对区域R i 进行分类：
区域R i 分类属于某个类别C的最大后验概率为：
Figure 2013103021240100002DEST_PATH_IMAGE003
……(Ⅰ)，
其中，，+1表示目标区域，-1表示背景区域；
在朴素贝叶斯假设的条件下，区域R i 中的每一个像素点对应的局部描述子f j 是相互独立的，故：
Figure 2013103021240100002DEST_PATH_IMAGE005
……(Ⅱ)；
将式(Ⅱ)中
Figure 943743DEST_PATH_IMAGE006
用高斯核密度估计进行近似，令k=1，可得：
Figure 2013103021240100002DEST_PATH_IMAGE007
……(Ⅲ)；
其中，NN C (f j )表示在类别C特征库中与f j 最相近的局部描述子，σ表示标准偏差；
       将式(Ⅲ)中的
Figure 47965DEST_PATH_IMAGE008
代入式(Ⅱ)，得：
Figure DEST_PATH_IMAGE009
……(Ⅳ)；
定义
Figure 735298DEST_PATH_IMAGE010
为区域R i 到某个类别C的距离；
       NBNN分类算法的目标就是找到使区域R i 到某个类别C距离最小的那个类别，即NBNN分类器对区域R i 的分类结果为：
Figure DEST_PATH_IMAGE011
  。
步骤（3）具体为：定义目标函数
Figure 689610DEST_PATH_IMAGE012
，对于图像I进行计算，选出的区域作为病灶检测的结果。
优选的，具体通过如下方式获得样本集A1的目标特征库F1和背景特征库F2：
（0.1）定义作为样本的医学病灶图像中的病灶区域为目标区域，病灶以外区域为背景区域，对目标区域内的每个像素点提取其局部描述子；对背景区域内的每个像素点同样提取其局部描述子，其中，像素点对应的局部描述子的提取过程如下所述： 
（0.1.1）以该像素点为圆心，分别按照半径R=﹛20、21、22、23、24、……、2n﹜像素作圆，其中n为自然数，得到一系列半径在对数尺度等间隔变化的同心圆，令最大的同心圆的半径
Figure 298446DEST_PATH_IMAGE014
像素，得到5个同心圆；
(0.1.2)以所述圆心为起点沿 8个方向画直线，8个方向的直线将圆面均分为八等份，并取这8条直线与所述的5个同心圆的相交的像素点，把圆心以及这些与同心圆相交的像素点以从圆心向外的顺序、并以逆时针方向排列成一个队列，然后将队列中，位于图像内的像素点分别替换为所对应的灰度值，位于图像外的像素点的灰度值分别赋予零，得到该像素点的灰度特征向量；
(0.1.3)获取该像素点在图像上的坐标(x, y)，对(x, y)进行归一化处理得到归一化的坐标
Figure DEST_PATH_IMAGE015
，即，其中，w和h分别为图像的宽和长，把步骤（0.1.2）中所得到的灰度特征向量与坐标
Figure DEST_PATH_IMAGE017
首尾连接排列成新的向量，得到该像素点的局部描述子；
（0.2）把步骤（0.1）中所产生的样本集A 1 中所有目标区域的局部描述子组成目标特征库F 1 ，同样，把所有背景区域的局部描述子组成背景特征库F 2 。
本发明的一种局灶性病灶的非参数自动检测方法，具有以下的有益效果：
    由于目标特征库F 1 和背景特征库F 2 都是直接从样本集A 1 中每幅图像的每个像素点提取局部描述子集合得到的，不需要经过局部描述子量化，从而完整地保留了局部描述子的鉴别能力。
    本发明使用非参数的NBNN分类方法来对图像进行分类，避免了参数模型的建立，以及参数的训练过程，从而增加了临床使用的便捷性。
    本发明以超像素（superpixel）为基本单元来判断区域的类别，相比于以像素点为单位的分类方法，更能获得稳定的类别估计。
故，本发明的一种医学图像中局灶性病灶的非参数自动检测方法，不需预先建立参数化模型，且不需要对图像局部描述子进行量化，应用灵活，图像局部描述子鉴别能力强，能精确检测局灶性病灶。
附图说明
利用附图对本发明作进一步的说明，但附图中的内容不构成对本发明的任何限制。
图1为本发明一种局灶性病灶的非参数自动检测方法的流程示意图；
图2为本发明方法用于肝部CT图像局灶性病灶检测的示例图，其中，编号为Ⅰ的轮廓线是用本发明的非参数模型得到的局灶性病灶的检测结果；编号为Ⅱ的轮廓线是局灶性病灶的手工勾画结果。
具体实施方式
下面结合具体的实施例对本发明进行详细描述。
实施例1。
本发明的一种局灶性病灶的非参数自动检测方法，如图1所示，预先以已经勾画好病灶的医学病灶图像作为样本构建样本集A1数据库,分别提取样本集A1中每一幅图像的目标区域和背景区域的局部描述子，得到样本集A1中所有目标区域局部描述子组成的目标特征库F1和样本集A1中所有背景区域局部描述子组成的背景特征库F2。
具体通过如下方式获得样本集A1的目标特征库F1和背景特征库F2：
（0.1）定义作为样本的医学病灶图像中的病灶区域为目标区域，病灶以外区域为背景区域，对目标区域内的每个像素点提取其局部描述子；对背景区域内的每个像素点同样提取其局部描述子。其中，像素点对应的局部描述子的提取过程如下所述： 
（0.1.1）以该像素点为圆心，分别按照半径R=﹛20、21、22、23、24、……、2n﹜像素作圆，其中n为自然数，得到一系列半径在对数尺度等间隔变化的同心圆，所取半径越大，相应的局部描述子的描述能力越强，但所产生的特征维数也越高，会导致计算量增加，为了平衡计算量和描述力，令最大的同心圆的半径
Figure 685882DEST_PATH_IMAGE018
像素，得到5个同心圆；
(0.1.2)以所述圆心为起点沿 8个方向画直线，8个方向的直线将圆面均分为八等份，并取这8条直线与（0.1.1）所述的5个同心圆的相交的像素点，把圆心以及这些与同心圆相交的像素点以从圆心向外的顺序、并以逆时针方向排列成一个队列，然后将队列中，位于图像内的像素点分别替换为所对应的灰度值，位于图像外的像素点的灰度值分别赋予零，得到该像素点的灰度特征向量；
(0.1.3)获取该像素点在图像上的坐标(x, y)，对(x, y)进行归一化处理得到归一化的坐标
Figure DEST_PATH_IMAGE019
，即
Figure 653838DEST_PATH_IMAGE020
，其中，w和h分别为图像的宽和长，把步骤（0.1.2）中所得到的灰度特征向量与坐标
Figure DEST_PATH_IMAGE021
首尾连接排列成新的向量，得到该像素点的局部描述子；
（0.2）把步骤（0.1）中所产生的样本集A 1 中所有目标区域的局部描述子组成目标特征库F 1 ，同样，把所有背景区域的局部描述子组成背景特征库F 2 。
处理时具体包括下列步骤：
（1）将待处理的医学病灶图像I分为多个子区域；
具体是将待处理的医学病灶图像I进行过分割，优选采用超像素算法将待处理的医学病灶图像I进行过分割 ，过分割即把I分成一系列大小几乎一致的，连续同质的小区域，第i个小区域记为R i ，一般来说，I中的病灶区域都包含1个或以上相连的小区域。
（2）用NBNN分类器对每个子区域进行分类。
步骤（2）用NBNN分类器对每个小区域进行分类，具体包括：
（2.1）对区域R i 内的每个像素点提取其局部描述子，记为f j ，用区域R i 内的所有局部描述子组成R i 的特征集
Figure 15549DEST_PATH_IMAGE022
， ，其中，n为区域R i 内局部描述子的个数；
（2.2）用NBNN分类器对区域R i 进行分类：
区域R i 分类属于某个类别C的最大后验概率为：
Figure 117366DEST_PATH_IMAGE024
……(Ⅰ)，
其中，
Figure DEST_PATH_IMAGE025
，+1表示目标区域，-1表示背景区域；
在朴素贝叶斯假设的条件下，区域R i 中的每一个像素点对应的局部描述子f j 是相互独立的，故：
Figure 779292DEST_PATH_IMAGE026
……(Ⅱ)；
将式(Ⅱ)中
Figure 918149DEST_PATH_IMAGE006
用高斯核密度估计进行近似，其中，当类别C中的样本数N c 足够大时，
Figure DEST_PATH_IMAGE027
可以进一步用k近邻来近似，在NBNN分类算法中，令k=1，可得：
Figure 908102DEST_PATH_IMAGE028
……(Ⅲ)；
其中，NN C (f j )表示在类别C特征库中与f j 最相近的局部描述子，σ表示标准偏差；
       将式(Ⅲ)中的
Figure DEST_PATH_IMAGE029
代入式(Ⅱ)，得：
……(Ⅳ)；
定义
Figure 674249DEST_PATH_IMAGE010
为区域R i 到某个类别C的距离；
       NBNN分类算法的目标就是找到使区域R i 到某个类别C距离最小的那个类别，分别计算区域R i 中所包含的局部描述子f j 到目标特征库F 1 和背景特征库F 2 中的各个局部描述子的距离，即计算
Figure DEST_PATH_IMAGE031
和
Figure 984008DEST_PATH_IMAGE032
。 
由于目标特征库F 1 和背景特征库F 2 所包含的局部描述子的数量是非常巨大的，如果对R i 中每个局部描述子都要分别找到其在F 1 和F 2 中最接近的局部描述子，是十分费时的，因此，采用随机KD-树来解决这个问题；
       NBNN分类算法的目标就是找到使区域R i 到某个类别C距离最小的那个类别，即NBNN分类器对区域R i 的分类结果为：
Figure DEST_PATH_IMAGE033
  。
（3）计算目标函数，得到病灶检测结果。
步骤（3）具体为：定义目标函数
Figure 209060DEST_PATH_IMAGE034
，对于图像I进行计算，选出
Figure DEST_PATH_IMAGE035
的区域作为病灶检测的结果。
本发明的一种局灶性病灶的非参数自动检测方法， 
（1）由于目标特征库F 1 和背景特征库F 2 都是直接从样本集A 1 中每幅图像的每个像素点提取局部描述子集合得到的，不需要经过局部描述子量化，从而完整地保留了局部描述子的鉴别能力；
（2）本发明使用非参数的NBNN分类方法来对图像进行分类，避免了参数模型的建立，以及参数的训练过程，从而增加了临床使用的便捷性；
（3）本发明以超像素（superpixel）为基本单元来判断区域的类别，相比于以像素点为单位的分类方法，更能获得稳定的类别估计。
故，本发明的一种医学图像中局灶性病灶的非参数自动检测方法，不需预先建立参数化模型，且不需要对图像局部描述子进行量化，应用灵活，图像局部描述子鉴别能力强，能精确检测局灶性病灶。
实施例2。
以一具体实施例对本发明的方法进行说明。
本实施例所使用的数据库存有458幅肝部CT图像样本集，其中包括：肝癌178幅、肝囊肿98幅、肝血管瘤182幅，且每张CT图像已人工勾画好了病灶轮廓。
以下参照图1详细描述通过所述数据库中CT图像的局灶性病灶的非参数自动检测方法。
读取所述数据库中的已勾画好病灶的CT图像，把所有图像分为2份，其中一份作为建立特征库的样本集，记为样本集A 1 ，剩下一份作为测试样本集，记为样本集A 2 。
对样本集A 1 中的每一幅CT图像进行以下处理：
首先对图像进行预处理，其过程为：由于在CT图像中，肝脏及其病灶的灰度值在[-50,250]HU的范围内，因此，可以用阈值化方法把图像中的空气、脂肪和骨头等部分去掉；把预处理后剩下的最大连通区域作为后续计算的感兴趣区域，经过预处理这一步骤，可以减少后续计算的计算量。
（0.1）定义CT图像中的病灶区域为目标区域，病灶以外、感兴趣区域以内的区域为背景区域，对目标区域内的每个像素点提取其局部描述子，同样，对背景区域内的每个像素点提取其局部描述子；其中，某个像素点对应的局部描述子的提取过程如下所述： 
（0.1.1）以该像素点为圆心，分别按照半径R=﹛20、21、22、23、24、……、2n﹜像素作圆，其中n为自然数，得到一系列半径在对数尺度等间隔变化的同心圆，所取半径越大，相应的局部描述子的描述能力越强，但所产生的特征维数也越高，会导致计算量增加，为了平衡计算量和描述力，我们令最大的同心圆的半径
Figure 793625DEST_PATH_IMAGE036
像素，可得到5个同心圆。
（0.1.2）以步骤（0.1.1）所述的圆心为起点，向东、西、南、北、东北、西北、西南、东南等8个方向画直线，并取这8条直线与步骤（0.1.1）所述的5个同心圆的相交的像素点，把圆心点以及这些与同心圆相交的像素点以从圆心向外的顺序，并以逆时针方向排列成一个队列，然后将队列中，位于图像内的像素点分别替换为所对应的灰度值，位于图像外的像素点的灰度值分别赋予零，从而得到该点的灰度特征向量。
（0.1.3）获取该像素点在图像上的坐标(x, y)，对(x, y)进行归一化处理得到归一化的坐标，即
Figure DEST_PATH_IMAGE037
，其中，w和h分别为图像的宽和长，把步骤（0.1.2）中所得到的灰度特征向量与坐标
Figure 520589DEST_PATH_IMAGE038
首尾连接排列成新的一个向量，便得到该像素点的局部描述子。
（0.2）把步骤（0.1）中所产生的样本集A 1 中所有目标区域的局部描述子组成目标特征库F 1 ，同样，把所有背景区域的局部描述子组成背景特征库F 2 。
对样本集A 2 的任一幅CT图像，用NBNN方法对病灶进行自动检测；具体包括：
首先把该CT图像记为I，对I进行预处理，用阈值化方法把图像中的空气、脂肪和骨头等部分去掉；把预处理后剩下的最大连通区域作为后续计算的感兴趣区域ROI。
（1）对图像I实行过分割，过分割可以采用超像素算法（superpixel algorithms），即把I分成一系列大小几乎一致的，连续同质的小区域，第i个小区域记为R i ，一般来说，I中的病灶区域都包含1个或以上相连的小区域。
（2）在此，我们只针对I中的ROI区域，对ROI区域中的每个小区域进行以下处理：
（2.1）对区域R i 内的每个像素点提取其局部描述子，记为f j ，用区域R i 内的所有局部描述子组成R i 的特征集， ，其中，n为区域R i 内局部描述子的个数；
（2.2）用NBNN分类器对区域R i 进行分类：
区域R i 分类属于某个类别C的最大后验概率为：
……(Ⅰ)，
其中，
Figure DEST_PATH_IMAGE041
，+1表示目标区域，-1表示背景区域；
在朴素贝叶斯假设的条件下，区域R i 中的每一个像素点对应的局部描述子f j 是相互独立的，故：
Figure 551365DEST_PATH_IMAGE042
……(Ⅱ)；
将式(Ⅱ)中
Figure DEST_PATH_IMAGE043
用高斯核密度估计进行近似，其中，当类别C中的样本数N c 足够大时，
Figure 202926DEST_PATH_IMAGE044
可以进一步用k近邻来近似，在NBNN分类算法中，令k=1，可得：
Figure DEST_PATH_IMAGE045
……(Ⅲ)；
其中，NN C (f j )表示在类别C特征库中与f j 最相近的局部描述子，σ表示标准偏差；
       将式(Ⅲ)中的代入式(Ⅱ)，得：
Figure 174610DEST_PATH_IMAGE046
……(Ⅳ)；
定义
Figure DEST_PATH_IMAGE047
为区域R i 到某个类别C的距离；
       NBNN分类算法的目标就是找到使区域R i 到某个类别C距离最小的那个类别，分别计算区域R i 中所包含的局部描述子f j 到目标特征库F 1 和背景特征库F 2 中的各个局部描述子的距离，即计算和
Figure DEST_PATH_IMAGE049
。 
       由于目标特征库F 1 和背景特征库F 2 所包含的局部描述子的数量是非常巨大的，如果对R i 中每个局部描述子都要分别找到其在F 1 和F 2 中最接近的局部描述子，是十分费时的，因此，采用随机KD-树来解决这个问题；
       NBNN分类算法的目标就是找到使区域R i 到某个类别C距离最小的那个类别，即NBNN分类器对区域R i 的分类结果为：
Figure 749128DEST_PATH_IMAGE050
  。
（3）计算目标函数，得到病灶检测结果。
步骤（3）具体为：定义目标函数
Figure DEST_PATH_IMAGE051
，对于图像I进行计算，选出
Figure 16162DEST_PATH_IMAGE052
的区域作为病灶检测的结果。
 图2是本发明用于肝部CT图像局灶性病灶检测的结果示例图，其中，编号为Ⅰ的轮廓线是用本发明所述的非参数模型得到的局灶性病灶的检测结果；编号为Ⅱ的轮廓线是局灶性病灶的手工勾画结果。从图2可以看出，使用本发明的方法能够比较精确地检测局灶性病灶。
 本发明的一种局灶性病灶的非参数自动检测方法， 
（1）由于目标特征库F 1 和背景特征库F 2 都是直接从样本集A 1 中每幅图像的每个像素点提取局部描述子集合得到的，不需要经过局部描述子量化，从而完整地保留了局部描述子的鉴别能力；
（2）本发明使用非参数的NBNN分类方法来对图像进行分类，避免了参数模型的建立，以及参数的训练过程，从而增加了临床使用的便捷性；
（3）本发明以超像素（superpixel）为基本单元来判断区域的类别，相比于以像素点为单位的分类方法，更能获得稳定的类别估计。
故，本发明的一种医学图像中局灶性病灶的非参数自动检测方法，不需预先建立参数化模型，且不需要对图像局部描述子进行量化，应用灵活，图像局部描述子鉴别能力强，能精确检测局灶性病灶。
  最后应当说明的是，以上实施例仅用以说明本发明的技术方案而非对本发明保护范围的限制，尽管参照较佳实施例对本发明作了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的实质和范围。

Claims (6)

1.一种局灶性病灶的非参数自动检测方法，其特征在于： 
预先以已经勾画好病灶的医学病灶图像作为样本构建样本集A1数据库,分别提取样本集A1中每一幅图像的目标区域和背景区域的局部描述子，得到样本集A1中所有目标区域局部描述子组成的目标特征库F1和样本集A1中所有背景区域局部描述子组成的背景特征库F2；
处理时具体包括下列步骤：
（1）将待处理的医学病灶图像I分为多个子区域；
（2）用NBNN分类器对每个子区域进行分类；
（3）计算目标函数，得到病灶检测结果。
2.根据权利要求1所述的局灶性病灶的非参数自动检测方法，其特征在于：步骤（1）具体是将待处理的医学病灶图像I进行过分割，第i个小区域记为R i 。
3.根据权利要求2所述的局灶性病灶的非参数自动检测方法，其特征在于：步骤（1）具体是采用超像素算法将待处理的医学病灶图像I进行过分割。
4. 根据权利要求3所述的局灶性病灶的非参数自动检测方法，其特征在于：所述步骤（2）用NBNN分类器对每个小区域进行分类，具体包括：
（2.1）对区域R i 内的每个像素点提取其局部描述子，记为f j ，用区域R i 内的所有局部描述子组成R i 的特征集                                                ， 
Figure 984528DEST_PATH_IMAGE002
，其中，n为区域R i 内局部描述子的个数；
（2.2）用NBNN分类器对区域R i 进行分类：
区域R i 分类属于某个类别C的最大后验概率为：
Figure 998621DEST_PATH_IMAGE003
……(Ⅰ)，
其中，，+1表示目标区域，-1表示背景区域；
在朴素贝叶斯假设的条件下，区域R i 中的每一个像素点对应的局部描述子f j 是相互独立的，故：
Figure 53481DEST_PATH_IMAGE005
……(Ⅱ)；
将式(Ⅱ)中
Figure 978712DEST_PATH_IMAGE006
用高斯核密度估计进行近似，令k=1，可得：
Figure 203020DEST_PATH_IMAGE007
……(Ⅲ)；
其中，NN C (f j )表示在类别C特征库中与f j 最相近的局部描述子，σ表示标准偏差；
       将式(Ⅲ)中的
Figure 356528DEST_PATH_IMAGE008
代入式(Ⅱ)，得：
Figure 640878DEST_PATH_IMAGE009
……(Ⅳ)；
定义为区域R i 到某个类别C的距离；
       NBNN分类算法的目标就是找到使区域R i 到某个类别C距离最小的那个类别，即NBNN分类器对区域R i 的分类结果为：
Figure 550246DEST_PATH_IMAGE011
  。
5.根据权利要求4所述的局灶性病灶的非参数自动检测方法，其特征在于：所述步骤（3）具体为：定义目标函数
Figure 731828DEST_PATH_IMAGE012
，对于图像I进行计算，选出
Figure 921501DEST_PATH_IMAGE013
的区域作为病灶检测的结果。
6.根据权利要求1至5任意一项所述的局灶性病灶的非参数自动检测方法，其特征在于：
具体通过如下方式获得样本集A1的目标特征库F1和背景特征库F2：
（0.1）定义作为样本的医学病灶图像中的病灶区域为目标区域，病灶以外区域为背景区域，对目标区域内的每个像素点提取其局部描述子；对背景区域内的每个像素点同样提取其局部描述子，其中，像素点对应的局部描述子的提取过程如下所述： 
  （0.1.1）以该像素点为圆心，分别按照半径R=﹛20、21、22、23、24、……、2n﹜像素作圆，其中n为自然数，得到一系列半径在对数尺度等间隔变化的同心圆，令最大的同心圆的半径
Figure 414799DEST_PATH_IMAGE014
像素，得到5个同心圆；
   (0.1.2)以所述圆心为起点沿 8个方向画直线，8个方向的直线将圆面均分为八等份，并取这8条直线与所述的5个同心圆的相交的像素点，把圆心以及这些与同心圆相交的像素点以从圆心向外的顺序、并以逆时针方向排列成一个队列，然后将队列中，位于图像内的像素点分别替换为所对应的灰度值，位于图像外的像素点的灰度值分别赋予零，得到该像素点的灰度特征向量；
   (0.1.3)获取该像素点在图像上的坐标(x, y)，对(x, y)进行归一化处理得到归一化的坐标，即
Figure 751420DEST_PATH_IMAGE016
，其中，w和h分别为图像的宽和长，把步骤（0.1.2）中所得到的灰度特征向量与坐标
Figure 111994DEST_PATH_IMAGE017
首尾连接排列成新的向量，得到该像素点的局部描述子；
（0.2）把步骤（0.1）中所产生的样本集A 1 中所有目标区域的局部描述子组成目标特征库F 1 ，同样，把所有背景区域的局部描述子组成背景特征库F 2 。
CN201310302124.0A 2013-07-19 2013-07-19 一种局灶性病灶的非参数自动检测方法 Active CN103366183B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310302124.0A CN103366183B (zh) 2013-07-19 2013-07-19 一种局灶性病灶的非参数自动检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310302124.0A CN103366183B (zh) 2013-07-19 2013-07-19 一种局灶性病灶的非参数自动检测方法

Publications (2)

Publication Number Publication Date
CN103366183A true CN103366183A (zh) 2013-10-23
CN103366183B CN103366183B (zh) 2017-05-10

Family

ID=49367483

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310302124.0A Active CN103366183B (zh) 2013-07-19 2013-07-19 一种局灶性病灶的非参数自动检测方法

Country Status (1)

Country Link
CN (1) CN103366183B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103886318A (zh) * 2014-03-31 2014-06-25 武汉天仁影像科技有限公司 尘肺病大体成像中病灶区域的提取与分析方法
CN106462974A (zh) * 2014-06-12 2017-02-22 皇家飞利浦有限公司 用于分割图像的参数优化
CN108447046A (zh) * 2018-02-05 2018-08-24 龙马智芯（珠海横琴）科技有限公司 病灶的检测方法和装置、设备、计算机可读存储介质
CN111369576A (zh) * 2020-05-28 2020-07-03 腾讯科技（深圳）有限公司 图像分割模型的训练方法、图像分割方法、装置及设备

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090016589A1 (en) * 2007-07-10 2009-01-15 Siemens Medical Solutions Usa, Inc. Computer-Assisted Detection of Colonic Polyps Using Convex Hull
CN101669828A (zh) * 2009-09-24 2010-03-17 复旦大学 基于pet/ct图像纹理特征的肺部恶性肿瘤与良性结节检测系统
CN101900737A (zh) * 2010-06-10 2010-12-01 上海理工大学 基于支持向量机的尿沉渣有形成分自动识别系统
US20110103656A1 (en) * 2009-04-17 2011-05-05 Gheorghe Iordanescu Quantification of Plaques in Neuroimages
CN102521618A (zh) * 2011-11-11 2012-06-27 北京大学 局部描述子的提取方法、图片检索方法及图像匹配方法
CN102945289A (zh) * 2012-11-30 2013-02-27 苏州搜客信息技术有限公司 基于cgci-sift局部特征的图像检索方法
CN103150611A (zh) * 2013-03-08 2013-06-12 北京理工大学 Ii型糖尿病发病概率分层预测方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090016589A1 (en) * 2007-07-10 2009-01-15 Siemens Medical Solutions Usa, Inc. Computer-Assisted Detection of Colonic Polyps Using Convex Hull
US20110103656A1 (en) * 2009-04-17 2011-05-05 Gheorghe Iordanescu Quantification of Plaques in Neuroimages
CN101669828A (zh) * 2009-09-24 2010-03-17 复旦大学 基于pet/ct图像纹理特征的肺部恶性肿瘤与良性结节检测系统
CN101900737A (zh) * 2010-06-10 2010-12-01 上海理工大学 基于支持向量机的尿沉渣有形成分自动识别系统
CN102521618A (zh) * 2011-11-11 2012-06-27 北京大学 局部描述子的提取方法、图片检索方法及图像匹配方法
CN102945289A (zh) * 2012-11-30 2013-02-27 苏州搜客信息技术有限公司 基于cgci-sift局部特征的图像检索方法
CN103150611A (zh) * 2013-03-08 2013-06-12 北京理工大学 Ii型糖尿病发病概率分层预测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
于梅等: "非张量积小波的肝脏CT图像检索", 《计算机工程与应用》 *
赵敏等: "K-means与朴素贝叶斯在商务智能中的应用", 《计算机技术与发展》 *
赵玉霞等: "贝叶斯方法在玉米叶部病害图像识别中的应用", 《计算机工程与应用》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103886318A (zh) * 2014-03-31 2014-06-25 武汉天仁影像科技有限公司 尘肺病大体成像中病灶区域的提取与分析方法
CN103886318B (zh) * 2014-03-31 2017-03-01 武汉天仁影像科技有限公司 尘肺病大体成像中病灶区域的提取与分析方法
CN106462974A (zh) * 2014-06-12 2017-02-22 皇家飞利浦有限公司 用于分割图像的参数优化
CN106462974B (zh) * 2014-06-12 2020-07-28 皇家飞利浦有限公司 用于分割图像的参数优化
CN108447046A (zh) * 2018-02-05 2018-08-24 龙马智芯（珠海横琴）科技有限公司 病灶的检测方法和装置、设备、计算机可读存储介质
CN111369576A (zh) * 2020-05-28 2020-07-03 腾讯科技（深圳）有限公司 图像分割模型的训练方法、图像分割方法、装置及设备

Also Published As

Publication number Publication date
CN103366183B (zh) 2017-05-10

Similar Documents

Publication Publication Date Title
US10096108B2 (en) Medical image segmentation method and apparatus
Yao et al. An improved K-means clustering algorithm for fish image segmentation
Wang et al. Hard exudate detection based on deep model learned information and multi-feature joint representation for diabetic retinopathy screening
Jia et al. Visual tracking via adaptive structural local sparse appearance model
Shao et al. A detection approach for solitary pulmonary nodules based on CT images
Zawbaa et al. An automatic flower classification approach using machine learning algorithms
Jing et al. Lung nodule classification combining rule-based and SVM
Deng et al. Graph cut based automatic aorta segmentation with an adaptive smoothness constraint in 3D abdominal CT images
CN108427969A (zh) 一种多尺度形态学结合卷积神经网络的纸张缺陷分类方法
Wang et al. Nucleus segmentation of cervical cytology images based on depth information
Yan et al. Improved mask R-CNN for lung nodule segmentation
Guo et al. Hair segmentation using heuristically-trained neural networks
CN103366183A (zh) 一种局灶性病灶的非参数自动检测方法
Kumar et al. An unsupervised approach for overlapping cervical cell cytoplasm segmentation
Liu et al. Extracting lungs from CT images via deep convolutional neural network based segmentation and two-pass contour refinement
CN110634118A (zh) 基于人工智能的乳腺影像识别系统及方法
CN110211117B (zh) 医学影像中识别线管状物的处理系统及优化分割的方法
Sarath et al. A two-stage multiple instance learning framework for the detection of breast cancer in mammograms
Akbar et al. Tumor localization in tissue microarrays using rotation invariant superpixel pyramids
Devasia et al. Automatic Optic Disc Localization and Segmentation using Swarm Intelligence.
Lyu et al. HRED-net: high-resolution encoder-decoder network for fine-grained image segmentation
Cheng et al. Improved faster RCNN for white blood cells detection in blood smear image
Kim et al. Quantification of hepatorenal index for computer-aided fatty liver classification with self-organizing map and fuzzy stretching from ultrasonography
de Brito Silva et al. Classification of breast masses in mammograms using geometric and topological feature maps and shape distribution
Gomez et al. Feature selection and classifier performance in computer-aided diagnosis for breast ultrasound

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant