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

Z时代
2024-01-10
分类：IT

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

运算：用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;
}