一种基于不变矩的图形畸变分析方法技术

一种基于不变矩的图形畸变分析方法，包括：１）选择合适的矩不变量：在动态形体模型中，使用点分布模型表示图像中的对象轮廓；２）模型畸变的主动检测方法：２．１）对动态形体模型中的ＰＤＭ，其中，每一个标注点都有相应的权重，根据权重将一个ＰＤＭ分解为子图形；２．２）计算训练集中每个子图形和模型在设定范围内进行变形后每一个子图形的边界矩不变量；然后计算这些边界矩不变量的聚类中心值；２．３）进行模型的畸变检测的步骤为：①将变形后的模型按照预先定义的分块策略进行分块；②计算每个子图形的边界矩不变量；③根据式（２）（３）（４）计算每个子图形的畸变概率；３）模型畸变的主动修正方法。本发明专利技术在搜索过程中，能有效检测并修正畸变，可靠性好。

【技术实现步骤摘要】

本专利技术涉及图形图像领域，尤其是一种图形畸变分析方法，该方法结合使用 动态形体模型方法、图像处理、计算机图形学、医学、数学、应用物理、统计学、 生物学等技术，主要用于分析图形拟合质量，尤其是医学图像中的形状分割和形 态畸变分析。技术背景可变形模型是20世纪80年代后期发展起来的一种有效图像分割、配准与识 别技术。可变形模型拥有形体变化能力，它能代表柔性物体。将可变形模型应用 于图像，可以通过变形寻找图像上待分割、配准或识别的物体。可变形模型方法 一般包含两个阶段，即训练阶段和实际匹配阶段。一般可变形模型的分类有两种方式， 一种是以参数化方式来分类，包括局部 参数化和全局参数化两种，局部参数模型又被称为分布式参数模型，是一种基于 物理变形的模型，全局参数模型又被称为参数模型，它一般是基于先验知识的模 型；另一种分类方式是按照是否是基于描述还是基于原型来进行区分的，基于描 述的模型的一个重要特征是它是通过抽样出形体来代表实际的物体，基于描述的 方式还可以分成固定的拓扑结构和自适应的拓扑结构两种，而基于原型的模型一 般是要建立一个简单的原始模型来进行分割。动态形体模型(Active Shape Model)是可变形模型中一种重要的方法，由英 国学者Cootesl995年首次提出，并用于医学图像处理和人脸识别中。动态形体模 型非常适合于形状比较复杂图像的分割，但是它也存在一些问题，首先，由于动 态形体模型搜索范围不是很广，如果初始搜索的位置不是很理想，就可能导致失 败的结果，其次，由于训练的过程需要专家的经验和手动标定对象的形状，势必要求进行多次训练来获取最佳的模型，但这样是很花时间和精力。因此很多人提 出了对动态形体模型的改进，包括Hicks提出了一种针对生物形体的自动标注的 方式。高原提出将主动轮廓模型的部分内容与统计模型相集成的方法，提出了一 种基于多层混合模型的图像分割方法。Rhodri提出了对动态形体模型的参数的一 系列优化。Hamameh将2D的动态形体模型和时间相结合，提出了分时形体模型 (TS-shapes)。蒋漫设计了两种基于Fisher变化的局部外观模型，提高了动态形体 模型在图像配准中的鲁棒性。胡海舟应用动态形体模型确定红眼区域中的瞳孔轮 廓，提出了一个自动检测和取出照片红眼的算法。除了对现有问题的一些改进和应用，动态形体模型现在最常见的有以下一些 改进在模型建立阶段，通常采用设置权重的方法；在搜索阶段，需要追求搜索 速度和准确度的平衡，通常采用图像分层的方法，这个方法首先将图像的像素降 低，然后用一个对应的模型來搜索，分割结束后提高像素继续搜索，直到恢复到 原来的像素，最终实现目标的准确和快速分割。动态形体模型通过训练的方法得到对象的统计知识，从而建立对象的统计模 型，其中包括对象的点分布模型X = f+加上式中，X表示新的形状，^表示从平均图形变换得到的初始形状，p则是经过主成分分析得到的主成分，矢量6则是每个主成分的权重。新的形状与初始 模型的相似程度完全由6决定。因此，为了对模型的变形进行约束以避免过度变形，将6限制在一个由训练集确定的范围之内。但是过于严格的限制减弱模型的变形能力，不能反映出实际对象的变化。现有的动态形体模型方法在搜索过程中，不能检测并修正畸变，可靠性差。
技术实现思路
为了克服已有动态形体模型方法在搜索过程中，不能检测并修正畸变，可靠性差的不足，本专利技术提供一种在搜索过程中，能有效检测并修正畸变，可靠性好 的基于不变矩的图形畸变分析方法。本专利技术解决其技术问题所采用的技术方案是，所述分析方法包括-1) 、选择合适的矩不变量在动态形体模型中，使用点分布模型表示图像中的对象轮廓，选择以下六个等式作为描述对象轮廓的矩不变量0 =(73。 -3 712)(%3 +72I) + (3"21 -%3)0721 +%3)XP073。 + 712)2 +7。3)2]豸=3^」—屮^)2 —+%)2W3/fc+"B)XR% +^)"《％(4 = ("40 —6"22 —7M)2 —《％ +&)2]+16(%<, —7m)('732, _《)(4 =(%—%)2 +化)2]——+%X% —6% +%)疾=%)M)4 _2线％ -" -" -《(1)上式中，《力=1...6)表示六个矩不变量，7,表示第，^阶中心矩；2) 、模型畸变的主动检测方法2.1) 、对动态形体模型中的点分布模型PDM，其中，每一个标注点都有相应 的权重，根据权重将一个PDM分解为子图形；每一个子图形都是由权重相近的 标注点组成，每个子图形至少包含三个标注点，每个子图形的标注点不超过6个；2.2) 、计算训练集中每个子图形和模型在设定范围内进行变形后每一个子图 形的边界矩不变量;然后计算这些边界矩不变量的聚类中心值；假设子图形的边 界矩不变量在每个聚类中心附近符合高斯分布，则每个聚类中心是高斯分布的均 值/V然后计算高斯分布的方差 ；其中/表示第/个子图形，j'表示第/个子图形 的第J'个聚类中心；对于任意一个子图形的边界矩不变量，首先采用式(2)来计算其相对于该子图 形所有聚类形成的高斯分布的概率<formula>formula see original document page 8</formula>然后，定义这个子图形相对每一个聚类的畸变概率如下 A,(X) = 1-W，力丫i，") (3)上式中，m是子图形的个数，"是每个子图形的聚类数目； 最后取";/中最小的值作为该子图形最终的畸变概率D,(X) = mi,W) (4)2.3)、进行模型的畸变检测的步骤为① 、将变形后的模型按照预先定义的分块策略进行分块；② 、计算每个子图形的边界矩不变量；③ 、根据式(2)(3)(4)计算每个子图形的畸变概率。作为优选的一种方案所述分析方法还包括3)、模型畸变的主动修正方法， 包括以下步骤3.1) 、修正策略的触发设第/次迭代时的第/个子图形的边界矩不变量为/i^，定义用于判断图形是 否发生突变的边界矩不变量的差分如下,=max—J+l -,)-气—!I)/H^,川-气|权-'l) (5) 当tffi^大于设定阈值时，启动修正策略；3.2) 、修正畸变模型① 、计算第/个子图形的畸变概率A;② 、按照式(6)计算畸变子图形区域中每个标注点相对于上一次迭代结果的距上式中，；，r力,是当前图形中第y个标注点的坐标，^^和^是上一次迭代结果中对应的第y'个标注点的坐标；③、然后根据对子图形中的每一点按照式(7)计算标注点的新的坐标；<formula>formula see original document page 9</formula>上式中，《"和^"是修正后的标注点坐标，并计算修正距离《"如下<formula>formula see original document page 9</formula>(8) 、对模型的其他点做以下修正<formula>formula see original document page 9</formula>上式中，不。和i;。是模型中第o个标注点在当前本文档来自技高网...

