C++博客-bobokids-文章分类-图像处理算法基础

滤波器类

晾衣竿 — Sat, 27 Sep 2008 19:57:00 GMT

BOOL   CFilter::Filter(const   int   nDefTemplate,BYTE   *pBuffer,DWORD   Size,DWORD   Width,   DWORD   Height)
  {

  SetTemplate(nDefTemplate);
  if(!pBuffer)//Buf为空
  {
  return   false;
  }

  if((Width+16)*Height<Size||(Width-16)*Height>Size) //非8位色彩的位图
  {
  //MYTODO:转化到8位
  //应该在CBitmap类中RGB->YUV,只取亮度信息
  AfxMessageBox("不支持非256色位图");
  //return   false;
  }
  Multiply(pBuffer,Size,Width,Height);
  return   true;

  }


  void   CFilter::SetTemplate(const   int   nDefTemplate)
  {
  if(m_pTemplate)
  {
  delete[]   m_pTemplate;
  }

  switch(nDefTemplate)
  {
  case   TEMPLATE_LOG:
  {
  m_TemplateWidth=5;
  m_fCoef=1;
  BYTE   btTemplateSize=m_TemplateWidth*m_TemplateWidth;
  m_pTemplate=new   char[btTemplateSize];
  char   btTemp[25]={
  -2,-4,-4,-4,-2,
  -4,   0,   8,   0,-4,
  -4,   8,24,   8,-4,
  -4,   0,   8,   0,-4,
  -2,-4,-4,-4,-2};
  ::memcpy(m_pTemplate,&btTemp,btTemplateSize);
  }
  break;
  case   TEMPLATE_LAPLACIAN:
  {
  m_TemplateWidth=3;
  m_fCoef=1;
  BYTE   btTemplateSize=m_TemplateWidth*m_TemplateWidth;
  m_pTemplate=new   char[btTemplateSize];
  char   btTemp[9]={-1,-1,-1,-1,9,-1,-1,-1,-1};

  ::memcpy(m_pTemplate,&btTemp,btTemplateSize);
  }
  break;
  case   TEMPLATE_GAUSS:
  {
  m_TemplateWidth=3;
  m_fCoef=1/16.0; //不能1/16,否则会变成0
  BYTE   btTemplateSize=m_TemplateWidth*m_TemplateWidth;
  m_pTemplate=new   char[btTemplateSize];
  char   btTemp[9]={1,2,1,2,4,2,1,2,1};

  ::memcpy(m_pTemplate,&btTemp,btTemplateSize);
  }
  break;
  default:
  TRACE("没有该默认的模板");
  }

  }

  void   CFilter::Multiply(BYTE   *pBuffer,   DWORD   Size,   DWORD   Width,   DWORD   Height)
  {
  int   nBorder=m_TemplateWidth/2; //需要空开的大小
  DWORD   dwOffset=0;
  float   fResult=0;
  int   nTempNum=0;
  BYTE*   pTempBuffer=new   BYTE[Size];
  for(DWORD   y=0+nBorder;y<Height-nBorder;y++) //因为窗口中心点要受模板大小影响
  for(DWORD   x=0+nBorder;x<Width-nBorder;x++)
  {
  dwOffset=Size-Width-y*Width+x;
  for(int   i=-nBorder;i<=nBorder;i++)
  for(int   j=-nBorder;j<=nBorder;j++)
  {

  nTempNum=m_pTemplate[i*m_TemplateWidth+j+nBorder*m_TemplateWidth+nBorder];
  fResult+=pBuffer[dwOffset+(-i)*Width+j]*nTempNum;
  m_lCount++;
  }


  fResult*=m_fCoef;
  pTempBuffer[dwOffset]=CutOverflow(fResult);
  fResult=0;
  }
  ::memcpy(pBuffer,pTempBuffer,Size);
  delete[]   pTempBuffer;

  }

  BYTE   CFilter::CutOverflow(float   fNum)
  {
  if(fNum>255.0)
  return   (BYTE)255;
  else   if(fNum<0.0)
  return   (BYTE)0;
  else
  return   (BYTE)fNum;
  }

