1 #include <string>
2 #include "20140318计算类的面积.cpp"
3
4 //////////////////////////////////////////////////////////////////////////
5 // 函数名称: histeq()
6 // 传入参数:
7 // BYTE*dstData 传入灰度图像内存空间的指针
8 // double *srcArray 模版的直方图累积，并进行归一化，大小为256
9 // int m_Width 匹配内存空间的宽度
10 // int m_Height 匹配内存空间的高度
11 // int m_pitch 匹配内存空间的每行所在内存大小
12 ////////////////////////////////////////////////////////////////////////
13
14 void histeq(int **dstData, double *srcArray,int m_Width,int m_Height,int m_pitch)
15 {
16 //void *memset(void *s,int c,size_t n)将已开辟内存空间 s 的首 n 个字节的值设为值 c
17 m_pitch=sizeof(data[0][0])*nXSize;
18 double dstHist[256];
19 memset(dstHist,0,256 * sizeof(double));
20 double dstArray[256];
21 memset(dstArray,0,256 * sizeof(double));
22 int i=0,j=0;
23 //统计直方图
24 for (i = 0;i < m_Height;i++)
25 {
26 for (j = 0;j < m_Width;j++)
27 {
28 dstHist[(int)dstData[i * m_pitch + j]]++;
29 }
30 }
31 //计算直方图累积
32 double m_Bytes = m_Width * m_Height;//m_Bytes初始化为像素总个数
33 dstArray[0] = dstHist[0];//第一个初始化为dstHist[0]=0
34 for (i = 1;i < 256;i++)
35 {
36 dstArray[i] = dstArray[i - 1] + dstHist[i];
37 }
38 //直方图累积归一化
39 for (i = 0;i < 256;i++)
40 {
41 dstArray[i] /= m_Bytes;
42 }
43 //直方图匹配
44 double m_diffA,m_diffB;
45 int k = 0;
46 float mapPixel[256];//每个颜色的value
47 memset(mapPixel,0,256 * sizeof(float));
48
49
50 for (i = 0;i < 256;i++)
51 {
52 m_diffB = 1;
53 for (j = k; j < 256;j++)
54 {
55 m_diffA = abs(dstArray[i] - srcArray[j]);
56 if (m_diffA - m_diffB < 1.0E-5)
57 {
58 m_diffB = m_diffA;
59 k = j;
60 }
61 else
62 {
63 k = j - 1;
64 break;
65 }
66 }
67 if (k == 255)
68 {
69 for (int l = i;l < 256;l++)
70 {
71 mapPixel[l] = (float) k;
72 }
73 break;
74 }
75 mapPixel[i] = (float) k;
76 }
77 //目标图像查找索引表
78 for (i = 0;i < m_Height;i++)
79 {
80 for (j = 0;j < m_Width;j++)
81 {
82 dstData[i * m_pitch + j] = mapPixel[(int)dstData[i * m_pitch + j]];
83 }
84 }
85 }
