① Sobel算子檢測(cè)邊緣
void CDemoView::Onsobel()
{
CDemoDoc *pDoc = GetDocument();
HDIB dib = pDoc->GetHDIB();
LPSTR lpDIB = (LPSTR)::GlobalLock((HGLOBAL)dib);
LPBITMAPINFOHEADER phead = (LPBITMAPINFOHEADER)lpDIB; //信息頭指針
int biByteCount = phead->biBitCount / 8; //每像素字節(jié)數(shù)
if (biByteCount == 3)
{
::GlobalUnlock(dib);
return;
} //彩色圖不做處理
int w = phead->biWidth;
int h = phead->biHeight;
int lineByte = (biByteCount * w + 3) / 4 * 4;
unsigned char *lpDIBBits = (unsigned char *)::FindDIBBits(lpDIB); //像素位置指針,unsigned char取值范圍是0~255(一個(gè)字節(jié))
unsigned char *buf = new unsigned char[w*h]; //新內(nèi)存存放運(yùn)算結(jié)果
unsigned char c[2]; //存2個(gè)模板的卷積結(jié)果
unsigned char a[8]; //取鄰域元素
int r; //存卷積結(jié)果
for (int i = 1; i < h - 1; i++) {
for (int j = 1; j < w - 1; j++) {
a[0] = *(lpDIBBits + (i + 1) * lineByte + j - 1);
a[1] = *(lpDIBBits + (i + 1) * lineByte + j);
a[2] = *(lpDIBBits + (i + 1) * lineByte + j + 1);
a[3] = *(lpDIBBits + i * lineByte + j + 1);
a[4] = *(lpDIBBits + (i - 1) * lineByte + j + 1);
a[5] = *(lpDIBBits + (i - 1) * lineByte + j);
a[6] = *(lpDIBBits + (i - 1) * lineByte + j - 1);
a[7] = *(lpDIBBits + i * lineByte + j - 1);
c[0] = abs(-a[0] + a[2] + a[3] * 2 + a[4] - a[6] - a[7] * 2);
c[1] = abs(a[0] + a[1] * 2 + a[2] - a[4] - a[5] * 2 - a[6]);
r = 2 * pow(c[0] * c[0] + c[1] * c[1], 0.5); //綜合垂直與水平兩個(gè)方向壹罚，乘2再放大一下翩剪，結(jié)果變亮一些
if (r > 255)
r = 255;
*(buf + i * w + j) = r;
}
}
for (int i = 1; i < h - 1; i++) {
for (int j = 1; j < w - 1; j++) {
*(lpDIBBits + i * lineByte + j) = *(buf + i * w + j);
}
}
delete[]buf;
::GlobalUnlock(dib);
Invalidate(); //區(qū)域重繪
}
② 直方圖均衡
void CDemoView::OnHistEqual()
{
CDemoDoc *pDoc = GetDocument();
HDIB dib = pDoc->GetHDIB();
LPSTR lpDIB = (LPSTR)::GlobalLock((HGLOBAL)dib);
LPBITMAPINFOHEADER phead = (LPBITMAPINFOHEADER)lpDIB; //信息頭指針
int biByteCount = phead->biBitCount / 8; //每像素字節(jié)數(shù)
if (biByteCount == 3)
{
::GlobalUnlock(dib);
return;
} //彩色圖不做處理
int w = phead->biWidth;
int h = phead->biHeight;
int lineByte = (biByteCount * w + 3) / 4 * 4;
unsigned char *lpDIBBits = (unsigned char *)::FindDIBBits(lpDIB); //像素位置指針,unsigned char取值范圍是0~255(一個(gè)字節(jié))
double p[256];
for (int i = 0; i < 256; i++) {
p[i] = 0;
}
unsigned char value[256];
//統(tǒng)計(jì)直方圖
unsigned char idx;
for (int i = 0; i < h; i++) {
for (int j = 0; j < w; j++) {
idx = *(lpDIBBits + i * lineByte + j);
p[idx] += 1.0 / (w*h);
}
}
//累加
for (int i = 1; i < 256; i++) {
p[i] += p[i - 1];
}
//計(jì)算新的灰度值
for (int i = 0; i < 256; i++) {
value[i] = (unsigned char)(p[i] * 255);
}
for (int i = 0; i < h; i++) {
for (int j = 0; j < w; j++) {
idx = *(lpDIBBits + i * lineByte + j);
*(lpDIBBits + i * lineByte + j) = value[idx];
}
}
::GlobalUnlock(dib);
Invalidate(); //區(qū)域重繪
}
③ 大津閾值
void CDemoView::Onotsu()
{
CDemoDoc *pDoc = GetDocument();
HDIB dib = pDoc->GetHDIB();
LPSTR lpDIB = (LPSTR)::GlobalLock((HGLOBAL)dib);
LPBITMAPINFOHEADER phead = (LPBITMAPINFOHEADER)lpDIB; //信息頭指針
int biByteCount = phead->biBitCount / 8; //每像素字節(jié)數(shù)
if (biByteCount == 3)
{
::GlobalUnlock(dib);
return;
} //彩色圖不做處理
int w = phead->biWidth;
int h = phead->biHeight;
int lineByte = (biByteCount * w + 3) / 4 * 4;
unsigned char *lpDIBBits = (unsigned char *)::FindDIBBits(lpDIB); //像素位置指針,unsigned char取值范圍是0~255(一個(gè)字節(jié))
//統(tǒng)計(jì)直方圖
float p[256];
for (int i = 0; i < 256; i++) {
p[i] = 0;
}
unsigned char idx;
for (int i = 0; i < h; i++) {
for (int j = 0; j < w; j++) {
idx = *(lpDIBBits + i * lineByte + j);
p[idx] += 1.0 / (w*h);
}
}
//計(jì)算全圖均值
float u = 0;
for (int i = 0; i < 256; i++) {
u += p[i] * i;
}
//統(tǒng)計(jì)類間最大方差diff(th)
float diff = 0;
float diffmax = 0;
int T = 0;
for (int th = 1; th < 256; th++) { //分為小于和大于等于兩部分
float w0 = 0;
float w1 = 0;
float u0 = 0;
float u1 = 0;
//w0 and u0
for (int i = 0; i < th; i++) {
w0 += p[i];
u0 += p[i] * i;
}
if (w0 == 0.0) //注意這里如果前半部分沒有像素，就不能當(dāng)作閾值
continue;
else
u0 = u0 / w0;
//w1 and u1
for (int i = th; i < 256; i++) {
u1 += p[i] * i;
}
w1 = 1.0 - w0;
if (w1 == 0.0)
continue;
else
u1 = u1 / w1;
//類間方差diff
diff = w0 * pow((u0 - u), 2) + w1 * pow((u1 - u), 2);
//diff = w0 * w1 * pow((u0 - u1), 2);
if (diffmax < diff) {
diffmax = diff;
T = th;
}
}
//閾值分割
for (int i = 0; i < h; i++) {
for (int j = 0; j < w; j++) {
if (*(lpDIBBits + i * lineByte + j) < T)
*(lpDIBBits + i * lineByte + j) = 0;
else
*(lpDIBBits + i * lineByte + j) = 255;
}
}
::GlobalUnlock(dib);
Invalidate(); //區(qū)域重繪
}
??：
（1）大津閾值算法中饶氏，注意遍歷的時(shí)候如果小于原圖灰度最小值或者大于灰度最大值，則不扔泄础疹启；
（2）大津閾值算法中，兩部分的均值計(jì)算要思考正確蔼卡；
（3）C語(yǔ)言給變量賦值的時(shí)候要考慮好變量的數(shù)據(jù)類型喊崖；
（4）在形如x += *的代碼中，要處理好x的初值；
（5）傅里葉反變換去噪時(shí)荤懂，實(shí)部和虛部都要做處理茁裙。
麻瓜寫于2018/11/30