晾衣竿 2008-09-28 03:57 发表评论

灰度拉伸

晾衣竿 — Sat, 27 Sep 2008 19:31:00 GMT

* 函数名称:
* GrayStrech()
*
* 参数:
* LPSTR  lpDIBBits  -指向原DIB图像的指针
* LONG  lWidth        -原图像的宽度(像素)
* LONG  lHeight  -原图像的高度(像素)
* BYTE  bX1  -灰度拉伸第一个点的X坐标
* BYTE  bY1  -灰度拉伸第一个点的Y坐标
* BYTE  bX2  -灰度拉伸第二个点的X坐标
* BYTE  bY2  -灰度拉伸第二个点的Y坐标
*
* 返回值:
* BOOL    -成功则返回TRUE,否则返回FALSE
*
* 说明:
* 该函数用来对图像进行灰度拉伸,改善输出图像
*
*************************************************************************/

BOOL WINAPI GrayStrech(LPSTR lpDIBBits,LONG lWidth,LONG lHeight,BYTE bX1,
        BYTE bY1,BYTE bX2,BYTE bY2)
{
unsigned char *lpSrc; //指向原图像的指针
LONG i,j; //循环变量
BYTE bMap[256];  //灰度映射表
LONG lLineBytes; //图像每行的字节数

lLineBytes = WIDTHBYTES(lWidth*8); //计算图像每行的字节数
//计算灰度映射表
for (i = 0;i <= bX1;i++){
  //判断bX1是否大于0(防止分母为0)
  if (bX1 > 0) {
  //线性变换
  bMap[i] = (BYTE) bY1*i/bX1;
  }
  else {
  //直接置0
  bMap[i] = 0;
  }

}
for (;i <= bX2;i++){
  //判断bX1是否等于bX2(防止分母为0)
  if (bX2 != bX1) {
  //线性变换
  bMap[i] = bY1 + (BYTE)((bY2 - bY1)*(i - bX1)/(bX2 - bX1));
  }
  else {
  //直接置为bY1
  bMap[1] = bY1;
  }
}
for (;i < 256;i++) {
  //判断bX2是否等于255(防止分母为0)
  if (bX2 != 255){
  //线性变换
  bMap[i] = bY2 + (BYTE)((255 - bY2)*(i - bX2)/(255 - bX2));
  }
  else {
  //直接置为255
  bMap[i] = 255;
  }
}
//每行
for (i = 0;i < lHeight;i++){
  //每列
  for (j = 0;j < lWidth;j++) {
  //指向DIB第i行第j个像素指针
  lpSrc = (unsigned char*)lpDIBBits + lLineBytes*(lHeight - 1 - i) + j;
  //计算新的灰度值
  *lpSrc = bMap[*lpSrc];
  }
}

return TRUE;
}

晾衣竿 2008-09-28 03:31 发表评论

直方图均衡子函数

晾衣竿 — Sat, 27 Sep 2008 19:21:00 GMT

/*************************************************************************
                                直方图均衡子函数
    ************************************************************************/
  BOOL   WINAPI   InteEqualize(LPSTR   lpDIBBits,   LONG   lWidth,   LONG   lHeight)
  {

  //   指向源图像的指针
  unsigned   char* lpSrc;

  //   临时变量
  LONG lTemp;

  //   循环变量
  LONG i;
  LONG j;

  //   灰度映射表
  BYTE bMap[256];

  //   灰度映射表
  LONG lCount[256];

  //   图像每行的字节数
  LONG lLineBytes;

  //   计算图像每行的字节数
  lLineBytes   =   WIDTHBYTES(lWidth   *   8);

  //   重置计数为0
  for   (i   =   0;   i   <   256;   i   ++)
  {
  //   清零
  lCount[i]   =   0;
  }

  //   计算各个灰度值的计数
  for   (i   =   0;   i   <   lHeight;   i   ++)
  {
  for   (j   =   0;   j   <   lWidth;   j   ++)
  {
  lpSrc   =   (unsigned   char   *)lpDIBBits   +   lLineBytes   *   i   +   j;

  //   计数加1
  lCount[*(lpSrc)]++;
  }
  }

  //   计算灰度映射表
  for   (i   =   0;   i   <   256;   i++)
  {
  //   初始为0
  lTemp   =   0;

  for   (j   =   0;   j   <=   i   ;   j++)
  {
  lTemp   +=   lCount[j];
  }

  //   计算对应的新灰度值
  bMap[i]   =   (BYTE)   (lTemp   *   255   /   lHeight   /   lWidth);
  }

  //   每行
  for(i   =   0;   i   <   lHeight;   i++)
  {
  //   每列
  for(j   =   0;   j   <   lWidth;   j++)
  {
  //   指向DIB第i行，第j个象素的指针
  lpSrc   =   (unsigned   char*)lpDIBBits   +   lLineBytes   *   (lHeight   -   1   -   i)   +   j;

  //   计算新的灰度值
  *lpSrc   =   bMap[*lpSrc];
  }
  }

  //   返回
  return   TRUE;
  }

