第三章低层预处理.ppt

无校正CCD图像效果示意图 不同? 值的效果示意图 不准确?值的效果示意图 相关与卷积运算 相关（加权平均） 给定图象f(x,y)大小N×N，模板T(i, j)大小m×m（m为奇数），常用的相关运算定义为: 使模板中心T((m-1)/2,(m-1)/2)与f(x,y)对应，则： 最优的边缘检测---Canny边缘检测方法 设二维高斯函数： 同态滤波器 同态滤波 同态滤波 同态滤波 同态滤波确定H(u,v) 三、逆滤波复原 逆滤波复原 用退化函数H除退化图像的傅立叶变换（G(u,v)）来计算原始图像的傅立叶变换估计 对上述式子求傅立叶逆变换就得到复原后的图像 逆滤波复原 * 例： 1 2 1 4 3 1 2 2 3 4 5 7 6 8 9 5 7 6 8 8 5 6 7 8 9 1 2 1 4 3 1 2 2 3 4 5 7 6 8 9 5 7 6 8 8 5 6 7 8 9 2 3 4 5 6 6 6 7 8 * * * 3.4 图像重建(Image Reconstruction) 一、引言 图像重建的定义 　　根据物体横剖面的投影数据，经过计算机处理后，得到物体横剖面的图像，这种从数据到图像的处理技术，称为图像重建。 图像重建的方法 直接法：代数迭代法； 间接法：傅立叶变换法、滤波逆投影法。 * 图像重建的发展 1895年，德国物理学家伦琴首先发现了X射线； 1917年奥地利数学家雷顿(Radon)提出了由投影重建图像的解析公式， Radon公式； 50年代，图像重建的数学技术开始独立地发展，并应用于许多领域的科研工作； 60年代后期，G.Hounsfield指出：可以用一定的数学方法由一系列X射线测量值来重建物体内部的结构； 1971年，第一台X射线CT扫描仪配置于英格兰的一家医院； 1975年5月，在加拿大蒙特利尔召开的第一届国际专题讨论会上，正式命名为“计算机断层摄影术”，即Computerrized Tomography(CT); 1975年，R.Ledley制成既可用于头部又可用于全身的CT机； 1979年10月11日，A.M.Cormack和G.Hounsfield共同获得诺贝尔生理医学浆； * 二、图像重建的基本原理 Lambert-Beer定理 X射线源 (x,y) 检测器 I0 I * 原理 μ11 μ12 . . . μ1n μ21 μ22 . . . μ2N . . . . . . . . . . . . . . . . . . . . . . . . μN1 μN2 . . . μNN X射线源 检测器 d * * 三、傅立叶变换重建 基本思想 一个三维（或二维）物体，它的二维（或一维）投影的傅立叶变换恰与此物体的傅立叶变换的主体部分相等。通过将投影进行旋转和傅立叶变换，可以首先构造整个傅立叶平面，然后，只需再通过傅立叶反变换，就可以得到重建后的物体。 * 雷顿变换 x s t y f(x,y) * =0时投影的傅立叶变换 u F(u, ) v * 时投影的傅立叶变换 * u w v * * * 2 T 1 T 算子名 特点 Roberts 边缘定位准，对噪声敏感 Prewitt 平均、微分对噪声有抑制作用 Sobel 加权平均边宽2象素 Isotropic Sobel 权值反比于邻点与中心点的距离，检测沿不同方向边缘时梯度幅度一致 * 图6-40 一阶微分算子的效果 （b）原图 （b）Robert算子 （c）Sobel算子 （d）Priwitt算子 * 二阶拉普拉斯算子 定义Laplacian算子为 Laplacian是各向同性（isotropic）的微分算子。 * 零交叉边缘检测 Marr和Hildreth提出的拉普拉斯边缘检 测算子被誉为最佳边缘检测器之一。 * * 在某一方向的二阶导数为： 令: 边缘方向为： 边缘强度为： * 高通滤波法 由于图象中的边缘及急剧变化部分与高频分量有关，所以利用高通滤波器衰减图象信号中的低频分量，就会相对强调其高频分量，从而加强了图象的边缘及急剧变化部分，达到图象锐化的目的。 1）理想高通滤波器 * * 图6-11 三种典型的高通滤波器，从上到下依次为理想高通滤波器、巴特沃思高通滤波器、指数滤波器。 * 高通滤波器效果1 图6-12 理想高通滤波效果，从左到右依次为D0=15，30，80。越小，振铃效应越明显 图6-13 BHPF效果，从左到右依次为D0=15，30，80。比IHPF的结果平滑得多 * 高通滤波器效果2