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图像处理（开运算）

猜你喜欢