完美滤波器

如下图所示,第j级为输入图像,其中第j-1级为第j级的尺寸减半的存在,直至为 1\times 1 的大小,这样的模式被称为图像金字塔

image-20231129132706007

设原图像像素点个数为 N^2,则图像金字塔的总像素个数为

N^{2}\left(1+\frac{1}{\left(4\right)^{1}}+\frac{1}{\left(4\right)^{2}}+\cdots+\frac{1}{\left(4\right)^{P}}\right){\leqslant}\frac{4}{3}N^{2}


对于图像金字塔建模,设第j级为图像降低分辨率后的近似图像,这可以视为由第j+1 级图像经过滤波操作和下采样实现后的存在,则第j+1级可以视为第j级经过上采样和插值操作后的存在,即如下图所示:

image-20231129134032195

其中对图像进行上采样操作,索引所对应的值为:

f_{2\uparrow}(n)=\begin{cases}f(n/2),&\quad n\text{为偶数}\\\quad0,&\quad\text{其他}\end{cases}


图像进行下采样操作,索引所对应的值为:

f_{2\downarrow}(n)=f(2n)


上采样可看成是在序列中的每一个样本后插人 0; 下采样可看成是每隔一个样本丢弃一个样本。

image-20231129140742261

设存在输入信号 f(n),其中 h_0(n)h_1(n) 分别为低通与高通滤波器,并输入信号一分为二,如下图所示

image-20231129141039757

并经过下采样得到 f_{lp}(n)f_{hp}(n)。然后经过上采样,与滤波 g_0(n)g_1(n)并将信号 f_{lp}(n)f_{hp}(n)合并得到信号 \hat f(n),若 \hat f(n)f(n)相等,可以称为采用了完美滤波。

存在一个Z变换:

X(z)=\sum_{-\infty}^\infty x(n)z^{-n}


其中离散傅里叶变换是Z变换的一个特殊形式,即 z=e^{jw}

若对z变换采用下采样,则可以得到:

x_{\mathrm{dowm}}(n)=x(2n)\iff X_{\mathrm{dowm}}(z)=\frac{1}{2}[X(z^{\frac{1}{2}})+X(-z^{-\frac{1}{2}})]

证明补充:

对信号 x(n)进行下采样处理,则可以得到 x_{\mathrm{dowm}}(n)=x(2n)

设信号 y(n)=x_{\mathrm{dowm}}(n),则存在Z变换:

\begin{aligned} Y(z) &=\sum_{-\infty}^\infty y(n)z^{-n}\\ &=\sum_{-\infty}^\infty x_{down}(n)z^{-n}\\ &=\sum_{-\infty}^\infty x(2n)z^{-n}\\ \end{aligned}


令t=2n ,则存在

\begin{aligned} Y(z) &=\sum_{-\infty}^\infty x(2n)z^{-n}\\ &=\sum_{-\infty}^\infty x(t)z^{-\frac{t}{2}} \end{aligned}


因为t=2n为偶数,但是此时t为所有整数,所以应该保证t为奇数是为0

所以构造

\begin{aligned} Y(z) &=\frac{1}{2}\left[\sum_{-\infty}^\infty x(t)z^{-\frac{t}{2}} +\sum_{-\infty}^\infty (-1)^{(t)} x(t)z^{-\frac{t}{2}}\right]\\ &=\frac{1}{2}\left[\sum_{-\infty}^\infty x(t)(z^{-\frac{1}{2}})^t +\sum_{-\infty}^\infty (-1)^{(t)} x(t)(z^{-\frac{1}{2}})^t\right]\\ &=\frac{1}{2}\left[\sum_{-\infty}^\infty x(t)(z^{-\frac{1}{2}})^t +\sum_{-\infty}^\infty x(t)((-z)^{-\frac{1}{2}})^t\right]\\ &=\frac{1}{2}[X(z^{\frac{1}{2}})+X(-z^{-\frac{1}{2}})] \end{aligned}



若对z变换采用上采样,则可以得到:

\left.x^\mathrm{up}(n)=\left\{\begin{matrix}x(n/2)&n=0,2,4,\cdots\\0&\text{其他}\end{matrix}\right.\right.\Longleftrightarrow X^\mathrm{up}(z)=X(z^2)


先对信号进行下采样然后进行上采样可得:

\hat{X}(z)=\frac{1}{2}[X(z)+X(-z)]


根据z变换的逆变换可得:

Z^{-1}[X(-z)]=(-1)^{n}x(n)


在根据完美滤波原理:

\begin{aligned} \hat{X}(z)&=\frac{1}{2}[H_0(z)X(z)+H_0(-z)X(-z)]G_0(z)\\ &+\frac{1}{2}[H_1(z)X(z)+H_1(-z)X(-z)]G_1(z)\\ &=\frac{1}{2}[H_0(z)G_0(z)+H_1(z)G_1(z)]X(z)\\ &+\frac{1}{2}[H_0(-z)G_0(z)+H_1(-z)G_1(-z)]X(-z) \end{aligned}


为了实现完美滤波则,应存在

\begin{aligned}H_0(-z)G_0(z)+H_1(-z)G_1(z)&=0\\\\H_0(z)G_0(z)+H_1(z)G_1(z)&=2\end{aligned}


即:

\begin{bmatrix}H_0(z)&H_1(z)\\H_0(-z)&H_1(-z)\end{bmatrix}\times\begin{bmatrix}G_0(z)\\G_1(z)\end{bmatrix}=\begin{bmatrix}2\\0\end{bmatrix}


利用克拉默法则可得

\left.\left[\begin{array}{c}{G_{0}(z)}\\{G_{1}(z)}\\\end{array}\right.\right]=\frac{2}{\det(\mathbf{H}_{m}(z))}\begin{bmatrix}{H_{1}(-z)}\\{-H_{0}(-z)}\\\end{bmatrix}


其中 \det(H_m(z))=\alpha z^{-(2k+1)},忽略时延,并令 \alpha=2,可得

即当 g_0(n)=(-1)^nh_1(n),g_1(n)=(-1)^{n+1}h_0(n)g_{0}(n)=\left(-1\right)^{n+1}h_{1}(n),g_{1}(n)=\left(-1\right)^{n}h_{0}(n)时,成立

其中 g_0g_1分别由 h_1h_0调制得到

因为 H_0(z)G_0(z)+H_1(z)G_1(z)=2

所以存在:

\sum_kg_0(k)h_0(n-k)+(-1)^n\sum_kg_0(k)h_0(n-k)=2\delta(n)


其中 \delta(n)为单位脉冲函数

当n为奇数时将会出现相消得情况,于是可以简化为

\sum_kg_0(k)h_0(2n-k)=<g_0(k),h_0(2n-k)>=\delta(n)


可以看成两个向量的内积,同理存在

\begin{aligned}&<g_1(k),h_1(2n-k)>=\delta(n)\\&<g_0(k),h_1(2n-k)>=0\\&<g_1(k),h_0(2n-k)>=0\end{aligned}


h_{0}(n),h_{1}(n),g_{0}(n)\text{和 }g_{1}(n)满足双正交

<h_i(2n-k),g_j(k)>=\delta(i-j)\delta(n)\quad i,j=\{0,1\}