概念:结构元素B(自己理解,百度很多)。
膨胀:就是先对B做关于原点的映射,在将映射连续的在A上移动而实现。
本文中结构元素采用:章毓晋老师的《图像分析》(第三版),285页
#include <stdio.h>
#include <stdlib.h>
#include <Windows.h>
int main(int* argc, char** argv)
{
FILE* fp = fopen("./threshold.bmp", "rb");
if (fp == 0)
return 0;
BITMAPFILEHEADER fileHead;
fread(&fileHead, sizeof(BITMAPFILEHEADER), 1, fp);
BITMAPINFOHEADER infoHead;
fread(&infoHead, sizeof(BITMAPINFOHEADER), 1, fp);
int width = infoHead.biWidth;
int height = infoHead.biHeight;
int biCount = infoHead.biBitCount;
int lineByte = (biCount*width / 8 + 3) / 4 * 4;
RGBQUAD* pColorTable;
pColorTable = new RGBQUAD[256];
fread(pColorTable, sizeof(RGBQUAD), 256, fp);
unsigned char* pBmpBuf;
pBmpBuf = new unsigned char[lineByte*height];
fread(pBmpBuf, lineByte*height, 1, fp);
fclose(fp);
// 新图
FILE* fop = fopen("dilate.bmp", "wb");
if (fop == 0)
return 0;
// 定义结构元素
// 腐蚀操作
int t = 0, d = 0, r = 0;
for (int i = 1; i < height; ++i){
for (int j = 0; j < width-1; ++j){
t = *(pBmpBuf + i*lineByte + j); // 当前点
d = *(pBmpBuf + (i - 1)*lineByte + j); // 下面点
r = *(pBmpBuf + i*lineByte + j + 1); // 右边点
if (t == 0 && d != 0){
*(pBmpBuf + (i - 1)*lineByte + j) = 0;//下边的置位1
}
if (t == 0 && r != 0){
*(pBmpBuf + i*lineByte + j+1) = 0;//右边的置位1
j = j + 1;
}
}
}
fwrite(&fileHead, sizeof(BITMAPFILEHEADER), 1, fop);
fwrite(&infoHead, sizeof(BITMAPINFOHEADER), 1, fop);
fwrite(pColorTable, sizeof(RGBQUAD), 255, fop);
fwrite(pBmpBuf, lineByte*height, 1, fop);
fclose(fop);
system("pause");
return 0;
}
实验结果:
二值图像(原图) 膨胀后
实验结果分析:
本实验中,判断的时候并不是对255灰度值判断,而是对0灰度判断(其实两者一样,只是看起来好看点),所以是对0灰度值进行膨胀。实验结果中,左眼旁边的头发部分、右眼旁边的头发和胡须部分都有所填充。
后续持续更新用C语言实现图像处理算法,敬请期待,欢迎关注。