医学图像的运算.pptVIP

3.1 医学图像运算概述 3.2 医学图像的灰度变换 (一)线性灰度变换 当图像成像时曝光不足或过度, 或由于成像设备的非线性和图像记录设备动态范围太窄等因素。都会产生对比度不足的弊病，使图像中的细节分辨不清。这时可将灰度范围线性扩展。 设原始图像f(x,y)灰度范围为[a,b]， 变换后图像g(x,y)灰度范围为[c,d], (二)非线性灰度变换 （1）对数变换 （2）指数变换 部分线性灰度变换 部分线性灰度变换是针对图像的某一感兴趣区域 灰度进行灰度扩展，其余部分进行灰度压缩。它 对感兴趣的灰度范围内的灰度变换是线性变换， 其余部分是非线性灰度变换，在医学图像处理 中，最常见的部分线性灰度变换就是高精度医学 图像的开窗显示，也称窗口技术。 部分线性灰度变换 由于临床诊断的需要，高精度医学图像具有较大 空间分辨率并且像素最高灰度值超过8bit。目前 国内医院中所用的普通电脑显示器不能直接显示 这种高精度医学图像，国外采用专门用于医学图 像的高分辨率数字显示器来解决高精度医学图像 的显示问题。 价格昂贵，不适合我国国情。 部分线性灰度变换 由于临床诊断的需要，高精度医学图像具有较大 空间分辨率并且像素最高灰度值超过8bit。目前 国内医院中所用的普通电脑显示器不能直接显示 这种高精度医学图像，国外采用专门用于医学图 像的高分辨率数字显示器来解决高精度医学图像 的显示问题。 价格昂贵，不适合我国国情。 3.3 医学图像的空间变换 图像的平移 图像的平移非常简单，所用到的是中学学过的直角坐标系的平移变换公式： 首先，在x方向上作线性插值，对上端的两个顶尖进行线性插值得： 水平镜像 垂直镜像 医学图像其它运算 在医学图像处理中，通常利用图像的加法运算把同一景物的多幅图像加起来求平均，以减少图像的随机噪声，也可以为突出某种增强效果的需要把多幅图像有目的的叠加在一起。 加法运算 设有一幅混入噪声的图像g (x, y) 是由原始图像 f (x, y)和噪声图像e (x, y)叠加而成，即： g (x, y) = f (x, y) + e (x, y) 假设各点的噪声是互不相关的，且具有零均值。可以将一系列的图像{gi (x, y) }相加来消除噪声。 图像平均 设M个图像相加求平均得到一幅新图像 即： 图像平均 考虑新图像和噪声图像各自均方差间的关系： 有： Addition: averaging for noise reduction 图像相减 通过减法运算对同一场景在不同时间拍摄的图像进行相减，用来消除前后影像中相同的背景内容，增强前后有所变化的部分。减法运算的这一特点在临床上具有重要的应用价值，数字减影血管造影（DSA）就是利用图像的相减。 二、图像的局部运算和全局运算 （一）图像的局部运算 如果输出图像在（x,y）点处的像素值不仅与输入图像在（ x,y ）点处的像素值有关，还与输入图像在（ x,y ）点处像素的周围像素的值有关，这种运算称为图像的局部运算。一般通过卷积运算获得。 MATLAB实现 imadjust函数来实现图像的对比度调整。 MATLAB实现 说明： 使用imadjust函数时，按照以下两个步骤 ① 绘制直方图，观察灰度范围； ② 将灰度范围转换为0.0~1.0之间的分数。 例： MATLAB实现 缩放变换函数imresize。 MATLAB实现 imtransform函数:对图像进行二维空间变换 MATLAB实现 maketform函数:产生空间变换参数结构体数据TFORM函数。 仿射变换 仿射变换 仿射变换 使用maketform函数可以创建TFORM结构体，参数为希望执行变换的类型和变换矩阵。 % 创建TFORM结构体 tform_translate = maketform(affine, xform); 执行变换 使用imtransform函数执行变换，参数为要变换的图像和TFORM结构体，函数将返回变换后的图像。 MATLAB实现 imrotate函数:对图像进行绕其中心的旋转变换 The basic principle of DSA 相减 DSA X-ray image obtained before and after injection of dye into arteries. Dye affects appearance of blood vessels but not surrounding tissues. S