基于基元结构定义与区域生长技术的虹膜分区及阳光放射沟提取方法

专利名称：基于基元结构定义与区域生长技术的虹膜分区及阳光放射沟提取方法
技术领域：
本发明涉及计算机图像处理技术领域，具体涉及一种利用虹膜进行疾病诊断的技 术。
背景技术：
目前，随着社会的进步和人民生活水平的日益提高，人们越来越关注自身的健康， 迫切需要一种方便可靠、无创伤非接触式的自我诊断手段。虹膜学与虹膜诊断正是这样一 种现代新型的诊病方法，它具有非接触性、无伤痛性、易于普及的特点，这也是未来医学与 预防医学诊断手段的重要发展趋势。中西医虹膜学认为人眼虹膜是人体中与全身各脏器密切相关的器官之一，各部 脏器在虹膜上有确切的映射区，其变化情况直接反映人体内部各脏器的健康状况。虹膜异 常主要体现在各个脏器映射区内，色调的变化、局部色斑的出现、虹膜纤维破洞的出现或虹 膜纤维疏密程度的变化等方面。通过观察虹膜不同区位的异常变化，可以对某些器官的疾 病做出诊断，如虹膜上方脑映射区内某种亮点表示脑神经有异常、肠胃映射区瞳孔周围出 现凹点说明有溃疡病等。利用虹膜诊病的关键所在是定位各部脏器在虹膜上的映射区。为此中西医虹膜学 均给出了虹膜图谱。从西方虹膜学之父匈牙利医生Ignatz Von Peczley到美国虹膜大师 Bernardjensen，西方虹膜学从实验观察和统计分析出发，其图谱定位准确，所以提出的诊 断标准比较明确和稳定。但西医对虹膜分区过于复杂，时常将一个完整病症被分割在几个 区域中，没有摆脱“人是无数零件组成的机器”的观点，因此在临床上无法把这些活体器官 之间的联系统一起来；另一方面西医分区在利用计算机辅助诊断方面的可操作性较差，因 为目前的虹膜图像采集设备和虹膜图像处理算法上还不能达到直接利用西医图谱诊病所 要求的精度。我国中医“眼诊学”认为眼虽是一个局部器官，但它与全身的脏腑、经络共同构成 了一个有机的整体在生理上相互协调，在病理上相互影响，有着非常密切的关系。所以中 医虹膜图谱在整体上给出了映射区在虹膜上的大致位置。因此造成了中医图谱精度较差， 并且目前还没有一个统一的、严格的和量化的标准。此外，中西医图谱还存在一个共同的问题病灶出现的位置会跨越相邻的两个分 区，出现这种情况包括两方面的情况一、相关脏器(如生殖泌尿系统中的各脏器)出现在 相邻的映射区中，从疾病的发生、发展来考察，他们是相互关联、相互影响的，其病灶也有可 能扩大到其它映射区中；另一方面，某个脏器病变较重(如重症冠心病，其心脏映射区病理 变化较大)。

