宫颈成像中的组织分类方法
专利名称：宫颈成像中的组织分类方法
技术领域：
本发明主要涉及医学成像，更具体地涉及利用宫颈组织识别和分类方法对子宫颈癌和癌前病变的计算机辅助检测和诊断。
背景技术：
正被公开的本发明与宫颈癌有关。子宫颈癌是全球妇女第二大最常见的癌症，每年有近 50 万新病例及超过 27 万的死亡(http://www- depdb. iarc. fr/globocan2002. htm.通过参考结合于此)。阴道镜检查是一种用于检测子宫颈癌和癌症前体的诊断方法(B. S. Apgar, Brotzman, G. L. and Spitzer, Μ. , Colposcopy: Principles and Practice, W. B. Saunders Company: Philadelphia, 2002,通过参考结合于此)。用于阴道镜检查的CAD (计算机辅助诊断)代表了医学图像处理的新应用。发明家已经开发出一种 CAD系统，其模仿或模拟阴道镜专家所用来评估异常的严重程度的诊断过程(Lange H. and Ferris, Daron G. ； Computer- Aided-Diagnosis (CAD) for colposcopy； SPIE Medical Imaging 2005; SPIE Proc. 5747，2005,通过参考结合于此)。阴道镜专家基于宫颈组织中不同的宏观上皮(与膜相关)特征形成阴道镜检查印象。阴道镜评级系统提供了结构化、系统化的方法，以精密的分析宫颈检查结果。一种系统化方法一旦被学会和掌握，就变为惯例，并且是非常有益的。阴道镜评级系统还容许阴道镜专家形成更准确的阴道镜检查印象。精心设计的阴道镜打分系统增强了阴道镜检查的重现性。当试图选择最合适的活检部位，尤其是遇到宫颈中大的、复杂病变(异常生长)时，阴道镜评级系统也是有用的。打分方案，像里德(Reid)阴道镜指数，是基于不同的特征作出阴道镜诊断的辅助方法(Reid Rj Scalzi P. Genital warts and cervical cancer. VII. An improved colposcopic index for differentiating benign papillomaviral infection from high-grade cervical intraepithelial neoplasia. Am J Obstet Gynecol 1985;153:611-618，通过参考结合于此；Reidj R.，Stanhope, C. R.，Herschmanj B. R.，Crumj C. P.，and Agronowj S. J.，Genital warts and cervical cancer. IV. A colposcopic index for differentiating subclinical papillomaviral infection from cervical intraepithelial neoplasia, Am. J. Obstet. Gynecol. 149 (8): 815-823. 1984，通过参考结合于此)，包括病变(异常生长)的边缘或边界，采用5%醋酸溶液后的病变颜色，病变中的血管特点，病变对采用卢戈氏(Lugol)碘液的反应。这些特点分别被评估并打分，然后把所有特征的分数相结合以得到对疾病严重程度分级的综合分数。Reid指数区分了低度宫颈疾病与高度宫颈疾病。因此，Reid阴道镜指数(RCI)并非设计为区分恶变前的与恶性的宫颈肿瘤。尽管如此，该指数提供了一种受欢迎的规范评价宫颈病变的方式。Rubin 禾口 Barbo (Rubin, Μ. Μ. and Barbo, D. Μ. , Ch. 9a: Rubin and Barbo Colposcopic Assessment System, in Colposcopy: Principles and practice, eds. Apgar, B. S. , Brotzman, G. L. , and Spitzer, Μ. , pp. 187-195. W. B. SaundersCompany, Philadelphia, 2002,通过参考结合于此)开发了一种评估方法，其保留了一些以前的阴道镜评级系统的最好的描述符，但去掉了数字，这可能很让人费解。此外，它将系统扩展为包括正常结果的描述符。更重要的是，它包括将临床医生的模式识别过程聚焦于微创的，或者坦白地说是侵入性疾病的可能性的描述符。例如，它不仅测量醋酸白色上皮变化的强度，还涉及其他色调的变化，例如红色，黄色，暗灰色，这些与侵入性癌症的存在更加相关。其他值得考虑的因素包括患者的年龄以及病变分布及大小。宫颈癌患者中很少见到年龄小于25岁的。高度鳞状上皮内疾病(特别地，CIN3)主要发现于年龄在四岁到32 岁之间的妇女。Cim病变倾向于相对小。Cim病变的平均长度大约是2. 8mm (毫米或千分之一米)。尽管如此，这些病变可达到11. 5mm的最大长度。相比之下，CIN2和CIN3病变更大。CIN2和CIN3的平均(平均的)长度分别为5. 8mm和7. 6mm。它们的最大长度分别为 18. 2mm和20. 6mm。低度病变可能只占据四分之一的宫颈，或占据外宫颈的表面区域的很小比例。cmi的分布随单病灶(病变发生在单一处)到多病灶(病变发生在多处)而不同。多个，不同的，小的，随机分布的病变是Cim的特点。CIN2总是出现在转移区(TZ，transformation zone)。TZ是宫颈的已被新化生的鳞状上皮取代的柱状上皮区域。相比于Cim(其可能被发现于TZ外)，CIN 2可能是多病灶， 但是单病灶病变更为常见。卫星状病变通常不是CIN2的代表。阴道镜专家通常沿鳞-柱交界处(SCJ)发现CIN2，要么位于外宫颈上要么位于宫颈内道内。CIN3病变倾向于是融合性的(汇至一起或混为一个)，比Cim或CIN2病变更长更宽。CIN3通常位于宫颈内的中央部分，在内弯道中朝向外口(os)。当出现在外宫颈时 CIN3很少是隐蔽的(眼睛看不到)。CIN3病变的线性长度，被定义为在尾部(位于或接近身体的后端)和头部(朝向身体的头部或前端)边缘之间的组织表面的距离，从2mm到22mm不等。平均线性长度范围为6mm到10mm。长的线性病变-那些长于IOmm的，特别是包含有宫颈内时，一直被怀疑为癌症。由于病变的表面区域增加至大于40平方毫米，所以应被怀疑为癌症。据报道，转换区TZ大小，病变(异常生长)大小，明显边缘，血管模式(血管模式)， 以及醋酸白颜色(应用醋酸后的颜色)与组织评级显著相关，并显示宫颈病变的大小可能具有丨临床重要性° (Kierkegaard, 0. , Byralsen, C, Hansen, K. C. , Frandsen, K. H., and Frydenberg, Μ. , Association between colposcopic findings and histology in cervical lesions: the significance of the size of the lesion, Gynecol. Oncol. 57(1) : 66-71. 1995，通过参考结合于此)。循证医学(EBM)被定义为“在对患者个体治疗作出决定时，对目前最好证据的认真、详细、明确的应用”。在EBM中，人类的感知和解释被认为对诊断说明非常有价值。然而， 人类的评估可能是主观的且具有低的观察者间的一致性(Ferris，D. G. and Litaker, Μ. , Interobserver agreement for colposcopy quality control using digitized colposcopic images during the ALTS trial, J. Low. Genit. Tract Dis. 9(1): 29-35，2005,通过参考结合于此；Jeronimo, J.，Massad, L. S.，Castle, P. Ε.， Wacholder, S. , and Schiffman, Μ. , Interobserver Agreement in the Evaluation of Digitized Cervical Images, Obstet. Gynecol. 110(4) : 833-840，2007,通过参考结合于此)。
计算机可被编程以执行推理任务，并且可用于解决诊断问题。阴道镜评级系统被阴道镜专家们用于获得临床诊断已超过40年。阴道镜专家们通过使用阴道镜的视觉检查定性地获得患者的阴道镜检查特征，确定疾病程度并基于他们的经验开始对患者诊治。这个过程在妇科医生的办公室需要花费10或15分钟。由于检查的主观性，阴道镜检查的准确度高度依赖于阴道镜专家的经验和专业知识。使用了量化的阴道镜特征的计算机实现的阴道镜打分算法是可取的以给医生提供建议的诊断决定。大部分已有的组织分类系统依赖于统计方法。一个例子就是Georgia Tech Vision (GTV)系统(Dickman, Ε. D. , Doll, Τ. J. , Chiu, C. K. , and Ferris, D. G., Identification of Cervical Neoplasia Using a Simulation of Human Vision, Journal of Lower Genital Tract Disease 5(3) : 144-152，2001，通过参考结合于此)，一种由人类视觉系统给予灵感的通用的计算机视觉系统识别正常和异常的宫颈特征。另一个例子是用于宫颈肿瘤的无创诊断的多模型高光谱成像(MHI)系统(Ferris， D. G.，Lawheadj R. A.，Dickman, E. D.，Holtzapplej N.，Miller, J. Α.，Groganj S. et al.，Multimodal hyperspectral imaging for the noninvasive diagnosis of cervical neoplasia, J. Low. Genit. Tract Dis. 5(2): 65-72，2001，通过参考结合于此)。最近，已经报道了用来区分高度上皮内病变与正常或低度上皮内病变的反射荧光系统。(Park, S. Y. , Fo 11 en, Μ. , Milbourne, Α. , Rhodes, H. , Malpica, Α., Mackinnon, N. et al. , Automated image analysis of digital colposcopy for the detection of cervical neoplasia, J. Biomed. Opt. 13(1) : 014029-1-014029-10，2008，通过参考结合于此)。此外，也已经报道了用于高度宫颈上皮内肿瘤的活体内(在有生命的体内)监测的成功的光学系统 LUMA (Huh, W. K. , Cestero, R. Μ. , Garcia, F. Α. , Gold, Μ. Α. , Guido, R. S. , Mclntyre—Seltman, K. et al. , Optical detection of high-grade cervical intraepithelial neoplasia in vivo: results of a 604- patient study, Am. J. Obstet. Gynecol. 190(5) : 1249-1257. 2004,通过参考结合于此；Kendrick, J. Ε. , Huh, W. K. , and Alvarez, R. D. , LUMA Cervical Imaging System, Expert Review of Medical Devices 4(2) : 121-129，2007，通过参考结合于此)。在 LUMA 系统中，在醋酸应用之后使用了 3种集成的光学测量执行整个人体宫颈的诊断扫描激光诱导荧光光谱，白光漫反射光谱，以及视频成像。进一步地，包括多变量分类算法，以利用特定光源与不同的宫颈组织类型相互作用中的光谱差异。最近报道了另一个称为DySIS (动态光谱成像系统)的活体内检测和量化的评级系统(Forthphotonics, DySISiDynamic Spectral Imaging System, http://www. forth-photonics. gr/dysis. php page=dysis_trials2008，通过参考结合于此；Soutterj W. P.， Diakomanolisj Ε.，Lyons, D.，and Haidopoulosj D.，Dynamic Spectral Imaging -In Vivo Detection and Quantitative Grading of Cervical Neoplasia, ASCCP 2008 Meeting Abstract. Journal of Lower Genital Tract Disease 12(2): 160， 2008，通过参考结合于此)。DySIS皿是第一个CE标记的(CE标记是在欧盟地区(EEA)的单一市场中投放的许多产品上的强制认证标记一CE标记证明产品符合欧盟的卫生，安全和环境要求， 这保证了消费者的安全)设备，以使用数字和动态光谱成像技术，用于帮助活体内自动(和用户无关)检测以及绘制宫颈瘤样病变(病变或肿瘤)。
本发明是一种基于阴道镜评级系统以准确评估及识别低度病变，高度病变和癌症的计算机辅助检测技术。它不同于现有技术中列出的统计学习方法。本发明是一种遵循阴道镜检查的基于规则的方法，所述方法使用上皮细胞和血管的形态学评价积累了疾病的证据。在统计学习中，诊断判断是基于经交叉验证从大量训练数据导出的统计特征。因此， 它是基于训练数据和统计理论的监督学习方法。目前的发明将现代阴道镜中的定性评级规则转移到定量计算机程序中。它不需要来自数据的任何训练。它是基于已有的医学知识的无监督方法。图1中示出了示出本发明概述的流程图。在应用本发明的组织分类算法之前，首先采用其他方法来提取信息，该信息用作组织分类算法的输入。其他方法包括解剖特征算法，醋酸白特征提取算法，镶嵌区域(mosaic)和点状区域(pimctation)检测算法(例如，W. Li and A. Poirson, Detection and characterization of abnormal vascular patterns in automated cervical image analysis, Lecture Notes in Computer Science - Advances in Visual Computing 4292, Second International Symposium, ISVC 2006, Lake Tahoe, NV, November 2006 Proceedings, Part II, 627-636 (Springer 2006],通过参考结合于此)以及非典型血管提取算法(如在共同
I^gih白勺、11 "Methods for Detection and Characterization of Atypical Vessels in Cervical Imagery"，2008年8月1日提交，序列号12/2213 的专利申请公开的，其通过参考结合于此)。这些方法的输出作为本文所描述的组织分类方法的输入。组织分类方法应用于在临床检查期间，应用醋酸之前和之后实时提取的宫颈的图像数据，以及对象的人口统计学数据，例如年龄。组织分类方法的输出包括具有以下检查结果的诊断判断病变类型(正常/NED (无疾病证据)，低度增生，高度增生，或癌症)；疾病位置或高度 (CIN2和CIN3)病变位置；以及判断的置信度。下面的专利和专利申请可认为与本发明的领域相关
Costa等的美国专利No. 7309867，其通过参考结合于此，公开了通过利用统计和非统计分类技术的结合，以及通过结合光谱数据和图像数据，确定组织样本的给定区域包含给定类型的组织，例如CIN，CIN 11/111，正常鳞状上皮，正常柱状以及组织转化的概率的方法。Zelenchuk的美国专利No. 731(^47，其通过参考结合于此，公开了一种系统和方法，用于在原位区分健康和发病组织。它包括直接将紫外线照射至组织标本，并收集荧光响应辐射的光纤探头。响应辐射在三个选定的波长被观察，其中一个对应于等吸光点。在一个例子中，等吸光点发生在大约431纳米。使用431nm强度将观察到的信号强度归一化。使用判别分析中的比例确定得分。根据判别分析的结果，基于疾病或健康的诊断确定被检查组织是否切除。Kaufman等的美国专利No. 7260M8，其通过参考结合于此，公开了与基于测量相似性的多个图像相关的方法。这些方法在分割组织的宫颈图像序列中有用。这些方法可以用于确定宫颈组织的醋酸变白测试中的组织特性。它公开了确定组织特性的方法，包括获得多个组织图像；确定这些图像中的两个或更多区域之间的关系；基于所述关系分割图像；并基于所述分割确定所述组织的特性。确定步骤包括将组织表征为正常，CIN I, CIN II，CIN III 或 CIN 11/111。Utzinger等的美国专利No. 6766184，其通过参考结合于此，公开了生成组织的多光谱图像的方法和装置。所述图像可用作宫颈癌检测。利用辐射源产生主要辐射，并过滤以选择第一波长和第一偏振。经过滤的主要辐射照射组织，生成次级辐射，其经过滤以选择第二波长和第二偏振。经过滤的次级辐射用探测器收集，并且根据第一、第二波长和偏振的不同组合，生成组织的多个多光谱图像。Roehrig等的美国专利No. 6198838，其通过参考结合于此，公开了通过利用独立计算的块(mass)和毛刺信息来检测数字乳房X线照片中的可疑部分的方法和系统。所述方法在计算机辅助诊断系统中使用，该系统设计为使纤维乳腺组织中的可疑或可能为癌的病变引起放射学家或其他医疗专业人士的关注。在优选实施例中，毛刺信息和块信息被独立计算，由于算出的毛刺信息不依赖于块信息计算的结果，从而具有更高的可靠性。Turner等的美国专利No. 6135965，其通过参考结合于此，公开了一种利用神经网络分析活体内荧光光谱测量来光谱检测宫颈组织中组织异常的装置和方法。Gaborski等的美国专利No. 5857030，其通过参考结合于此，公开了一种对放射图像进行数字图像处理的自动化的方法和系统，包括过滤的预处理阶段，对分割预选阶段， 以及包括神经网络分类的模式分类阶段。Clarke等的美国专利No. 5982917，其通过参考结合于此，公开了一种计算机辅助诊断(CAD)方法和装置，用于对数字X射线图像中的可疑区域检测和增强。Balas的美国专利公开No. 2006/0141633，其通过参考结合于此，公开了一种用于活体内检测并绘制(map)在癌症发展期间上皮组织的生化和/或功能特性引起的变化的方法和装置。Keller等的美国专利申请No. 2006/0184040，其通过参考结合于此，公开了检测组织异常的方法和设备，据此，所述方法包括从光源向组织照射光；将从组织反射的照射光通过光学器件引至多波长成像光学子系统；并分离感兴趣的一个或多个波长或波带；将感兴趣的一个或多个波长或波带引至一个或多个成像设备，并使用所述设备记录感兴趣的一个或多个波长或波带的图像；将图像数据从图像传输至计算系统；并对组织异常相关的一个或多个光谱模式分析所述图像。Gu等的美国专利公开No. 2008/0058593，其通过参考结合于此，公开了一种在用内窥镜检查期间从器官的视频数据提供计算机辅助诊断的过程，包括分析和增强来自视频的图像帧，并在检查期间实时检测和诊断图像帧中的任何病变。进一步地，图像数据可被用于生成器官的三维重建。
发明内容
这里以及下文更充分描述的本发明包括一种用于将宫颈组织分类的基于规则的无监督过程，其包括对宫颈组织应用包含一组分类器规则的分类器，通过串行地应用选定的分类器规则来获得最终预测，其中，采用选自由确定以下组成的组的分类器规则执行应用步骤纹理区域的大小、不透明度参数、醋酸白区域大小、粗大和细小点状区域的数目、粗大和细小镶嵌区域的数目、非典型血管大小以及人口统计学数据，由此，可将宫颈组织分类为无疾病证据、低度异型增生、高度异型增生或癌症。如果纹理区域大于宫颈面积的约1%，以及如果不透明度参数小于像素值的约4%， 或大于像素值的约12%，则将获得无疾病证据的最终预测。
如果纹理区域大于宫颈面积的约1%，所述不透明度参数介于像素值的约4%至像素值的约1 之间，醋酸白的区域大小大于宫颈面积的7%，以及如果口附近的醋酸白区域大小大于宫颈面积的0%，则将获得高度异型增生的最终预测。如果纹理区域大于宫颈面积的约1%，如果所述不透明度参数介于像素值的约4% 至像素值的约1 之间，醋酸白的区域大小小于宫颈面积的7%，并且口附近没有醋酸白区域，则将获得低度异型增生的最终预测。如果纹理区域小于宫颈面积的约1%，并且如果不透明度参数小于像素值的5%，粗大点状区域的数目小于2，且没有粗大的镶嵌区域血管，则将获得低度异型增生的最终预测。如果纹理区域小于宫颈面积的约1%，并且如果不透明度参数大于像素值的10%， 口附近的醋酸白区域大小大于15%，则将获得癌症的预测。如果纹理区域小于宫颈面积的约1%，并且不透明度参数大于8%，口附近的醋酸白区域的大小大于0. 5%，并且粗大点状区域的数目大于100，并且没有粗大的镶嵌区域血管， 并且异常血管大小大于0. 5毫米，则将得到癌症的预测。如果纹理区域小于宫颈面积的1%，并且如果不透明度参数大于像素值的5%，口附近的醋酸白区域的大小小于像素值的0. 5%，粗大点状区域血管的数目小于100，不存在粗大镶嵌区域血管，并且非典型血管的大小小于0. 5%，则将得到高度异型增生的判断。
图1展示了本发明的概述的流程图。
图2是用于组织分类方法的输入参数的列表。
图3是组织诊断方案的第一部分。
图4是组织诊断方案的第二部分。
图5是NED (没有疾病证据)的对象的示例的照片。
图6 (a)是低度疾病对象的示例的照片。
图6 (b)是图6 (a)的对象的对应低度疾病位置的图像处理过的照片。
图7 (a)是癌症对象的示例的照片。
图7 (b)是图7 (a)的对象的对应高度疾病位置的图像处理过的照片。
图7 (c)是图7 (a)中的对象的对应低度疾病位置的图像处理过的照片。
图8 (a)是高度疾病对象的示例的照片。
图8 (b)是图8 (a)中的对象的对应的高度疾病位置的图像处理过的照片
图8 (c)是图8 (a)中对象的对应的低度疾病位置的图像处理过的照片。
具体实施例方式
1.预方法应用
本文所述的本发明目前优选的实施例优选使用几个预方法应用的输出(输出为如图2 中所描述的某些参数)作为组织分类方法的输入。很多参数或变量被描述为图像中特定特征的面积相比于宫颈的总面积的百分比。所述的点状区域特征被定义为图像中存在的点状区域的数目。非典型性血管的大小优选的以毫米(mm)来计算。
该过程首先收集宫颈的数字彩色图像(使用红色，绿色，蓝色通道或RGB)，——一个在施加醋酸之前，一个在应用醋酸之后。优选地，图像被正交极化，以抑制镜面反射(闪烁)并被校准。接着，采用几种预方法应用，包括解剖特征提取，醋酸白特征提取，镶嵌区域和点状区域检测，以及非典型血管提取。解剖特征提取优选地包括检测子宫颈，除去阴道壁(水平和竖直)，口检测，柱状上皮检测以及纹理区域(异常/化生)检测。醋酸白特征提取优选地包括应用醋酸到感兴趣的解剖区域，以产生醋酸白区域； 拍摄至少一张醋酸白区域的后醋酸图像；检测感兴趣的具有上皮细胞的解剖区域；从后醋酸图像提取出纹理信息，以定义纹理区域；从后醋酸图像中提取色彩信息以定义色彩区域； 并将纹理区域和色彩区域相结合生成候选的醋酸白区域。优选地，它还包括获得至少一个上皮细胞的前醋酸图像；在前醋酸图像之一中和在后醋酸图像之一中记录醋酸白区域， 以形成醋酸白区域的记录图像；从醋酸白区域的记录图像中提取不透明度参数，如共同待决共同受让的2008年8月4日提交的序列号为12/221645的题为“Computerized Image Analysis for Acetic Acid Induced Cervical Intraepithelial Neoplasia，，的专利申请，其通过参考结合于此。该申请公开了通过聚类在记录的前醋酸图像和后图像的时间变化中的差异，来创建的不透明度参数。所述不透明度参数对病变诊断提供了高预测值—— 正常和低度病变具有比高度病变和癌症情形低得多的不透明度。其他的不透明度参数提取方法也可被使用，包括在不同色彩空间中的比率成像。在本发明的一优选实施例中，在CIE-Lab色彩空间中的a通道中的色彩差异被用于不透明度参数提取。在本发明另一优选实施例中，RGB色彩空间中的g通道的色彩差异被用于不透明度参数提取。由于不成熟的化生和柱状上皮组织在醋酸施加后短暂变为白色， 但并未表现出发育异常的组织变化，他们不应被排除在感兴趣的醋酸白区域之外。这些组织区域通常表现出轻微的不透明度变化。因此本发明在色彩差异特征空间内采用了两步平均移位聚类算法。第一步是分割主要的不透明度变化，并消除轻微不透明度变化。第二步是从第一步获得的前景区域中分割最不透明的变化。计算不透明度参数，作为最不透明区域的平均色差。最不透明区域被定义为具有最大平均色差的区域。不透明度参数通过以下公式来表示
权利要求
1.一种用于将宫颈组织分类的基于规则的无监督方法，包括对所述宫颈组织应用分类器，所述分类器包括一组分类器规则；并且通过串行地应用在所述组中选定的分类器规则来获取最终预测，其中使用选自由确定以下组成的组的分类器规则来执行所述应用步骤纹理区域的大小， 不透明度参数，醋酸白区域的大小，粗大和细小点状区域的数目，粗大和细小镶嵌区域的大小，非典型血管的大小以及人口统计学数据；由此，所述宫颈组织可被分类为无疾病证据、低度异型增生、高度异型增生或癌症。
2.根据权利要求1所述的方法，其中如果所述纹理区域大于所述宫颈面积的约1%，并且所述不透明度参数小于像素值的约 4%或大于像素值的约12%，获取无疾病证据的最终预测。
3.根据权利要求1所述方法，其中如果所述纹理区域大于所述宫颈面积的约1%，所述不透明度参数介于像素值的约4% 和像素值的约1 之间，醋酸白区域的大小大于所述宫颈面积的7%，并且口附近的醋酸白区域的大小大于所述宫颈面积的0%，则获取高度异型增生的最终预测。
4.根据权利要求1所述方法，其中，如果所述纹理区域大于所述宫颈面积的约1%，所述不透明度参数介于像素值的约4% 和约12%之间，醋酸白区域的大小小于宫颈面积的7%，并且在口附近没有醋酸白区域，则获取低度异型增生的最终预测。
5.根据权利要求1所述方法，其中，如果所述纹理区域小于所述宫颈区域的约1%，所述不透明度参数小于像素值的5%，粗大点状区域的数目小于2，并且没有粗大镶嵌区域血管，则获取低度异型增生的最终预测。
6.根据权利要求1所述方法，其中如果所述纹理区域小于所述宫颈面积的约1%，所述不透明度参数大于像素值的10%， 并且口附近的所述醋酸白区域的大小大于所述宫颈面积的15%，则获取癌症的最终预测。
7.根据权利要求1所述方法，其中如果所述纹理区域小于所述宫颈面积的约1%，所述不透明度参数大于像素值的8%，口附近的所述醋酸白区域的大小大于宫颈面积的0. 5%，粗大点状区域的数目大于100，没有粗大镶嵌区域血管，并且异常血管的大小大于0. 5mm，则获取癌症的最终预测。
8.根据权利要求1所述方法，其中，如果所述纹理区域小于所述宫颈面积的约1%，所述不透明度参数大于像素值的5%，口附近的醋酸白区域的大小小于所述宫颈面积的0. 5%，粗大点状区域血管的数目小于100， 没有粗大镶嵌区域血管，并且非典型血管的大小为小于0. 5mm,则获取高度异型增生的最终预测。
9.根据权利要求1-8任一项所述方法，进一步包括如果患者年龄低于25岁或高于 60岁，则将最终的癌症预测中的置信度降低约10%-约20%之间。
10.一种用于将宫颈组织分类的基于规则的无监督方法，包括对所述宫颈组织应用分类器，所述分类器包括一组分类器规则；并且通过串行地应用在所述系列中选定的分类器规则来获取最终预测，其中使用选自由确定以下组成的组的分类器规则来执行所述应用步骤纹理区域的大小，不透明度参数，醋酸白区域的大小，粗大和细小点状区域的数目，粗大和细小镶嵌区域的大小，非典型血管的大小以及人口统计学数据；由此，所述宫颈组织可被分为无疾病证据、低度异型增生、高度异型增生或癌症； 其中如果所述纹理区域大于所述宫颈面积的约1%，并且所述不透明度参数小于像素值的约4%或大于像素值的约12%，则获取无疾病证据的最终预测；如果所述纹理区域大于所述宫颈面积的约1%，所述不透明度参数介于像素值的约4% 和像素值的约1 之间，醋酸白区域的大小大于所述宫颈面积的7%，并且如果口附近的醋酸白区域的大小大于所述宫颈面积的0%，则获取高度异型增生的最终预测；其中如果所述纹理区域小于所述宫颈面积的约1%，所述不透明度参数大于像素值的 5%，口附近的所述醋酸白区域的大小小于宫颈面积的0. 5%，粗大点状区域血管的数目小于 100，没有粗大镶嵌区域血管，并且异常血管的大小小于0. 5mm,则获取高度异型增生的最终预测；其中如果所述纹理区域大于所述宫颈区域的约1%，所述不透明度参数介于像素值的约 4%和像素值的约1 之间，醋酸白区域的大小小于宫颈面积的7%，并且在口附近没有醋酸白区域，则获取低度异型增生的最终预测；其中如果所述纹理区域小于所述宫颈区域的约1%，所述不透明度参数小于像素值的 5%，粗大点状区域的数目小于2，并且没有粗大镶嵌区域血管，则获取低度异型增生的最终预测；其中如果所述的纹理区域小于所述宫颈面积的约1%，所述不透明度参数大于像素值的 10%，口附近的所述醋酸白区域的大小大于所述宫颈面积的15%，则获取癌症的预测；其中如果所述的纹理区域小于所述宫颈面积的约1%，所述不透明度参数大于像素值的 8%，口附近的所述醋酸白区域的大小大于宫颈面积的0. 5%，粗大点状区域的数目大于100， 没有粗大镶嵌区域血管，并且异常血管的大小大于0. 5mm，则获取癌症的最终预测。
11.根据权利要求10所述方法，进一步包括如果患者年龄低于25岁或高于60岁，则将最终的癌症预测中的置信度降低约10%-约 20%之间。
12.一种用于将宫颈组织分类的基于规则的无监督方法，包括 应用分类器，所述分类器包括确定不透明度参数的分类器规则；由此，所述宫颈组织可被分为无疾病证据、低度异型增生、高度异型增生或癌症。
全文摘要
一种用于将宫颈组织分类的基于规则的无监督方法，所述分类通过串行地应用选自由确定以下组成的组的分类器纹理区域的大小、不透明度参数、醋酸白区域的大小、粗大和细小点状区域的数目、粗大和细小镶嵌区域的大小、非典型血管的大小以及人口统计学数据，使得所述宫颈组织可被分类为无疾病证据、低度异型增生、高度异型增生或癌症。
文档编号G06K9/00GK102282569SQ200980139997
公开日2011年12月14日 申请日期2009年10月9日 优先权日2008年10月10日
发明者李 W. 申请人:国际科学技术医疗系统有限责任公司
技术领域：
本发明主要涉及医学成像，更具体地涉及利用宫颈组织识别和分类方法对子宫颈癌和癌前病变的计算机辅助检测和诊断。
背景技术：
正被公开的本发明与宫颈癌有关。子宫颈癌是全球妇女第二大最常见的癌症，每年有近 50 万新病例及超过 27 万的死亡(http://www- depdb. iarc. fr/globocan2002. htm.通过参考结合于此)。阴道镜检查是一种用于检测子宫颈癌和癌症前体的诊断方法(B. S. Apgar, Brotzman, G. L. and Spitzer, Μ. , Colposcopy: Principles and Practice, W. B. Saunders Company: Philadelphia, 2002,通过参考结合于此)。用于阴道镜检查的CAD (计算机辅助诊断)代表了医学图像处理的新应用。发明家已经开发出一种 CAD系统，其模仿或模拟阴道镜专家所用来评估异常的严重程度的诊断过程(Lange H. and Ferris, Daron G. ； Computer- Aided-Diagnosis (CAD) for colposcopy； SPIE Medical Imaging 2005; SPIE Proc. 5747，2005,通过参考结合于此)。阴道镜专家基于宫颈组织中不同的宏观上皮(与膜相关)特征形成阴道镜检查印象。阴道镜评级系统提供了结构化、系统化的方法，以精密的分析宫颈检查结果。一种系统化方法一旦被学会和掌握，就变为惯例，并且是非常有益的。阴道镜评级系统还容许阴道镜专家形成更准确的阴道镜检查印象。精心设计的阴道镜打分系统增强了阴道镜检查的重现性。当试图选择最合适的活检部位，尤其是遇到宫颈中大的、复杂病变(异常生长)时，阴道镜评级系统也是有用的。打分方案，像里德(Reid)阴道镜指数，是基于不同的特征作出阴道镜诊断的辅助方法(Reid Rj Scalzi P. Genital warts and cervical cancer. VII. An improved colposcopic index for differentiating benign papillomaviral infection from high-grade cervical intraepithelial neoplasia. Am J Obstet Gynecol 1985;153:611-618，通过参考结合于此；Reidj R.，Stanhope, C. R.，Herschmanj B. R.，Crumj C. P.，and Agronowj S. J.，Genital warts and cervical cancer. IV. A colposcopic index for differentiating subclinical papillomaviral infection from cervical intraepithelial neoplasia, Am. J. Obstet. Gynecol. 149 (8): 815-823. 1984，通过参考结合于此)，包括病变(异常生长)的边缘或边界，采用5%醋酸溶液后的病变颜色，病变中的血管特点，病变对采用卢戈氏(Lugol)碘液的反应。这些特点分别被评估并打分，然后把所有特征的分数相结合以得到对疾病严重程度分级的综合分数。Reid指数区分了低度宫颈疾病与高度宫颈疾病。因此，Reid阴道镜指数(RCI)并非设计为区分恶变前的与恶性的宫颈肿瘤。尽管如此，该指数提供了一种受欢迎的规范评价宫颈病变的方式。Rubin 禾口 Barbo (Rubin, Μ. Μ. and Barbo, D. Μ. , Ch. 9a: Rubin and Barbo Colposcopic Assessment System, in Colposcopy: Principles and practice, eds. Apgar, B. S. , Brotzman, G. L. , and Spitzer, Μ. , pp. 187-195. W. B. SaundersCompany, Philadelphia, 2002,通过参考结合于此)开发了一种评估方法，其保留了一些以前的阴道镜评级系统的最好的描述符，但去掉了数字，这可能很让人费解。此外，它将系统扩展为包括正常结果的描述符。更重要的是，它包括将临床医生的模式识别过程聚焦于微创的，或者坦白地说是侵入性疾病的可能性的描述符。例如，它不仅测量醋酸白色上皮变化的强度，还涉及其他色调的变化，例如红色，黄色，暗灰色，这些与侵入性癌症的存在更加相关。其他值得考虑的因素包括患者的年龄以及病变分布及大小。宫颈癌患者中很少见到年龄小于25岁的。高度鳞状上皮内疾病(特别地，CIN3)主要发现于年龄在四岁到32 岁之间的妇女。Cim病变倾向于相对小。Cim病变的平均长度大约是2. 8mm (毫米或千分之一米)。尽管如此，这些病变可达到11. 5mm的最大长度。相比之下，CIN2和CIN3病变更大。CIN2和CIN3的平均(平均的)长度分别为5. 8mm和7. 6mm。它们的最大长度分别为 18. 2mm和20. 6mm。低度病变可能只占据四分之一的宫颈，或占据外宫颈的表面区域的很小比例。cmi的分布随单病灶(病变发生在单一处)到多病灶(病变发生在多处)而不同。多个，不同的，小的，随机分布的病变是Cim的特点。CIN2总是出现在转移区(TZ，transformation zone)。TZ是宫颈的已被新化生的鳞状上皮取代的柱状上皮区域。相比于Cim(其可能被发现于TZ外)，CIN 2可能是多病灶， 但是单病灶病变更为常见。卫星状病变通常不是CIN2的代表。阴道镜专家通常沿鳞-柱交界处(SCJ)发现CIN2，要么位于外宫颈上要么位于宫颈内道内。CIN3病变倾向于是融合性的(汇至一起或混为一个)，比Cim或CIN2病变更长更宽。CIN3通常位于宫颈内的中央部分，在内弯道中朝向外口(os)。当出现在外宫颈时 CIN3很少是隐蔽的(眼睛看不到)。CIN3病变的线性长度，被定义为在尾部(位于或接近身体的后端)和头部(朝向身体的头部或前端)边缘之间的组织表面的距离，从2mm到22mm不等。平均线性长度范围为6mm到10mm。长的线性病变-那些长于IOmm的，特别是包含有宫颈内时，一直被怀疑为癌症。由于病变的表面区域增加至大于40平方毫米，所以应被怀疑为癌症。据报道，转换区TZ大小，病变(异常生长)大小，明显边缘，血管模式(血管模式)， 以及醋酸白颜色(应用醋酸后的颜色)与组织评级显著相关，并显示宫颈病变的大小可能具有丨临床重要性° (Kierkegaard, 0. , Byralsen, C, Hansen, K. C. , Frandsen, K. H., and Frydenberg, Μ. , Association between colposcopic findings and histology in cervical lesions: the significance of the size of the lesion, Gynecol. Oncol. 57(1) : 66-71. 1995，通过参考结合于此)。循证医学(EBM)被定义为“在对患者个体治疗作出决定时，对目前最好证据的认真、详细、明确的应用”。在EBM中，人类的感知和解释被认为对诊断说明非常有价值。然而， 人类的评估可能是主观的且具有低的观察者间的一致性(Ferris，D. G. and Litaker, Μ. , Interobserver agreement for colposcopy quality control using digitized colposcopic images during the ALTS trial, J. Low. Genit. Tract Dis. 9(1): 29-35，2005,通过参考结合于此；Jeronimo, J.，Massad, L. S.，Castle, P. Ε.， Wacholder, S. , and Schiffman, Μ. , Interobserver Agreement in the Evaluation of Digitized Cervical Images, Obstet. Gynecol. 110(4) : 833-840，2007,通过参考结合于此)。
计算机可被编程以执行推理任务，并且可用于解决诊断问题。阴道镜评级系统被阴道镜专家们用于获得临床诊断已超过40年。阴道镜专家们通过使用阴道镜的视觉检查定性地获得患者的阴道镜检查特征，确定疾病程度并基于他们的经验开始对患者诊治。这个过程在妇科医生的办公室需要花费10或15分钟。由于检查的主观性，阴道镜检查的准确度高度依赖于阴道镜专家的经验和专业知识。使用了量化的阴道镜特征的计算机实现的阴道镜打分算法是可取的以给医生提供建议的诊断决定。大部分已有的组织分类系统依赖于统计方法。一个例子就是Georgia Tech Vision (GTV)系统(Dickman, Ε. D. , Doll, Τ. J. , Chiu, C. K. , and Ferris, D. G., Identification of Cervical Neoplasia Using a Simulation of Human Vision, Journal of Lower Genital Tract Disease 5(3) : 144-152，2001，通过参考结合于此)，一种由人类视觉系统给予灵感的通用的计算机视觉系统识别正常和异常的宫颈特征。另一个例子是用于宫颈肿瘤的无创诊断的多模型高光谱成像(MHI)系统(Ferris， D. G.，Lawheadj R. A.，Dickman, E. D.，Holtzapplej N.，Miller, J. Α.，Groganj S. et al.，Multimodal hyperspectral imaging for the noninvasive diagnosis of cervical neoplasia, J. Low. Genit. Tract Dis. 5(2): 65-72，2001，通过参考结合于此)。最近，已经报道了用来区分高度上皮内病变与正常或低度上皮内病变的反射荧光系统。(Park, S. Y. , Fo 11 en, Μ. , Milbourne, Α. , Rhodes, H. , Malpica, Α., Mackinnon, N. et al. , Automated image analysis of digital colposcopy for the detection of cervical neoplasia, J. Biomed. Opt. 13(1) : 014029-1-014029-10，2008，通过参考结合于此)。此外，也已经报道了用于高度宫颈上皮内肿瘤的活体内(在有生命的体内)监测的成功的光学系统 LUMA (Huh, W. K. , Cestero, R. Μ. , Garcia, F. Α. , Gold, Μ. Α. , Guido, R. S. , Mclntyre—Seltman, K. et al. , Optical detection of high-grade cervical intraepithelial neoplasia in vivo: results of a 604- patient study, Am. J. Obstet. Gynecol. 190(5) : 1249-1257. 2004,通过参考结合于此；Kendrick, J. Ε. , Huh, W. K. , and Alvarez, R. D. , LUMA Cervical Imaging System, Expert Review of Medical Devices 4(2) : 121-129，2007，通过参考结合于此)。在 LUMA 系统中，在醋酸应用之后使用了 3种集成的光学测量执行整个人体宫颈的诊断扫描激光诱导荧光光谱，白光漫反射光谱，以及视频成像。进一步地，包括多变量分类算法，以利用特定光源与不同的宫颈组织类型相互作用中的光谱差异。最近报道了另一个称为DySIS (动态光谱成像系统)的活体内检测和量化的评级系统(Forthphotonics, DySISiDynamic Spectral Imaging System, http://www. forth-photonics. gr/dysis. php page=dysis_trials2008，通过参考结合于此；Soutterj W. P.， Diakomanolisj Ε.，Lyons, D.，and Haidopoulosj D.，Dynamic Spectral Imaging -In Vivo Detection and Quantitative Grading of Cervical Neoplasia, ASCCP 2008 Meeting Abstract. Journal of Lower Genital Tract Disease 12(2): 160， 2008，通过参考结合于此)。DySIS皿是第一个CE标记的(CE标记是在欧盟地区(EEA)的单一市场中投放的许多产品上的强制认证标记一CE标记证明产品符合欧盟的卫生，安全和环境要求， 这保证了消费者的安全)设备，以使用数字和动态光谱成像技术，用于帮助活体内自动(和用户无关)检测以及绘制宫颈瘤样病变(病变或肿瘤)。
本发明是一种基于阴道镜评级系统以准确评估及识别低度病变，高度病变和癌症的计算机辅助检测技术。它不同于现有技术中列出的统计学习方法。本发明是一种遵循阴道镜检查的基于规则的方法，所述方法使用上皮细胞和血管的形态学评价积累了疾病的证据。在统计学习中，诊断判断是基于经交叉验证从大量训练数据导出的统计特征。因此， 它是基于训练数据和统计理论的监督学习方法。目前的发明将现代阴道镜中的定性评级规则转移到定量计算机程序中。它不需要来自数据的任何训练。它是基于已有的医学知识的无监督方法。图1中示出了示出本发明概述的流程图。在应用本发明的组织分类算法之前，首先采用其他方法来提取信息，该信息用作组织分类算法的输入。其他方法包括解剖特征算法，醋酸白特征提取算法，镶嵌区域(mosaic)和点状区域(pimctation)检测算法(例如，W. Li and A. Poirson, Detection and characterization of abnormal vascular patterns in automated cervical image analysis, Lecture Notes in Computer Science - Advances in Visual Computing 4292, Second International Symposium, ISVC 2006, Lake Tahoe, NV, November 2006 Proceedings, Part II, 627-636 (Springer 2006],通过参考结合于此)以及非典型血管提取算法(如在共同
I^gih白勺、11 "Methods for Detection and Characterization of Atypical Vessels in Cervical Imagery"，2008年8月1日提交，序列号12/2213 的专利申请公开的，其通过参考结合于此)。这些方法的输出作为本文所描述的组织分类方法的输入。组织分类方法应用于在临床检查期间，应用醋酸之前和之后实时提取的宫颈的图像数据，以及对象的人口统计学数据，例如年龄。组织分类方法的输出包括具有以下检查结果的诊断判断病变类型(正常/NED (无疾病证据)，低度增生，高度增生，或癌症)；疾病位置或高度 (CIN2和CIN3)病变位置；以及判断的置信度。下面的专利和专利申请可认为与本发明的领域相关
Costa等的美国专利No. 7309867，其通过参考结合于此，公开了通过利用统计和非统计分类技术的结合，以及通过结合光谱数据和图像数据，确定组织样本的给定区域包含给定类型的组织，例如CIN，CIN 11/111，正常鳞状上皮，正常柱状以及组织转化的概率的方法。Zelenchuk的美国专利No. 731(^47，其通过参考结合于此，公开了一种系统和方法，用于在原位区分健康和发病组织。它包括直接将紫外线照射至组织标本，并收集荧光响应辐射的光纤探头。响应辐射在三个选定的波长被观察，其中一个对应于等吸光点。在一个例子中，等吸光点发生在大约431纳米。使用431nm强度将观察到的信号强度归一化。使用判别分析中的比例确定得分。根据判别分析的结果，基于疾病或健康的诊断确定被检查组织是否切除。Kaufman等的美国专利No. 7260M8，其通过参考结合于此，公开了与基于测量相似性的多个图像相关的方法。这些方法在分割组织的宫颈图像序列中有用。这些方法可以用于确定宫颈组织的醋酸变白测试中的组织特性。它公开了确定组织特性的方法，包括获得多个组织图像；确定这些图像中的两个或更多区域之间的关系；基于所述关系分割图像；并基于所述分割确定所述组织的特性。确定步骤包括将组织表征为正常，CIN I, CIN II，CIN III 或 CIN 11/111。Utzinger等的美国专利No. 6766184，其通过参考结合于此，公开了生成组织的多光谱图像的方法和装置。所述图像可用作宫颈癌检测。利用辐射源产生主要辐射，并过滤以选择第一波长和第一偏振。经过滤的主要辐射照射组织，生成次级辐射，其经过滤以选择第二波长和第二偏振。经过滤的次级辐射用探测器收集，并且根据第一、第二波长和偏振的不同组合，生成组织的多个多光谱图像。Roehrig等的美国专利No. 6198838，其通过参考结合于此，公开了通过利用独立计算的块(mass)和毛刺信息来检测数字乳房X线照片中的可疑部分的方法和系统。所述方法在计算机辅助诊断系统中使用，该系统设计为使纤维乳腺组织中的可疑或可能为癌的病变引起放射学家或其他医疗专业人士的关注。在优选实施例中，毛刺信息和块信息被独立计算，由于算出的毛刺信息不依赖于块信息计算的结果，从而具有更高的可靠性。Turner等的美国专利No. 6135965，其通过参考结合于此，公开了一种利用神经网络分析活体内荧光光谱测量来光谱检测宫颈组织中组织异常的装置和方法。Gaborski等的美国专利No. 5857030，其通过参考结合于此，公开了一种对放射图像进行数字图像处理的自动化的方法和系统，包括过滤的预处理阶段，对分割预选阶段， 以及包括神经网络分类的模式分类阶段。Clarke等的美国专利No. 5982917，其通过参考结合于此，公开了一种计算机辅助诊断(CAD)方法和装置，用于对数字X射线图像中的可疑区域检测和增强。Balas的美国专利公开No. 2006/0141633，其通过参考结合于此，公开了一种用于活体内检测并绘制(map)在癌症发展期间上皮组织的生化和/或功能特性引起的变化的方法和装置。Keller等的美国专利申请No. 2006/0184040，其通过参考结合于此，公开了检测组织异常的方法和设备，据此，所述方法包括从光源向组织照射光；将从组织反射的照射光通过光学器件引至多波长成像光学子系统；并分离感兴趣的一个或多个波长或波带；将感兴趣的一个或多个波长或波带引至一个或多个成像设备，并使用所述设备记录感兴趣的一个或多个波长或波带的图像；将图像数据从图像传输至计算系统；并对组织异常相关的一个或多个光谱模式分析所述图像。Gu等的美国专利公开No. 2008/0058593，其通过参考结合于此，公开了一种在用内窥镜检查期间从器官的视频数据提供计算机辅助诊断的过程，包括分析和增强来自视频的图像帧，并在检查期间实时检测和诊断图像帧中的任何病变。进一步地，图像数据可被用于生成器官的三维重建。
发明内容
这里以及下文更充分描述的本发明包括一种用于将宫颈组织分类的基于规则的无监督过程，其包括对宫颈组织应用包含一组分类器规则的分类器，通过串行地应用选定的分类器规则来获得最终预测，其中，采用选自由确定以下组成的组的分类器规则执行应用步骤纹理区域的大小、不透明度参数、醋酸白区域大小、粗大和细小点状区域的数目、粗大和细小镶嵌区域的数目、非典型血管大小以及人口统计学数据，由此，可将宫颈组织分类为无疾病证据、低度异型增生、高度异型增生或癌症。如果纹理区域大于宫颈面积的约1%，以及如果不透明度参数小于像素值的约4%， 或大于像素值的约12%，则将获得无疾病证据的最终预测。
如果纹理区域大于宫颈面积的约1%，所述不透明度参数介于像素值的约4%至像素值的约1 之间，醋酸白的区域大小大于宫颈面积的7%，以及如果口附近的醋酸白区域大小大于宫颈面积的0%，则将获得高度异型增生的最终预测。如果纹理区域大于宫颈面积的约1%，如果所述不透明度参数介于像素值的约4% 至像素值的约1 之间，醋酸白的区域大小小于宫颈面积的7%，并且口附近没有醋酸白区域，则将获得低度异型增生的最终预测。如果纹理区域小于宫颈面积的约1%，并且如果不透明度参数小于像素值的5%，粗大点状区域的数目小于2，且没有粗大的镶嵌区域血管，则将获得低度异型增生的最终预测。如果纹理区域小于宫颈面积的约1%，并且如果不透明度参数大于像素值的10%， 口附近的醋酸白区域大小大于15%，则将获得癌症的预测。如果纹理区域小于宫颈面积的约1%，并且不透明度参数大于8%，口附近的醋酸白区域的大小大于0. 5%，并且粗大点状区域的数目大于100，并且没有粗大的镶嵌区域血管， 并且异常血管大小大于0. 5毫米，则将得到癌症的预测。如果纹理区域小于宫颈面积的1%，并且如果不透明度参数大于像素值的5%，口附近的醋酸白区域的大小小于像素值的0. 5%，粗大点状区域血管的数目小于100，不存在粗大镶嵌区域血管，并且非典型血管的大小小于0. 5%，则将得到高度异型增生的判断。
图1展示了本发明的概述的流程图。
图2是用于组织分类方法的输入参数的列表。
图3是组织诊断方案的第一部分。
图4是组织诊断方案的第二部分。
图5是NED (没有疾病证据)的对象的示例的照片。
图6 (a)是低度疾病对象的示例的照片。
图6 (b)是图6 (a)的对象的对应低度疾病位置的图像处理过的照片。
图7 (a)是癌症对象的示例的照片。
图7 (b)是图7 (a)的对象的对应高度疾病位置的图像处理过的照片。
图7 (c)是图7 (a)中的对象的对应低度疾病位置的图像处理过的照片。
图8 (a)是高度疾病对象的示例的照片。
图8 (b)是图8 (a)中的对象的对应的高度疾病位置的图像处理过的照片
图8 (c)是图8 (a)中对象的对应的低度疾病位置的图像处理过的照片。
具体实施例方式
1.预方法应用
本文所述的本发明目前优选的实施例优选使用几个预方法应用的输出(输出为如图2 中所描述的某些参数)作为组织分类方法的输入。很多参数或变量被描述为图像中特定特征的面积相比于宫颈的总面积的百分比。所述的点状区域特征被定义为图像中存在的点状区域的数目。非典型性血管的大小优选的以毫米(mm)来计算。
该过程首先收集宫颈的数字彩色图像(使用红色，绿色，蓝色通道或RGB)，——一个在施加醋酸之前，一个在应用醋酸之后。优选地，图像被正交极化，以抑制镜面反射(闪烁)并被校准。接着，采用几种预方法应用，包括解剖特征提取，醋酸白特征提取，镶嵌区域和点状区域检测，以及非典型血管提取。解剖特征提取优选地包括检测子宫颈，除去阴道壁(水平和竖直)，口检测，柱状上皮检测以及纹理区域(异常/化生)检测。醋酸白特征提取优选地包括应用醋酸到感兴趣的解剖区域，以产生醋酸白区域； 拍摄至少一张醋酸白区域的后醋酸图像；检测感兴趣的具有上皮细胞的解剖区域；从后醋酸图像提取出纹理信息，以定义纹理区域；从后醋酸图像中提取色彩信息以定义色彩区域； 并将纹理区域和色彩区域相结合生成候选的醋酸白区域。优选地，它还包括获得至少一个上皮细胞的前醋酸图像；在前醋酸图像之一中和在后醋酸图像之一中记录醋酸白区域， 以形成醋酸白区域的记录图像；从醋酸白区域的记录图像中提取不透明度参数，如共同待决共同受让的2008年8月4日提交的序列号为12/221645的题为“Computerized Image Analysis for Acetic Acid Induced Cervical Intraepithelial Neoplasia，，的专利申请，其通过参考结合于此。该申请公开了通过聚类在记录的前醋酸图像和后图像的时间变化中的差异，来创建的不透明度参数。所述不透明度参数对病变诊断提供了高预测值—— 正常和低度病变具有比高度病变和癌症情形低得多的不透明度。其他的不透明度参数提取方法也可被使用，包括在不同色彩空间中的比率成像。在本发明的一优选实施例中，在CIE-Lab色彩空间中的a通道中的色彩差异被用于不透明度参数提取。在本发明另一优选实施例中，RGB色彩空间中的g通道的色彩差异被用于不透明度参数提取。由于不成熟的化生和柱状上皮组织在醋酸施加后短暂变为白色， 但并未表现出发育异常的组织变化，他们不应被排除在感兴趣的醋酸白区域之外。这些组织区域通常表现出轻微的不透明度变化。因此本发明在色彩差异特征空间内采用了两步平均移位聚类算法。第一步是分割主要的不透明度变化，并消除轻微不透明度变化。第二步是从第一步获得的前景区域中分割最不透明的变化。计算不透明度参数，作为最不透明区域的平均色差。最不透明区域被定义为具有最大平均色差的区域。不透明度参数通过以下公式来表示
权利要求
1.一种用于将宫颈组织分类的基于规则的无监督方法，包括对所述宫颈组织应用分类器，所述分类器包括一组分类器规则；并且通过串行地应用在所述组中选定的分类器规则来获取最终预测，其中使用选自由确定以下组成的组的分类器规则来执行所述应用步骤纹理区域的大小， 不透明度参数，醋酸白区域的大小，粗大和细小点状区域的数目，粗大和细小镶嵌区域的大小，非典型血管的大小以及人口统计学数据；由此，所述宫颈组织可被分类为无疾病证据、低度异型增生、高度异型增生或癌症。
2.根据权利要求1所述的方法，其中如果所述纹理区域大于所述宫颈面积的约1%，并且所述不透明度参数小于像素值的约 4%或大于像素值的约12%，获取无疾病证据的最终预测。
3.根据权利要求1所述方法，其中如果所述纹理区域大于所述宫颈面积的约1%，所述不透明度参数介于像素值的约4% 和像素值的约1 之间，醋酸白区域的大小大于所述宫颈面积的7%，并且口附近的醋酸白区域的大小大于所述宫颈面积的0%，则获取高度异型增生的最终预测。
4.根据权利要求1所述方法，其中，如果所述纹理区域大于所述宫颈面积的约1%，所述不透明度参数介于像素值的约4% 和约12%之间，醋酸白区域的大小小于宫颈面积的7%，并且在口附近没有醋酸白区域，则获取低度异型增生的最终预测。
5.根据权利要求1所述方法，其中，如果所述纹理区域小于所述宫颈区域的约1%，所述不透明度参数小于像素值的5%，粗大点状区域的数目小于2，并且没有粗大镶嵌区域血管，则获取低度异型增生的最终预测。
6.根据权利要求1所述方法，其中如果所述纹理区域小于所述宫颈面积的约1%，所述不透明度参数大于像素值的10%， 并且口附近的所述醋酸白区域的大小大于所述宫颈面积的15%，则获取癌症的最终预测。
7.根据权利要求1所述方法，其中如果所述纹理区域小于所述宫颈面积的约1%，所述不透明度参数大于像素值的8%，口附近的所述醋酸白区域的大小大于宫颈面积的0. 5%，粗大点状区域的数目大于100，没有粗大镶嵌区域血管，并且异常血管的大小大于0. 5mm，则获取癌症的最终预测。
8.根据权利要求1所述方法，其中，如果所述纹理区域小于所述宫颈面积的约1%，所述不透明度参数大于像素值的5%，口附近的醋酸白区域的大小小于所述宫颈面积的0. 5%，粗大点状区域血管的数目小于100， 没有粗大镶嵌区域血管，并且非典型血管的大小为小于0. 5mm,则获取高度异型增生的最终预测。
9.根据权利要求1-8任一项所述方法，进一步包括如果患者年龄低于25岁或高于 60岁，则将最终的癌症预测中的置信度降低约10%-约20%之间。
10.一种用于将宫颈组织分类的基于规则的无监督方法，包括对所述宫颈组织应用分类器，所述分类器包括一组分类器规则；并且通过串行地应用在所述系列中选定的分类器规则来获取最终预测，其中使用选自由确定以下组成的组的分类器规则来执行所述应用步骤纹理区域的大小，不透明度参数，醋酸白区域的大小，粗大和细小点状区域的数目，粗大和细小镶嵌区域的大小，非典型血管的大小以及人口统计学数据；由此，所述宫颈组织可被分为无疾病证据、低度异型增生、高度异型增生或癌症； 其中如果所述纹理区域大于所述宫颈面积的约1%，并且所述不透明度参数小于像素值的约4%或大于像素值的约12%，则获取无疾病证据的最终预测；如果所述纹理区域大于所述宫颈面积的约1%，所述不透明度参数介于像素值的约4% 和像素值的约1 之间，醋酸白区域的大小大于所述宫颈面积的7%，并且如果口附近的醋酸白区域的大小大于所述宫颈面积的0%，则获取高度异型增生的最终预测；其中如果所述纹理区域小于所述宫颈面积的约1%，所述不透明度参数大于像素值的 5%，口附近的所述醋酸白区域的大小小于宫颈面积的0. 5%，粗大点状区域血管的数目小于 100，没有粗大镶嵌区域血管，并且异常血管的大小小于0. 5mm,则获取高度异型增生的最终预测；其中如果所述纹理区域大于所述宫颈区域的约1%，所述不透明度参数介于像素值的约 4%和像素值的约1 之间，醋酸白区域的大小小于宫颈面积的7%，并且在口附近没有醋酸白区域，则获取低度异型增生的最终预测；其中如果所述纹理区域小于所述宫颈区域的约1%，所述不透明度参数小于像素值的 5%，粗大点状区域的数目小于2，并且没有粗大镶嵌区域血管，则获取低度异型增生的最终预测；其中如果所述的纹理区域小于所述宫颈面积的约1%，所述不透明度参数大于像素值的 10%，口附近的所述醋酸白区域的大小大于所述宫颈面积的15%，则获取癌症的预测；其中如果所述的纹理区域小于所述宫颈面积的约1%，所述不透明度参数大于像素值的 8%，口附近的所述醋酸白区域的大小大于宫颈面积的0. 5%，粗大点状区域的数目大于100， 没有粗大镶嵌区域血管，并且异常血管的大小大于0. 5mm，则获取癌症的最终预测。
11.根据权利要求10所述方法，进一步包括如果患者年龄低于25岁或高于60岁，则将最终的癌症预测中的置信度降低约10%-约 20%之间。
12.一种用于将宫颈组织分类的基于规则的无监督方法，包括 应用分类器，所述分类器包括确定不透明度参数的分类器规则；由此，所述宫颈组织可被分为无疾病证据、低度异型增生、高度异型增生或癌症。
全文摘要
一种用于将宫颈组织分类的基于规则的无监督方法，所述分类通过串行地应用选自由确定以下组成的组的分类器纹理区域的大小、不透明度参数、醋酸白区域的大小、粗大和细小点状区域的数目、粗大和细小镶嵌区域的大小、非典型血管的大小以及人口统计学数据，使得所述宫颈组织可被分类为无疾病证据、低度异型增生、高度异型增生或癌症。
文档编号G06K9/00GK102282569SQ200980139997
公开日2011年12月14日 申请日期2009年10月9日 优先权日2008年10月10日
发明者李 W. 申请人:国际科学技术医疗系统有限责任公司
最新回复(0)