专利名称：基于目标跟踪的计算机辅助胃癌诊断方法
技术领域：
本发明属于图像处理技术领域，更进一步涉及一种计算机视觉领域中处理胃癌诊断中的淋巴结转移识别领域中的基于目标跟踪的计算机辅助胃癌诊断方法。本发明可用于胃癌诊断中淋巴结转移的判断、识别过程，完成淋巴结转移的判断、识别、检测，更好的辅助医学研究中的临床诊断。
背景技术：
目前，胃癌影像研究主要是基于传统的检测方法将被检查的样本通过一系列诸如稀释、染色等化学处理之后制作成目标涂片；有着丰富经验的兵力医生在显微镜下自习观察涂片中目标核、目标浆的形貌、颜色特点，根据相关的医学标准并结合自身的经验，对待检查的样本作出诊断结论，在普查中，待检查人群分布地区广、数量大，要有效地进行健康普查存在以下几个问题1)工作量非常大，医生容易疲劳而引起误判；幻工作效率低下， 容易延误诊断结果；3)检查结果往往容易过多加入医生的主观因素，诊断结果因人而异； 4)我国中小城市一些医生缺乏临床经验和专业知识，诊断水平仍令人堪忧；因此，为了更好的开展防癌普查工作，尽量做到早期诊断、早期治疗、检查能有效花费，必须研制适合我国国情的辅助诊断系统。希森美康株式会社申请的专利“胃癌淋巴结转移的判断方法”(申请号 200810000445. 4，公开号CN101260428)公开一种胃癌淋巴结转移的判断方法，包括1)定量步骤，定量用疑有胃癌转移的淋巴结组织制备的检测试样中的角蛋白19mRNAJ)判断步骤，根据所得所述mRNA的定量结果判断胃癌的淋巴结转移，当定量值低于阈值时，判断胃癌淋巴结转移阴性，当定量值超过阈值时，判断胃癌淋巴结转移阳性。该方法存在的不足是，对于阈值的选取会因为实验条件、环境因素、专业素养等因素而有所差异，会影响诊断结果的有效性。威里德克斯有限责任公司申请的专利“对胃癌淋巴结转移的检测”(申请号 200710100645. 2，公开号CN101(^9342)公开的胃癌淋巴结转移检测方法1)从患者获取生物学样品；幻测量样品中的标记水平，其中该水平高于预定截止值水平值的认为是胃癌淋巴结转移。该方法存在的不足是，人工设定截止水平值限定了这种方法的有效性，而且截止水平值会因为选取的目标样本不同而有所差异，最终导致结果不够准确有效。
发明内容
本发明的目的在于克服上述现有技术中的不足，提出一种基于目标跟踪的计算机辅助胃癌诊断方法。该方法运用卡尔曼框架下特征匹配的目标跟踪方法，用于检测胃部切片图当中的淋巴结，从而辅助胃癌的诊断。实现本发明目的的技术思路是，对胃部切片图序列中每幅图进行感兴趣区域提取、分割，去除其他不相关区域，得到的切片图中仅存在淋巴结、血管和杂点，运用目标跟踪的方法，对连续切片图中的目标进行逐个跟踪，得到每个目标的运动轨迹段，前后帧面积变化较大且其中心位置变化小于一定阈值的目标判定为淋巴结，将前后帧面积变化很小且中心位置变化较大的目标判定为血管，其余判定为杂点，从而实现淋巴结转移的判断。本发明的具体实现步骤包括如下1、一种基于目标跟踪的计算机辅助胃癌诊断方法，其具体实现步骤如下(1)分割图像la)采用KSVD字典学习的方法提取胃壁周围感兴趣区域；lb)运用字典学习的稀疏表示方法检测每帧切片图中的目标。(2)提取被跟踪目标2a)采用四邻域连通标记算法，对所有待跟踪的胃部切片图像中的目标区域依次进行中心标记；2b)提取每个目标区域中心标记，得到每个目标的中心位置，将每一帧中所有目标的中心位置确定为被跟踪目标的中心位置。⑶预测跟踪3a)按照中心标记的先后顺序依次从切片图中提取被跟踪目标；3b)采用卡尔曼预测方法，对被跟踪目标进行预测，得到该目标在下一帧胃部切片图像的预测坐标点，如果预测坐标点在下一帧切片图某个区域的内部，则将此区域目标的坐标作为被跟踪目标在下一帧切片图中的坐标，记录到目标轨迹段中；如果预测坐标点不在下一帧切片图某个区域的内部，则返回步骤3a)。(4)判断是否全部目标预测跟踪完毕，如果全部目标预测跟踪完毕，则进行下一步骤；否则，返回步骤(3)。(5)特征匹配5a)将记录的每个目标轨迹段中的起始和终止坐标分别存于起始坐标集合和终止坐标集合中；5b)将被跟踪目标在当前胃部切片序列图像中的坐标作为中心，在下一帧胃部切片序列图像中建立一个矩形区域，将该矩阵区域作为被跟踪目标的邻域匹配区域；将当前胃部切片序列图像中轨迹终止坐标代表的目标作为待匹配目标；5c)利用相似度计算公式计算被跟踪目标与每个待匹配目标的相似度，将相似度最大的待匹配目标与被跟踪目标确定为同一个目标，再将其跟踪轨迹连接到被跟踪目标的跟踪轨迹后；5d)重复执行步骤5a)、步骤5b)、步骤5c)，直到最后一帧图像，完成对所有目标的特征匹配，将属于同一目标的跟踪轨迹，按顺序连接，形成完整的跟踪轨迹。(6)淋巴结识别6a)依次对每个目标轨迹中目标前后帧面积差求其平均值，如果平均值大于10， 则将此目标标记为候选淋巴结，若小于5，则标记为候选血管，其余判定为杂点目标；6b)对候选淋巴结依次计算其跟踪轨迹段中心的前后帧欧氏距离，计算每个目标中心欧氏距离的平均值，如果平均值小于3，则判定为淋巴，否则，判定为血管；6c)对候选血管目标查看其对应的跟踪轨迹，如果其出现的帧数大于10帧则判定为血管，否则，判定为杂点。本发明与现有的技术相比具有以下优点
1、由于本发明采用了基于目标跟踪的方法用于检测胃部切片图的淋巴结，加入计算机辅助诊断之后，较之传统的人工检测淋巴结的方法提高了检测速度及正确率；2、由于本发明首次使用计算机辅助检测识别的方法用于检测胃部切片图的淋巴结，可以很好的区分切片图中的血管和淋巴结，解决了现有技术中存在的血管和淋巴难以区分的问题。
图1为本发明的流程图；图2为未经处理的胃部切片图；图3为经本发明图像分割步骤后获得的切片图；图4为经本发明预测跟踪步骤后获得的切片图。
具体实施例方式下面结合附图对本发明做进一步的详细描述。参照图1，本发明的具体实施步骤如下步骤1、分割图像首先，采用KSVD字典学习的方法提取胃壁周围感兴趣区域第一步采用KSVD字典学习的方法对脂肪组织和非脂肪组织进行学习，生成脂肪组织Dl和非脂肪组织D2两个字典；第二步以原始图像的每个像素点为中心取5X5的邻域得到像素块Q，对像素块Q 提取15维的梯度特征值和25维的灰度特征值，组成一个40维的特征向量v，用字典D1和字典D2分别逼近特征向量v，得到逼近误差ei和％，如果ei < e2则认为该像素点为感兴趣像素点，标记类标为1，如果ei ^ e2则该像素点为不感兴趣点，标记类标为0 ；对于类标为0 的像素点，将原始图像中该像素点的灰度值赋值为0，对于类标为1的像素点，保持该像素点的灰度值为原始图像中的灰度值，得到脂肪图；第三步检测脂肪图中人体组织的最外层边缘，以脂肪图中心点处为基准，将最外层边缘向中心点缩进30个像素的半径，得到感兴趣脂肪图；第四步以感兴趣脂肪图的每个像素点为中心，将0°到360°方向均分为32个方向，如果32个方向中有超过观个方向上存在灰度值非0的像素点，则判定该像素点为感兴趣点，将该像素点灰度值赋值为原始图像中对应的灰度值，得到感兴趣区域检测结果图。然后，运用字典学习的稀疏表示方法检测每帧切片图中的目标第一步，在感兴趣区域检测结果图中分别为疑似淋巴结类和非目标类(包括脂肪区域和脂肪边界区域)选取5*5大小的图像块作为训练样本；第二步，对训练样本提取纹理和灰度特征；使用字典学习的稀疏表示(SRC)方法对每一幅图像中的图像块进行测试，标记目标类；第三步，删去一像素点个数小于10的非目标区域并输出最终分割图像。步骤2、提取被跟踪目标首先，采用四邻域连通标记算法，对所有待跟踪的切片图像中的目标区域依次进行中心标记。四邻域标记算法按照如下方法实现判断此点四邻域中的最左，最上有没有点，如果都没有点，则表示一个新的区域的开始；如果此点四邻域中的最左有点，最上没有 点，则标记此点为最左点的值，如果此点四邻域中的最左没有点，最上有点，则标记此点为 最上点的值；如果此点四邻域中的最左有点，最上都有点，则标记此点为这两个中的最小标 记点，并修改大标记为小标记；然后，提取每个目标区域中心标记，得到每个目标的中心位置，将每一帧中所有目 标的中心位置确定为被跟踪目标的中心位置。步骤3、预测跟踪首先，按照中心标记的先后顺序依次从切片图中提取被跟踪目标，然后采用卡尔 曼预测方法，对被跟踪的目标进行预测，得到该目标在下一帧切片图像的预测坐标点。本发明采用先基于卡尔曼滤波进行整体跟踪，然后采用特征匹配来连接未跟踪上 的片段，这样正确率得到了很大的提升。具体步骤如下第一歩根据被跟踪目标的运动特点，为每个目标建立勻速运动动力学模型，并对 模型中的状态向量和误差协方差进行初始化，状态向量定义为乂= (x y ~^)\量测向量定 义为Z = (x y)T，x, y分别为目标中心横坐标和纵坐标，xv, yv分别为目标在x轴方向和y 轴方向的速度，状态转移矩阵和量测矩阵分别定义为
权利要求
1.一种基于目标跟踪的计算机辅助胃癌诊断方法，其具体实现步骤如下(1)分割图像la)采用KSVD字典学习的方法提取胃壁周围感兴趣区域；lb)运用字典学习的稀疏表示方法检测每帧切片图中的目标；(2)提取被跟踪目标2a)采用四邻域连通标记算法，对所有待跟踪的胃部切片图像中的目标区域依次进行中心标记；2b)提取每个目标区域中心标记，得到每个目标的中心位置，将每一帧中所有目标的中心位置确定为被跟踪目标的中心位置(3)预测跟踪3a)按照中心标记的先后顺序依次从切片图中提取被跟踪目标；3b)采用卡尔曼预测方法，对被跟踪目标进行预测，得到该目标在下一帧胃部切片图像的预测坐标点，如果预测坐标点在下一帧切片图某个区域的内部，则将此区域目标的坐标作为被跟踪目标在下一帧切片图中的坐标，记录到目标轨迹段中；如果预测坐标点不在下一帧切片图某个区域的内部，则返回步骤3a)；(4)判断是否全部目标预测跟踪完毕，如果全部目标预测跟踪完毕，则进行下一步骤； 否则，返回步骤⑶；(5)特征匹配5a)将记录的每个目标轨迹段中的起始和终止坐标分别存于起始坐标集合和终止坐标集合中；5b)将被跟踪目标在当前胃部切片序列图像中的坐标作为中心，在下一帧胃部切片序列图像中建立一个矩形区域，将该矩阵区域作为被跟踪目标的邻域匹配区域；将当前胃部切片序列图像中轨迹终止坐标代表的目标作为待匹配目标；5c)利用相似度计算公式计算被跟踪目标与每个待匹配目标的相似度，将相似度最大的待匹配目标与被跟踪目标确定为同一个目标，再将其跟踪轨迹连接到被跟踪目标的跟踪轨迹后；5d)重复执行步骤5a)、步骤5b)、步骤5c)，直到最后一帧图像，完成对所有目标的特征匹配，将属于同一目标的跟踪轨迹，按顺序连接，形成完整的跟踪轨迹；(6)淋巴结识别6a)依次对每个目标轨迹中目标前后帧面积差求其平均值，如果平均值大于10，则将此目标标记为候选淋巴结，若小于5，则标记为候选血管，其余判定为杂点目标；6b)对候选淋巴结依次计算其跟踪轨迹段中心的前后帧欧氏距离，计算每个目标中心欧氏距离的平均值，如果平均值小于3，则判定为淋巴，否则，判定为血管；6c)对候选血管目标查看其对应的跟踪轨迹，如果其出现的帧数大于10帧则判定为血管，否则，判定为杂点。
2.根据权利要求1所述的基于目标跟踪的计算机辅助胃癌诊断方法，其特征在于步骤la)所诉的KSVD字典学习的方法是按照以下步骤求感兴趣区域第一步采用KSVD字典学习的方法对脂肪组织和非脂肪组织进行学习，生成脂肪组织 Dl和非脂肪组织D2两个字典；第二步以原始图像的每个像素点为中心取5X5的邻域得到像素块Q，对像素块Q提取15维的梯度特征值和25维的灰度特征值，组成一个40维的特征向量v，用字典D1和字典D2分别逼近特征向量v，得到逼近误差ei和％，如果ei < e2则认为该像素点为感兴趣像素点，标记类标为1，如果ei ^ e2则该像素点为不感兴趣点，标记类标为0 ；对于类标为0的像素点，将原始图像中该像素点的灰度值赋值为0，对于类标为1的像素点，保持该像素点的灰度值为原始图像中的灰度值，得到脂肪图；第三步检测脂肪图中人体组织的最外层边缘，以脂肪图中心点处为基准，将最外层边缘向中心点缩进30个像素的半径，得到感兴趣脂肪图；第四步以感兴趣脂肪图的每个像素点为中心，将0°到360°方向均分为32个方向， 如果32个方向中有超过观个方向上存在灰度值非0的像素点，则判定该像素点为感兴趣点，将该像素点灰度值赋值为原始图像中对应的灰度值，得到感兴趣区域检测结果图。
3.根据权利要求1所述的基于目标跟踪的计算机辅助胃癌诊断方法，其特征在于步骤lb)所述字典学习的稀疏表示方法是按照下面步骤进行分割第一步，在感兴趣区域检测结果图中分别为疑似淋巴结类和非目标类(包括脂肪区域和脂肪边界区域)选取5*5大小的图像块作为训练样本；第二步，对训练样本提取纹理和灰度特征；使用字典学习的稀疏表示(SRC)方法对每一幅图像中的图像块进行测试，标记目标类；第三步，删去一像素点个数小于10的非目标区域并输出最终分割图像。
4.根据权利要求1所述的基于目标跟踪的计算机辅助胃癌诊断方法，其特征在于步骤3b)所述卡尔曼预测方法的具体步骤如下步骤1，根据被跟踪目标的运动特点，为每个目标建立勻速运动动力学模型； 步骤2，建立目标运动模型； Xi (k+1) = FiXi (^ + TiWi (k)其中，Xi(k+1)是第i个目标在下一帧中的状态向量，k是帧数,Fi是第i个目标的状态转移矩阵，Xi (k)是第i个目标在当前帧中的状态向量，Γ i是第i个目标的噪声矩阵，WiGO 是与第i个目标关联的系统噪声； 步骤3，建立量测向量模型 ZiGO = HiXiGO+ViGO其中，Zi (k)是第i个目标的量测向量，Hi是第i个目标的量测矩阵，Xi (k)是第i个目标的状态向量，ViGO是与第i个目标关联的量测噪声，k是帧数，； 步骤4，对状态向量及其误差协方差矩阵进行初始化； 步骤5，按照步骤3中的量测向量模型求预测值。
5.根据权利要求1所述的基于目标跟踪的计算机辅助胃癌诊断方法，其特征在于所述步骤5c)计算被跟踪目标与每个待匹配目标的相似度按照下式计算求得sim(L j)=---aD(iJ)+PA(iJ)+YE(iJ)其中，sim(i，j)表示相似度，i表示被跟踪目标，j表示待匹配目标，α表示距离特征， D(i,j)表示待匹配目标与被跟踪目标的距离相似度，β表示面积特征，A(i，j)表示待匹配目标与被跟踪目标的面积相似度，Y表示偏心率特征的权值，E(i，j)表示待匹配目标与被跟踪目标的偏心率相似度。
全文摘要
本发明公开了一种基于目标跟踪的计算机辅助胃癌诊断方法，主要解决医学影像领域中胃癌诊断过程中存在的淋巴结转移判定问题问题。其实现步骤为(1)分割图像；(2)提取被跟踪目标；(3)预测跟踪；(4)判断是否全部目标预测跟踪完毕，如果全部目标预测跟踪完毕，则进行下一步骤；否则，返回步骤(3)；(5)特征匹配；(6)淋巴结识别。本发明是将计算机辅助医学诊断的方法用于胃癌淋巴结转移的检测，可以在很短的时间内处理大量胃部切片图，较之传统的胃癌淋巴结转移的检测提高了检测速率和正确率，实现胃部切片图淋巴结的自动识别，极大的缩短临床诊断时间，为患者争取更好的治疗时机。
文档编号G06F19/00GK102436551SQ20111035608
公开日2012年5月2日 申请日期2011年11月10日 优先权日2011年11月10日
发明者公茂果, 刘芳, 周治国, 焦李成, 王爽, 高婷婷 申请人:西安电子科技大学