[发明专利]一种医学X光影像滤波处理的二级修正方法

权利要求书

1.一种医学X光影像滤波处理的二级修正方法，其特征在于该方法具体是：

步骤1、对源影像P(i，j)做超低通滤波，制作超低频分量图像F(i，j)；

步骤2、用超低频分量图像F(i，j)对源影像P(i，j)做完全减影运算，得到超低频高通图像H(i，j)；

步骤3、非线性转换超低频高通图像H(i，j)像素值，压抑弱信号，保留强信号，获得超低频高通非线性像素图L₂(i，j)；

步骤4、在超低频高通非线性像素图L₂(i，j)的引导下对影像做适宜性低通滤波的二级修正；

步骤5、计算频率滤波二级修正的较高频分量H₂(i，j)；

步骤6、对源影像P(i，j)用超低频分量图像F(i，j)做部分减影，再利用超低频高通非线性像素图L₂(i，j)和较高频分量H₂(i，j)做二级修正，获取最终处理后影像P_f(i，j)；

所述的步骤4具体是：

首先，找出超低频高通非线性像素图L₂(i，j)绝对值的最大值h_max；

则在h_max处的低通小波的应用次数为

其中m为较高频分量二级修正深度，取5-20；n₀是一个不为零的基数，n是低通小波的应用次数；

在超低频高通非线性像素图L₂(i，j)值为零处的低通小波的应用次数为

n_min＝n

这里，n_max＜n_min；

然后，用线性插值法决定在超低频高通非线性像素图L₂(i，j)∈[-h_max，h_max]时的低通小波的应用次数n’(i，j)；

在得到n’(i，j)矩阵后，用该矩阵控制低通小波对源影像做适宜性低通滤波，具体是：

(i)做一变量N，初始赋值为1；做两个图像存储空间，用指针F’和F”分别指向这两个存储空间；将源影像P整体拷贝到F’；

(ii)对所有n’(i，j) ≥ N的像素，把F’的相应位置的像素作为输入做低通小波L滤波，将结果赋值于F”(i，j)；对所有n’(i，j)＜N的像素位置，直接把F’(i，j)的像素值拷贝到F”(i，j)；

(iii)N加1，F’和F”存储空间的指针互换；如果N≤n_min，重复步骤(ii)；当N＞n_min时，流程结束，F’是适宜性超低频低通滤波的结果，是对F做了二级修正的含有较高频成分的低通图像。

2.根据权利要求1所述的一种医学X光影像滤波处理的二级修正方法，其特征在于：所述的步骤1具体是：采用小波低通冲击反应函数反复应用于源影像来获得超低频分量图像。

3.根据权利要求1所述的一种医学X光影像滤波处理的二级修正方法，其特征在于：所述的步骤6具体是：

对源影像P(i，j)用超低频分量图像F(i，j)做部分减影，获取的是没有二级修正的处理后影像P₀(i，j)；

P₀(i，j)＝P(i，j)-r F(i，j)

其中，r是图像处理强度，r∈[0，1)；

对影像P₀(i，j)做二级修正，得到最终处理后影像P_f(i，j)；

P_f(i，j)＝P₀(i，j)+r c{L₂(i，j)-c₂H₂(i，j)}

其中，c为二级修正强度，c∈[0，1]；c₂是二级修正区较高频分量修正强度，c₂∈[0，1]。