数字图像处理A-第7章 小波和多分辨率处理
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
(7.1.10)
或
(7.1.11)
更有普遍意义的表达式
hi (2n k), g j (k) (i j) (n) i, j {0,1}
(7.1.12)
满足该条件的滤波器组称为具有双正交
《数字图像处理A》
7.1.2 子带编码
表7.0 完美重建滤波器族
(正交镜像滤波器)
(共轭正交滤波器)
在双正交的基础上进一步要求
的滤波器h5(n)由下式定义,称之为h1(n) 的调 制形式：
h5(n)=(-1)nh1(n)
《数字图像处理A》
(7.1.8)
7.1.2 子带编码
最后图7.5（f）显示的序列是h1(n)的顺序 反转形式，它被调制了：
h6(n)=(-1)nh1(K-1-n)
(7.1.9)
综上6个序列说明了这样一个事实，即在规 定两个滤波器之间的关系时，符号反转、顺 序反转和调制有时时合并在一起的。
在多分辨率分析( MRA )中，尺度函数被用 于建立某一函数或图像的一系列近似值，相 邻两近似值之间的近似度相差2倍。
被称为小波的附加函数用于对相邻近似值之 间的差异进行编码。
《数字图像处理A》
7.2.1 级数展开
信号或函数常常可以被很好地分解为一系列 展开函数的线性组合
f (x) kk (x)
《数字图像处理A》
7.1.2 子带编码
表7.0中的一维滤波器也可用于图像处理的二维可分 离滤波器。如图7.7所示。可分离滤波器首先应用于 某一维(如垂直向)，再应用于另一维(如水平向)。
行（沿m）
列（沿n） 列
列 行
列
图7.7 子带图像编码的一个二维4带宽滤波器组 《数字图像处理A》
7.1.2 子带编码
7.1.3 哈尔变换
例如，N=4时，k，q和p值如下：
4×4哈尔变换矩阵H4
1
H4
1
2
4 2
0
《数字图像处理A》
1 1 1
1
1
1
2 0 0
0 2 2
(7.1.19)
7.1.3 哈尔变换
• 例7.3 离散小波变换的哈尔函数
64×64 128×128
256×256
图7.10 (a) 用H2哈尔基函数的离散小波变换，并显示了
图7.2 (b) 创建近似和预测残差金字塔的一个简单系统 《数字图像处理A》
7.1.1 图像金字塔
例7.1 高斯和拉普拉斯金字塔
图7.3显示了图7.1中花瓶的一种可能的近似值和预 测残差金字塔。
（a）
(b)
图7.3 两种图像金字塔及它们的直方图：（a）近似金字塔；（
《数字图像处理A》
b）预测残差金字塔
图7.2 (b) 创建近似和预测残差金字塔的一个简单系统 《数字图像处理A》
7.1.1 图像金字塔
传递由3个连续步骤组成：
2.对上一步的输出进行内插—因子仍为2—并进行 过滤。这将生成与输入等分辨率的预测图像。由于 在步骤1的输出像素之间进行插值运算，插入滤波 器决定了预测值与步骤1的输入之间的近似程度。 如果插入滤波器被忽略了，预测值将是步骤1输出 的内插形式，复制像素的块效应将变得很明显。
《数字图像处理A》
主要内容
7.1 背景 7.2 多分辨率展开 7.3 一维小波变换 7.4 快速小波变换 7.5 二维小波变换 7.6 小波包
《数字图像处理A》
主要内容
7.1 背景
图像金字塔 子带编码 哈尔变换
7.2 多分辨率展开 7.3 一维小波变换 7.4 快速小波变换 7.5 二维小波变换 7.6 小波包
哈尔变换的变换矩阵H包含哈尔基函数hk(z) ，它们定义在连续闭区间z∈[0,1]， k=0,1,2…,N-1，这里N=2n。
为生成H矩阵，定义整数k，即k=2p+q-1(这 里0≤p≤n-1, p=0时，q=0或1，p≠0时， 0≤q≤2p)。
可得哈尔基函数为：
h0 (z) h00 (z)
图7.2 (b) 创建近似和预测残差金字塔的一个简单系统 《数字图像处理A》
7.1.1 图像金字塔
传递由3个连续步骤组成：
3.计算步骤2的预测值和步骤1的输入之间的差异。
以j 级预测残差进行标识的这个差异将用于原始图
像的重建(见例7.1)。在没有量化差异的情况下，预 测残差金字塔可以用于生成相应的近似金字塔，包 括原始图像，而没有误差。
数字图像处理A-第7章 小波和多分辨率 处理
第七章：小波和多分辨率处理
小波变换是基于具有变化的频率和有限持续时 间的小型波进行的。它是多分辨率理论的分析 基础。
多分辨率理论将多种学科的技术有效地统一在 一起，其优势很明显—某种分辨率下所无法发 现的特性在另一种分辨率下将很容易被发现。
本章将从多分辨率的角度解释小波变换。
局部直方图的变化；（b）~（d）几种不同的近似 《数字图像处理A》
主要内容
7.1 背景 7.2 多分辨率展开
级数展开 尺度函数 小波函数
7.3 一维小波变换 7.4 快速小波变换 7.5 二维小波变换 7.6 小波包
《数字图像处理A》
7.2 多分辨率展开
图像金字塔、子带编码和哈尔变换，在数学 理论多分辨率分析中扮演了重要角色。
《数字图像处理A》
7.1.3 哈尔变换
哈尔变换(Haar)是与多分辨率分析有关的图 像处理手段之一。
哈尔变换本身是可分离的，也是对称的，可 以用下述矩阵形式表达：
T=HFHT
(7.1.15)
其中，F是一个N×N图像矩阵，H是N×N 变换矩阵，T是N×N变换的结果
《数字图像处理A》
7.1.3 哈尔变换
1 N
z [0,1]
(7.1.16)
《数字图像处理A》
7.1.3 哈尔变换
且
hk (z) hpq (z)
2p 2
1 2 p 2
N
0
(q-1)/2 p z (q 0.5) / 2 p (q-0.5)/2 p z q / 2 p
otherwise, z [0,1]
(7.1.17)
N×N哈尔变换矩阵的第i行包含了元素hi(z) ，其中z=0/N，1/N，2/N，…，(N-1)/N
《数字图像处理A》
7.1.2 子带编码
图7.6(a)显示了两段子带编译码系统的基本部分
a
f
综合滤波器组 分析滤波器组
b
低频带
高频带
图7.6 （a）一个二带子带编码和解码系统；（
《数字图像处理A》
b）频谱分裂属性
7.1.2 子带编码
因此，FIR综合滤波器是分析滤波器的交叉调制的 副本，有且仅有一个符号相反。为完美重构，综合 滤波器和分析滤波器的冲激响应必须按如下两种方 式之一联系起来：
第j级的大小为2j×2j
（0j J）
共有J+1级，但是通常 第J-1级 我们截短到P＋1级，其
第J级（底部） 中1 PJ
图7.2 (a) 一个图像金字塔 《数字图像处理A》
7.1.1 图像金字塔
如图7.2(b) 框图所表明的，近似值和预测残差 金字塔都是以一种迭代的方式进行计算。
下采样器（行和列） 近似滤波器
k
展开集合的闭合跨度，表示为：
(7.2.1)
V Spkan{k (x)}
(7.2.2)
k ~k (x), f (x) ~k*(x) f (x)dx
《数字图像处理A》
(7.2.3)
7.2.1 级数展开
由于展开集合的正交性，该计算可以是3种 可能形式中的一种。
情况1：如果该展开函数构成了V的一个正
第0级（顶点） 第1级 第2级
第J-1级
第J级（底部）
图7.2 (a) 一个图像金字塔 《数字图像处理A》
7.1.1 图像金字塔
金字塔的底部是待处理图像的高分辨率表示， 而顶部是低分辨率的近似。当向金字塔的上层 移动时，尺寸和分辨率就降低。
第0级（顶点） 第1级 第2级
基础级J的大小为N×N
（J=log2N） 顶点级0的大小为1×1
图7.5显示了6个功能上相关的滤波器的冲激
响应。
h1(n)
h2(n)=-h1(n)
h3(n)=h1(-n)
h4.(n)=h1(K-1-n)
h5(n)=(-1)nh1(n)
h6(n)=(-1)nh1(K-1-n)
图7.5 6个功能上相关的滤波器的冲激响应:(a)参考响应;（b） 符号的反转;（c~d）顺序反转; (e)调制;（f）顺序反转和调制
第j-1级近似
上采样器（行和列）
插值滤波器
预测
第j级输入 图像
第j级预测 残差
图7.2 (b) 创建近似和预测残差金字塔的一个简单系统
《数字图像处理A》
7.1.1 图像金字塔
传递由3个连续步骤组成：
1.计算输入图像减少的分辨率近似值。通过对输入 进行滤波并以2为步长进行抽样实现(即子抽样)。没 有滤波器，在金字塔的上一层混淆变得显著，子抽 样点对所采取的区域没有很好的代表性。
。
《数字图像处理A》
7.1.3 哈尔变换
2×2哈尔变换矩阵H4
H2
1 1 2 1
1 1
(7.1.18)
它的基函数仅定义了2抽头FIR滤波器族
，可满足表7.1中第一行第一列的QMF
滤波器原型的规范。
相应QMF分析滤波器h0(n)和h1(n)的系数 分别是矩阵H2的第一行和第二行的元素。
《数字图像处理A》
《数字图像处理A》
7.1 背景
从数学的观点看，图像是一个亮度值的二维矩阵，像 边界和对比强烈区域那样的突变特性的不同组合会产生 统计值的局部变化。如图7.1所示。
图7.1 一幅自然图像和它的局部直方图变化 《数字图像处理A》
7.1.1 图像金字塔
图像金字塔是以多分辨率来解释图像的一种有 效但概念简单的结构。
子带来完成
《数字图像处理A》
7.1.2 子带编码
因为分解和重建是借助数字滤波器实现的，所以 我们的讨论从数字信号处理（DSP）和数字信号 滤波的简介开始。图7.4为简单数字滤波器。
f(n)
f(n-0)
f(n-1)
单位延迟
单位延迟
f(n-2)
…
f(n-K+1) 单位延迟
h(0)
h(1)
h(2)
h(K-1)
交基，即：
j (x),k (x)
jk
0 1
jk
(7.2.4)
jk
基与它的对偶相等。即：k (x) ~k (x)
7.1.2 子带编码
子带编码也是多分辨率相关的重要图像技术 在子带编码中，一幅图像被分解为一系列限带
分量的几何，称为子Fra Baidu bibliotek。
子带可以重组在一起无失真地重建原始图象。 每个子带通过对输入进行带通滤波而得到。 子带带宽小于原始图像带宽，子带可以进行无信息
损失的抽样 原始图象的重建可以通过内插、滤波、和叠加单个
《数字图像处理A》
7.1.2 子带编码
图7.5(c)和(d)中的滤波器h3(n)和h4(n)的顺序 反转形式：
h3(n)=h1(-n) h4.(n)=h1(K-1-n)
(7.1.6) (7.1.7)
滤波器h3(n)是h1(n)关于垂直轴的映像；滤 波 器 h4(n) 是 h1(n) 的 映 像 和 平 移 形 式 。 忽 略 平移两个滤波器的响应相同。图7.5（e）中
例7.2 图7.1中花瓶的4带宽子带编码
图7.8显示了一个8抽头正交滤波器的冲激响应。
《数字图像处理A》
7.1.2 子带编码
使用7.7所示的子带编码系统对图7.1中的花瓶进行 4个子带分离得到图7.9
a
b
c
d
图7.9 (a)近似子带；（b）水平细节子带；（c）垂直细 节子带；（d）对角线细节子带
…
+
+
+
h(0)f(n)+h(1)f(n-1) h(0)f(n)+h(1)f(n-1)+h(2)f(n-2)
《数字图像处理A》
图7.4 （a）数字滤波器
7.1.2 子带编码
数字滤波器由延迟单元、乘法器和加法器组 成。从滤波器的顶部开始，延迟单元依次连 接建立输入法延迟形式。
如图7.4(a)中的注释所指出的，输入序列 f(n)的延迟单元输出的K-1延迟序列分别与常 数h(0),h(1),h(k-1)相乘后求和，可产生滤波 后的输出序列：
gi (n), g j (n 2m) (i j) (m) i, j {0,1}
(7.1.13)
这对可完美重建的滤波器族定义了正交性。
《数字图像处理A》
7.1.2 子带编码
除式（7.1.13）外，可以证明正交滤波器满 足如下两个条件：
(7.1.14)
其中，Keven的下标指出滤波器系数值必须是能被2整 除的数。正如7.14指出的那样，综合滤波器g1通过顺 序反转和调制与g0建立联系。
(7.1.3)
其中，* 表示卷积，乘数K 称为滤波系数。
《数字图像处理A》
7.1.2 子带编码
如果输入到图7.4(a)的滤波器是图7.4（b）和 4.2.3节中的离散单位冲激，则式（7.1.3）变 为 (7.1.4)
图7.4 （b）单位离散冲激响应;（c）滤波器的冲激响应 《数字图像处理A》
7.1.2 子带编码
或
(7.1.11)
更有普遍意义的表达式
hi (2n k), g j (k) (i j) (n) i, j {0,1}
(7.1.12)
满足该条件的滤波器组称为具有双正交
《数字图像处理A》
7.1.2 子带编码
表7.0 完美重建滤波器族
(正交镜像滤波器)
(共轭正交滤波器)
在双正交的基础上进一步要求
的滤波器h5(n)由下式定义,称之为h1(n) 的调 制形式：
h5(n)=(-1)nh1(n)
《数字图像处理A》
(7.1.8)
7.1.2 子带编码
最后图7.5（f）显示的序列是h1(n)的顺序 反转形式，它被调制了：
h6(n)=(-1)nh1(K-1-n)
(7.1.9)
综上6个序列说明了这样一个事实，即在规 定两个滤波器之间的关系时，符号反转、顺 序反转和调制有时时合并在一起的。
在多分辨率分析( MRA )中，尺度函数被用 于建立某一函数或图像的一系列近似值，相 邻两近似值之间的近似度相差2倍。
被称为小波的附加函数用于对相邻近似值之 间的差异进行编码。
《数字图像处理A》
7.2.1 级数展开
信号或函数常常可以被很好地分解为一系列 展开函数的线性组合
f (x) kk (x)
《数字图像处理A》
7.1.2 子带编码
表7.0中的一维滤波器也可用于图像处理的二维可分 离滤波器。如图7.7所示。可分离滤波器首先应用于 某一维(如垂直向)，再应用于另一维(如水平向)。
行（沿m）
列（沿n） 列
列 行
列
图7.7 子带图像编码的一个二维4带宽滤波器组 《数字图像处理A》
7.1.2 子带编码
7.1.3 哈尔变换
例如，N=4时，k，q和p值如下：
4×4哈尔变换矩阵H4
1
H4
1
2
4 2
0
《数字图像处理A》
1 1 1
1
1
1
2 0 0
0 2 2
(7.1.19)
7.1.3 哈尔变换
• 例7.3 离散小波变换的哈尔函数
64×64 128×128
256×256
图7.10 (a) 用H2哈尔基函数的离散小波变换，并显示了
图7.2 (b) 创建近似和预测残差金字塔的一个简单系统 《数字图像处理A》
7.1.1 图像金字塔
例7.1 高斯和拉普拉斯金字塔
图7.3显示了图7.1中花瓶的一种可能的近似值和预 测残差金字塔。
（a）
(b)
图7.3 两种图像金字塔及它们的直方图：（a）近似金字塔；（
《数字图像处理A》
b）预测残差金字塔
图7.2 (b) 创建近似和预测残差金字塔的一个简单系统 《数字图像处理A》
7.1.1 图像金字塔
传递由3个连续步骤组成：
2.对上一步的输出进行内插—因子仍为2—并进行 过滤。这将生成与输入等分辨率的预测图像。由于 在步骤1的输出像素之间进行插值运算，插入滤波 器决定了预测值与步骤1的输入之间的近似程度。 如果插入滤波器被忽略了，预测值将是步骤1输出 的内插形式，复制像素的块效应将变得很明显。
《数字图像处理A》
主要内容
7.1 背景 7.2 多分辨率展开 7.3 一维小波变换 7.4 快速小波变换 7.5 二维小波变换 7.6 小波包
《数字图像处理A》
主要内容
7.1 背景
图像金字塔 子带编码 哈尔变换
7.2 多分辨率展开 7.3 一维小波变换 7.4 快速小波变换 7.5 二维小波变换 7.6 小波包
哈尔变换的变换矩阵H包含哈尔基函数hk(z) ，它们定义在连续闭区间z∈[0,1]， k=0,1,2…,N-1，这里N=2n。
为生成H矩阵，定义整数k，即k=2p+q-1(这 里0≤p≤n-1, p=0时，q=0或1，p≠0时， 0≤q≤2p)。
可得哈尔基函数为：
h0 (z) h00 (z)
图7.2 (b) 创建近似和预测残差金字塔的一个简单系统 《数字图像处理A》
7.1.1 图像金字塔
传递由3个连续步骤组成：
3.计算步骤2的预测值和步骤1的输入之间的差异。
以j 级预测残差进行标识的这个差异将用于原始图
像的重建(见例7.1)。在没有量化差异的情况下，预 测残差金字塔可以用于生成相应的近似金字塔，包 括原始图像，而没有误差。
数字图像处理A-第7章 小波和多分辨率 处理
第七章：小波和多分辨率处理
小波变换是基于具有变化的频率和有限持续时 间的小型波进行的。它是多分辨率理论的分析 基础。
多分辨率理论将多种学科的技术有效地统一在 一起，其优势很明显—某种分辨率下所无法发 现的特性在另一种分辨率下将很容易被发现。
本章将从多分辨率的角度解释小波变换。
局部直方图的变化；（b）~（d）几种不同的近似 《数字图像处理A》
主要内容
7.1 背景 7.2 多分辨率展开
级数展开 尺度函数 小波函数
7.3 一维小波变换 7.4 快速小波变换 7.5 二维小波变换 7.6 小波包
《数字图像处理A》
7.2 多分辨率展开
图像金字塔、子带编码和哈尔变换，在数学 理论多分辨率分析中扮演了重要角色。
《数字图像处理A》
7.1.3 哈尔变换
哈尔变换(Haar)是与多分辨率分析有关的图 像处理手段之一。
哈尔变换本身是可分离的，也是对称的，可 以用下述矩阵形式表达：
T=HFHT
(7.1.15)
其中，F是一个N×N图像矩阵，H是N×N 变换矩阵，T是N×N变换的结果
《数字图像处理A》
7.1.3 哈尔变换
1 N
z [0,1]
(7.1.16)
《数字图像处理A》
7.1.3 哈尔变换
且
hk (z) hpq (z)
2p 2
1 2 p 2
N
0
(q-1)/2 p z (q 0.5) / 2 p (q-0.5)/2 p z q / 2 p
otherwise, z [0,1]
(7.1.17)
N×N哈尔变换矩阵的第i行包含了元素hi(z) ，其中z=0/N，1/N，2/N，…，(N-1)/N
《数字图像处理A》
7.1.2 子带编码
图7.6(a)显示了两段子带编译码系统的基本部分
a
f
综合滤波器组 分析滤波器组
b
低频带
高频带
图7.6 （a）一个二带子带编码和解码系统；（
《数字图像处理A》
b）频谱分裂属性
7.1.2 子带编码
因此，FIR综合滤波器是分析滤波器的交叉调制的 副本，有且仅有一个符号相反。为完美重构，综合 滤波器和分析滤波器的冲激响应必须按如下两种方 式之一联系起来：
第j级的大小为2j×2j
（0j J）
共有J+1级，但是通常 第J-1级 我们截短到P＋1级，其
第J级（底部） 中1 PJ
图7.2 (a) 一个图像金字塔 《数字图像处理A》
7.1.1 图像金字塔
如图7.2(b) 框图所表明的，近似值和预测残差 金字塔都是以一种迭代的方式进行计算。
下采样器（行和列） 近似滤波器
k
展开集合的闭合跨度，表示为：
(7.2.1)
V Spkan{k (x)}
(7.2.2)
k ~k (x), f (x) ~k*(x) f (x)dx
《数字图像处理A》
(7.2.3)
7.2.1 级数展开
由于展开集合的正交性，该计算可以是3种 可能形式中的一种。
情况1：如果该展开函数构成了V的一个正
第0级（顶点） 第1级 第2级
第J-1级
第J级（底部）
图7.2 (a) 一个图像金字塔 《数字图像处理A》
7.1.1 图像金字塔
金字塔的底部是待处理图像的高分辨率表示， 而顶部是低分辨率的近似。当向金字塔的上层 移动时，尺寸和分辨率就降低。
第0级（顶点） 第1级 第2级
基础级J的大小为N×N
（J=log2N） 顶点级0的大小为1×1
图7.5显示了6个功能上相关的滤波器的冲激
响应。
h1(n)
h2(n)=-h1(n)
h3(n)=h1(-n)
h4.(n)=h1(K-1-n)
h5(n)=(-1)nh1(n)
h6(n)=(-1)nh1(K-1-n)
图7.5 6个功能上相关的滤波器的冲激响应:(a)参考响应;（b） 符号的反转;（c~d）顺序反转; (e)调制;（f）顺序反转和调制
第j-1级近似
上采样器（行和列）
插值滤波器
预测
第j级输入 图像
第j级预测 残差
图7.2 (b) 创建近似和预测残差金字塔的一个简单系统
《数字图像处理A》
7.1.1 图像金字塔
传递由3个连续步骤组成：
1.计算输入图像减少的分辨率近似值。通过对输入 进行滤波并以2为步长进行抽样实现(即子抽样)。没 有滤波器，在金字塔的上一层混淆变得显著，子抽 样点对所采取的区域没有很好的代表性。
。
《数字图像处理A》
7.1.3 哈尔变换
2×2哈尔变换矩阵H4
H2
1 1 2 1
1 1
(7.1.18)
它的基函数仅定义了2抽头FIR滤波器族
，可满足表7.1中第一行第一列的QMF
滤波器原型的规范。
相应QMF分析滤波器h0(n)和h1(n)的系数 分别是矩阵H2的第一行和第二行的元素。
《数字图像处理A》
《数字图像处理A》
7.1 背景
从数学的观点看，图像是一个亮度值的二维矩阵，像 边界和对比强烈区域那样的突变特性的不同组合会产生 统计值的局部变化。如图7.1所示。
图7.1 一幅自然图像和它的局部直方图变化 《数字图像处理A》
7.1.1 图像金字塔
图像金字塔是以多分辨率来解释图像的一种有 效但概念简单的结构。
子带来完成
《数字图像处理A》
7.1.2 子带编码
因为分解和重建是借助数字滤波器实现的，所以 我们的讨论从数字信号处理（DSP）和数字信号 滤波的简介开始。图7.4为简单数字滤波器。
f(n)
f(n-0)
f(n-1)
单位延迟
单位延迟
f(n-2)
…
f(n-K+1) 单位延迟
h(0)
h(1)
h(2)
h(K-1)
交基，即：
j (x),k (x)
jk
0 1
jk
(7.2.4)
jk
基与它的对偶相等。即：k (x) ~k (x)
7.1.2 子带编码
子带编码也是多分辨率相关的重要图像技术 在子带编码中，一幅图像被分解为一系列限带
分量的几何，称为子Fra Baidu bibliotek。
子带可以重组在一起无失真地重建原始图象。 每个子带通过对输入进行带通滤波而得到。 子带带宽小于原始图像带宽，子带可以进行无信息
损失的抽样 原始图象的重建可以通过内插、滤波、和叠加单个
《数字图像处理A》
7.1.2 子带编码
图7.5(c)和(d)中的滤波器h3(n)和h4(n)的顺序 反转形式：
h3(n)=h1(-n) h4.(n)=h1(K-1-n)
(7.1.6) (7.1.7)
滤波器h3(n)是h1(n)关于垂直轴的映像；滤 波 器 h4(n) 是 h1(n) 的 映 像 和 平 移 形 式 。 忽 略 平移两个滤波器的响应相同。图7.5（e）中
例7.2 图7.1中花瓶的4带宽子带编码
图7.8显示了一个8抽头正交滤波器的冲激响应。
《数字图像处理A》
7.1.2 子带编码
使用7.7所示的子带编码系统对图7.1中的花瓶进行 4个子带分离得到图7.9
a
b
c
d
图7.9 (a)近似子带；（b）水平细节子带；（c）垂直细 节子带；（d）对角线细节子带
…
+
+
+
h(0)f(n)+h(1)f(n-1) h(0)f(n)+h(1)f(n-1)+h(2)f(n-2)
《数字图像处理A》
图7.4 （a）数字滤波器
7.1.2 子带编码
数字滤波器由延迟单元、乘法器和加法器组 成。从滤波器的顶部开始，延迟单元依次连 接建立输入法延迟形式。
如图7.4(a)中的注释所指出的，输入序列 f(n)的延迟单元输出的K-1延迟序列分别与常 数h(0),h(1),h(k-1)相乘后求和，可产生滤波 后的输出序列：
gi (n), g j (n 2m) (i j) (m) i, j {0,1}
(7.1.13)
这对可完美重建的滤波器族定义了正交性。
《数字图像处理A》
7.1.2 子带编码
除式（7.1.13）外，可以证明正交滤波器满 足如下两个条件：
(7.1.14)
其中，Keven的下标指出滤波器系数值必须是能被2整 除的数。正如7.14指出的那样，综合滤波器g1通过顺 序反转和调制与g0建立联系。
(7.1.3)
其中，* 表示卷积，乘数K 称为滤波系数。
《数字图像处理A》
7.1.2 子带编码
如果输入到图7.4(a)的滤波器是图7.4（b）和 4.2.3节中的离散单位冲激，则式（7.1.3）变 为 (7.1.4)
图7.4 （b）单位离散冲激响应;（c）滤波器的冲激响应 《数字图像处理A》
7.1.2 子带编码