第五章-医学图像描述new.pptVIP

医学影像图像处理 生物医学工程系 5.1 图像描述 2. 内容： 1）内部特性：灰度、颜色、纹理 2）外部特性：形状 3）不同区域间关系 4）目标物特征：灰度、纹理、几何形状特性等 3. 标准：精确、简练、受干扰小 像素间的基本关系 利用连通概念对图像中对象间关系的描述 5.2图像的几何特征 2. 方向 也就是说，要找出一条直线，使下式定义的E值最小： (1) 当把图像中的像素看作单位面积小方块时，则图像中的区域和背景均由小方块组成。区域的周长即为区域和背景缝隙的长度和，此时边界用隙码表示。因此，求周长就是计算隙码的长度。 像素计数面积 最简单的(未校准的)面积计算方法是统计边界内部(也包括边界上)的像素的数目。求出域边界内像素点的总和即可，计算公式如下： 5.3 图像的形状特征 4.不变矩 不变矩及其组合具备了好的形状特征应具有的某些性质， 但是要区别相似的形体需要一个很大的特征集。这样所产生的高维分类器对噪声和类内变化十分敏感。 （一）链码定义 1. 链码定义：一系列具有特定长度和方向的相连的直线段。 2. 常用两类 3. 起点归一化链码 (1) 归一化链码 若起点不同，则同一边界的链码就不同，为此引 入起点归一化链码。 方法：任选一起点S 求出原链码，将链码看作 由各方向数构成的自然数，将该码按一个方向循 环，使其构成的自然数的值最小，此时形成（起点 ）归一化链码，也称规格化链码。 （2）差分码和归一化差分码 目标旋转则链码会发生变化。可利用链码的一阶差分来重新构造一个表示原链码各段之间方向变化的新序列，这相当于把链码进行旋转归一化 ① 差分码(微分码) 通过计算相邻两个元素方向变化（逆时针方向）的数字得到。即前一个数字变化到后一个数字需要经过的步数 ② 归一化差分码：对差分码进行起点归一化(值最小的差分码) ③ 特点：具有旋转归一化（不变性），也具有平移不变性。 二、边界的形状数表示 1. 定义：形状数是基于链码的一种边界形状描述 符。它定义为归一化差分码，即取值最小的差分码。 2. 举例（8 链码表示） 三、边缘线（边界）的链码跟踪——以8 链码为例 1. 方法：从起点开始，按顺时钟方向向左看 跟踪，可得到与起点8 连通的最长边缘曲线；若 最终能回到起点，则为一闭合曲线，即目标物边 界。对4 链码（4连通）也类似。 2. 举例： 4、边界链码跟踪中面积的求法 5.4 图像纹理描述 5.4.1统计分析方法 空间灰度共生矩阵 灰度共生矩阵就是从N*N的图像f(x,y)的灰度为i的像素出发，统计与距离为 灰度为 j的像素同时出现的概率 用数学表达式则为： 1. 0o方向灰度共生矩阵 当 θ=45° 时，dx=1, dy=0，由于所给图像中只有4个灰度级，因此所求得的灰度共生矩阵的大小为 。 2. 45o方向灰度共生矩阵 当 θ=45° 时，dx=1, dy=-1 3. 90o方向灰度共生矩阵 当 θ=90°时，dx=0,dy=-1 4. 135o方向灰度共生矩阵 当θ=135° 时， dx=-1,dy=-1 灰度共生矩阵计算结果 5.4.3频谱法 该方法主要基于滤波器理论。包含了傅里叶变换法，Gabor变换法，小波变换法。 傅里叶变换法：只能完成图像的频率分解，因而获得的信息不是很充分。 Gabor变换法： Gabor函数可以捕捉到相当多的纹理信息，且具有极佳的空间/频域联合分辨率因此在实际中获得了较广泛的应用。 小波变换法：该方法大体分金字塔形小波变换法和树形小波变换法。 Gabor滤波 Gabor变换属于加窗傅立叶变换，Gabor函数可以在频域不同尺度、不同方向上提取相关的特征。 Gabor函数与人眼的生物作用相仿，所以经常用作纹理识别上，并取得了较好的效果 基于Gabor滤波器的纹理特征提取 将输入图像分成小的图像块； 选择频率与方向（ v, u），组成滤波器组； 滤波器组与每个图像块做空域卷积； 提取特征（均值与方差）。 Gabor滤波 5.4.4图像纹理分析的应用 超声图像纹理分析判断HIFU凝固性坏死 对边界的离散傅立叶变换表达，可以作为定量描述边界形状的基础。采用傅立叶描述的一个优点是将二维问题简化为一维问题。 将x-y平面中的曲线段转化为一维函数f(r)(在r-f(r)平面上)， 也可将x-y平面中的曲线段转化为复平面上的一个序列。 具体就是将x-y平