晾衣竿 2008-09-28 03:21 发表评论

图象的平滑(去噪声)、锐化

晾衣竿 — Sat, 27 Sep 2008 18:14:00 GMT

摘要: 平滑先举个例子说明一下什么是平滑(smoothing)，如下面两幅图所示：可以看到，图3.2比图3.1柔和一些(也模糊一些)。是不是觉得很神奇？其实实现起来很简单。我们将原图中的每一点的灰度和它周围八个点的灰度相加，然后除以9，作为新图中对应点的灰度，就能实现上面的效果。 ... 阅读全文

晾衣竿 2008-09-28 02:14 发表评论

图象的几何变换

晾衣竿 — Sat, 27 Sep 2008 18:09:00 GMT

摘要: 几何变换，包括图象的平移、旋转、镜象变换、转置、放缩等。如果你熟悉矩阵运算，你将发现，实现这些变换是非常容易的。平移平移(translation)变换大概是几何变换中最简单的一种了。如图2.1所示，初始坐标为(x0,y0)的点经过平移(tx,ty)(以向右，向下为正方向)后，坐标变为(x1,y1)。这两点之间的关系是x1=x0+tx ，y1=y0+ty。图2.1 &... 阅读全文

晾衣竿 2008-09-28 02:09 发表评论

位图和调色板的概念

晾衣竿 — Sat, 27 Sep 2008 17:56:00 GMT

摘要: 普通的显示器屏幕是由许许多多点构成的，我们称之为象素。显示时采用扫描的方法：电子枪每次从左到右扫描一行，为每个象素着色，然后从上到下这样扫描若干行，就扫过了一屏。为了防止闪烁，每秒要重... 阅读全文

晾衣竿 2008-09-28 01:56 发表评论

对比度扩展

晾衣竿 — Sat, 27 Sep 2008 17:45:00 GMT

摘要: 假设有一幅图，由于成象时光照不足，使得整幅图偏暗(例如，灰度范围从0到63)；或者成象时光照过强，使得整幅图偏亮(例如，灰度范围从200到255)，我们称这些情况为低对比度，即灰度都挤在一起，没有拉开。灰度扩展的意思就是把你所感性趣的灰度范围拉开，使得该范围内的象素，亮的越亮，暗的越暗，从而达到了增强对比度的目的。我们可以用图5.5来说明对比度扩展(contrast stretching)的原理。... 阅读全文

晾衣竿 2008-09-28 01:45 发表评论

灰度直方图均衡化

晾衣竿 — Sat, 27 Sep 2008 17:34:00 GMT

摘要: 在介绍灰度直方图均衡化(histogram equalization)之前，先讲讲直方图修正。所谓直方图修正，就是通过一个灰度映射函数Gnew=F(Gold)，将原灰度直方图改造成你所希望的直方图。所以，直方图修正的关键就是灰度映射函数。我们刚才介绍的阈值化、削波、灰度窗口变换等等，都是灰度映射函数。直方图均衡化是一种最常用的直方图修正。它是把给定图象的直方图分布改造成均匀直方图分布。由信息学... 阅读全文

晾衣竿 2008-09-28 01:34 发表评论