基于纳维-斯托克斯方程的图像修补模型的实现.pdfVIP

1O6 福 建 电 脑 2008年第7期 基于纳维一斯托克斯方程的图像修补模型的实现 黄陈思 (福州大学数学与计算机科学学院 福建福州 350108) 【摘 要】：图像修补是对图像损毁区域进行适当修补以满足人眼视觉要求的一种技术。Bertalmio等人证明了流体力学 中的纳堆一斯托克斯方程可以应用于图像修补．提出了基于纳堆一斯托克斯方程的修补模型。本文对此模型的设计思想以及 数值实现进行了详细的阐述．并分析了实验结果。 【关键字】：图像修补；纳维一斯托克斯方程；数值实现 1、引言 无散向量。因此有V·u=O。对于二维流体，可以引入一个流函数 图像修补作为一项特殊的图像恢复任务．其目的是利用待 来表示速度 V =l‘。在流体中。流速向量“的旋度定义为流体 修补区域周围的信息修复图像中的破损区域或去除图像中的文 微团的涡量，以∞表示即w=-~Txu。在(3)式中消去P，并在等式两 字和物体以满足人眼视觉系统的要求。这项技术已逐渐成为研 端取旋度．可得 究的热点。大多数图像修补技术都需要经过人工干预事先选定 + 。V∞=uArn (4) 需要修复的区域。然后利用修补技术自动修复缺失信息的区域。 在粘度系数为0的情况下。(4)的稳定状态为： 目前．图像修补技术中主要有基于偏微分方程的修补模型以及 · V∞=V 、y．V△ ≈0 (5) 基于纹理的修补模型。其中．最初出现”image inpainting”(图像修 注意到(5)与(1)的相似性。对应可知流函数 在图像中对应为亮 补1词语的是Bertalmio等人提出的BSCB模型【l1。该模型依据艺 度函数l。因此，流速向量∞：V×u则对应了等照度线方向．而涡 术家修复图片的方法．利用等照度线的传输将图像平滑度填入 量“=V =V 则对应了图像的Laplacian算子，粘度系数则对 目标区域。由于等照度线之间没有信息交互．因此需要强制加入 应着各向异性扩散系数。通过以上的分析。可以建立起流体动力 各向异性扩散以保证模型的稳定性。模型缺乏数学依据嘲。 学与图像处理之间的联系。由f41可得基于纳维一斯托克斯方程的 Chn等人提出了一系列基于变分问题的图像修补模型．总 图像修补模型为： 变分frⅥ模型田通过求有界变分空间中总体变分最小值来修复 q+l|·V∞=vV·(g(1w1)w)， (6) 图像．由于其模型中最终曲率为零，因此TV模型的扩散结果在 其中，g为一递减函数，在P-M各向异性扩散【6】中可取g为以下 图像中形成直线．且不满足连接性准则 。针对TV模型的问题． 形式： ． Chn等人又提出了基于曲率驱动的CDD(curvature driven dittu． g )=]／0+ ／ )， 仃、 sio 模 [31和修复曲线的欧拉弹性(Eulers Elastica)模型1．1。基于 本文算法的实现取K=2。 变分问题的修补模型具有较好的理论基础．但是修复的图像不 4、数值实现 自然．选定目标区域轮廓明显l 。 4．1算法 Bertalmio等人在闭中证明了无粘性流体力学中的纳维一斯托 基于纳维一斯托克斯方程的修补模型为方程f61，但是修补的 克斯方程与BSCB修补模型的等价性。由此．为BSCB模型提供 最终目的是求得待修补区域中的像素灰度值．因此同时要求解 了有力的理论依据。 的还有方程∞=XT!。因此先假设图像中待修补区域 为一矩形区 2、BSCB修补模型 域。 BSCB修补模型将图像平滑度沿等照度线方向推进到待修 4．2数值差分形式 补区域中『I1。其修复模型是： 将(6)展开后可得： = V I·V