数字图像处理第10章图像复原

（8）携带遥感仪器的飞机或 在研究图像原因的基础上，以图像为依据，根据一定的先验知识，建立一个模型，然后用相反的运算，以恢复原始 n(x,f(x, g输入图像f(x,y)经过一个 算子H(x,y)后产 图像g(x,y)可以表示为：g(x,y)=H[f(x, g(x,y)=H[f(x,y)]+n(x, 的图像是由成像系统的加上额若已知H(x,y)和n(x,y)，图像复原是在退化图像的基础上，作逆运算，得到f(x,y)由于存在可能导致图像复原的 入图像f(x,y)经 后的输出为g(x,y)：g(x,y)H[f(x,y)]H

f(,)H[(x,yf(,)h(x,y 又被称 g(x,y) f(,)h(x,y)ddn(x, f(x,y)*h(x,y)n(x,

G(u,v)F(u,v)H(u,v)N(u, G(u,v)、F(u,v)、N(u,v)分别是g(x,y)、f(x,y)、n(x,y) H(u,v)是h(x,y) fe(x,y)

f(x,

0xA1且0yB1 h(x,y)h(x,

0xC1且0yD10.1.210.1.2 M Nge(x,y) fe(m,n)he(xm,y

HH H H

H H H

2HH

H0

H3 HM- H0 he(j, he(j,N he(j,N H

he(j,N he(j,he(he(j,

he(j,N

he(j,N ge(x,y) fe(m,n)he(xm,yn)ne(x, g=Hf 图像g(x,y)、 噪声分布n(x,y)，就可以得到原始图像f的估计。 1,ifLih(i) 其他 h(i,j)K

i2j 这是由于相机聚焦确引起的，虽然不聚焦由许多参数决定，如相究中为了简单起见，我们用下列函数表示聚焦引起的模糊h(i,j)R

i2j2其他 1,ifLi,jh(i,j) 其他 模糊信噪比(BSNRtheBlurredSignal-to-NoiseRatio [g(i,j)g(i,j)]2BSNR10log10

i,2 2n2叠加噪声n(ij)的方差g(i,j)y(i,jn(i,j)ng(i,jE(g)是g(i,j)ISNR（theImprovementinISNR10

[f(i,j)y(i,j)]2i, 10[f(i,j)fˆ(i,j)]2 i, f(ij)－－原始图像y(i,j)－－降质图像ˆ(i11.2 n=当对n的统计特性不确定时，希望对原始图像f的估计满足这样的条件，使Hˆ在最小二乘意义上近似于g。希望找到一个ˆ使得噪声项的范数||n||2nT

ˆ(HTH)1HTgH1 ˆ(u,v)G(u, H(u,对传递函数进行修正，以降低由于传递函数而造成的恢复不稳2222J(ˆ)

Q

gHˆ

n 日乘子,f的最佳估ˆ(HTHQTQ)1HT

Rf=E{ff Rn= ˆ(HTHR1R)1HT

ˆ(u,v)

|H(u,v)|2[S(u,v)/ (u,v)]G(u,v)

|H(u,v)

|H(u,v)2|[S(u,v)/ (u,v)] , J=deconvwnr(IPSFNSR)J=deconvwnr(I,PSF,NCORR,为噪声信号功率比，默认值为0；NCORR,PSF=

使用较“长”的PSF复原（c）使用较“陡峭”的 g(x,y)S[b(x,y)]n(x,b(x,y)h(x,;y,)f(,通过指数变换得到复原图像(x,y)g(x,y)=i(x,y) logg(x,y)=logi(x,y)+log ˆ(x, =10l(x,y)*[logi(x,y)+log 1 |L(u,v)||H(u,v)

S(u,v)/S(u,v) 。 [J,PSF]=deconvblind(I,INITPSF,NUMIT,DAMPAR,WEIGHT,。 这类畸变一