VC++编程实现对曝光不足照片的修复
2010-06-15 20:39:54 来源:WEB开发网 闂傚倸鍊搁崐鎼佸磹閹间礁纾归柟闂寸绾惧綊鏌熼梻瀵割槮缁炬儳缍婇弻鐔兼⒒鐎靛壊妲紒鐐劤缂嶅﹪寮婚悢鍏尖拻閻庨潧澹婂Σ顔剧磼閹冣挃闁硅櫕鎹囬垾鏃堝礃椤忎礁浜鹃柨婵嗙凹缁ㄧ粯銇勯幒瀣仾闁靛洤瀚伴獮鍥敍濮f寧鎹囬弻鐔哥瑹閸喖顬堝銈庡亝缁挸鐣烽崡鐐嶆棃鍩€椤掑嫮宓佸┑鐘插绾句粙鏌涚仦鎹愬闁逞屽墰閹虫捇锝炲┑瀣╅柍杞拌兌閻ゅ懐绱撴担鍓插剱妞ゆ垶鐟╁畷銉р偓锝庡枟閻撴洘銇勯幇闈涗簼缂佽埖姘ㄧ槐鎾诲礃閳哄倻顦板┑顔硷工椤嘲鐣烽幒鎴旀瀻闁规惌鍘借ⅵ濠电姷鏁告慨顓㈠磻閹剧粯鈷戞い鎺嗗亾缂佸鏁婚獮鍡涙倷閸濆嫮顔愬┑鐑囩秵閸撴瑦淇婇懖鈺冪<闁归偊鍙庡▓婊堟煛鐏炵硶鍋撻幇浣告倯闁硅偐琛ラ埀顒冨皺閺佹牕鈹戦悙鏉戠仸闁圭ǹ鎽滅划鏃堟偨缁嬭锕傛煕閺囥劌鐏犻柛鎰ㄥ亾婵$偑鍊栭崝锕€顭块埀顒佺箾瀹€濠侀偗婵﹨娅g槐鎺懳熺拠鑼舵暱闂備胶枪濞寸兘寮拠宸殨濠电姵纰嶉弲鎻掝熆鐠虹尨宸ョ€规挸妫濆铏圭磼濡搫顫嶇紓浣风劍閹稿啿鐣烽幋锕€绠婚悹鍥у级瀹撳秴顪冮妶鍡樺鞍缂佸鍨剁粋宥夋倷椤掍礁寮垮┑鈽嗗灣閸樠勭妤e啯鍊垫慨妯煎亾鐎氾拷

照片修复程序的实现
根据上述算法就可以对照片进行均衡化处理了。对于照片一般是将其当作DIB(设备无关位图)来处理的,在将其装载到内存并获取到其HDIB 型句柄 hDIB后,需要预先获取该位图的一些相关参数以便后续使用:
// 锁定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);
根据算法设计,首先需要对原始位图进行灰度统计并对灰度分布密度进行计算,对于24位真彩位图,一个象素是由三个字节(R、G、B三颜色分量)组成的,因此需要分别对各颜色分量单独进行统计,同样,在灰度均衡处理时也是各自独立进行的。
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);
}
接下来就要根据前面计算得到的各颜色分量的灰度分布密度来计算R、G、B三分量各自灰度级的累计直方图分布,并对其进行取整以得出从源到目标图像的灰度映射关系表。这一部分实际也是本文的核心,根据这一步计算得到的映射关系表就可以确定出源图象的某个灰度级在目的图象中对应于哪个新的灰度:
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);
}
更多精彩
赞助商链接