[发明专利]一种切伦科夫荧光图像冲击噪声去除方法及系统

权利要求书

1.一种切伦科夫荧光图像冲击噪声去除方法，其特征在于，所述切伦科夫荧光图像冲击噪声去除方法包括：复制光学探测器采集到的原始切伦科夫荧光图像边缘，对图像进行扩展；以扩展图像的元素为中心点，取3×3窗口内的9个值，从小到大进行排序，得到中间值，中间值构成的图像为中值图像；计算原始切伦科夫荧光图像与中值图像的差值图像；计算差值图像中各个像素的模糊值；根据模糊值，计算得到滤除冲击噪声后的切伦科夫荧光图像。

2.如权利要求1所述的切伦科夫荧光图像冲击噪声去除方法，其特征在于，所述切伦科夫荧光图像冲击噪声去除方法包括以下步骤：

步骤一，复制光学探测器采集到的具有K行N列像素的原始切伦科夫荧光图像X的上下边缘也即第1行和最后1行，扩展得到大小为K+2行N列的行扩展图像

步骤二，复制行扩展图像的左右边缘也即第1列和最后1列，得到大小为K+2行N+2列的扩展图像Y；

步骤三，以扩展图像Y的第2行第2列的元素y_2,2为中心点，取3×3窗口内的9个元素，记为{y_1,1,y_1,2,y_1,3,y_2,1,y_2,2,y_2,3,y_3,1,y_3,2,y_3,3}，对9个元素从小到大进行排序，得到中间值，记为m_1,1；

步骤四，从左至右，从上至下滑动(3)中的3×3窗口，计算扩展图像Y中除边缘也即第1行、第K+2行、第1列和第N+2列以外的所有像素的中值，得到中值图像M；

步骤五，计算原始切伦科夫荧光图像与中值图像的差值图像V＝X-M；

步骤六，对步骤五得到的差值图像，计算各个像素的模糊值；

步骤七，根据步骤六得到的模糊值，计算得到滤除冲击噪声后的切伦科夫荧光图像。

3.如权利要求2所述的切伦科夫荧光图像冲击噪声去除方法，其特征在于，所述步骤六中计算各个像素的模糊值为：

式中，v_i,j为差值图像V的第i行第j列像素，且1≤i≤K，1≤j≤N，T₁和T₂为两个阈值。

4.如权利要求2所述的切伦科夫荧光图像冲击噪声去除方法，其特征在于，所述步骤七中计算得到滤除噪声后的切伦科夫荧光图像为：

z_i,j＝r_i,j×m_i,j+(1-r_i,j)×x_i,j；

式中，m_i,j和x_i,j分别为中值图像M和原始切伦科夫荧光图像X的第i行第j列像素。

5.一种实现权利要求1所述切伦科夫荧光图像冲击噪声去除方法的切伦科夫荧光图像冲击噪声去除系统，其特征在于，所述切伦科夫荧光图像冲击噪声去除系统包括：

采集模块，用于在光密闭环境中采集原始切伦科夫荧光图像；

扩展模块，用于原始切伦科夫荧光图像边缘扩展；

排序模块，用于以图像的元素为中心点，取3×3窗口内的9个元素，从小到大进行排序，得到中间值；

中值图像计算模块，用于从左至右，从上至下滑动3×3窗口，计算图像中像素的中值，得到中值图像；

差值图像计算模块，用于计算原始切伦科夫荧光图像与中值图像的差值图像；

模糊值计算模块，用于根据得到的差值图像计算各个像素的模糊值；

图像恢复模块，用于根据模糊值，计算得到滤除冲击噪声后的切伦科夫荧光图像。

6.一种应用权利要求1～4任意一项所述切伦科夫荧光图像冲击噪声去除方法的信息处理终端。