李梦婕等:
基于
犔
1
范数和曲波系数双约束的稀疏角度微分相位衬度计算机层析成像重建方法
θ
=
∫
δ
(
狓
,
狔
)
狔
r
d
狓
r
=
∫
δ
(
狓
,
狔
)
d
狓
r
狔
r
.
(
1
)
同时,根据图
1
所示的几何关系可得
狔
=
狓
r
sin
+
狔
r
cos
狓
=
狓
r
cos
-
狔
r
sin
烅
烄
烆
.
(
2
)
将(
2
)式代入(
1
)式可得
∫
δ
(
狓
,
狔
)
狔
d
狓
r
=
sin
∫
δ
(
狓
,
狔
)
狓
r
d
狓
r
+
cos
∫
δ
(
狓
,
狔
)
狔
r
d
狓
r
∫
δ
(
狓
,
狔
)
狓
d
狓
r
=
cos
∫
δ
(
狓
,
狔
)
狓
r
d
狓
r
-
sin
∫
δ
(
狓
,
狔
)
狔
r
d
狓
烅
烄
烆
r
, (
3
)
式中
δ
(
狓
,
狔
)/
狓
r
表示相位项沿
狓
r
方向的梯度,如
果将空气的折射系数认为是
0
,则
∫
δ
(
狓
,
狔
)/
狓
r
d
狓
r
≡
0.
(
4
)
将(
4
)式代入(
3
)式化简后,可得
θ
cos
=
∫
δ
(
狓
,
狔
)
狔
d
狓
r
θ
sin
=-
∫
δ
(
狓
,
狔
)
狓
d
狓
烅
烄
烆
r
.
(
5
)
以相位项梯度
δ
为目标的重建算法是以(
5
)式
为基础的。将相位项水平方向的梯度
δ
(
狓
,
狔
)/
狓
的
相反数作为待重建图像
狌
,将折射角
θ
sin
作为测量
值
狆
,二维平面离散化的待重建物体,如图
1
(
b
)所示,
根据(
5
)式中水平方向的梯度公式可知,第
犻
条射线
穿过物体后的折射角
θ
犻
满足
θ
犻
sin
犻
=
∑
犜
狋
狉
犻
,
狋
狌
狋
,
犻
=
1
,
2
,…,
犐
, (
6
)
式中
犻
为射线与
狓
轴的夹角;
狌
狋
为待重建图像的第
狋
个像素;
犐
为射线总个数;
犜
为像素值总个数;如果射
线
犻
穿过像素
狋
,
狉
犻
,
狋
为
1
,否则为
0
。由于在实际运用
中,测量误差和噪声无法避免,故 在考虑 存在误 差
犲
的情况下,(
6
)式可用矩阵表示为
狆
=
犚
·
狌
+
犲
, (
7
)
式中
犲
为误差向量,
犚
为投影矩阵,定义为
犚
=
狉
1
,
1
狉
1
,
2
…
狉
1
,
犜
狉
2
,
1
狉
2
,
2
…
狉
2
,
犜
狉
犐
,
1
狉
犐
,
2
…
狉
犐
,
熿
燀
燄
燅
犜
.
(
8
)
3
算法原理
3.1
算法原理
本文将(
7
)式的解集表示为
犆
0
=
狌
狘
‖
狆
-
犚
·
狌
‖
2
2
<
ξ
{ }
0
, (
9
)
式中
ξ
0
为反 映 误差 的 门限 参 数。 在 稀 疏角 度 条件
下,由于数据量的缺失,观测数据个数往往远小于待
重建图像的像素个数。此时,集合
犆
0
中符合要求的
解往往有无穷多个,而 待 重建 图 像只 是 这些 解 当中
的一个。若要从这 些 解中 找 出待 重 建图 像,则需 要
增加约束条件。
与光学图 像不同,物体 的断层 图像的 像素取 决
于该位置的材质,而相 同 材质 的 物质 总 是聚 集 在一
起的,因此,如果用
δ
代表断层图像的边沿,则
δ
具有良 好 的稀 疏 性。众 所周 知,
犔
1
范数 是 求 解 稀 疏
信号时常见 的 约 束 条 件。本文 用 集 合 的 形 式,将 该
约束条件表示为
犆
1
=
{
狌
狘
‖
狌
‖
1
<
ξ
1
}, (
10
)
式中
ξ
1
为
犔
1
范数的门限参数。
δ
通常表现为若干曲线,但
犔
1
范数只 能保证
重建出的图像是稀 疏 的
,不 能 保证 非 零像 素 点可 以
连接成曲线。曲波变换是一种对描述图像中曲线有
着相当强表达能力的变换方法,与小波变换相比,曲
波变换中的曲波基 具 有方 向 性,对于 表达 具 有方 向
性的曲线有着无可比拟的稀疏优越性。第一代曲波
变换利用了脊波变 换 的相 关 知识,而 第 二代 曲 波变
换主要在频域进行,在频域中通 过快速 算法
[
11
]
来得
到曲波 系 数 的 频 域 表 达 式。
2005
年,
Candes
等
[
11
]
提出了两种快速离 散 曲波 变 换算 法,即非 等 距快 速
傅里叶变换(
USFFT
)算法 和 频域 装 配 (
FW
)算 法,
并证明了频域装配算法在计算速度等方面优于非等
距快速傅里叶变换算法,因此,本文选用频域装配算
法。
空域曲 波 基
φ
犼
,
犾
(
犽
1
,
犽
2
)都 是 由 尺 寸 参 数
犼
、表
示方向的角度参 数
犾
和 空 间 坐 标(
犽
1
,
犽
2
)这 三 个 参
量来标定的。图
2
(
a
)为一个在固定尺度
犼
和方向角
度参 数
犾
范 围 内 的 曲 波 基
φ
犼
,
犾
(
犽
1
,
犽
2
)的 示 意 图,
图
2
(
b
)为傅里 叶 域中 的
φ
犼
,
犾
(
ω
1
,
ω
2
)。由图
2
(
a
)可
见,离散曲波基
φ
犼
,
犾
(
犽
1
,
犽
2
)具有明显的方向性:当曲
波基
φ
犼
,
犾
(
犽
1
,
犽
2
)与曲线
狇
大体重合时,如图
2
(
c
)所
示,快 速算法 得出的 曲波系 数
犮
犼
.
犾
(
犽
1
,
犽
2
)是一 个幅
值比较大的数;否则,如 图
2
(
d
)所 示,系 数
犮
犼
.
犾
(
犽
1
,
犽
2
)是一个接近于
0
的数。因此,用非常少的曲波系
01110033