发明名称 一种辐射成像中图像信息增强的方法
摘要 一种辐射成像中图像信息增强的方法,涉及辐射成像中的图像信息处理技术。本发明方法的主要步骤为:在获取的被测物体的辐射图像A(x,y)中查找像素灰度最大值A<SUB>max</SUB>;将图像A(x,y)归一化处理后的图像B(x,y)进行均值滤波处理得到图像C(x,y);将图像B(x,y)减去图像C(x,y)得到的差值与像素点的边缘域值比较得到图像D(x,y),并对其进行非线性对比度增强变换得到图像E(x,y);将确定的增强比例系数α(x,y)和E(x,y)、C(x,y)相乘得到边缘增强图像F(x,y);最终增强图像G(x,y)=(F(x,y)+C(x,y))*A<SUB>max</SUB>。同现有技术相比,本发明方法对图像信息增强处理速度快、算法简单、成像清晰,消除了图像的噪声,满足扫描图像对于暗区域的对比度需要增强更多的要求。
申请公布号 CN100405400C 申请公布日期 2008.07.23
申请号 CN200510077797.6 申请日期 2005.06.27
申请人 同方威视技术股份有限公司;清华大学 发明人 张丽;陈志强;邢宇翔;李元景;刘以农;赵自然;肖永顺
分类号 G06T5/00(2006.01) 主分类号 G06T5/00(2006.01)
代理机构 代理人
主权项 1.一种辐射成像中图像信息增强的方法,它是对包括使用辐射源、机械控制装置、探测装置、数据采集装置及计算机成像装置的辐射成像系统所获得的信号进行处理,其主要步骤为:1)通过辐射成像系统扫描获取被测物体的辐射图像A(x,y);2)在图像A(x,y)中查找像素灰度最大值Amax;3)对图像A(x,y)进行归一化处理,使处理后图像B(x,y)中的像素灰度值分布在(0,1)之间,归一化处理的公式为<math><mrow><mi>B</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>A</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><msub><mi>A</mi><mi>max</mi></msub></mfrac><mo>;</mo></mrow></math> 4)对图像B(x,y)进行均值滤波处理得到图像C(x,y);5)计算图像中每个像素点的相对标准偏差,相对标准偏差的3倍值即为图像中每个像素点的边缘域值,标准偏差的计算公式为<math><mrow><mi>&sigma;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><msqrt><mo>{</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>=</mo><mo>-</mo><mfrac><mi>N</mi><mn>2</mn></mfrac></mrow><mfrac><mi>N</mi><mn>2</mn></mfrac></munderover><msup><mrow><mo>[</mo><mi>A</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><msub><mi>y</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>M</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup><mo>}</mo><mo>/</mo><mrow><mo>(</mo><msup><mi>N</mi><mn>2</mn></msup><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msqrt><mo>,</mo></mrow></math>公式中的M(x,y)为像素(x,y)的像素邻域均值;6)将图像B(x,y)减去图像C(x,y),得到的差值与像素点的边缘域值比较,当差值大于或等于边缘域值时,保存该差值到图像D(x,y)中,否则将D(x,y)对应的像素值设为0;7)对图像D(x,y)进行非线性对比度增强变换得到图像E(x,y);8)确定一个大于1的常数c,得到增强比例系数<math><mrow><mi>&alpha;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>c</mi><msqrt><mfrac><mrow><mi>A</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow></mrow><msub><mi>A</mi><mi>max</mi></msub></mfrac></msqrt></mfrac><mo>;</mo></mrow></math> 9)将α(x,y)、E(x,y)和C(x,y)相乘,得到边缘增强图像F(x,y)=α(x,y)*E(x,y)*C(x,y);10)最终增强图像G(x,y)=(F(x,y)+C(x,y))*Amax。
地址 100083北京市清华同方科技广场A座2907