专利名称：：基于变形垂足曲线的人脸轮廓提取方法
技术领域：
：本发明涉及一种计算机视觉
技术领域：
，特别涉及一种基于变形垂足曲线的人脸轮廓提取方法。
背景技术：
：人脸轮廓是一种描述人脸形状和结构的特征，作为人脸特征的重要组成部分，是后续其它人脸特征检测提取和人脸识别的基础，对人脸信息处理和分析有着重要的作用。由于人脸形状的多样性和复杂性，在使用刚性模型提取轮廓时遇到了很大困难，而常规的边缘检测算子得到的边缘是不连续的，曲线进化方法是20世纪80年代后期发展起来的一种图像分割方法，特别适用于建模和提取任意形状的变形轮廓，因此也必然适用于人脸这样一个高度非刚体目标的轮廓提取。水平集方法是求解曲线进化的一种新颖方法，它最初由Osher和Sethian提出，其基本思想是将平面闭合曲线隐含地表达为二维曲面函数的水平集，即具有相同函数值的点集，通过水平集函数曲面的进化隐含地求解曲线的运动。这种方法的数学基础是Crandall和Lions建立的Hamilton-Jacobi方程的粘性解理论，它使象水平集方程这样的几何偏微分方程，能够通过计算机模拟精确地找到Lipschitz连续的唯一解，从而可以避免陷入局部极值和对初始位置敏感的缺陷。此外，这种曲线进化方式的最大优点是，即使隐含在水平集函数中的闭合曲线发生了拓扑结构变化(合并或分裂)，水平集函数仍保持为一个有效的函数，即它可以自然地处理曲线的拓扑变化。而且，曲线内在的几何特征(如单位法向向量和曲率等)也可以利用水平集函数直接计算出来。由于水平集方法的良好特性，已经引起了越来越多的关注，在图像处理和计算机视觉的许多领域中得到了广泛的应用。垂足曲线是微分几何中的一种重要的数学曲线，它可以很方便地同水平集方法结合以实现曲线变形的能力，从而可以用于描述非刚体目标的形状，并且可以同时处理形状的全局变形和局部变形。
发明内容本发明是针对现有人脸轮廓提取困难的问题，提出了一种基于变形垂足曲线的人脸轮廓检测方法，借助于变形垂足曲线能同时表示全局形状和局部形变的能力，将人脸形状的先验知识嵌入到该模型中，同时利用水平集方法做数值求解，从而能自然地处理曲线的拓扑变化，而且该方法在检测出人脸轮廓的同时可以粗略估计出人脸的大小位置和水平旋转角度等信息。本发明的技术方案为一种基于变形垂足曲线的人脸轮廓提取方法，其特征在于，方法包括如下步骤1)人脸全局形状参数估计人脸的全局形状由垂足曲线的发生器控制，人脸椭圆全局形状参数{(Xo，y。)，0，a,b}利用垂足曲线的一阶和二阶矩来快速估计，其中(X(l，y0)为椭圆中心，a，b分别为椭圆的长短轴半径，水平旋转角度0，设垂足曲线用pe(t)表示，10而二咎式中则椭圆中心点坐标(XQ，yQ)可以通过pe(t)的重心估计而=为pe(t)的i+j阶原点矩；椭圆的倾斜角e可以由Pe(t)的中心矩估计》式中口。为Pe(t)的i+j阶中心矩；椭圆的长短轴a和b可以由p6(t)的最大、最小惯性矩估计’其中，最大惯性矩Imax和最小惯性矩Imin分别定义如下2)人脸局部形状变形人脸形状的局部变形通过垂足曲线的进化来实现，采用Chan-Vese模型来控制垂足曲线的进化，同时利用水平集方法做数值求解，垂足点位置由曲线P(t)控制，垂足曲线用pe(t)表示，用向量形式表示为巩其设垂足曲线P6(t)是二维曲面函数扒的零水平集，扒的进印2，式化方程釆用Chm-Vese改进模型，即0为固定常数，K为水平集函数曲率，参数计算如下'其中冲Ho…为Heaviside函数；垂足点控制曲线p(t)通过对二维曲面函数牝进行进化来求解，牝的进化方程通过P(t)%pe(t)的关系以及小2的进化方程推导得到牝的进化方程为具体实施例方式本发明提出的基于变形垂足曲线的人脸轮廓提取方法，包括以下步骤(1)人脸全局形状参数估计人脸的全局形状由垂足曲线的发生器控制。考虑到人脸本质上是一个具有一定形变的椭圆，即可以用椭圆曲线C*描述为其中，(x0,y0)为椭圆中心，a,b分别为椭圆的长短轴半径。由于图像中人脸的姿态未知，故引入水平旋转角度e，人脸椭圆模型如附图1所示。垂足曲线的发生器a(t)可以用椭圆曲线方程表示如下\(t)可表示为为简单起见，人脸椭圆全局形状参数{(&，％)，e，a，b}利用垂足曲线的一阶和二阶矩来快速估计。设垂足曲线用Pe(t)表示，则椭圆中心点坐标(X(l，y(l)可以通过pe(t)的重心估计式中为pe(t)的i+j阶原点矩。椭圆的倾斜角9可以由pe(t)的中心矩估计垂足曲线的发生器a(t)f因此a(t)的法线方向Ja(5)3)根据^的进化方程迭代计算U+1，由得到垂足点控制曲线P(t)，然后通过P(t)对给定平面参数曲线做垂足操作得到垂足曲线pe(t)，由pe(t)得到ctu/，可设置最大迭代次数，或两次迭代误差小于设定常数为迭代终止条件，满足终止条件后，垂足曲线Pe(t)即为最终得到的人脸轮廓提取结果。本发明的有益效果在于本发明基于变形垂足曲线的人脸轮廓检测方法，可以提取出噪声图像及平面旋转的人脸图像中的完整连续人脸轮廓，而且同时可以粗略估计出人脸的大小位置和水平旋转角度等信息。图1为本发明基于变形垂足曲线的人脸轮廓检测方法中人脸椭圆形状轴图；图2为本发明基于变形垂足曲线的人脸轮廓检测方法中垂足曲线示意图；图3为本发明基于变形垂足曲线的人脸轮廓检测方法中变形垂足曲线示意图；图4为本发明基于变形垂足曲线的人脸轮廓检测方法中径向相关策略示意图。图5为本发明p(t)、pe(t)、小”小2的关系图。式中ii为pe(t)的i+j阶中心矩。椭圆的长短轴a和b可以由p6(t)的最大、最小惯性矩估计其中，最大惯性矩Imax和最小惯性矩Imin分别定义如下(2)人脸局部形状变形人脸形状的局部变形通过垂足曲线的进化来实现，采用Chan-Vese模型来控制垂足曲线的进化。实际计算时，并不需要对垂足曲线直接进行进化，而是先根据垂足曲线和垂足点控制曲线之间的关系，推导得到垂足点控制曲线的进化方程，通过对该曲线进化后，再对其进行垂足操作得到垂足曲线。垂足曲线是微分几何中的一种重要曲线，给定平面参数曲线a(s)和定点p(x，y),sG为曲线参数，则点P在a(s)的动切线上的投影轨迹3(s)称为曲线a(s)关于定点P的垂足曲线，如附图2所示，这里曲线a(s)称为垂足曲线3(s)的发生器，p点称为垂足点。垂足曲线可以用于描述物体的形状，并且可以同时处理形状的全局变形和局部变形，其局部变形由垂足点的位置控制，全局变形则由发生器控制。当垂足点P动态变化即由一条曲线控制时，垂足曲线可以用于描述更多的物体形状，由此产生的垂足曲线称为变形垂足曲线，其示意图如附图3所示，其中1为变形垂足曲线，2为垂足点控制曲线，3为发生器。设垂足点位置由曲线p(t)控制，p(t)与发生器a(t)之间的对应关系采用径向相关法，即对发生器上的每一点a(Si)，该点与其中心o的连线与垂足点控制曲线的交点P(Si)即为对应的垂足点，如附图4所示，4为垂足点控制曲线，5为发生器。垂足曲线用Pe(t)表示，为使其具有更大的变形能力，垂足曲线的定义可做如下修改这种修改使垂足曲线进一步具有膨胀和收缩的能力当垂足点位于发生器内部时，垂足曲线具有收缩的能力；当垂足曲线位于发生器外部时，垂足曲线具有膨胀的能力。用向量形式表示为Pe=JAP"b设垂足曲线P。(t)是二维曲面函数中，的零水平集，中，的进化方程采用than—Vese改进模型，即其中为Heaviside函数。设垂足点控制曲线P(t)通过对二维曲面函数中、进行进化来求解，中、的进化方程通过P(七)与P。(t)的关系以及中，的进化方程推导得到。实际计算时，并不需要对中，进行直接进化来求解P。(t)，而是通过对中、进化求解P(t)，获得P(t)后，再对P(t)进行垂足操作来得到P。(t)。P(t)、P。(t)、中、、中，的关系如附图5所示。根据P。(t)和P(t)的关系以及中，的进化方程可以推导出中、的进化方程如下式中由于上述方程的等号右边仅含有V中、和V中，，二者之间存在如下关系V中l一丁c’V中或V中一J信’V中l(12)由此可得中、的进化方程为(3)根据上述方程迭代计算中㈠、n+l，由中㈠、n／‘得到垂足点控制曲线P(t)，然后通过P(t)对o(t)做垂足操作得到垂足曲线P。(t)，由P。(t)得到中，』”。可设置最大迭代次数，或两次迭代误差小于某个很小的常数为迭代终止条件，满足终止条件后，垂足曲线P。(t)即为最终得到的人脸轮廓提取结果。实例(1)设迭代次数n—l，初始化垂足点控制曲线P(t)和水平集函数中㈠、”，给定初始椭圆发生器o(t)，由P(t)对o(t)做垂足操作得到垂足曲线P。(t)，根据P。(t)初始化水平集函数I,:。(2)根据式(13)迭代计算(^，un+1，由小，ijn+1得到垂足点控制曲线P(t)，然后通过p(t)对a(t)做垂足操作得到垂足曲线pe(t)，由pjt)得到小2,^。方程(13)采用有限差分方法来求解，其等号右边第一项为与曲率有关的项，须用中心差分近似；第二项为曲线在其法线方向上的常量演化，必须用逆风有限差分方法来近似，因此其最终数值计算公式如下式中其中Di/x，Di^和Di/x分别表示JdtV小！在x方向上的前向、后向和中心差分，Di/y，Di厂和Di/y分别表示JDT▽小i在y方向上的前向、后向和中心差分。中的参数Cl，c2根据小2，J计算得到。设Jdt可表示为牝可写为因此Ke的计算公式如下式中各阶偏导用相应的中心差分来逼近。(3)若满足迭代终止条件则停止，否则根据垂足曲线pe(t)更新椭圆发生器参数，具体计算公式见⑷(7)。令n=n+l，转⑵。权利要求一种基于变形垂足曲线的人脸轮廓提取方法，其特征在于，方法包括如下步骤1)人脸全局形状参数估计人脸的全局形状由垂足曲线的发生器控制，人脸椭圆全局形状参数{(x0，y0)，θ，a，b}利用垂足曲线的一阶和二阶矩来快速估计，其中(x0，y0)为椭圆中心，a，b分别为椭圆的长短轴半径，水平旋转角度θ，设垂足曲线用pe(t)表示，则椭圆中心点坐标(x0，y0)可以通过pe(t)的重心估计<mrow><msub><mi>x</mi><mn>0</mn></msub><mo>=</mo><mfrac><msub><mi>M</mi><mn>10</mn></msub><msub><mi>M</mi><mn>00</mn></msub></mfrac><mo>,</mo></mrow><mrow><msub><mi>y</mi><mn>0</mn></msub><mo>=</mo><mfrac><msub><mi>M</mi><mn>01</mn></msub><msub><mi>M</mi><mn>00</mn></msub></mfrac></mrow>式中Mij为pe(t)的i+j阶原点矩；椭圆的倾斜角θ可以由pe(t)的中心矩估计<mrow><mi>θ</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>arctan</mi><mrow><mo>(</mo><msub><mrow><mn>2</mn><mi>μ</mi></mrow><mn>11</mn></msub><mo>/</mo><mrow><mo>(</mo><msub><mi>μ</mi><mn>20</mn></msub><mo>-</mo><msub><mi>μ</mi><mn>02</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow></mrow>式中μij为pe(t)的i+j阶中心矩；椭圆的长短轴a和b可以由pe(t)的最大、最小惯性矩估计<mrow><mi>a</mi><mo>=</mo><msup><mrow><mo>(</mo><mn>4</mn><mo>/</mo><mi>π</mi><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>4</mn></mrow></msup><msup><mrow><mo>(</mo><msubsup><mi>I</mi><mi>max</mi><mn>3</mn></msubsup><mo>/</mo><msub><mi>I</mi><mi>min</mi></msub><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>,</mo></mrow><mrow><mi>b</mi><mo>=</mo><msup><mrow><mo>(</mo><mn>4</mn><mo>/</mo><mi>π</mi><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>4</mn></mrow></msup><msup><mrow><mo>(</mo><msubsup><mi>I</mi><mi>min</mi><mn>3</mn></msubsup><mo>/</mo><msub><mi>I</mi><mi>min</mi></msub><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><mo>,</mo></mrow>其中，最大惯性矩Imax和最小惯性矩Imin分别定义如下<mrow><msub><mi>I</mi><mi>min</mi></msub><mo>=</mo><munder><mi>Σ</mi><mrow><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>∈</mo><msub><mi>p</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></munder><msup><mrow><mo>[</mo><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>cos</mi><mi>θ</mi><mo>-</mo><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>sin</mi><mi>θ</mi><mo>]</mo></mrow><mn>2</mn></msup><mo>,</mo></mrow><mrow><msub><mi>I</mi><mi>max</mi></msub><mo>=</mo><munder><mi>Σ</mi><mrow><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>∈</mo><msub><mi>p</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow></munder><msup><mrow><mo>[</mo><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>sin</mi><mi>θ</mi><mo>-</mo><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>cos</mi><mi>θ</mi><mo>]</mo></mrow><mn>2</mn></msup><mo>;</mo></mrow>2)人脸局部形状变形人脸形状的局部变形通过垂足曲线的进化来实现，采用ChanVese模型来控制垂足曲线的进化，同时利用水平集方法做数值求解，垂足点位置由曲线p(t)控制，垂足曲线用pe(t)表示，用向量形式表示为pe＝JApb其中<mrow><msub><mi>J</mi><mi>A</mi></msub><mo>=</mo><mfencedopen='['close=']'><mtable><mtr><mtd><msub><mi>J</mi><mrow><mi>A</mi><mn>11</mn></mrow></msub></mtd><mtd><msub><mi>J</mi><mrow><mi>A</mi><mn>12</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>J</mi><mrow><mi>A</mi><mn>21</mn></mrow></msub></mtd><mtd><msub><mi>J</mi><mrow><mi>A</mi><mn>22</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfencedopen='['close=']'><mtable><mtr><mtd><mn>1</mn><mo>+</mo><mfrac><msup><msub><mi>n</mi><mn>1</mn></msub><mn>2</mn></msup><mrow><msup><msub><mi>n</mi><mn>1</mn></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>n</mi><mn>2</mn></msub><mn>2</mn></msup></mrow></mfrac></mtd><mtd><mfrac><mrow><msub><mi>n</mi><mn>1</mn></msub><msub><mi>n</mi><mn>2</mn></msub></mrow><mrow><msup><msub><mi>n</mi><mn>1</mn></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>n</mi><mn>2</mn></msub><mn>2</mn></msup></mrow></mfrac></mtd></mtr><mtr><mtd><mfrac><mrow><msub><mi>n</mi><mn>1</mn></msub><msub><mi>n</mi><mn>2</mn></msub></mrow><mrow><msup><msub><mi>n</mi><mn>1</mn></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>n</mi><mn>2</mn></msub><mn>2</mn></msup></mrow></mfrac></mtd><mtd><mn>1</mn><mo>+</mo><mfrac><msup><msub><mi>n</mi><mn>2</mn></msub><mn>2</mn></msup><mrow><msup><msub><mi>n</mi><mn>1</mn></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>n</mi><mn>2</mn></msub><mn>2</mn></msup></mrow></mfrac></mtd></mtr></mtable></mfenced><mo>,</mo></mrow><mrow><mi>b</mi><mo>=</mo><mfencedopen='['close=']'><mtable><mtr><mtd><msub><mi>b</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><msub><mi>b</mi><mi>y</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfrac><mrow><msub><mi>α</mi><mn>1</mn></msub><msub><mi>n</mi><mn>1</mn></msub><mo>+</mo><msub><mi>α</mi><mn>2</mn></msub><msub><mi>n</mi><mn>2</mn></msub></mrow><mrow><msup><msub><mi>n</mi><mn>1</mn></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>n</mi><mn>2</mn></msub><mn>2</mn></msup></mrow></mfrac><mfencedopen='['close=']'><mtable><mtr><mtd><msub><mi>n</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>n</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced></mrow>设垂足曲线pe(t)是二维曲面函数φ2的零水平集，φ2的进化方程采用ChanVese改进模型，即<mrow><mfrac><msub><mrow><mo>∂</mo><mi>φ</mi></mrow><mn>2</mn></msub><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>=</mo><mo>|</mo><mo>▿</mo><msub><mi>φ</mi><mn>2</mn></msub><mo>|</mo><mo>[</mo><mi>μκ</mi><mo>-</mo><mi>v</mi><mo>-</mo><msub><mi>λ</mi><mn>1</mn></msub><msup><mrow><mo>(</mo><mi>I</mi><mo>-</mo><msub><mi>c</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msub><mi>λ</mi><mn>2</mn></msub><msup><mrow><mo>(</mo><mi>I</mi><mo>-</mo><msub><mi>c</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo><mo>,</mo></mrow>式中μ≥0，v≥0，λ1，λ2＞0为固定常数，κ为水平集函数曲率，参数c1，c2计算如下<mrow><msub><mi>c</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><msub><mo>∫</mo><mi>Ω</mi></msub><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mi>H</mi><mrow><mo>(</mo><mi>φ</mi><mo>)</mo></mrow><mi>dxdy</mi></mrow><mrow><msub><mo>∫</mo><mi>Ω</mi></msub><mi>H</mi><mrow><mo>(</mo><mi>φ</mi><mo>)</mo></mrow><mi>dxdy</mi></mrow></mfrac><mo>,</mo></mrow><mrow><msub><mi>c</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><msub><mo>∫</mo><mi>Ω</mi></msub><mi>I</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>H</mi><mrow><mo>(</mo><mi>φ</mi><mo>)</mo></mrow><mo>)</mo></mrow><mi>dxdy</mi></mrow><mrow><msub><mo>∫</mo><mi>Ω</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>H</mi><mrow><mo>(</mo><mi>φ</mi><mo>)</mo></mrow><mo>)</mo></mrow><mi>dxdy</mi></mrow></mfrac><mo>,</mo></mrow>其中<mrow><mi>H</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><mfencedopen='['close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>z</mi><mo>></mo><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>z</mi><mo>≤</mo><mn>0</mn></mtd></mtr></mtable></mfenced></mrow>为Heaviside函数；垂足点控制曲线p(t)通过对二维曲面函数φ1进行进化来求解，φ1的进化方程通过p(t)与pe(t)的关系以及φ2的进化方程推导得到φ1的进化方程为<mrow><mfrac><mrow><mo>∂</mo><msub><mi>φ</mi><mn>1</mn></msub></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>=</mo><mo>|</mo><msup><msub><mi>J</mi><mi>D</mi></msub><mi>T</mi></msup><msub><mrow><mo>▿</mo><mi>φ</mi></mrow><mn>1</mn></msub><mo>|</mo><mo>[</mo><mi>μ</mi><mo>·</mo><msub><mi>κ</mi><mi>e</mi></msub><mo>+</mo><msub><mi>F</mi><mi>CV</mi></msub><mo>]</mo><mo>,</mo></mrow>其中<mrow><msub><mi>κ</mi><mi>e</mi></msub><mi>div</mi><mrow><mo>(</mo><mfrac><mrow><msup><msub><mi>J</mi><mi>D</mi></msub><mi>T</mi></msup><msub><mrow><mo>▿</mo><mi>φ</mi></mrow><mn>1</mn></msub></mrow><mrow><mo>|</mo><msup><msub><mi>J</mi><mi>D</mi></msub><mi>T</mi></msup><msub><mrow><mo>▿</mo><mi>φ</mi></mrow><mn>1</mn></msub><mo>|</mo></mrow></mfrac><mo>)</mo></mrow><mo>,</mo></mrow>FCV＝vλ1(Ic1)2+λ2(Ic2)2，JD＝JC1，<mrow><msub><mi>J</mi><mi>C</mi></msub><mo>=</mo><mfencedopen='['close=']'><mtable><mtr><mtd><mfrac><msub><mrow><mo>∂</mo><mi>p</mi></mrow><mi>ex</mi></msub><msub><mrow><mo>∂</mo><mi>p</mi></mrow><mi>x</mi></msub></mfrac></mtd><mtd><mfrac><msub><mrow><mo>∂</mo><mi>p</mi></mrow><mi>ex</mi></msub><msub><mrow><mo>∂</mo><mi>p</mi></mrow><mi>y</mi></msub></mfrac></mtd></mtr><mtr><mtd><mfrac><msub><mrow><mo>∂</mo><mi>p</mi></mrow><mi>ey</mi></msub><mrow><mo>∂</mo><msub><mi>p</mi><mi>x</mi></msub><mo></mo></mrow></mfrac></mtd><mtd><mfrac><msub><mrow><mo>∂</mo><mi>p</mi></mrow><mi>ey</mi></msub><mrow><mo>∂</mo><msub><mi>p</mi><mi>y</mi></msub><mo></mo></mrow></mfrac></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>3)根据φ1的进化方程迭代计算φ1，ijn+1，由φ1，ijn+1得到垂足点控制曲线p(t)，然后通过p(t)对给定平面参数曲线做垂足操作得到垂足曲线pe(t)，由pe(t)得到φ2，ijn，可设置最大迭代次数，或两次迭代误差小于设定常数为迭代终止条件，满足终止条件后，垂足曲线pe(t)即为最终得到的人脸轮廓提取结果。全文摘要本发明涉及一种基于变形垂足曲线的人脸轮廓提取方法，借助于变形垂足曲线能同时表示全局形状和局部形变的能力，采用人脸的先验椭圆形状作为垂足曲线的发生器，以此控制人脸的全局形状；人脸形状的局部变形通过垂足曲线的进化来实现，采用Chan-Vese模型来控制垂足曲线的进化，利用Chan-Vese模型强大的检测模糊目标边缘的能力来达到精确定位人脸轮廓的目的。由于利用水平集方法做数值求解，从而能自然地处理曲线的拓扑变化。该方法对噪声图像及平面旋转的人脸图像有一定的处理能力，可以提取出噪声图像及平面旋转的人脸图像中的完整连续人脸轮廓，而且同时可以粗略估计出人脸的大小位置和水平旋转角度等信息。文档编号G06T7/60GK101930612SQ20091019618公开日2010年12月29日申请日期2009年9月23日优先权日2009年9月23日发明者吕学勤,黄福珍申请人:上海电力学院