基于空域运算的快速股骨头X线片图象增强方法

1、第十四章 基于空域运算的快速股骨头X线片图象增强方法14.1 股骨头X线片图象高通滤波实验分析(a) 股骨头X线片原图(b) 拉普拉斯锐化算子锐化结果14.1 股骨头X线片图象高通滤波实验分析(c) Sobel锐化算子1 的锐化结果(d) Sobel锐化算子2 的锐化结果14.1 股骨头X线片图象高通滤波实验分析(e)Butterworth滤波结果(f)直方图均衡化之后的Butterworth滤波结果 14.2 新的滤波算法空域中有许多著名的锐化算子，例如用以下的脉冲响应阵列在空域中利用卷积实现高通滤波。注意到这些脉冲响应阵列中元值之和为1，起到归一化的作用。0101510101H111191

2、1112H1212521213H12, 122, 121, 122, 121 , 1212, 12, 122, 11 , 112, 12, 11, 12, 11 , 1.) 12(11.rrrrrrrrrrrrrrrrrhhhhhhhrhhhhhhhH 14.2 新的滤波算法这一脉冲响应阵列如果能够有效地增强骨小梁检测能力，那么它的速度的提升不仅可以由于是空域运算而加快，甚至可以不进行卷积，因为hij 是相同的，所以图象与脉冲响应阵列的卷积可以用以下的快速算法来实现：BAH121121,1, 1rkrllkrljrkiijhab21,1,1,1,11,211()()(146)(21)ri ji

3、 ji ji ji k rj ri k rj rkbbaaaar 14.2 新的滤波算法1( )( )( )( )( )2f xf xxf xg xb0( )1(0)20 xbxg xbxbbbxb sin(2)sin(2)( )( )( )( ) 122bsbsF sF sF sF sbsbssin(2)( )12bsH sbs 14.2 新的滤波算法注意到H阵列的各元素之和为0，没有归一化，求得图象 之后，我们对所有的灰度值加上128：上式的运算只是对图象B的表示空间的转移，以利于视觉检测，如果对 H作以下改变：,128i ji jbb1,11,21,11,21,211,11,21,21,

4、21221,121,221,121,221,21.1.2.(21).rrrrrrrrrrrrrrrrrhhhhhhhhhHrhhhhh 14.3 锐化结果及分析高通滤波后的图象由于相对去除了低频成分，所以整体较为灰暗，我们对它使用直方图均衡化处理以提高视觉检测能力。(a)股骨头X线片原图(b) r=75时，由式(14-4)与原图象卷积的结果14.3 锐化结果及分析(c)b中的结果进行直方图均衡后的结果 (d)对图14.4(a)采用 卷积增强的结果（r=75）14.4 快速增强源代码void CDibView:OnFastEn() 习 题 十 四1.分析本章的空域运算算子为什么是快速的。2.分析本章的空域算子为什么是高通算子。 14.2 新的滤波算法这一脉冲响应阵列如果能够有效地增强骨小梁检测能力，那么它的速度的提升不仅可以由于是空域运算而加快，甚至可以不进行卷积，因为hij 是相同的，所以图象与脉冲响应阵列的卷积可以用以下的快速算法来实现：BAH121121,1, 1rkrllkrljrkiijhab21,1,1,1,11,211()()(146)(21)ri ji ji ji ji k rj ri k rj rkbbaaaar 14.3 锐化结果及分析高通滤波后的图象由于相对去除了