图像抠图的closed form算法(A Closed-Form Solution to Natural Image Matting)

2024-01-01 20:30

本文主要是介绍图像抠图的closed form算法(A Closed-Form Solution to Natural Image Matting),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

关于图像抠图算法,Levin等人在2007年基于图像的局部光滑假设,利用代数的方法推导出了alpha matte矩阵闭合解的形式。原文名称是”A Closed Form Solution to Natural Image Matting”。
在抠图问题中,假设第i个像素点的值 Ii I i 是由前景点 Fi F i 和背景点 Bi B i 按下式加权合成的:

Ii=αiFi+(1αi)Bi1 I i = α i F i + ( 1 − α i ) B i 式 1

其中 αi α i 称为前景透明度,也就是要求解的”alpha matte”矩阵,指定了原图中对应像素点的类型——属于背景( αi=0 α i = 0 )、前景( αi=1 α i = 1 )或是未知(0< αi α i <1)。鉴于前景图F和背景图B都是未知,要想求解出alpha矩阵,需要其他条件——用户输入提示和图像分布假设。

用户输入方式

不同的抠图系统所接受的用户输入一般有两种形式——trimap或scribbles。两种方式的目的都是要给出alpha matte矩阵的初始估计。trimap的方式需要用户指出所有的混合像素点,哪怕以掺入大量的前景点和背景点为代价,如(a)所示;而scribbles的方法只需用户在确定的前景和背景区域分别画出一些条状区域即可,如下图(b)所示,在转化成alpha matte矩阵时,α的值同样也是前景处为1,背景处为0,scribble的交互方式显然比trimap的方式要简单,尤其当对于图像前景轮廓比较复杂的时候。
这里写图片描述

灰度图

为了使式1可以被求解,必须对F和B施加一些假设,该算法即假定F和B是局部光滑的,也就是在一个小的图像窗口w中(如3x3),F和B固定不变,那么式(1)就可以写成:

αiαIi+b,iw2 α i ≈ α I i + b , ∀ i ∈ w 式 2

也就说α与 Ii I i 是线性相关(α= 1FB 1 F − B , b= BFB − B F − B 是常数)。那么就可以简单的定义要优化的目标函数为α的估计值与实际值的差:
这里写图片描述 式3
从表达式可以看出,子窗口存在大量重叠,正是该性质使得像素信息可以进行传播。该算法的巧妙之处在于通过一些代数推导,可以将a和b从式3中消去,转换为只有一个变量α的优化问题。
这里写图片描述 式4
其中,α是Nx1维的列向量(N是像素总数),L是一个NxN的对称正定矩阵:
这里写图片描述
只有当i=j时, δij δ i j 才等于1,其他情况下为0, uk u k σ2k σ k 2 分别是子窗口内像素的均值和方差。
在抠图算法中,L是一个很重要的矩阵,人们称之为matting laplace矩阵,它除了是一个对称的正定矩阵之外,还是一个稀疏矩阵,只有满足|i−j|<| wk w k |−1的像素处L才是非0。所以,虽然L的维度很高,但仍可以采用高效的方式求解,尤其是利用MATLAB自带的一些稀疏矩阵函数。
推导过程也比较简单,主要就是先将式3平方和累加的形式转换为列向量的模,恰好是一个最小二乘法的形式,再根据最小二乘法解的公式求出使得J最小时a,b的值,即可消去a,b,只留下α。在抠图时,除了希望得到alpha matte矩阵,一般还希望可以得到对应的前景图和背景图的值,所以在这里详细推导一遍,同时给出a,b与α之间的关系。
用列向量模的形式表示式3:
这里写图片描述
其中,
这里写图片描述
这里写图片描述
这是一个典型的最小二乘问题,求解结果就是:
这里写图片描述
带入J(α)的表达式中即可整理出L。

彩色图

如果是rgb图,因为有了3个通道,所以需要对每个通道重复上述过程,但原文却采取了另外一种假设模型——color line model。与灰度图的局部光滑假设不同,color line 模型假定子窗口内的前景和背景像素值不是固定的,而是位于一条直线上:
这里写图片描述
将此假设带入到式1中,经过推导(原文中给出了证明),可以得到彩色图像的线性模型:
这里写图片描述
ac a c 和b是仅与F,B,I有关的变量。将式5带入灰度图像的推导过程,同样可以将目标函数化简成式4,只是L略有不同:
这里写图片描述
此时 k ∑ k 变成了窗口内所有像素的协方差矩阵, uk u k 是三维列向量,分别对应三个通道的像素均值, I3 I 3 是3x3的单位阵。L又被称为matting laplace矩阵。

用户输入约束下的优化问题

施加了像素分布的假设后,还需要考虑用户指定的前景和背景的先验分布。当用户按图(b)的方式在原始图像上画出一系列区域时,α中相应像素点被赋值1(前景)或0(背景),那么式4的优化问题就变成了:
这里写图片描述
S是用户所指定的像素点的集合, si s i 取0或1。
为求解式5,首先对L中的像素位置进行重组,再次强调之前所说的,L是对称矩阵:
这里写图片描述
LM L M 是已知像素点的集合, LU L U 是未知像素点的集合,也就是将图像中所有标记出的像素排在前面,要求解的未知像素排在后面,当然为保证变换前后等式关系一致,α也需要相应调整:
这里写图片描述
带入式4:
这里写图片描述
要求J(α)的最小值,需要进行矩阵求导,矩阵求导的通用公式如下:
这里写图片描述
这里写图片描述
因为 αM α M 已知,所以上式中的第一项导数为0:
这里写图片描述
LU L U 是对称的,所以最后求出的解为:
这里写图片描述

