WPF：图像处理（五）迭代法

作者：Splash

转自：http://blog.csdn.net/jhqin/article/details/7472219

[csharp] view plain copy

1. using System;  
2.   
3. namespace Splash.Imaging  
4. {  
5.     /// <summary>  
6.     /// 图像处理：迭代法二值化阈值计算方法  
7.     /// </summary>  
8.     public static partial class Binarize  
9.     {  
10.         /// <summary>  
11.         /// 迭代法计算阈值  
12.         /// </summary>  
13.         /// <param name="grayArray">灰度数组</param>  
14.         /// <returns>二值化阈值</returns>   
15.         public static Int32 IterativeThreshold(Byte[,] grayArray)  
16.         {   // 建立统计直方图  
17.             Int32[] Histogram = new Int32[256];  
18.             Array.Clear(Histogram, 0, 256);     // 初始化  
19.             foreach (Byte b in grayArray)  
20.             {  
21.                 Histogram[b]++;                 // 统计直方图  
22.             }  
23.   
24.             // 总的质量矩和图像点数  
25.             Int32 SumC = grayArray.Length;    // 总的图像点数  
26.             Int32 SumU = 0;  
27.             for (Int32 i = 1; i < 256; i++)  
28.             {  
29.                 SumU += i * Histogram[i];     // 总的质量矩                  
30.             }  
31.   
32.             // 确定初始阈值  
33.             Int32 MinGrayLevel = Array.FindIndex(Histogram, NonZero);       // 最小灰度值  
34.             Int32 MaxGrayLevel = Array.FindLastIndex(Histogram, NonZero);   // 最大灰度值  
35.             Int32 T0 = (MinGrayLevel + MaxGrayLevel) >> 1;  
36.             if (MinGrayLevel != MaxGrayLevel)  
37.             {  
38.                 for (Int32 Iteration = 0; Iteration < 100; Iteration++)  
39.                 {   // 计算目标的质量矩和点数  
40.                     Int32 U0 = 0;  
41.                     Int32 C0 = 0;  
42.                     for (Int32 i = MinGrayLevel; i <= T0; i++)  
43.                     {   // 目标的质量矩和点数                  
44.                         U0 += i * Histogram[i];  
45.                         C0 += Histogram[i];  
46.                     }  
47.   
48.                     // 目标的平均灰度值和背景的平均灰度值的中心值  
49.                     Int32 T1 = (U0 / C0 + (SumU - U0) / (SumC - C0)) >> 1;  
50.                     if (T0 == T1) break; else T0 = T1;  
51.                 }  
52.             }  
53.   
54.             // 返回最佳阈值  
55.             return T0;  
56.         }          
57.     }  
58. }