发明内容
本发明的目的在于提供一种有机结合中、西医的虹膜学成果一依据西医虹膜学中较准确定位肌体脏器与虹膜映射区域的关系、根据中医理论对虹膜映射区形成的阐述， 给出基于基元结构定义与区域生长技术的虹膜分区及阳光放射沟提取方法。本发明的目的是这样实现的1、基于中西医虹膜学的计算机辅助诊断虹膜图谱(CADIC computer aided diagnosis irischar)本发明的虹膜图谱定义如下对于一幅虹膜图像，经过前期预处理与虹膜内外边 缘精确定位后，将虹膜图像置于极坐标系(Θ，r)中，把虹膜分为以瞳孔圆心为极点、以瞳 孔边缘为起始边界的五个环、12个区的虹膜分区图谱。设一幅虹膜图像用A表示，Rt为瞳孔半径，R1为虹膜半径，Rd为图像中卷缩轮的理 论半径，R' D为卷缩轮的实际半径。极坐标的极轴方向取为时钟的零点方向，极角按照顺 时针方向定义。将A划分为五个环12个区后，其中的任意一个区用Aik(e，r)来表示。其中k = 1，2，· · ·，12 ;i = 1，2，· · ·，5，每个分区的定义如下虹膜分区是按时钟顺时针方向从12点标记开始每30度为一个区；在虹膜环形的 半径方向上5等分，构成5个环。根据虹膜学理论，虹膜中卷缩轮大小占虹膜的三分之一左 右，所以最里边的两个环被合并表达卷缩轮部分。再根据虹膜的另外三个环和分区，将虹膜 分成各个诊断区域。本发明根据制定的分区规则和医学虹膜图谱进行了比照，得到了一一 对应关系，表1-8给出了左右眼虹膜CADIC图谱中各个分区对应的肌体脏器在虹膜上的映 射部位。2、基于基元定义和区域生长技术的虹膜处理方法本发明提出的动态自适应虹膜分区(图谱覆盖)方法针对虹膜图像纹理复杂、区 域变化不明显特性，在基于串行区域生长与合并技术基础上，通过对特征先验知识，定义区 域生长与合并准则，以“基元”为区域生长点，具体地方案包含如下步骤1.虹膜图像的采集；2.虹膜图像的预处理；3.虹膜定位；4. CADIC虹膜图谱覆盖；5.基于基元定义和区域生长的阳光放射沟提取；6.虹膜自适应分区。各个处理步骤的细节如下1.虹膜图像的采集采样方法是对视频图像抓拍；采样者 部固定并正视；采样时所采用的光源产生的光斑打在瞳孔和虹膜的交界处；采样图像为768*576的24位位图。每个采样者采3到5 组虹膜图像。2.虹膜图像的预处理在对虹膜图像进行分区分析之前，要进行预处理工作。虹膜图像的预处理主要包 括彩色图像灰度转换、噪声处理、几何校正、图像增强。
2. 1)图像灰度化
对采集到的彩色图像，应用如下公式将其变换为灰度图像。 Y = O. 299 X R+0. 587 X G+0. 114 X B (2-2)式中R——真彩图中每个像素的红分量；G——真彩图中每个像素的绿分量；B——真彩图中每个像素的蓝分量；Y——该像素所对应的灰度2.2)虹膜图像灰度归一化通过将原始虹膜图像中感兴趣的区域映射到相同的尺寸和相同的对应区域中进 行虹膜归一化，从而消除平移、旋转对纹理检测的影响。2. 3)光斑填充采用数字图像的自适应插值法来对光斑进行填充。零阶插值算法，即令输出像素 的灰度值等于离它所映射的位置最近的输入像素的灰度值。3.虹膜定位除去那些于虹膜无关的信息(如眼睑Eyelid，睫毛Eyelash，瞳孔Pupil等)，也就 是从虹膜图像中定位出虹膜的内边界和外边界。3. 1)虹膜粗定位通过对虹膜图像的分析，得出以下的先验知识瞳孔的直径最小也要大于120个 像素；瞳孔是眼睛图像中面积最大的低灰度区域，而其圆心位置应该在上下左右灰度和值 最小的区域附近，利用这些先验知识，可以减小搜索的空间提高精度。本文采用传统的灰度 投影法对虹膜内、外边缘粗定位，得到虹膜内径r、外径rl和瞳孔中心(x。，y。)；3. 2)虹膜精定位本发明采用Daugman算法来对虹膜进行精定位，其核心在于利用一个有效的积分 微分算子来计算圆参数。虹膜具有良好的环状特性，所以可采用圆形探测算子对虹膜边界 进行精定位。在虹膜粗定位瞳孔中心(x。，y。)和半径r的基础上，采用圆探测算子实现虹膜精定 位；所述虹膜精定位方法为基于圆探测算子的定位方法a、对于中心矩阵内的每一个元素，分别计算半径为r同一圆周上的像素点灰度和 sum，其中中心矩阵的范围是(x。-5，y。_5)，半径变化范围为r_5 r+5 ；b、对于r的变化值d，计算出对应的高斯函数值
8 c、计算灰度梯度值向量det = g(d+l)-g(d) I X sum；d、依次计算所有元素对应的det向量值；e、在所有的det向量中取元素最大值，即为精定位的结果。4. CADIC虹膜图谱覆盖在虹膜内外边缘精确定位后，将虹膜部分所在区域的直角坐标(X，y)形式转换为 极坐标(r，θ)形式。根据虹膜图谱和r、θ的值来将虹膜部分分成若干个区域。这是第一 次分区，也就是理论分区。5.基于基元定义和区域生长的阳光放射沟提取本发明选取虹膜图像中疾病形成的阳光放射沟为处理对象，首先定义纹理基元并 讨论基于基元的纹理表达方法。设纹理基元为h(x，y)，排列规则为r(x，y)，则纹理t(x，y)定义为t{x,y) = h{x,y) r{x,y)(4-1)设xm，ym为脉冲函数的位置坐标，则r(x, y) =Σ δ (x-xm, y-ym) (4-2)根据卷积定理，在频域有T (U，ν) = H (U，ν) R (U，ν)(4-3)R (U，ν) = T (U，ν) H (U，νΓ1 (4-4)那么对于给定纹理基元h(x，y)的描述，即可得到反卷积滤波器H (u，ν) Λ将这个 滤波器用于待处理纹理图像，每个脉冲都在纹理基元的中心。对整幅图像中不同纹理基元 的分布进行统计可获得图像的全面纹理信息。通过纹理基元确定对了对待分割的区域找到一个种子子图一即基元脉冲中心， 然后将种子像素周围邻域中与种子像素有相同或相似性质的像素(通过定义的相同或相 似规则)合并到种子像素所在区域中。在将新像素作为种子继续扩展，直到再没有满足条 件的像素可被包括，这样的一个区域就长成。对于一个给定的目标图像X和一个结构基元S,以符合S的种子区域ΧΘ5为中心 在图像上进行迭代生长。对于种子Z φ S的每一个邻近像素X，在生长过程中应始终保持if S(x) = TRUE(4-5) S卩，如果区域邻近像素χ能符合结构基元定义，则将像素χ合并到特征区域，使之 成为新种子的一部分。如果对于种子1十^的所有邻近像素X，均有S (χ) = FALSE(4-7)则生长停止。对于阳光放射沟特征，通过前期对临床病理的研究分析，从获得的先验统计知识 定义像素生长规则如下。规则1 沟长度为L个像素，5彡L彡R虹膜；
规则2 沟宽度为W个像素，3彡W彡5 ；规则3 沟的延长线与瞳孔交点的切线相垂直；规则4 根据虹膜图像灰度直方图确定沟区域灰度最优阈值Te ；规则5 沟的起点为卷缩轮边界点。根据CADIC虹膜分区图，对每个感兴趣的诊断候选区域，执行如下阳光放射沟检 测算法。1)计算阳光放射沟灰度阈值对整个虹膜统计灰度直方图，定义灰度阈值Te，使得灰度小于Te的像素数目为整 个虹膜像素数目的10%。根据不同的虹膜图像采集设备，此比例数应根据人工标注阳光放 射沟情况的先验统计数据做出相应的调整。2)对任一个候选区Q，执行以下3) 10)。3)沟起点选取对于任意像素χ e Q，如果有符合规则5的像素χ的灰度值f (x) ( tG,则认为χ是 一个沟起点。把找到的所有沟起点记为dn，η = 1,2,3...N04)如果N = 0，即没有找到任何沟起点，则说明没有阳光放射沟存在，转10)结束 本算法。5)沟起点筛选与合并，形成种子候选对沟底点记dn，η = 1,2,3... N，先后执行生长与合并操作，合并成种子候选d’ n， η = 1，2，3· · · N，，N，< N。6)种子候选筛选，形成种子对每个种子候选d’ 定义它在平行于虹膜内边缘方向上的像素跨度数为其宽度 W(d’。，删除所有^ n) > 5的种子候选d’ n。所有剩余的种子候选d’ n成为种子，记为zm，m = 1，2，3. . . M。在不造成混淆的情 况下，也用Zm代表本种子的中心像素。7)如果M = 0，即没有找到任何种子，则说明没有阳光放射沟存在，转10)结束本算法。8)区域生长，形成阳光放射沟候选在每个种子zm，m = 1,2,3. .. M周围邻域中搜索同种子像素具有相同或相似性质 的像素，合并到种子像素所在区域中；在将新像素作为种子继续迭代扩展，直到再没有满足 条件的像素可被包括，这样一个区域就长成了。8. 1)定义每个种子(中心)Zm的邻域为NB(Zm)，NB(Zm)应符合规则1，2和3 ；8. 2)定义基元结构为S，此处的S定义为典型阳光放射沟的横截线上的像素集合。8. 3)对于任意像素χ e NB (Zm)，根据公式(4_5)、(4-6)所定义的生长规则进行区 域生长。9)阳光放射沟筛选对获得的阳光放射沟候选，根据规则1和2进行筛选。10)输出结果，结束算法。通过上述方法实验，可以获得满足裂缝基元定义的线条，即为阳光放射沟。6.虹膜动态自适应分区
具体做法体现在对病理特征(此处为阳光放射沟，也可以推广到其他虹膜特征 上)在分区边界上的处理当分区时病理特征被分割开来，不能被完整地表达在同一个区 域中时，进一步从保持病理特征完整性的角度，重新确定该特征所覆盖区域。6. 1)首先，通过估计病理特征的区域进行病理特征筛选对病理特征边缘形状的提取采用估算的方法。求出组成特征的像素点在极坐标中 的角度和极径的最大和最小值，这样就大概确定出了纹理所覆盖的区域的面积。从前面特 征提取的结果可以看到，只有一部分重要的纹理特征(阳光放射沟)具有诊病价值。我们 根据所确定的每个特征区域的面积，选取最重要的若干个特征，比如5个；删除其余的不重 要特征。6. 2)结合图15，进行CADIC图谱覆盖，方法如前述。6. 3)判断是否有跨区域的病理特征对每一个虹膜图谱上的分区注意考察，如果发现其边界同重要特征有交叠，则记 录这个交叠。6.4)进行区域动态调整，完成自适应分区如果一个纹理特征覆盖多个区域，进行如下处理对每个虹膜图谱分区，统计区域 外内特征的像素点数目M ；之后对其每处交叠边界，考察此区域外该特征的像素点并进行 统计数目M’。如果区域外像素点中只有极少部分Μ’ << M(例如M小于5%的M)，则对这些少 数像素点所分布的区域不予考虑；如果区域外像素点过于重要M’> M(说明该特征不是主要处于本区域)，则同样不 予考虑，并删除本区域中的特征；此外，说明区域外特征比较重要；则将此边界邻接的区域同本区域合并，实现了重 新的区域划分，也就是进行了自适应分区。本发明提出了一种融合中西医虹膜学虹膜分区的新虹膜分区图谱(CADIC)；并发 明了应用此图谱所必需的基于基元定义的区域生长技术用于典型虹膜特征一阳光放射 沟的提取，并发明了应用此图谱所必需的基于虹膜图像动态分区技术和虹膜图谱的自适应 覆盖方法。本发明在保证了肌体各部脏器在虹膜上的准确定位同时，突出了虹膜图谱的实 用性，实现了医学意义下完整虹膜特征的动态分割技术及其在虹膜图谱上的定位技术。通 过实验证明，此技术能够大大提高了虹膜识别和诊病可靠性。本发明涉及到的技术领域属 于计算机图像处理、模式识别以及中西医结合与现代化领域，有机结合中、西医的虹膜学成 果，根据西医虹膜学叫准确定位肌体脏器与病症区域、根据中医理论解释、并处理病症的不 同表现，使中医诊疗方法、技术更好地为大众服务，也使中医理论精萃得以传承和发扬，利 用计算机强大的信息处理能力，解决中医诊疗在现代化和信息化建设上面临的困难，这是 本发明的主要意义。采用本发明，可以在保证了肌体各部脏器在虹膜上的准确定位，同时有 效检测虹膜特征，从而大大提高了虹膜识别和诊病的可靠性与实用性。