算法运行效果

原图和用户输入的scribble图如下:
这里写图片描述
程序输出的结果如下:
这里写图片描述
可以看出,尽管该算法比较简单,但即使是与一些采用繁复的迭代的算法相比,其效果也是不错的。
当然,最终结果的好坏还是取决于两个因素——图像是否满足color line model的假设以及用户输入的先验信息,据原文的说法,一般的自然图像都可以较好的满足color line model。

文章来源

转载地址:https://blog.csdn.net/edesignerj/article/details/53349663

这篇关于图像抠图的closed form算法(A Closed-Form Solution to Natural Image Matting)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/560423

相关文章

SpringBoot实现MD5加盐算法的示例代码

《SpringBoot实现MD5加盐算法的示例代码》加盐算法是一种用于增强密码安全性的技术,本文主要介绍了SpringBoot实现MD5加盐算法的示例代码,文中通过示例代码介绍的非常详细,对大家的学习... 目录一、什么是加盐算法二、如何实现加盐算法2.1 加盐算法代码实现2.2 注册页面中进行密码加盐2.

Java时间轮调度算法的代码实现

《Java时间轮调度算法的代码实现》时间轮是一种高效的定时调度算法,主要用于管理延时任务或周期性任务,它通过一个环形数组(时间轮)和指针来实现,将大量定时任务分摊到固定的时间槽中,极大地降低了时间复杂... 目录1、简述2、时间轮的原理3. 时间轮的实现步骤3.1 定义时间槽3.2 定义时间轮3.3 使用时

如何通过Golang的container/list实现LRU缓存算法

《如何通过Golang的container/list实现LRU缓存算法》文章介绍了Go语言中container/list包实现的双向链表,并探讨了如何使用链表实现LRU缓存,LRU缓存通过维护一个双向... 目录力扣:146. LRU 缓存主要结构 List 和 Element常用方法1. 初始化链表2.

使用Python开发一个图像标注与OCR识别工具

《使用Python开发一个图像标注与OCR识别工具》:本文主要介绍一个使用Python开发的工具,允许用户在图像上进行矩形标注,使用OCR对标注区域进行文本识别,并将结果保存为Excel文件,感兴... 目录项目简介1. 图像加载与显示2. 矩形标注3. OCR识别4. 标注的保存与加载5. 裁剪与重置图像

golang字符串匹配算法解读

《golang字符串匹配算法解读》文章介绍了字符串匹配算法的原理,特别是Knuth-Morris-Pratt(KMP)算法,该算法通过构建模式串的前缀表来减少匹配时的不必要的字符比较,从而提高效率,在... 目录简介KMP实现代码总结简介字符串匹配算法主要用于在一个较长的文本串中查找一个较短的字符串(称为

通俗易懂的Java常见限流算法具体实现

《通俗易懂的Java常见限流算法具体实现》:本文主要介绍Java常见限流算法具体实现的相关资料,包括漏桶算法、令牌桶算法、Nginx限流和Redis+Lua限流的实现原理和具体步骤,并比较了它们的... 目录一、漏桶算法1.漏桶算法的思想和原理2.具体实现二、令牌桶算法1.令牌桶算法流程:2.具体实现2.1

基于WinForm+Halcon实现图像缩放与交互功能

《基于WinForm+Halcon实现图像缩放与交互功能》本文主要讲述在WinForm中结合Halcon实现图像缩放、平移及实时显示灰度值等交互功能,包括初始化窗口的不同方式,以及通过特定事件添加相应... 目录前言初始化窗口添加图像缩放功能添加图像平移功能添加实时显示灰度值功能示例代码总结最后前言本文将

Python中的随机森林算法与实战

《Python中的随机森林算法与实战》本文详细介绍了随机森林算法,包括其原理、实现步骤、分类和回归案例,并讨论了其优点和缺点,通过面向对象编程实现了一个简单的随机森林模型,并应用于鸢尾花分类和波士顿房... 目录1、随机森林算法概述2、随机森林的原理3、实现步骤4、分类案例:使用随机森林预测鸢尾花品种4.1

不懂推荐算法也能设计推荐系统

本文以商业化应用推荐为例,告诉我们不懂推荐算法的产品,也能从产品侧出发, 设计出一款不错的推荐系统。 相信很多新手产品,看到算法二字,多是懵圈的。 什么排序算法、最短路径等都是相对传统的算法(注:传统是指科班出身的产品都会接触过)。但对于推荐算法,多数产品对着网上搜到的资源,都会无从下手。特别当某些推荐算法 和 “AI”扯上关系后,更是加大了理解的难度。 但,不了解推荐算法,就无法做推荐系

基于人工智能的图像分类系统

目录 引言项目背景环境准备 硬件要求软件安装与配置系统设计 系统架构关键技术代码示例 数据预处理模型训练模型预测应用场景结论 1. 引言 图像分类是计算机视觉中的一个重要任务,目标是自动识别图像中的对象类别。通过卷积神经网络(CNN)等深度学习技术,我们可以构建高效的图像分类系统,广泛应用于自动驾驶、医疗影像诊断、监控分析等领域。本文将介绍如何构建一个基于人工智能的图像分类系统,包括环境