[发明专利]一种螺旋断层放疗图像质量提高方法

摘要

一种螺旋断层放疗图像质量提高方法，它是涉及一种基于视网膜‑大脑皮层理论的螺旋断层放疗图像质量提高方法，它有四大步骤步骤一计算机在MATLAB环境下读取螺旋断层放疗图像；步骤二对图像进行双边滤波去噪；步骤三对图像进行基于视网膜‑大脑皮层理论的对比度提升；步骤四使用高斯‑赛格尔迭代法对图像边缘进行增强。本发明解决了原螺旋断层放疗图像噪声较多、对比度差和边缘不清晰的问题，取得了较好的质量增强结果，在螺旋断层治疗领域里具有广阔的应用前景。

主权项

一种螺旋断层放疗图像质量提高方法，它是涉及一种基于视网膜‑大脑皮层理论的螺旋断层放疗图像质量提高方法，其特征在于：该方法的具体步骤如下：步骤一：首先通过螺旋断层放疗仪器输出计算机数字图像，然后使用Matlab语言中的imread函数读取该图像，将其信息转为Matlab矩阵形式，使用Matlab语言对其进行处理；螺旋断层放疗图像为512像素*512像素*R通道的数字图像，即读入的矩阵数据为512*512*R维度；其中符号R表示该图像中包含的断层数量，每一断层为一幅灰度图像，R通道数目的断层图像堆叠成完整的螺旋断层放疗图像；以下使用大写字母I表示图像被Matlab读入后的矩阵数据其中的某一通道，即512*512维度的矩阵；上标数字表明不同步骤的中间图像矩阵数据，如I0表示原始图像矩阵数据；使用下标表示矩阵中某一元素，即对应图像中的某个坐标上的像素，即Ix,y表示矩阵I中位于坐标(x,y)的元素；步骤二：在Matlab软件中，对上一步骤读入的矩阵，即I0进行如下运算：Ix,y1=ΣΩIi,j0w(x,y,i,j)ΣΩw(x,y,i,j)]]>其中，Ω表示(x,y)坐标的邻域范围，坐标(i,j)满足(i,j)∈Ω；权重系数w的取值为值域权重与空间域权重的乘积；该步骤为图像的双边滤波操作，在保持图像边缘的同时有良好的去噪效果；步骤三：在Matlab软件中，对上一步骤的矩阵进行如下操作：Ix,y2=(11+exp(-a*Rx,y+b)+c)c+1]]>其中，a,b,c为超参数，由手动设置，并有：Rx,y=Ix,y1Ix,y1*G]]>这里，符号*代表卷积，符号G表示高斯核，即分母为分子高斯核模糊后的图像；该步骤即为改进的Retinex算法，它没有采用以往常见的Retinex方法中采用的对数形式，而是使用了新的曲线函数形式；步骤四：首先，使用Matlab中的索贝尔算子函数，对矩阵I2进行操作，获得的边缘图像矩阵标记为Gra；然后，对边缘图像矩阵Gra进行逐像素遍历，寻找邻域内的边缘最大值，得到矩阵GraM来标记各个元素是否为矩阵Gra中的邻域最大值，再通过：Grax,y1=Grax,y·fx,y]]>更新矩阵Gra，其中fx,y为一以当前点(x,y)距离最近边缘最大值点的距离为自变量的自定义函数，若当前点(x,y)与最近边缘最大值点距离较小时，该函数值大于1；若当前点(x,y)与最近边缘最大值点距离较大时，该函数值小于1；这样使新的Gra1相对应的边缘图像中的边缘的分布变窄、峰值变高；最后，使用高斯‑赛德尔迭代方法，利用矩阵Gra1反推出边缘增强的图像，步骤如下：步骤4.1：对矩阵Gra1再次进行差分运算，得到目标图像的二阶导图像矩阵Lap；步骤4.2：初始化矩阵dst0＝I2步骤4.3：反复进行如下迭代，直到迭代次数达到5次，其中n为迭代次数，n＝0,1,2…：1)2)经过以上步骤迭代，获得最终图像，完成对螺旋断层放疗图像的质量提高。