621常见的退化函数模型
(a)原始图像 (b)运动造成的图像模糊 (c)Wiener滤波去图像模糊 维纳滤波复原实例 I = imread('lena.jpg'); figure(1),imshow(I); LEN = 21; THETA = 11; PSF = fspecial('motion', LEN, THETA); blurred = imfilter(I, PSF, 'conv', 'circular'); figure(2), imshow(blurred) estimated_nsr = noise_var / var(I(:)); wnr= deconvwnr(blurred_noisy, PSF, estimated_nsr); figure(3), imshow(wnr) (a)原始图像 (b)运动模糊图像 (c)维纳滤波复原 6.6.2 约束最小二乘滤波复原 用deconvreg函数实现对模糊图像的约束最小二乘复原的调用格式为: J=deconvreg(I,PSF) J=deconvreg(I,PSF,NP) J=deconvreg(I,PSP,NP,LRANGE) J=deconvreg(I,PSF,NP,LRANGE,REGOP) [J,LAGRA]=deconvreg(I,PSP,…) 使用该函数时应注意:输出图像J可能会受到算法中使用的离散傅里叶变换而出现振铃现象。 为了减少振铃现象的影响,在使用deconvreg 函数前先使用edgetaper 函数处理图像。 例如,I=edgetaper(I,PSF)。 几种复原法的比较 令 则有 这是已知退化模型的傅里叶变换式。 下面以实例介绍图像模糊的模型。 例如,设原图像f(x,y)只在x方向以给定的速度做匀速直线运动,则有 和 当t=T时,图像f(x,y)在水平x方向的移动距离为a。则 若y分量也发生变化 ,则退化函数变为 6.5 图像的几何校正 图像在获取或显示生成过程中,由于成像系统本身具有的几何非线性以及由于视像管摄像机及阴极射线管显示器的扫描偏转系统有一定的非线性,或者摄像时视角的不同,都会使生成的图像产生几何失真或几何畸变。 图像的几何失真实质上也是一种图像退化的过程。 解决图像的几何失真校正的办法包括如下两个步骤: (1)空间变换:对图像平面上的像素进行重新排列以恢复原空间关系。 (2)灰度插值:对空间变换后的像素赋予相应的灰度值以恢复原位置的灰度值。 几何变换是图像处理中一种基本的、常用的图像预处理方法,其主要用途是: 1.实现数字图像的放大、缩小及旋转; 2.实现畸变(畸变原因可以多种多样,如摄影系统或镜头畸变)图像的校正; 3.实现不同来源图像(如航空摄影,卫星遥感,合成孔径雷达等不同来源)的配准; 4.显示和打印图像时的一种图像排版工具; 5.可以使处理后的图像具有多种不同的特殊效果。 6.5.1 几何畸变的描述 图像的几何畸变是指在成像过程中所产生的图像像元的几何位置相对于参照系统(地面实际位置或地形图)发生的挤压、伸展、偏移和扭曲等变形,使图像的几何位置、尺寸、形状、方位等发生改变。 图像中所产生几何畸变大致分为两大类:(1)内部畸变;(2)外部畸变。 比例尺 歪斜 中心移动 扫描非线性 辐射状畸变 扭曲 图像的几何畸变示例 (a)内部畸变 倾斜引起的 高度变化引起的 地形起状引起 地球曲率引起 投影畸变 比例尺误差 的畸变 的畸变 (b)外部畸变 图像的几何畸变实例 (a)扭曲 (b)扫描非线性 6.5.2 图像空间变换 图像空间几何坐标变换以及像素点灰度值的确定这两部分内容是几何校正的基础,几何校正需要两个独立的算法。 一个算法是几何空间变换本身,用它描述每个像素如何从其初始位置移动到终止位置,即每个像素的运动。 同时,还需要另一个算法用于灰度级的插值。