数字图象处理-Read.ppt

第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 第5章 图象恢复和重建 * * 5.1 概述 图象增强：改善图象的视觉效果 图象恢复：分析图象退化的原因，建立模型，恢复图象 图象重建：根据投影数据重建图象中物体内部结构 图象恢复分类： 从技术上：有约束、无约束恢复 从策略上：自动处理、交互处理 从方法上：空间域、代数域、频域 5.2 退化模型和对角化 一．退化模型： 1．简单的通用图象退化模型： g(x,y)=H[f(x,y)]+n(x,y) 5.2 退化模型和对角化 一．退化模型： 1．简单的通用图象退化模型： 假设n(x,y)=0 ，H的性质： ⑴假设为线性系统： H[k1f1(x,y)+k2f2(x,y)]=k1H[f1(x,y)]+ k2H[f2(x,y)] ⑵k1= k2时， H[f1(x,y)+ f2(x,y)]=H[f1(x,y)]+H[f2(x,y)] ⑶一致性：若f2(x,y)＝0，则H[k1f1(x,y)]=k1H[f1(x,y)] ⑷位置不变性：H[f(x-a,y-b)]=g(x-a,y-b) 2．常见的图象变质： ⑴记录显示时的非线性 ⑵聚焦模糊 5.2 退化模型和对角化 一．退化模型： 2．常见的图象变质： ⑶运动模糊 ⑷随机噪声的叠加 3．建立退化模型 变质图象g(x,y)=f(x,y)*h(x,y)+n(x,y) ⑴1D离散情况： 函数：f(x): x=0…A-1 h(x): x=0…B-1 扩展后：fe(x):x=0…M-1 he(x):x=0…M-1 5.2 退化模型和对角化 一．退化模型： 3．建立退化模型 ⑴1D离散情况： 不考虑噪声时：ge(x)= fe(m)he(x-m) x=0…M-1 写成矩阵形式：g=Hf H= H为循环矩阵 5.2 退化模型和对角化 一．退化模型： 3．建立退化模型 ⑵ 2D离散情况： 扩展后：fe(x,y):0≤x≤M-1,0≤y≤N-1 he(x,y):0≤x≤M-1,0≤y≤N-1 考虑噪声时：ge(x,y)= fe(m,n)he(x-m,y-n)+ne(x,y) 写成矩阵形式：g=Hf+n 块循环矩阵H为： H= 5.2 退化模型和对角化 一．退化模型： 3．建立退化模型 ⑵ 2D离散情况： 其中 二．循环矩阵的对角化 1．循环矩阵的对角化 循环矩阵H的本征矢量： w(k)=[1 exp(j2 k/M) ……exp(j2 (M-1)k/M)]T k=0…M-1 5.2 退化模型和对角化 二．循环矩阵的对角化 1．循环矩阵的对角化 H的本征值： (k)=he(0)+he(M-1)exp(j2 k/M)+……+ he(1)exp(j2 (M-1)k/M) W=[w(0) w(1)……w(M-1)] 其中w(k,l)= exp(j2 kl/M) w-1(k,l)= exp(-j2 kl/M)/M H=WDW-1 其中D(k,k)= (k) 2．块循环矩阵的对角化 w(i,m)= exp(j2 im/M)WN w-1(i,m)= exp(-j2 im/M)WN-1/M wN(k,n)= exp(j2 kn/N) wN-1(k,n)= exp(-j2 kn/N)/N H=WDW-1 5.2 退化模型和对角化 二．循环矩阵的对角化 3．退化模型对角化的效果 ⑴1D情况： g=Hf W-1g=DW-1f G(k)=MH(k)F(k) 即矩阵相乘转换为傅氏变换求取 ⑵2D情况： g=Hf+n W-1g=DW-1f+W-1n G(u,v)=H(u,v)F(u,v