[发明专利]一种血管内超声图像序列中钙化斑块帧的自动检索方法

权利要求书

1.一种血管内超声图像序列中钙化斑块帧的自动检索方法，其特征是，所述方法根据IVUS图像径向灰度曲线的斜率值来判断斑块是否存在，具体步骤如下：

a、对各帧IVUS图像进行各向异性扩散滤波，减少噪声和无用信息，同时保留、增强图像的边缘信息；

b、对各向异性扩散滤波后的各帧IVUS图像进行极坐标变换，将其变换到以图像中心为坐标原点的极坐标系中，得到极坐标视图；

c、根据极坐标视图求得各个角度的径向灰度变化曲线；

d、初步检索含钙化斑块的图像：

将每一径向灰度变化曲线的斜率最大值、即极径上的灰度跳变极值，与预设的参考阈值进行对比，如果某方向上的灰度变化曲线的斜率最大值超过所设阈值，则初步认定该方向属于钙化区域，否则不属于；

e、精细检索含钙化斑块的图像：

如果初步认定极坐标视图中的某列像素属于钙化区域，那么判断其邻域内的连续n列像素，如果都属于初步认定的钙化区域，则该区域属于钙化区域，否则该区域不是钙化斑块，最后将最终检测到的钙化斑块的具体位置显示在图像中。

2.根据权利要求1所述的一种血管内超声图像序列中钙化斑块帧的自动检索方法，其特征是，对各帧IVUS图像进行各向异性扩散滤波的具体方法如下：

遍历整幅图像，对当前像素(x, y)，首先计算其在8-邻域方向上的灰度梯度即                                                的值（k=1,2,…,8），然后，将原始图像的灰度值I₀(x,y)作为各向异性扩散滤波器零时刻的输入f(x,y;0)，求得滤波后图像：

 ，

其中f(x,y;t)是各向异性扩散滤波器在t时刻的输出，是梯度算子，c_k(x,y;t)是t时刻当前像素在第k个8-邻域方向上的扩散系数，其构造方法是：

，

其中λ是边缘幅值的门限参数。

3.根据权利要求1或2所述的一种血管内超声图像序列中钙化斑块帧的自动检索方法，其特征是，对各向异性扩散滤波后的各帧IVUS图像进行极坐标变换的具体方法如下：

设原直角坐标图像为f(x, y)，那么利用直角坐标(x, y)与极坐标(ρ, θ)的转换关系：

即可得到其极坐标视图

其中ρ是极径，θ是极角。