本发明涉及数字图像处理技术领域，更具体的说是涉及一种数字病理全切片图像空白区域自动分割方法。
背景技术
数字病理全切片图像，是将传统的玻璃病理切片通过全自动显微镜或光学放大系统扫描采集得到的大尺寸高分辨数字图像，是病理医生在诊断时的重要依据。近年来，随着病理学和计算机技术的快速发展，数字病理全切片图像的数量迅速增长，通过图像算法或机器学习方法对其进行处理和分析，对辅助病理医生快速诊断以及计算机自动诊断有很大参考价值。
但是，数字病理全切片图像在扫描时除细胞组织等内容外，图像外围通常为没有实际内容的空白区域，在处理数字病理全切片图像时如果不加区分一同处理空白区域，不仅会增加处理时间，也会影响算法精度。此外，由于采集环境不同，导致不同的数字病理全切片图像切片的亮度和色调有差异，因此使用固定灰度阈值对所有数字病理全切片进行空白区域分割可能造成结果不准确。
因此，如何提供一种结果准确的数字病理全切片图像空白区域的分割方法，是本领域技术人员亟需解决的问题。
技术实现要素：
有鉴于此，本发明提供的一种数字病理全切片图像空白区域自动分割方法，通过颜色空间变换和自适应阈值方法对数字病理全切片进行二值化，然后进行形态学后处理，最终使计算机能够准确地对数字病理全切片空白区域进行自动分割，结果更加准确。
为了实现上述目的，本发明采用如下技术方案：
一种数字病理全切片图像空白区域自动分割方法，包括以下操作步骤：
s1、将数字病理全切片图像由rgb颜色空间转换到hsv颜色空间；
s2、求出所述hsv颜色空间中饱和度通道和亮度通道的阈值；
s3、使用所述饱和度通道阈值和所述亮度通道阈值对所述数字病理全切片图像进行二值化操作，得到初步分割结果；
s4、对所述初步分割结果进行形态学后处理，得到最终分割结果。
优选地，步骤s2中，所述饱和度通道阈值计算和所述亮度通道阈值计算包括全局均值法或最大类间方差法。图像的二值化，就是将图像上的像素点的灰度值设置为0或255，也就是将整个图像呈现出明显的只有黑和白的视觉效果。其中，全局均值法是对整幅图像选取单一的阈值来进行二值化；最大类间方差法是指是一种自适应的阈值确定的方法，又叫大津法，它是按图像的灰度特性，将图像分成背景和目标两部分。背景和目标之间的类间方差越大，说明构成图像的两部分的差别越大，当部分目标错分为背景或部分背景错分为目标都会导致两部分差别变小。因此，使类间方差法最大的分割意味着错分概率最小。
优选地，步骤s4中，所述形态学后处理包括开运算、闭运算或去除小面积噪点。开运算，就是先腐蚀后膨胀的过程，可以用来消除小物体，在纤细点处分离物体，并且在平滑较大物体的边界的同时不明显改变其面积。闭运算，就是先膨胀后腐蚀的过程，闭运算可以用来排除小型黑洞，即黑色区域。小面积噪点就是指面积小于某个阈值的连通区域，去除指让这种小面积噪点在二值图中的值变得与周围相同。
经由上述的技术方案可知，与现有技术相比，本发明公开提供了一种数字病理全切片图像空白区域自动分割方法，无需手动设置阈值，即可实现任何全切片区域自动分割，根据分割进行后续处理时，可以将空白区域剔除，只处理包含病理组织部分，不仅能够提高方法的准确度，还能大幅度减少计算量；同时，本发明得到的空白区域的分割结果对后续基于机器学习或图像处理的数字病理全切片分析方法具有很大的提升效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
图1附图为本发明提供的自动分割方法的流程图。
图2附图为本发明提供的自动分割方法的流程中的各个状态图，其中，
图(a)为数字病理全切片的原始状态图；
图(b)为数字病理全切片的初步分割结果的状态图；
图(c)为数字病理全切片的最终分割结果的状态图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
参照图1-2，一种数字病理全切片图像空白区域自动分割方法，应用于任意全切片，具体包括以下步骤：
s1、将数字病理全切片图像由rgb颜色空间转换到hsv颜色空间；
rgb，即红绿蓝颜色空间是依据人眼识别的颜色定义出的空间，可表示大部分颜色，但是它将色调、饱和度和亮度三个量放在一起表示，难以进行图像处理，因此需要转换到hsv，即色调、饱和度和亮度颜色空间。
s2、求出数字病理全切片图像中s饱和度通道阈值和v亮度通道的阈值；
设数字病理全切片像素宽度和高度分别为w和h，位置坐标为(i，j)的像素的s通道和v通道数值分别为s(i，j)和v(i，j)，数字病理全切片在s通道和v通道的阈值分别为st和vt。
st和vt的计算采用全局均值计算方法时，阈值计算公式为：
st和vt的计算采用最大类间方差法时，在饱和度通道中小于st的像素个数及其占整幅图像的比例分别记为n0和ω0，其平均灰度为μ0；不小于st的像素个数及其占整幅图像的比例分别记为n1和ω1，其平均灰度为μ1。图像的总平均灰度记为μ,两种像素的类间方差记为g。则有：
g＝ω0(μ0-μ)2+ω1(μ1-μ)2＝ω0ω1(μ0-μ1)2
阈值st计算公式为：
vt的计算与st相似。
s3、使用步骤s2中得到的饱和度通道阈值和亮度通道阈值对数字病理全切片图像进行二值化操作，得到初步分割结果；
数字病理全切片在hsv颜色空间中，如果一个像素点的饱和度通道值小于其阈值且亮度通道值大于其阈值，则标记为空白，否则标记为有内容。如图2(b)所示，设数字病理全切片的初步分割结果为二值图像b，b中位置坐标为(i，j)的值为b(i，j)，其余符号与步骤s2中相同，则计算公式为
s4、对步骤s3中得到的初步分割结果进行形态学后处理，得到最终分割结果。
如图2(c)所示，为了抑制初步分割结果中的噪声点与不平滑轮廓，可以使用形态学操作进行后处理，方法中包括开运算、闭运算或去除小面积噪点。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。