程序: // 锁定DIB并返回指向DIB的指针LPSTR lpDIB = (LPSTR) ::GlobalLock((HGLOBAL) hDIB);// 找到DIB图像象素起始位置并返回指向DIB象素指针LPSTR lpDIBBits = m_clsDIB.FindDIBBits(lpDIB);// 获取DIB的宽度LONG lWidth = m_clsDIB.DIBWidth(lpDIB);// 获取DIB的高度LONG lHeight = m_clsDIB.DIBHeight(lpDIB);for (i = 0; i < lHeight; i ++) // 对各像素进行灰度转换{ for (j = 0; j < lWidth * 3; j ++) { // 对各像素进行灰度统计 unsigned char R = *((unsigned char *)lpDIBBits + lWidth * 3 * i + j); nNs_R[R]++; j++; unsigned char G = *((unsigned char *)lpDIBBits + lWidth * 3 * i + j); nNs_G[G]++; j++; unsigned char B = *((unsigned char *)lpDIBBits + lWidth * 3 * i + j); nNs_B[B]++; }}for(i=0;i<256;i++) // 计算灰度分布密度{ fPs_R[i] = nNs_R[i] / (lHeight * lWidth * 1.0f); fPs_G[i] = nNs_G[i] / (lHeight * lWidth * 1.0f); fPs_B[i] = nNs_B[i] / (lHeight * lWidth * 1.0f);}for(i = 0; i < 256; i++){ //计算累计直方图分布 if(i == 0) { temp_r[0] = fPs_R[0]; temp_g[0] = fPs_G[0]; temp_b[0] = fPs_B[0]; } else { temp_r[i] = temp_r[i-1] + fPs_R[i]; temp_g[i] = temp_g[i-1] + fPs_G[i]; temp_b[i] = temp_b[i-1] + fPs_B[i]; } //累计分布取整,nNs_R[]、nNs_G[]、nNs_B[]保存有计算出来的灰度映射关系 nNs_R[i] = (int)(255.0f * temp_r[i] + 0.5f); nNs_G[i] = (int)(255.0f * temp_g[i] + 0.5f); nNs_B[i] = (int)(255.0f * temp_b[i] + 0.5f);}for (i = 0; i < lHeight; i ++){ for (j = 0; j < lWidth * 3; j ++) { //对R分量进行灰度映射(均衡化) unsigned char R = *((unsigned char *)lpDIBBits + lWidth * 3 * i + j); *((unsigned char *)lpDIBBits + lWidth * 3 * i + j) = nNs_R[R]; j++; //对G分量进行灰度映射(均衡化) unsigned char G = *((unsigned char *)lpDIBBits + lWidth * 3 * i + j); *((unsigned char *)lpDIBBits + lWidth * 3 * i + j) = nNs_G[G]; j++; //对B分量进行灰度映射(均衡化) unsigned char B = *((unsigned char *)lpDIBBits + lWidth * 3 * i + j); *((unsigned char *)lpDIBBits + lWidth * 3 * i + j) = nNs_B[B]; }}
Opencv13(直方图和均衡化)