发明名称 舌头彩色数码照片的舌苔的分割提取方法
摘要 本发明属于医用诊断用彩色图像处理技术领域,其特征在于依次含有的下步骤:预处理,去除反光点;对去除反光点后的图像进行分裂合并计算对颜色较相近的部位用同一颜色代替;用Ostu阈值法对分类合并后的图像进行分割,区分舌质和舌苔部分;在色调空间,对舌苔部分用Ostu阈值法进行分割,分割出厚苔,其余部分为薄苔;用不同的颜色分别标示舌质、厚苔和薄苔后列出。本发明对舌苔进行二次分割,具有分割细致,结果兼顾医生主观判断的优点。
申请公布号 CN100405402C 申请公布日期 2008.07.23
申请号 CN200610112441.6 申请日期 2006.08.18
申请人 清华大学 发明人 白净;张永红;吴佳;史云迪
分类号 G06T5/00(2006.01) 主分类号 G06T5/00(2006.01)
代理机构 代理人
主权项 1.舌头彩色数码照片的舌苔的分割提取方法,其特征在于,所述方法在一台计算机中依次按以下步骤进行:步骤(1),把舌头彩色数码照片输入计算机,对舌头以外的部分用白色填充,保留舌头部分;步骤(2),预处理,按以下步骤去除反光点,并对去除的反光点计数储存:把步骤(1)得到的舌头彩色图映射为灰度图,并按下列准则判断舌头上的区域内反光点:f(x,y)>threshold_Flash;其中,(x,y)为坐标,f(x,y)是该灰度图在坐标(x,y)点像素的灰度值,取值范围是从0到255;threshold_Flash取得范围为160到180,f(x,y)>threshold_Flash者为反光点;步骤(3),把去除反光点后的舌头彩色图像分成小块A,每一块为8×8或者16×16个像素;步骤(4),先把步骤(3)得到的小块图像A分开成四个子区域,用A1,A2,A3,A4标志,子区域的大小为4×4或者8×8个像素,然后,再按照以下分裂准则和合并准则,判断所述四个子区域应该按现状分裂,还是应合并为一小块图像A:若:小块图像A满足下列任一个分裂准则,则把原来用同一颜色表达的A块分成四个颜色不同的小区域A1,A2,A3,A4: |R(Ak)-R(A)|>thresholdR “或” |G(Ak)-G(A)|>thresholdG “或” |B(Ak)-B(A)|>thresholdB 其中,R(A),G(A),B(A)代表区域A块中红色、绿色、蓝色的平均颜色,取值范围是从0到255,R(Ak),G(Ak),B(Ak)分别代表子区域Ak块中的红色、绿色、蓝色平均颜色,k=1,2,3,4,threshodR,threshodG,threshodB,依次分别表示为红色、绿色、蓝色的阈值,取值范围均为10到25;若:区域A块不需要分裂则按以下合并准则判断该A块相邻的四个包括A块在内的都不需要分裂的小区域X1,X2,X3,X4是否应该合并成一个大区域X,大小为16×16或者32×32: |R(Xk)-R(X)|<thresholdR “且” |G(Xk)-G(X)|<thresholdG “且” |B(Xk)-B(X)|<thresholdB 其中,R(Xk),G(Xk),B(Xk)依次分别代表Xk块中的红色、绿色、蓝色平均颜色,k=1,2,3,4;R(X),G(X),B(X)依次代表X块中红色、绿色、蓝色平均颜色,取值范围是从0到255;thresholdR,thresholdG,thresholdB同上述,若每一块Xk的平均颜色和四块综合后的块X的平均颜色之差低于设定阈值时,则把这四块区域合并,做同一块区域处理;步骤(5),把分裂合并后的舌头彩色图像,映射到灰度空间内,用Ostu阈值法对该图像进行分割,区分舌质和舌苔:若:f(x,y)>threshold_Coating,则(x,y)为舌质,若:f(x,y)≤threshold_Coating,则(x,y)为舌苔,其中,阈值threshold_Coating用Ostu法根据图像自动求得,且可在±10的范围内调整;步骤(6),把步骤(5)得到的舌头彩色图像转入到色调空间,对舌苔进行第二次分割,从舌苔中分割出厚苔,其分割公式如下:从RGB空间转换到色调空间的表达式为:<math><mrow><mi>H</mi><mo>=</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><msup><mi>cos</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>[</mo><mfrac><mrow><mrow><mo>(</mo><mi>R</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mi>G</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mi>R</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mi>B</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow><mrow><mn>2</mn><msqrt><msup><mrow><mo>(</mo><mi>R</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mi>G</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mrow><mo>(</mo><mi>R</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mi>B</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mi>G</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mi>B</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>)</mo></mrow></msqrt></mrow></mfrac><mo>]</mo><mo>,</mo></mrow></math> H(x,y)为点(x,y)的色调值,R(x,y),G(x,y),B(x,y)分别为该点的红色,绿色,蓝色分量;若:H(x,y)>thresholdH,则(x,y)为厚苔,若:H(x,y)≤thresholdH,则(x,y)为薄苔,阈值thresholdH用Ostu阈值法根据色调图像自动求得,且可在±10的范围内调整;步骤(7),计算机用不同颜色分别表示舌质,薄苔,厚苔部分,列出结果。
地址 100084北京市100084-82信箱
您可能感兴趣的专利