Kass
,
Witkin
等
[
1
]
首 次 提 出 了 基 于 能 量 最 小
化模型的参 数 曲 线 演 化 方 法
snake.Snake
是 在 图
像域内定义可变形 曲 线,通过对其能量函数的最小
化及变形和调整从 而 产 生 连 续 平 滑 的 轮 廓,当能量
达到最小 时 曲 线 收 敛 到 目 标 边 界
.Caselles
等
[
2
]
和
Malladi
等
[
3
]
提 出 基 于 几 何 形 变 模 型 (
g
eometric
activecontour
)的曲线演化方法,将曲线的形状变化
用曲率、法向量和梯度等几何度量表示,作为曲线演
化的速度 函 数
.Caselles
等
[
4
]
根 据 能 量 最 小 化 模 型
提出了测地线(
g
eodesicactivecontours
)演化方 程,
通过引入约束项使曲线即使跨过弱边界也可以被拉
回来
.
文献[
15
]中的方法一般依据图像梯度信息进
行演化计算,在灰度值变化剧烈的边缘容易收敛,但
在边界模糊 或 渐 变 区 域 会 导 致 演 化 抖 动 或 越 过 边
缘
.
另外,图像噪声 可 以 使 模 型 陷 入 局 部 最 优,干扰
对图像轮廓的准确 提 取
.
虽然可以对图像进行滤波
去噪后再提取,但由于边缘信息被弱化,使得后续处
理变得较为 困难
.
近年来
Mumford
和
Shah
等
[
6
]
提
出的基于最优划分的
MS
图 像 分 割 模 型 以 及
Chan
和
Vese
等
[
7
]
提出的
CV
模型成为研究的热点,一些
有关的新方 法 和 新 观 点 也 相 继 被 提 出
.
例 如,文 献
[
8
]对
CV
方法的欧拉
拉格朗日偏 微分方程 求解做
了改进,提 出 了 快 速 构 造
SDF
的 方 法;文 献 [
9
]对
CV
模型进行 改 进,提 出 了 曲 线 分 层 演 化 模 型 和 多
水平集方程求解方 法
.Han
等提出 了边界提取与配
准结合的边界匹配算 法
[
10
]
,利用
MS
模型及 其自动
耦合函数进行边界 特 征 表 示,将曲线演化划分为线
性与非线性两部分,分别用有限元和梯度下降法迭
代计算
.Du
等改进了
MS
模型
[
11
]
,用高频滤 波分量
作为曲线演化的耦合函数,并用单变量
PDE
和高频
卷积滤波加速演化函数的求解
.
本文将改进的
MS
模型作为图像轮廓提取的曲
线演化模型,在模型 中 引 入 模 糊 方 法 定 义 轮 廓 拟 合
能量函数,同时改进 基 于 全 局 信 息 的 轮 廓 演 化 能 量
函数,降低由于边缘 模 糊 和 灰 度 不 均 匀 等 因 素 所 产
生的轮廓线演化抖 动;对模型的数值处理方法进行
研究,推导出最小化能量模型的水平集演化方程,给
出了基于加性
乘 性 算 子 分 裂 的 水 平 集 函 数 迭 代 计
算方法
.
2
犕狌犿犳狅狉犱犛犺犪犺
和
犆犺犪狀犞犲狊犲
模型
MumfordShah
模型(简称
MS
)是
20
世纪
80
年
代提出的、具有曲线 演 化 拓 扑 自 适 应 能 力 的 分 析 模
型,它在理论上是一 类 具 有 体 积 能 量 和 低 维 测 度 的
变分问题
.MS
模 型 将 局 部 边 缘 信 息 与 全 局 同 质 区
域信息相结合,可以 实 现 对 轮 廓 曲 线 的 比 较 精 确 的
定位,并对噪声具有一定的抑制作用
.
MumfordShah
模型的一般形式为
犉
(
狌
,
犆
)
=
α
∫
Ω
\
犆
|
狌
|
2
d
狓
d
狔
+
β
∫
Ω
|
狌
-
狌
0
|
d
狓
d
狔
+
γ
·
犾犲狀
犵
狋犺
(
犆
) (
1
)
式中,
狌
0
为开集
Ω
∈
犚
2
上的噪声图像,闭集
犆
∈
犚
2
为
Ω
上不连续集合,表征图像上不连续的突变点集合,
如边缘、角 点 等;
狌
为 近 似 于
狌
0
且 仅 与
犆
相 关 的 分
段光滑函数
.MS
模 型 的 求 解 可 归 为 寻 找 合 适 的 数
对(
狌
,
犆
)使
犉
(
狌
,
犆
)逼近极小值,即
{
狌
,
犆
}
=inf
狌
,
犆
犉
(
狌
.
犆
)
.
由于
狌
和
犆
分属不同维数空间集合,且
犉
(
狌
,
犆
)
非凸,函数容易陷入局部极小值,因而在实际应用中
求数值解或数值 逼 近 比 较 困 难
.
在
CV
模 型 中 通 过
水平集曲 线演化方 法,将
MS
模型转 化为简化 的隐
式模型,其能量泛函为
犉
(
犆
,
犮
1
,
犮
2
)
=
μ
·
犾犲狀
犵
狋犺
(
犆
)
+
ν
·
犪狉犲犪
(
犻狀狊犻犱犲犆
)
+
λ
1
∫
犻狀狊犻犱犲
(
犆
)
|
狌
0
-
犮
1
|
2
d
狓
d
狔
+
λ
2
∫
狅狌狋狊犻犱犲
(
犆
)
|
狌
0
-
犮
2
|
2
d
狓
d
狔
(
2
)
式中,
μ
,
ν
0
,
λ
1
,
λ
2
>
0
为各项的权重系数;图像
狌
0
被闭合曲线
犆
分割为 内 外 两 个 区 域,
犮
1
和
犮
2
分别为
这两个区域 的 灰 度 均 值
.
当
犆
没 有 位 于
狌
0
的 边 界
时,能量泛函无法 达 到 极 小 值,曲 线 会 向
狌
0
的 区 域
边界处演化
.
只有
犆
逼近
狌
0
中同质区域的边界时,
犆
的演化逐渐趋于稳定,
犉
(
犆
,
犮
1
,
犮
2
)趋近极小值
.
在曲
线演化
PDE
的每次迭代过程中,
犮
1
和
犮
2
的计算函数
对图像具有光滑效 果,其作用虽然可以在一定程度
上去噪,但同时轮廓提取的准确性也降低了
.MR
图
像边缘往往存在着边界模糊和灰度渐进等现象
.MS
和
CV
模型虽然 可 以 使 区 域 分 割 逼 近 全 局 最 优,但
在处理 渐 进 和 模 糊 边 界 时,仍 难 以 精 确 刻 画 边 界
的轮廓
.
另外,不 同 的 初 始 曲 线 对 边 缘 提 取 结 果 也
有影响
.
3
轮廓提取模型
3.1
图像轮廓提取模型
用曲 线演化进 行 轮 廓 提 取 的 过 程 可 解 释 为 在
给定初 始 轮 廓 线
犆
的 条 件 下,在 演 化 能 量 的 作 用
下,通 过
犆
的 演 化 使 得 轮 廓 演 化 模 型 的 能 量 泛 函
962
2
期 原
达等:基于
MumfordShah
模型的高精度
MR
图像轮廓提取算法