C语言实现BMP图像开运算处理

开运算可以把比结构元素小的特定图像细节出去，同时保证不产生全局的几何失真。滤掉比结构元素小的突刺，切断细长搭接而起到分离作用。

运算：用B开启A就是选出了A中某些与B相匹配的点，这些点可由完全包含在A中的结构元素B的平移得到。也就是先腐蚀后加膨胀。

#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("open.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;
 }
}
}
// 膨胀操作
// 初始化
unsigned char* pBmpBuf2;
pBmpBuf2 = new unsigned char[lineByte*height];
for (int i = 0; i < height; ++i){
for (int j = 0; j < width; ++j){
 *(pBmpBuf2 + i*lineByte + j) = 255;
}
}

// 腐蚀操作
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 && r == 0){
  *(pBmpBuf2 + i*lineByte + j) = 0;       // 当前点
 }
}
}
// 结构元素向上反转180度，对最下面一排处理
for (int j = 0; j < width - 1; ++j){
t = *(pBmpBuf + j);            // 当前点
d = *(pBmpBuf + lineByte + j); // 上面点
r = *(pBmpBuf + j + 1);        // 右边点
if (t == 0 && d == 0 && r == 0){
 *(pBmpBuf2 + j) = 0;       // 当前点
}
}
// 结构元素向右反转，对最右边一列处理
for (int i = 1; i < height; ++i){
t = *(pBmpBuf + i*lineByte + width - 1);
d = *(pBmpBuf + (i - 1)*lineByte + width - 1);
r = *(pBmpBuf + i*lineByte + width - 2);
if (t == 0 && d == 0 && r == 0){
 *(pBmpBuf2 + i*lineByte + width - 1) = 0;       // 当前点
}
}

fwrite(&fileHead, sizeof(BITMAPFILEHEADER), 1, fop);
fwrite(&infoHead, sizeof(BITMAPINFOHEADER), 1, fop);
fwrite(pColorTable, sizeof(RGBQUAD), 255, fop);
fwrite(pBmpBuf2, lineByte*height, 1, fop);
fclose(fop);

system("pause");
return 0;
}

实验结果：

C语言实现BMP图像开运算处理