【技术保护点】
一种基于不变矩的图形畸变分析方法，其特征在于：所述分析方法包括：　　　　１）、选择合适的矩不变量：在动态形体模型中，使用点分布模型表示图像中的对象轮廓，选择以下六个等式作为描述对象轮廓的矩不变量：　　　　φ↓［１］＝（η↓［３０］－３η↓［１２］）（η↓［０３］＋η↓［２１］）［（η↓［３０］＋η↓［１２］）↑［２］－３（η↓［０３］＋η↓［２１］）↑［２］］＋（３η↓［２１］－η↓［０３］）（η↓［２１］＋η↓［０３］）×［３（η↓［３０］＋η↓［１２］）↑［２］－（η↓［２１］＋η↓［０３］）↑［２］］　　　　φ↓［３］＝３（η↓［２１］－η↓［０３］）（η↓［３０］＋η↓［１２］）［（η↓［３０］＋η↓［１２］）↑［２］－３（η↓［０３］＋η↓［２１］）↑［２］］＋（３η↓［１２］－η↓［３０］）（η↓［２１］＋η↓［０３］）×［３（η↓［３０］＋η↓［１２］）↑［２］－（η↓［２１］＋η↓［０３］）↑［２］］　　　　φ↓［４］＝（η↓［４０］－６η↓［２２］＋η↓［０４］）［（η↓［４０］－η↓［０４］）↑［２］－４（η↓［３１］＋η↓［１３］）↑［２］］＋１６（η↓［４０］－η↓［０４］）（η↓［３１］↑［２］－η↓［１３］↑［２］）　　　　φ↓［５］＝（η↓［３１］－η↓［１３］）［（η↓［４０］－η↓［０４］）↑［２］－４（η↓［３１］＋η↓［１３］）↑［２］］－（η↓［４０］－η↓［０４］）（η↓［３１］＋η↓［１３］）（η↓［４０］－６η↓［２２］＋η↓［０４］）　　　　φ↓［６］＝η↓［４０］η↓［２２］η↓［０４］－２η↓［３１］η↓［２２］η↓［１３］－η↓［４０］η↓［１３］↑［２］－η↓［０４］η↓［３１］↑［２］－η↓［２２］↑［３］　　（１）　　　　上式中，φ↓［ｉ］（ｉ＝１…６）表示六个矩不变量，η↓［ｊｋ］表示第ｊ＋ｋ阶中心矩；　　　　２）、模型畸变的主动检测方法：　　　　２．１）、对动态形体模型中的点分布模型ＰＤＭ，其中，每一个标注点都有相应的权重，根据权重将一个ＰＤＭ分解为子图形；每一个子图形都是由权重相近的标注点组成，每个子图形至少包含三个标注点，每个子图形的标注点不超过６个；　　　　２．２）、计算训练集中每个子图形和模型在设定范围内进行变形后每一个子图形的边界矩不变量；然后计算这些边界矩不变量的聚类中心值；假设子图形的边界矩不变量在每个聚类中心附近符合高斯分布，则每个...

【技术特征摘要】

【专利技术属性】
技术研发人员：陈胜勇，张剑华，许艺强，陈敏，刘盛，
申请(专利权)人：浙江工业大学，
类型：发明
国别省市：86[中国|杭州]