主权项 |
1.一种基于尺度空间分解与重构的X射线图像均衡显示处理方法,其特征在于,该方法具体的步骤是:a.归一化处理将输入的X射线图像通过归一化处理及对数变换得到F0(x,y);b.尺度空间分解定义算子Mt为滤波算子,其尺度参数为t,则用尺度为t0的滤波算子Mt0对F0(x,y)进行滤波,则将滤波后的图像F1(x,y)表示为:<math-cwu><![CDATA[<math> <mrow> <msub> <mi>F</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>M</mi> <msub> <mi>t</mi> <mn>0</mn> </msub> </msub> <msub> <mi>F</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow></math>]]></math-cwu>D0(x,y)=F0(x,y)-F1(x,y)重复上述步骤,用尺度为tN-1的滤波算子MtN-1对DN-2(x,y)进行平滑,得到FN(x,y),DN-2(x,y)与FN(x,y)之差即DN-1(x,y),如下所示:<math-cwu><![CDATA[<math> <mrow> <msub> <mi>F</mi> <mi>N</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>M</mi> <msub> <mi>t</mi> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </msub> <msub> <mi>D</mi> <mrow> <mi>N</mi> <mo>-</mo> <mn>2</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow></math>]]></math-cwu>DN-1(x,y)=DN-2(x,y)-FN(x,y)因此有:F0(x,y)=D0(x,y)+F1(x,y)D0(x,y)=D1(x,y)+F2(x,y)D1(x,y)=D2(x,y)+F3(x,y)DN-2(x,y)=DN-1(x,y)+FN(x,y)所以,输入的F0(x,y)表示为:<math-cwu><![CDATA[<math> <mrow> <msub> <mi>F</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>F</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>D</mi> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow></math>]]></math-cwu>c.重建每一次分解过程都看作将待分解图像Di-1(x,y)分解为Di-1(x,y)的大尺度部分Di(x,y)与小尺度部分Fi+1(x,y)的过程,保留每次分解获得的Fi+1(x,y),以及DN-1(x,y),即可恢复原始图像F0(x,y),恢复得到的图像表示为F0′(x,y),因此,加上均衡效果的重建过程表示为:<math-cwu><![CDATA[<math> <mrow> <msubsup> <mi>F</mi> <mn>0</mn> <mo>′</mo> </msubsup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>F</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <msub> <mrow> <mo>×</mo> <mi>α</mi> </mrow> <mi>i</mi> </msub> <mo>+</mo> <msub> <mi>D</mi> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <msub> <mrow> <mo>×</mo> <mi>α</mi> </mrow> <mrow> <mi>N</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> </mrow></math>]]></math-cwu>其中,αi为均衡系数,i=1,2,…,N+1;d.反归一化及反对数变换最后将F0′(x,y)通过反归一化及反对数变换,及灰度范围重映射得到均衡处理后的图像。 |