发明名称 一种基于不变矩的图形畸变分析方法
摘要 一种基于不变矩的图形畸变分析方法,包括:1)选择合适的矩不变量:在动态形体模型中,使用点分布模型表示图像中的对象轮廓;2)模型畸变的主动检测方法:2.1)对动态形体模型中的PDM,其中,每一个标注点都有相应的权重,根据权重将一个PDM分解为子图形;2.2)计算训练集中每个子图形和模型在设定范围内进行变形后每一个子图形的边界矩不变量;然后计算这些边界矩不变量的聚类中心值;2.3)进行模型的畸变检测的步骤为:①将变形后的模型按照预先定义的分块策略进行分块;②计算每个子图形的边界矩不变量;③根据式(2)(3)(4)计算每个子图形的畸变概率;3)模型畸变的主动修正方法。本发明在搜索过程中,能有效检测并修正畸变,可靠性好。
申请公布号 CN101256627B 申请公布日期 2011.12.28
申请号 CN200810059520.4 申请日期 2008.01.25
申请人 浙江工业大学 发明人 陈胜勇;张剑华;许艺强;陈敏;刘盛
分类号 G06K9/00(2006.01)I;G06K9/64(2006.01)I;G06T11/00(2006.01)I 主分类号 G06K9/00(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;王利强
主权项 一种基于不变矩的图形畸变分析方法,其特征在于:所述分析方法包括:1)、选择合适的矩不变量:在动态形体模型中,使用点分布模型表示图像中的对象轮廓,选择以下五个等式作为描述对象轮廓的矩不变量:φ1=(η30‑3η12)(η03+η21)[(η30+η12)2‑3(η03+η21)2]+(3η21‑η03)(η21+η03)×[3(η30+η12)2‑(η21+η03)2]φ3=3(η21‑η03)(η30+η12)[(η30+η12)2‑3(η03+η21)2]+(3η12‑η30)(η21+η03)×[3(η30+η12)2‑(η21+η03)2] <mrow> <msub> <mi>&phi;</mi> <mn>4</mn> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>&eta;</mi> <mn>40</mn> </msub> <mo>-</mo> <mn>6</mn> <msub> <mi>&eta;</mi> <mn>22</mn> </msub> <mo>+</mo> <msub> <mi>&eta;</mi> <mn>04</mn> </msub> <mo>)</mo> </mrow> <mo>[</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>&eta;</mi> <mn>40</mn> </msub> <mo>-</mo> <msub> <mi>&eta;</mi> <mn>04</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mn>4</mn> <msup> <mrow> <mo>(</mo> <msub> <mi>&eta;</mi> <mn>31</mn> </msub> <mo>+</mo> <msub> <mi>&eta;</mi> <mn>13</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>]</mo> <mo>+</mo> <mn>16</mn> <mrow> <mo>(</mo> <msub> <mi>&eta;</mi> <mn>40</mn> </msub> <mo>-</mo> <msub> <mi>&eta;</mi> <mn>40</mn> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msubsup> <mi>&eta;</mi> <mn>31</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>&eta;</mi> <mn>13</mn> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> </mrow>φ5=(η31‑η13)[(η40‑η04)2‑4(η31+η13)2]‑(η40‑η04)(η31+η13)(η40‑6η22+η04) <mrow> <msub> <mi>&phi;</mi> <mn>6</mn> </msub> <mo>=</mo> <msub> <mi>&eta;</mi> <mn>40</mn> </msub> <msub> <mi>&eta;</mi> <mn>22</mn> </msub> <msub> <mi>&eta;</mi> <mn>04</mn> </msub> <mo>-</mo> <mn>2</mn> <msub> <mi>&eta;</mi> <mn>31</mn> </msub> <msub> <mi>&eta;</mi> <mn>22</mn> </msub> <msub> <mi>&eta;</mi> <mn>13</mn> </msub> <mo>-</mo> <msub> <mi>&eta;</mi> <mn>40</mn> </msub> <msubsup> <mi>&eta;</mi> <mn>13</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msub> <mi>&eta;</mi> <mn>04</mn> </msub> <msubsup> <mi>&eta;</mi> <mn>31</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>&eta;</mi> <mn>22</mn> <mn>3</mn> </msubsup> </mrow>                                                                   (1)上式中,φi(i=1...6,且i≠2)表示五个矩不变量,ηjk表示第j+k阶中心矩;2)、模型畸变的主动检测方法:2.1)、对动态形体模型中的点分布模型PDM,其中,每一个标注点都有相应的权重,根据权重将一个PDM分解为子图形;每一个子图形都是由权重相近的标注点组成,每个子图形至少包含三个标注点,每个子图形的标注点不超过6个;采用c‑均值聚类算法将所有标注点的权重进行聚类,保证每一类的标注点的数目在3到6之间;2.2)、计算训练集中每个子图形和模型在设定范围内进行变形后每一个子图形的边界矩不变量;然后计算这些边界矩不变量的聚类中心值;假设子图形的边界矩不变量在每个聚类中心附近符合高斯分布,则每个聚类中心是高斯分布的均值μij,然后计算高斯分布的方差σij;其中i表示第i个子图形,j表示第i个子图形的第j个聚类中心;对于任意一个子图形的边界矩不变量,首先采用式(2)来计算其相对于该子图形所有聚类形成的高斯分布的概率: <mrow> <msub> <mi>f</mi> <mi>ij</mi> </msub> <mrow> <mo>(</mo> <mi>X</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <msub> <mi>&sigma;</mi> <mi>ij</mi> </msub> <msqrt> <mn>2</mn> <mi>&pi;</mi> </msqrt> </mrow> </mfrac> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mi>X</mi> <mo>-</mo> <msub> <mi>&mu;</mi> <mi>ij</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>/</mo> <mrow> <mo>(</mo> <mn>2</mn> <msubsup> <mi>&sigma;</mi> <mi>ij</mi> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>然后,定义这个子图形相对每一个聚类的畸变概率如下:Dij(X)=1‑fij(X)i(1,m),j(1,n)    (3)上式中,m是子图形的个数,n是每个子图形的聚类数目;最后取Dij中最小的值作为该子图形最终的畸变概率:Di(X)=min(Dij(X))i(1,m),j(1,n)  (4)2.3)、进行模型的畸变检测的步骤为:①、将变形后的模型按照预先定义的分块策略进行分块;②、计算每个子图形的边界矩不变量;③、根据式(2)(3)(4)计算每个子图形的畸变概率。
地址 310014 浙江省杭州市下城区朝晖六区