主权项 |
一种基于二维模糊熵的粘连物分割方法,其特征在于,该方法包括以下步骤:S1:采集粘连物图像并提取目标粘连物图像;S2:在步骤S1的基础上获得所述粘连物图像的初始分割区域;S3:计算初始分割区域中的过分割区域和正确分割区域之间的最佳分割阈值,根据所述最佳分割阈值得到过分割区域和正确分割区域;S4:计算过分割区域合并后的区域特征;S5:将所述区域特征与标准样本特征模型进行比较,选择与标准样本特征模型距离最小的所述区域特征对应的区域作为合并区域,进而得到粘连物的分割结果;其中,所述步骤S3具体包括:S31:确定粘连物图像的隶属度函数;S32:根据所述隶属度函数对粘连物图像进行二维划分;S33:通过所述二维划分得到二维模糊熵;S34:对所述二维模糊熵做最大化处理,得到初始分割区域的过分割区域和正确分割区域的最佳分割阈值,得到过分割区域和正确分割区域;所述最佳分割阈值的计算公式为:(a,b,c,d)=Argmax<sub>M×N</sub>(H)其中:a为过分割区域像素最小分割阈值;b为过分割区域像素最大分割阈值;c为过分割区域相对高度最小分割阈值;d为过分割区域相对高度最大分割阈值;H为二维模糊熵;M为集水盆中像素个数;N为分水岭和集水盆最低处的高度差的等级个数;二维模糊熵步骤具体为:H=H(R<sub>B</sub>)+H(R<sub>O</sub>)H为二维模糊熵;H(R<sub>B</sub>)为正确分割的熵;H(R<sub>O</sub>)为过分割的熵;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><msub><mi>R</mi><mi>B</mi></msub><mo>)</mo></mrow><mo>=</mo><mo>-</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>μ</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mfrac><mrow><mi>c</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mrow><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>c</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></mfrac><mi>log</mi><mrow><mo>(</mo><mfrac><mrow><mi>c</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mrow><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>c</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000627693420000021.GIF" wi="1319" he="341" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>H</mi><mrow><mo>(</mo><msub><mi>R</mi><mi>O</mi></msub><mo>)</mo></mrow><mo>=</mo><mo>-</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>μ</mi><mi>O</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mfrac><mrow><mi>c</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mrow><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>c</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></mfrac><mi>log</mi><mrow><mo>(</mo><mfrac><mrow><mi>c</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow><mrow><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>c</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000627693420000022.GIF" wi="1324" he="345" /></maths>μ<sub>B</sub>(i,j)为正确分割的隶属度函数;c(i,j)表示集水盆的像素个数为i,分水岭和集水盆最低处的相对高度等级为j的参数对的发生次数;μ<sub>O</sub>(i,j)为过分割的隶属度函数。 |