Efficient Graph-Based Image Segmentation

本文相关

原文链接

基础知识

  • 一张图是由不同的像素点构成的,本文的计算和构建都是基于像素点的运算,即(RGB)
  • 高斯模糊/拉普拉斯变换:用于转换图像,减少图像噪声的平滑算法
  • 最小生成树(Minimum Spanning Tree | MST)指的是,在图中建立一个连通图并且没有回路是生成树,而最小生成树指的是构成结果权值最小
  • 不同像素点之间的差:即RGB值之间的欧氏距离
    • \sqrt{(R_1-R_2)^2+(G_1-G_2)^2+(B_1-B_2)^2}
  • 并查集算法(union find set)以及克鲁斯卡尔算法(Kruskal),使用边建立并查集,并且使用kruskal进行搜索合并

早期的分割方法

  • Zahn提出了一种基于图的最小生成树(MST)的分割方法,用来进行点聚类以及图像分割,前者权值是点间距离,后者权值是像素差异。
  • 不足:根据阈值不同,会导致高可变性(大约是色彩对比强的一个区域)区域划分为多个区域;将ramp和constant region合并到一起。
  • Urquhart提出用边相连的点中边权值最小的进行归一化,找周围相似的。
  • 根据各个区域是否符合某种均匀性标准来分割,找均匀强度或梯度的区域,不适用于某个变化很大的区域。
  • 使用特征空间聚类:通过平滑数据——给定半径的超球面对各个点扩张其连通分量,找到簇,来保持该区域的边界,并对数据进行转换。

基于图的分割

定义

  • G:将图像由像素点转化为图
  • V:每一个像素点都是图中的点
  • E:任意两个相邻像素点之间边
  • C:被划分的Segmentation,一个C中有至少1个像素点
  • Int(C):区域内最小生成树权值最大的边,表示的是,记为
    • Int(C) = \max{w(e)} , e∈MST(C,E)
  • Dif(C1,C2):表示C1和C2之间的距离,记为
    • Dif(C1,C2) = \min{w(vi,vj)} ,vi∈C1,vj∈C2,(vi,vj)∈E
  • 最后要形成的分组要求是(表示了所有区域之间的最小距离都比区域内的最大距离和权值的和要大)
    • D(C1,C2)=\left\{ \begin{aligned} true, & & \ if Dif(C1,C2)>MInt(C1,C2) \\ false, & & \ otherwise \end{aligned} \right.
  • 其中MInt(C1,C2)的值为:
    • MInt(C1,C2) = (Int(C1)+τ(C1),Int(C2)+τ(C2))
  • 阈值设定的原因是为了在开始时,因为只有单个像素点,那么点内的距离为0,而点之间的距离还存在,那么导致无法合并。所以加入阈值。
    • τ(C) = k/|C|

分割算法(与克鲁斯卡尔算法构建最小生成树有密切关系。)

  • 输入是一个有n个节点和m条边的图G,输出是一系列区域。步骤如下:
  • 0.将边按照权重值以非递减方式排序
  • 1.最初的分割记为S(0),即每一个节点属于一个区域。
  • 2.按照以下的方式由S(q-1)构造S(q):记第q条边连接的两个节点为vi和vj,如果在S(q-1)中vi和vj是分别属于两个区域并且第q条边的权重小于两个区域的区域内间距,则合并两个区域。否则令S(q) = S(q-1)。
  • 3.从q=1到q=m,重复步骤2。
  • 4.返回S(m)即为所求分割区域集合。

补充

高斯滤波器

  • 高斯变换就是用高斯函数对图像进行卷积,高斯滤波器是一种线性滤波器,能够有效抑制噪声,并平滑图像。其实质是取滤波器窗口内像素的均值作为输出。
  • 高斯函数公式如下:
    • f(x) = \frac{1}{\sigma\sqrt{2\pi}}e^{ -\frac{(x-\mu)^2}{2\sigma^2}}
      其中,ux的均值,σ是方差。
  • 由一维函数,我们可以推导出二维函数的公式如下:
    • f(x,y) = \frac{1}{2\pi \sigma^2} e^{-\frac{(x^2+y^2)}{2\sigma^2} }
  • 高斯函数在图像处理中的使用,实际上就是对每个像素点的周边像素取平均值,从而达到平滑的效果,在取值(周边半径)时,周围像素点的半径越大,则图像的模糊度就越强。在实际计算时,利用高斯模糊按正态曲线分配周边像素的权重,从而求中心点的加权平均值。
  • 高斯模糊的具体计算方式如下:
    • 1.将中心点周围的八个点带入到高斯函数中,从而得到权重矩阵A1;
    • 2.为使归一化,将矩阵A1中的各个点除以所有点(9个点)的权重和,得到归一化后的权重矩阵A2;
    • 3.图片原始的像素矩阵分别乘以A2中各自的权重值,将得到的所有点的值加起来求平均,便得到中心点的高斯模糊值。图像中其余点相同求法。
    • 注:1.彩色图片,可对RGB三通道分别作高斯模糊。
  • 2.σ代表数据的离散程度,σ越大,中心系数越小,图像越平滑;反之,反之。

拉普拉斯变换:是为解决傅立叶变换等幅振荡的缺点。

  • 首先了解一下傅立叶变换:傅立叶变换是一种物理上探究频谱的方法,三角公式是:
    • f(t) = \sum_{n=1}^\infty A_ncos(nw_0t+\varphi_n)+B
    • 其中,w0表示基波。
  • 由欧拉公式:
    • \left\{ \begin{aligned} e^{ix} =cosx+isinx, \\ e^{-ix} =cosx -isinx, \end{aligned} \right.
  • 将傅立叶三角形式公式中的正余弦函数用指数函数表示,改写为用复指数表示的公式,如下:
    • f(t) = \sum_{-\infty}^\infty F(nw_0)e^{jw_0t}
  • 将上述公式改为积分形式,即得到复指数形式公式为:
    • F(w) =\int_{-\infty}^\infty f(t)e^{-jwt}dt
  • 但由于傅立叶变换是等幅振荡的正弦波,故当f(t)不断趋向无穷时,此时函数将不再收敛,这时候便不再适合使用傅立叶变换。于是,我们引入一个衰减因子,对其作变换。对函数y=f(t)乘上一个e^{\sigma t},其中,σ>0。
    • F(w) =\int_{-\infty}^\infty f(t)e^{-\sigma t}e^{-jwt}dt
  • 对上式进行合并同类项,可得到F(w) =\int_{-\infty}^\infty f(t)e^{-t(\sigma+jw)}dt
  • 我们将指数中的σ+jw最初的分割记为S,于是得到拉普拉斯公式:
    • \Rightarrow  F(w) =\int_{-\infty}^\infty f(t)e^{-st}dt
  • 由上式推导,很清楚的知道,当s=jw时,拉普拉斯函数就变成了傅立叶函数,也就相当于拉氏不再具有衰减功能。
  • 又由上述公式可以很直观地看到当取值\sigma_0刚好收敛时,则\sigma>\sigma_0的区域全都收敛。
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 158,736评论 4 362
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 67,167评论 1 291
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 108,442评论 0 243
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 43,902评论 0 204
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 52,302评论 3 287
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 40,573评论 1 216
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 31,847评论 2 312
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 30,562评论 0 197
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 34,260评论 1 241
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 30,531评论 2 245
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 32,021评论 1 258
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 28,367评论 2 253
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 33,016评论 3 235
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 26,068评论 0 8
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 26,827评论 0 194
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 35,610评论 2 274
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 35,514评论 2 269

推荐阅读更多精彩内容