图1为中医虹膜图谱；图2为Jensen虹膜图谱；
图3-4为CADIC图谱示意图；图5为具体实施方式
流程图；图6为彩色虹膜图像；图7为灰度化后的图像；图8-9为光斑填充前后对比图像；图10为虹膜粗定位图；图11为内、外边缘精定位图；图12为理论分区图；图13为图像区域分布图；图14为使用基元定义的裂缝识别结果图；图15为CADIC图谱覆盖结果图；图16为虹膜动态自适应分区结果示意图一阳光放射沟的自适应分区结果图。
具体实施例方式下面结合附图对本发明作进一步的说明1 基于中西医虹膜学的计算机辅助诊断虹膜图谱(CADIC computer aided diagnosis irischar)计算机辅助诊断虹膜图谱分区方法，既要避免西医虹膜图谱的复杂与实用性差的 弊端，又要克服中医图谱精度不足和定位不准的弱点；其分区方法应能够保证各大脏器在 映射区的完整性和定位的精确性，还要利于计算机实现又方便样本图像的分类存储。根据图1，图2，本实施例的虹膜图谱定义如下对于一幅虹膜图像，经过前期预处 理与虹膜内外边缘精确定位后，将虹膜图像置于极坐标系(Θ，r)中。把虹膜分为以瞳孔圆 心为极点、以瞳孔边缘为起始边界的五个环、12个区的虹膜分区图谱，如图3所示。设一幅虹膜图像用A表示，Rt为瞳孔半径，R1为虹膜半径，Rd为图像中卷缩轮的理 论半径，R' D为卷缩轮的实际半径。极坐标的极轴方向取为时钟的零点方向，极角按照顺 时针方向定义。将A划分为五个环12个区后，其中的任意一个区用AikO，r)来表示。其中k =
1 ，2， · · · ， 12；i= 1 ,2，...，5，每个分区的定义如下 虹膜分区是按时钟顺时针方向从12点标记开始每30度为一个区；在虹膜环形的 半径方向上5等分，构成5个环。根据虹膜学理论，虹膜中卷缩轮大小占虹膜的三分之一左 右，所以最里边的两个环被合并表达卷缩轮部分。再根据虹膜的另外三个环和分区，将虹膜 分成各个诊断区域。我们根据制定的分区规则和医学虹膜图谱进行了比照，得到了一一对应关系，可见，我们所提出的实用化虹膜图谱同传统的中西医虹膜图谱是十分吻合的。见图 4，表1-8给出了左右眼虹膜CADIC图谱中各个分区对应的肌体脏器在虹膜上的映射部位。表1右眼虹膜CADIC图谱 表2右眼虹膜CADIC图谱 表3右眼虹膜CADIC图谱
13 表4右眼虹膜CADIC图谱
冈 表5左眼虹膜CADIC图谱 表6左眼虹膜CADIC图谱
表7左眼虹膜CADIC图谱 表8左眼虹膜CADIC图谱 2、基于基元定义和区域生长技术的虹膜处理方法理论分区是通过对正常人虹膜图像学习而制定的，如果人体有病变发生，在对应 的虹膜区域就会出现明显的纹理变化。而随着病情的发展和治疗手段的引入，纹理变化会 更加明显，甚至出现变化扩散到其它的区域的现象。为了更准确的通过虹膜反映人体病变 情况，就要针对虹膜的纹理变化进行自适应分区。从虹膜解剖学我们知道，正常的卷缩论范围在虹膜图像的1/3。虹膜诊断学理论认
16为，卷缩轮范围变化与消化系统疾病密切相关。如卷缩轮紧缩表示胃部疾病，卷缩轮膨胀， 反映着肠道系统疾病。又如人眼正常虹膜不会出现裂缝，在慢性疾病病理中，会出现这种 长长的裂缝。本发明提出的基于基元的串行区域生长与合并的虹膜图像动态分区技术，较 好地解决了病理情况下虹膜病变特征的定位与测量问题。串行区域生长与合并技术是分析图像中区域结构和相似性的模式识别方法，是建 立在集合代数基础上，用集合论与概率论方法定量描述图像特征的技术。根据分割过程中 处理策略的不同，分割技术又可分为并行技术与串行技术，虽然串行技术所需要计算时间 比并行技术略长，但其抗噪声能力通常较强。我们所提出的动态自适应虹膜分区(图谱覆盖)方法针对虹膜图像纹理复杂、区 域变化不明显特性，在基于串行区域生长与合并技术基础上，通过对特征先验知识，定义区 域生长与合并准则，以“基元”为区域生长点，这样在理论上满足医学图像处理的统计学要 求，又可以便于计算机处理，通过实验得到了很好的结果。具体地方案实施如图5所示，包 含如下步骤1.虹膜图像的采集；2.虹膜图像的预处理；3.虹膜定位；4. CADIC虹膜图谱覆盖；5.基于基元定义和区域生长的阳光放射沟提取；6.虹膜自适应分区。各个处理步骤的细节如下(1)虹膜图像的采集虹膜具有对光的生理反应的特性；毋须物理接触，就可以在一定的距离获得被试 的纹理图样。采样方法是对视频图像抓拍；采样者头部固定并正视；采样时所采用的光源 产生的光斑打在瞳孔和虹膜的交界处；采样图像为768*576的24位位图。每个采样者采3 到5组虹膜图像。(2)虹膜图像的预处理由于虹膜图像的采集受到光照及仪器差异的影响，包含了大量的噪声，同时眼睑 和眼睫毛可能会掩盖一部分的虹膜，而睫状体对于虹膜诊病也是无用的信息，所以在对虹 膜图像进行分区分析之前，要进行预处理工作。虹膜图像的预处理主要包括彩色图像灰度 转换、噪声处理、几何校正、图像增强。(2. 1)彩色图像灰度化结合图6-图7，目前，因为我们只对虹膜图像的纹理特征感兴趣，为了减少计算 量，我们对得到的彩色图像进行灰度化。对采集到的彩色图像，应用如下公式将其变换为灰 度图像。
Y = O. 299 X R+0. 587 X G+0. 114 X B(2-2)式中R——真彩图中每个像素的红分量；
G—真彩图中每个像素的绿分量；B——真彩图中每个像素的蓝分量；Y——该像素所对应的灰度(2·2)虹膜图像灰度归一化由于受到光照、拍照角度等因素的影响，使得图像中虹膜的尺寸和旋转的角度产 生较大差异，这将极大的降低模板匹配的精度。虹膜归一化的目的是将原始虹膜图像中感 兴趣的区域映射到相同的尺寸和相同的对应区域中，从而消除平移、旋转对纹理检测的影 响。(2. 3)光斑填充使用带有双光源的CCD采集设备采集到的虹膜图像，来自光源的入射光在舌面形 成镜面反射，引起虹膜表面颜色的突变而形成光斑，且光斑颜色与光源颜色一致。因而虹膜 图像上属于反光点的像素与虹膜上其他像素相比具有最高的亮度值和最低的饱和度值。这 些光斑区域这会影响到定位结果，需要进行光斑填充。结合图8-9，采用数字图像的自适应插值法来对光斑进行填充。零阶插值算法，即 令输出像素的灰度值等于离它所映射的位置最近的输入像素的灰度值。(3)虹膜定位虹膜定位的目的就是除去那些于虹膜无关的信息(如眼睑Eyelid，睫毛Eyelash， 瞳孔Pupil等)，也就是从虹膜图像中定位出虹膜的内边界和外边界。准确的虹膜定位是实 现有效特征提取的前提。(3. 1)虹膜粗定位结合图10，粗定位的目的是估算出虹膜内、外边缘对应圆的圆心和半径。通过对虹膜图像的分析，得出以下的先验知识瞳孔的直径最小也要大于120个 像素；瞳孔是眼睛图像中面积最大的低灰度区域，而其圆心位置应该在上下左右灰度和值 最小的区域附近，利用这些先验知识，可以减小搜索的空间提高精度。本文采用传统的灰度 投影法对虹膜内、外边缘粗定位，得到虹膜内径r、外径rl和瞳孔中心(x。，y。)；(3. 2)虹膜精定位结合图11，本发明采用Daugman算法来对虹膜进行精定位，其核心在于利用一个 有效的积分微分算子来计算圆参数。虹膜具有良好的环状特性，所以可采用圆形如下算子 对虹膜边界进行精定位。在虹膜粗定位瞳孔中心(x。，y。)和半径r的基础上，采用圆探测算子实现虹膜精定 位；所述虹膜精定位方法为基于圆探测算子的定位方法a、对于中心矩阵内的每一个元素，分别计算半径为r同一圆周上的像素点灰度和 sum，其中中心矩阵的范围是(x。-5，y。_5)，半径变化范围为r-5 r+5 ；b、对于r的变化值d，计算出对应的高斯函数值 c、计算灰度梯度值向量det = g(d+l)-g(d) I X sum；d、依次计算所有元素对应的det向量值；e、在所有的det向量中取元素最大值，即为精定位的结果。
(4) CADIC虹膜图谱覆盖结合图12，在虹膜内外边缘精确定位后，将虹膜部分所在区域的直角坐标(X，y) 形式转换为极坐标(r，θ)形式。根据虹膜图谱和r、θ的值来将虹膜部分分成若干个区 域。这是第一次分区，也就是理论分区。(5)基于基元定义和区域生长的阳光放射沟提取结合图13，分析虹膜主要是分析其结构纹理特征，虹膜诊病更是要将病灶特征独 立而又完整地提取出来。本发明选取虹膜图像中疾病形成的阳光放射沟为处理对象，首先 定义纹理基元并讨论基于基元的纹理表达方法。设纹理基元为h(x，y)，排列规则为r(x，y)，则纹理t(x，y)定义为/0,少)二办(太,少)(8)厂0,力(4-1)设xm，ym为脉冲函数的位置坐标，则r(x, y) =Σ δ (x-xm, y-ym)(4—2)根据卷积定理，在频域有T (U，ν) = H (U，ν) R (U，ν)(4-3)R (U，ν) = T (U，ν) H (U，νΓ1(4-4)那么对于给定纹理基元h(x，y)的描述，即可得到反卷积滤波器H (u，ν) Λ将这个 滤波器用于待处理纹理图像，每个脉冲都在纹理基元的中心。对整幅图像中不同纹理基元 的分布进行统计可获得图像的全面纹理信息。通过纹理基元确定对了对待分割的区域找到一个种子子图一即基元脉冲中心， 然后将种子像素周围邻域中与种子像素有相同或相似性质的像素(通过定义的相同或相 似规则)合并到种子像素所在区域中。在将新像素作为种子继续扩展，直到再没有满足条 件的像素可被包括，这样的一个区域就长成。对于一个给定的目标图像X和一个结构基元S,以符合S的种子区域j十S为中心 在图像上进行迭代生长。对于种子Z φ ^的每一个邻近像素X，在生长过程中应始终保持if S(x) = TRUE(4-5)沩十十(4-6)S卩，如果区域邻近像素χ能符合结构基元定义，则将像素χ合并到特征区域，使之 成为新种子的一部分。如果对于种子Z Θ S的所有邻近像素X，均有S (χ) = FALSE(4-7)则生长停止。对于阳光放射沟特征，通过前期对临床病理的研究分析，从获得的先验统计知识 定义像素生长规则如下。规则1 沟长度为L个像素，5彡L彡R虹膜；规则2 沟宽度为W个像素，3彡W彡5 ；规则3 沟的延长线与瞳孔交点的切线相垂直；规则4 根据虹膜图像灰度直方图确定沟区域灰度最优阈值Te ；规则5 沟的起点为卷缩轮边界点。
结合图14，根据CADIC虹膜分区图，对每个感兴趣的诊断候选区域，执行如下阳光 放射沟检测算法。1)计算阳光放射沟灰度阈值对整个虹膜统计灰度直方图，定义灰度阈值Te，使得灰度小于Te的像素数目为整 个虹膜像素数目的10%。根据不同的虹膜图像采集设备，此比例数应根据人工标注阳光放 射沟情况的先验统计数据做出相应的调整。2)对任一个候选区Q，执行以下3) 10)。3)沟起点选取对于任意像素χ e Q，如果有符合规则5的像素χ的灰度值f (x) ( Tg,则认为χ是 一个沟起点。把找到的所有沟起点记为dn，η = 1,2,3...N04)如果N = 0，即没有找到任何沟起点，则说明没有阳光放射沟存在，转10)结束 本算法。5)沟起点筛选与合并，形成种子候选对沟底点记dn，η = 1,2,3... N，先后执行生长与合并操作，合并成种子候选d’ n， η = 1，2，3· · · Ν，，N，彡 N。6)种子候选筛选，形成种子对每个种子候选d’ 定义它在平行于虹膜内边缘方向上的像素跨度数为其宽度 W(d’。，删除所有奴^ n) > 5的种子候选d’ n。所有剩余的种子候选(!^成为种子，记为zm，m= 1，2，3...M。在不造成混淆的情 况下，也用Zm代表本种子的中心像素。7)如果M = 0，即没有找到任何种子，则说明没有阳光放射沟存在，转10)结束本算法。8)区域生长，形成阳光放射沟候选在每个种子zm，m = 1,2,3. .. M周围邻域中搜索同种子像素具有相同或相似性质 的像素，合并到种子像素所在区域中；在将新像素作为种子继续迭代扩展，直到再没有满足 条件的像素可被包括，这样一个区域就长成了。(8. 1)定义每个种子(中心)Zm的邻域为NB(Zm)，NB(Zm)应符合规则1，2和3 ；(8. 2)定义基元结构为S,此处的S定义为典型阳光放射沟的横截线上的像素集
I=I O(8. 3)对于任意像素χ e NB (Zm)，根据公式(4_5)、(4-6)所定义的生长规则进行 区域生长。9)阳光放射沟筛选对获得的阳光放射沟候选，根据规则1和2进行筛选。10)输出结果，结束算法。通过上述方法实验，可以获得满足裂缝基元定义的线条，即为阳光放射沟。(6)虹膜动态自适应分区具体做法体现在对病理特征(此处为阳光放射沟，也可以推广到其他虹膜特征 上)在分区边界上的处理当分区时病理特征被分割开来，不能被完整地表达在同一个区 域中时，进一步从保持病理特征完整性的角度，重新确定该特征所覆盖区域。
6. 1)首先，通过估计病理特征的区域进行病理特征筛选对病理特征边缘形状的提取采用估算的方法。求出组成特征的像素点在极坐标中 的角度和极径的最大和最小值，这样就大概确定出了纹理所覆盖的区域的面积。从前面特 征提取的结果可以看到，只有一部分重要的纹理特征(阳光放射沟)具有诊病价值。我们 根据所确定的每个特征区域的面积，选取最重要的若干个特征，比如5个；删除其余的不重 要特征。6. 2)结合图15，进行CADIC图谱覆盖，方法如前述。6. 3)判断是否有跨区域的病理特征对每一个虹膜图谱上的分区注意考察，如果发现其边界同重要特征有交叠，则记 录这个交叠。6. 4)进行区域动态调整，完成自适应分区如果一个纹理特征覆盖多个区域，进行如下处理对每个虹膜图谱分区，统计区域 外内特征的像素点数目M;之后对其每处交叠边界，考察此区域外该特征的像素点并进行 统计数目M’。如果区域外像素点中只有极少部分M’ << M(例如M’小于5%的M)，则对这些少 数像素点所分布的区域不予考虑；如果区域外像素点过于重要M’>M(说明该特征不是主要处于本区域)，则同样不 予考虑，并删除本区域中的特征；此外，说明区域外特征比较重要；则将此边界邻接的区域同本区域合并，实现了重 新的区域划分，也就是进行了自适应分区。例如在图16中跨区纹理特征进行自适应分区在(1-1-2)、(1-5-2、1-5-3)与 (1-7-2、1-7-3、1-7-4)区域中出现长的连续纹理特征一阳光放射沟，(1-1-2)区中特征没 有出现跨区域现象；而在5区与7区中，把两个连续特征分割在几个小的区域中，使诊断很 困难。应用上述发明技术，使图谱在覆盖虹膜同时，对检测到的完整特征实现自适应分区， 确保特征信息的完整性，有利于诊病判断。
权利要求
一种基于基元结构定义与区域生长技术的虹膜分区及阳光放射沟提取方法，其特征在于在它是在基于串行区域生长与合并技术基础上，通过对特征先验知识，定义区域生长与合并准则，以基元为区域生长点，具体地方案包含如下步骤(1)虹膜图像的采集；(2)虹膜图像的预处理；(3)虹膜定位；(4.)CADIC虹膜图谱覆盖；(5)基于基元定义和区域生长的阳光放射沟提取；(6)虹膜自适应分区。
2.根据权利要求1所述的基于基元结构定义与区域生长技术的虹膜分区及阳光放射 沟提取方法，其特征在于虹膜图谱定义为对于一幅虹膜图像，经过前期预处理与虹膜内外 边缘精确定位后，将虹膜图像置于极坐标系(9，r)中，把虹膜分为以瞳孔圆心为极点、以 瞳孔边缘为起始边界的五个环、12个区的虹膜分区图谱，设一幅虹膜图像用A表示，Rt为瞳 孔半径，民为虹膜半径，RD为图像中卷缩轮的理论半径，R'。为卷缩轮的实际半径。极坐 标的极轴方向取为时钟的零点方向，极角按照顺时针方向定义，将A划分为五个环12个区 后，其中的任意一个区用A、9，r)来表示。其中:k = 1,2, ... ,12 ；i = 1，2，…，5，每个 分区的定义如下
3.根据权利要求2所述的基于基元结构定义与区域生长技术的虹膜分区及阳光放射 沟提取方法，其特征在于虹膜图像的采集方法是对视频图像抓拍，采样者头部固定并正视， 采样时所采用的光源产生的光斑打在瞳孔和虹膜的交界处，采样图像为768*576的24位位 图，每个采样者采3到5组虹膜图像。
4.根据权利要求3所述的基于基元结构定义与区域生长技术的虹膜分区及阳光放射 沟提取方法，其特征在于虹膜图像的预处理方法是在对虹膜图像进行分区分析之前进行预 处理工作，虹膜图像的预处理方法包括彩色图像灰度转换、噪声处理、几何校正、图像增强 方法，图像灰度转换方法为对采集到的彩色图像，应用如下公式将其变换为灰度图像 Y = 0. 299XR+0. 587XG+0. 114XB 式中R——真彩图中每个像素的红分量； G——真彩图中每个像素的绿分量； B——真彩图中每个像素的蓝分量； Y——该像素所对应的灰度；再进行虹膜图像灰度归一化，通过将原始虹膜图像中感兴趣的区域映射到相同的尺寸 和相同的对应区域中进行虹膜归一化，然后进行光斑填充，采用数字图像的自适应插值法 来对光斑进行填充。零阶插值算法，即令输出像素的灰度值等于离它所映射的位置最近的输入像素的灰度值。
5.根据权利要求4所述的基于基元结构定义与区域生长技术的虹膜分区及阳光放射 沟提取方法，其特征在于虹膜定位方法是除去那些于虹膜无关的信息，也就是从虹膜图像 中定位出虹膜的内边界和外边界。包括(1)虹膜粗定位采用灰度投影法对虹膜内、外边缘 粗定位，得到虹膜内径r、外径rl和瞳孔中心(x。，yc) ； (2)虹膜精定位采用Daugman算法 来对虹膜进行精定位，利用一个有效的积分微分算子来计算圆参数，采用圆形探测算子对 虹膜边界进行精定位，在虹膜粗定位瞳孔中心(x。，yc)和半径r的基础上，采用圆探测算子 实现虹膜精定位；所述虹膜精定位方法为基于圆探测算子的定位方法a、对于中心矩阵内 的每一个元素，分别计算半径为r同一圆周上的像素点灰度和sum，其中中心矩阵的范围是 (xc-5, y。-5)，半径变化范围为r-5 r+5 ；b、对于r的变化值d，计算出对应的高斯函数值g⑷=(/(0.\x (V2^)))e(("(d"6)2)/(2X°'l2)); c、计算灰度梯度值向量det = | g (d+1) -g (d) | X sum ；d、依次计算所有元素对应的det向量值；e、在所有的det向量中取元素最大值，即为精定位 的结果。
6.根据权利要求5所述的基于基元结构定义与区域生长技术的虹膜分区及阳光放射 沟提取方法，其特征在于CADIC虹膜图谱覆盖方法是在虹膜内外边缘精确定位后，将虹膜 部分所在区域的直角坐标(x，y)形式转换为极坐标(r，0)形式，根据虹膜图谱和r、0的 值来将虹膜部分分成若干个区域。
7.根据权利要求6所述的基于基元结构定义与区域生长技术的虹膜分区及阳光放射 沟提取方法，其特征在于基于基元定义和区域生长的阳光放射沟提取方法是选取虹膜图像 中疾病形成的阳光放射沟为处理对象，首先定义纹理基元并讨论基于基元的纹理表达方 法设纹理基元为h(x，y)，排列规则为r(x，y)，则纹理t(x，y)定义为 设xm，yffl为脉冲函数的位置坐标，则 r (x, y) = E 6 (x-xm, y-ym) 根据卷积定理，在频域有 对于给定纹理基元h(x，y)的描述，可得到反卷积滤波器H(u，V)—1，将这个滤波器用于 待处理纹理图像，每个脉冲都在纹理基元的中心，对整幅图像中不同纹理基元的分布进行 统计可获得图像的全面纹理信息；通过纹理基元确定对了对待分割的区域找到一个种子子图一即基元脉冲中心，然后 将种子像素周围邻域中与种子像素有相同或相似性质的像素合并到种子像素所在区域中； 在将新像素作为种子继续扩展，直到再没有满足条件的像素可被包括，这样的一个区域就 长成；对于一个给定的目标图像X和一个结构基元S,以符合S的种子区域Z十S为中心在图 像上进行迭代生长；对于种子的每一个邻近像素x，在生长过程中应始终保持 if :S(x) = TRUEthen: X S <r- X S\j{x)即，如果区域邻近像素X能符合结构基元定义，则将像素X合并到特征区域，使之成为 新种子的一部分；如果对于种子Z ^的所有邻近像素X，均有 S (x) = FALSE 则生长停止；对于阳光放射沟特征，通过前期对临床病理的研究分析，从获得的先验统计知识定义 像素生长规则如下规则1 沟长度为L个像素，5彡L彡R虹膜;规则2 沟宽度为W个像素，3^ff^5；规则3 沟的延长线与瞳孔交点的切线相垂直；规则4 根据虹膜图像灰度直方图确定沟区域灰度最优阈值Te ；规则5 沟的起点为卷缩轮边界点。
8.根据权利要求7所述的基于基元结构定义与区域生长技术的虹膜分区及阳光放射 沟提取方法，其特征在于根据CADIC虹膜分区图，对每个感兴趣的诊断候选区域，执行如下 阳光放射沟检测算法1)计算阳光放射沟灰度阈值对整个虹膜统计灰度直方图，定义灰度阈值L，使得灰度小于L的像素数目为整个虹 膜像素数目的10% ；根据不同的虹膜图像采集设备，此比例数应根据人工标注阳光放射沟 情况的先验统计数据做出相应的调整；2)对任一个候选区Q，执行以下3) 10)；3)沟起点选取对于任意像素xG Q，如果有符合规则5的像素x的灰度值f (x)彡Tg, 则认为x是一个沟起点。把找到的所有沟起点记为dn，n = 1，2，3. . . N ；4)如果N= 0，即没有找到任何沟起点，则说明没有阳光放射沟存在，转10)结束本算法；5)沟起点筛选与合并，形成种子候选对沟底点记dn，n= 1，2，3. . . N，先后执行生长与 合并操作，合并成种子候选d’n，n= 1，2，3...N’，N’彡N;6)种子候选筛选，形成种子对每个种子候选d’n定义它在平行于虹膜内边缘方向上 的像素跨度数为其宽度W(d’ n)，删除所有W(d’ n) > 5的种子候选d’ n;所有剩余的种子候 选d’ 成为种子，记为zm，m= 1，2，3...M。在不造成混淆的情况下，也用zm代表本种子的 中心像素；7)如果M= 0，即没有找到任何种子，则说明没有阳光放射沟存在，转10)结束本算法；8)区域生长，形成阳光放射沟候选在每个种子zm，m=1，2，3... M周围邻域中搜索同 种子像素具有相同或相似性质的像素，合并到种子像素所在区域中；在将新像素作为种子 继续迭代扩展，直到再没有满足条件的像素可被包括，这样一个区域就长成了 ；8. 1)定义每个种子(中心)的邻域为NB(zm)，NB(zffl)应符合规则1，2和3 ； 8. 2)定义基元结构为S，此处的S定义为典型阳光放射沟的横截线上的像素集合； 8. 3)对于任意像素x G NB(zm)，根据上述定义的生长规则进行区域生长；9)阳光放射沟筛选对获得的阳光放射沟候选，根据规则1和2进行筛选；10)输出结果，结束算法；通过上述方法可得满足裂缝基元定义的线条，即为阳光放射沟。
9.根据权利要求8所述的基于基元结构定义与区域生长技术的虹膜分区及阳光放射 沟提取方法，其特征在于虹膜动态自适应分区方法是在对病理特征在分区边界上的处理 当分区时病理特征被分割开来，不能被完整地表达在同一个区域中时，进一步从保持病理 特征完整性的角度，重新确定该特征所覆盖区域1)首先，通过估计病理特征的区域进行病理特征筛选对病理特征边缘形状的提取采用估算的方法，求出组成特征的像素点在极坐标中的角 度和极径的最大和最小值，确定出了纹理所覆盖的区域的面积，从前面特征提取的结果可 以看到，只有一部分重要的纹理特征阳光放射沟具有诊病价值，根据所确定的每个特征区 域的面积，选取重要特征，删除其余的不重要特征；2)进行CADIC图谱覆盖，方法如前述；3)判断是否有跨区域的病理特征对每一个虹膜图谱上的分区注意考察，如果发现其 边界同重要特征有交叠，则记录这个交叠；4)进行区域动态调整，完成自适应分区如果一个纹理特征覆盖多个区域，进行如下 处理对每个虹膜图谱分区，统计区域外内特征的像素点数目M ；之后对其每处交叠边界， 考察此区域外该特征的像素点并进行统计数目M’ ；如果区域外像素点中只有极少部分M’ << M，则对这些少数像素点所分布的区域不予 考虑；如果区域外像素点过于重要M’ >M，则同样不予考虑，并删除本区域中的特征；此外，说明区域外特征比较重要；则将此边界邻接的区域同本区域合并，实现重新的区 域划分，也就是进行了自适应分区。
全文摘要
本发明提出了一种基于基元结构定义与区域生长技术的虹膜分区及阳光放射沟提取方法，包括一种融合中西医虹膜学虹膜分区的新虹膜分区图谱(CADIC)；并发明了应用此图谱所必需的基于基元定义的区域生长技术用于典型虹膜特征——阳光放射沟的提取，并发明了应用此图谱所必需的基于虹膜图像动态分区技术和虹膜图谱的自适应覆盖方法。本发明在保证了肌体各部脏器在虹膜上的准确定位同时，突出了虹膜图谱的实用性，实现了医学意义下完整虹膜特征的动态分割技术及其在虹膜图谱上的定位技术。通过实验证明，此技术能够大大提高了虹膜识别和诊病可靠性。
文档编号G06K9/34GK101882222SQ20091007237
公开日2010年11月10日 申请日期2009年6月26日 优先权日2009年6月26日
发明者张大鹏, 李乃民, 王宽全, 王钧慧, 马琳 申请人:哈尔滨工业大学