VC++数字图像处理源程序及详解之6（灰度均衡）.docVIP

9.6? 灰度均衡 本节介绍灰度均衡的原理和实现方法，它是增强图像的有效手段之一。 9.6.1? 基本原理 灰度均衡是以累计分布函数变换为基础的直方图修正法，它可以产生一幅灰度级分布概率均匀的图像。也就是说，经 过灰度均衡后的图像在每一级灰度上都具有相同数量的像素点，对应灰度直方图的每一级灰度具有相同的高度。灰度均衡同样也属于改进图像的方法，灰度均衡后的 图像具有最大的信息量。 图9-33显示了高亮度和低对比度图像经过灰度均 衡后的效果，可以看到灰度均衡对图像效果进行了重要的改进。从变换后图像的直方图来看，灰度分布更加均匀。 下面进行灰度均衡变换函数 的推导。 设转换前图像的密度函数为pr(r)， 其中0≤r≤1；转化后图像的密度函数为ps(s)，同样有0≤s≤1； 灰度均衡变换函数为s=T(r)。从概率理论可以得到如下公式： 图9-33? 利用灰度均衡处理图像 转化后图像灰度均匀分布，有ps(s)=1， 故： 两边取积分有： 这就是图像的累计分布函数。对于图像而言，密度函 数为： 其中x表 示灰度值，nx表示灰度级为x的像素个数，n表示图像总像素个数。 通过上面的公式就能推导出基于离散型的灰度均衡公 式： 其中y的取值范围是0至255。 9.6.2? 编程实现 1．灰度均衡函数 下面的GrayEqualize函数实现了灰度 均衡的功能： /*************************************************************************** *?? 作用: 灰度均衡函数 *?? 参数:??????? *?????? pixel: ???? 原始像素数组 *?????? tempPixel: 保存变换后图像的像素数组 *?????? width: ???? 原始图像宽度 ***************************************************************************/ void GrayEqualize(BYTE* pixel, BYTE* tempPixel, UINT width, UINT height) { ??? // 灰度映射表 ??? BYTE map[256]; ??? long lCounts[256]; ??? memset(lCounts, 0, sizeof(long) * 256); ??? // 计算各灰度值个数 ??? for (UINT i = 0; i width * height; i++) ??? { ??????? int x = pixel[i * 4]; ??????? lCounts[x]++; ??? } ??? // 保存运算中的临时值 ??? long lTemp; ??? for (int i = 0; i 256; i++) ??? { ??????? lTemp = 0; ??????? for (int j = 0; j = i; j++) ??????????? lTemp += lCounts[j]; ??????? map[i] = (BYTE)(lTemp * 255.0f / width / height); ??? } ??? // 变换后的值直接在映射表中查找 ??? for (UINT i = 0; i width * height; i++) ??? { ??????? int x = pixel[i * 4]; ??????? tempPixel[i*4] = tempPixel[i*4+1] = tempPixel[i*4+2] = map[x]; ??????? tempPixel[i*4+3] = 255; ??? } } 2．使用灰度均衡函数 打开“菜单栏”资源，在“点运算”菜单下添加一个 新项“灰度均衡”，其ID为ID_POINT_EQUALIZE，如图9-34所示。 图9-34? 在菜单栏资源上添加新项 最后添加响应函数，内容如下： void CMagicHouseView::OnPointEqualize() { ??? if (!m_bIsEditMode || m_nPos == -1) ??? { ??????? MessageBox(L请先打开图像文件，然后选择编辑模式！); ??????? return; ??? } ??? GrayEqualize(m_pImageBuffer, m_pImageTempBuffer, m_nPicWidth, m_nPicHeight); ??? // 创建预览对话框 ??? CPreviewDlg dlg; ??? dlg.m_nType = 2;??? // 对话框在关闭时提