图像工程chapv.pptVIP

采用傅立叶变换来计算余弦变换的步骤： 1)把f(x)延拓成 长度为2N 2) 求 的2N点的FFT 3）对u各项乘上对应的因子 4) 取实部，并乘上因子 Matlab函数 1) 余弦变换函数：g=dct2(f) 2)逆变换函数：f=idct2(g); 例如： f1 = imread(f:\hdl\images_ch03\fig0318.tif); f=dct2(f1); imshow(f); figure,imshow(fftshift(abs(f))); g=idct2(f); figure,imshow(g,[0 256]); 4.3 低通和高通滤波器 {低通滤波和高通滤波的目的或功能相反，但低通滤波器和高通滤波器又具有对偶性} 4.3.1 低通滤波器 4.3.2 高通滤波器 4.3.1 低通滤波器 低通滤波是要保留图象中的低频分量而除去高频分量 图象中的边缘和噪声都对应图象傅里叶变换中的高频部分，所以如要在频域中消弱其影响就要设法减弱这部分频率的分量 根据频域增强技术的原理，需要选择一个合适的H(u, v)以得到消弱F(u, v)高频分量的G(u, v) 4.3.1 低通滤波器 1. 理想低通滤波器 理想是指小于（截断/截止频率）D0的频率可以完全不受影响地通过滤波器，而大于D0的频率则完全通不过 4.3.1 低通滤波器 1. 理想低通滤波器 H(u, v)：转移 / 滤波函数 D0：截断频率（非负整数） D(u, v)是从点(u, v)到频率平面原点的距离 D(u, v) = (u2 +v2)1/2 4.3.1 低通滤波器 3. 巴特沃斯低通滤波器 减少振铃效应，高低频率间的过渡比较光滑 阶为 n 物理上可实现（理想低通滤波器在数学上定义得很清楚，在计算机模拟中也可实现，但在截断频率处直上直下的理想低通滤波器是不能用实际的电子器件实现的） 4.3.1 低通滤波器 3. 巴特沃斯低通滤波器 截断频率 使H最大值降到 某个百分比（常取 50％）的频率 在D(u, v) = D0时 H(u, v) = 1/2 H(u, v) = 1/21/2 4.3.1 低通滤波器 4. 其他低通滤波器（图4.3.9） 梯形 指数 n=1 4.3.2 高通滤波器 1. 理想高通滤波器 形状与低通滤波器的形状正好相反 4.3.2 高通滤波器 2. 巴特沃斯高通滤波器 形状与巴特沃斯低通滤波器的形状正好相反 截断频率 使H值上升到最大值 某个百分比（常取50％） 的频率 H(u, v) = 1/2 4.3.2 高通滤波器 3. 高频增强滤波器 傅里叶变换：G(u, v) = H(u, v)F(u,v) 高频增强转移函数：He(u, v) = k?H(u, v) + c 高频增强输出图的傅里叶变换： Ge(u, v) = k ? G(u, v) + c ? F(u, v) 反变换回去： ge(x, y) = k ? g(x, y) + c ? f(x, y) 高通 原图 4.3.2 高通滤波器 4. 高频提升滤波器 用原始图减去低通图得到高通滤波器的效果 把原始图乘以一个放大系数A再减去低通图就可构成高频提升（high-boost）滤波器 高通滤波器：A = 1 高频增强滤波器：A = ？ 4.3.2 高通滤波器 5. 其他高通滤波器 梯形 指数 4.4 带通和带阻滤波器 1. 带阻滤波器 阻止一定频率范围 （允许其它频率范围） 如果这个频率范围的下限是0（上限不为?），则带 阻滤波器成为高通滤波器 如果这个频率范围的上限是?（下限不为0），则带 阻滤波器成为低通滤波器 4.4 带通和带阻滤波器 1. 带阻滤波器 （以原点为中心） 巴特沃斯带阻滤波器 4.4 带通和带阻滤波器 2. 带通滤波器 巴特沃斯带通滤波器 4.4 带通和带阻滤波器 陷波滤波器 陷波滤波器必须两两对称地工作 4.5 同态滤波器 {同态滤波基于一个简单的亮度成象模型} 4.5.1 亮度成象模型 4.5.2 同态滤波增强 4.5.1 亮度成象模型 图象成象模型 2-D亮度函数：f (x, y) 亮度是能量的量度，一定不为零且为有限值 (1) 入射到可见场景上的光量 (2) 场景中目标对入射光反射的比率 照度成分 i(x, y) 和反射成分 r(x, y) 4.5.1 亮度成象模型 图象成象模型 f(x, y)与 i(x, y)和 r(x, y)都成正比 i(x, y)的值是由光源决定的 r(x, y)的值是由