邻域平均法

 /****************************************
******************************
*
* 函数名称：
*     AvgTemplate(int inputH, int inputW, int
inputX0, int inputY0, float 
* 
* pTemplate, float fCoef)
*
* 参数：
*     int inputH-模板的高度
*     int inputW-模板的宽度
*     int inputX0-模板的中心元素X坐标
*     int inputY0-模板的中心元素Y坐标
*     float *pTemplate-指向模板数组的指针
*     float f-模板系数
*
* 返回值：
*     void
*
* 说明：
*     该函数用指定的模板（任意大小）来对图像进行操作，
参数inputH指定模板
*     的高度，参数inputW指定模板的宽度，
参数inputX0和iinputY0指定模板的中心
*     元素坐标，参数pTemplate指定模板元素，f指定系数
*
********************************************
**************************/
void CImgEnhance::AvgTemplate(int inputH,
int inputW, int inputX0, int inputY0, 
float *pTemplate, float f)
{
unsigned char *pSrc, *pDst;   
int i,j,k,l;
float value;
if(m_pImgDataOut != NULL)
{
delete []m_pImgDataOut;
m_pImgDataOut = NULL;
}
 int lineByte = (m_imgWidth *
m_nBitCount / 8 + 3) / 4 * 4;

if(m_nBitCount != 8)
{
AfxMessageBox("只能处理8位灰度图像!");
return ;
}
//创建要复制的图像区域
m_nBitCountOut = m_nBitCount;
int lineByteOut = (m_imgWidth*m_
nBitCountOut / 8 + 3) / 4 * 4;
if (!m_pImgDataOut)
{
m_pImgDataOut = new unsigned char
[lineByteOut*m_imgHeight];
}

int pixelByte = m_nBitCountOut / 8;
for(i = 0; i < m_imgHeight; i++){
for(j = 0; j < m_imgWidth * pixelByte; j++)
*(m_pImgDataOut+i*lineByteOut+j)= *
(m_pImgData+i*lineByteOut+j);
}

//行处理（去掉边缘几行）
for (i = inputY0; i < m_imgHeight - 
inputH + inputY0 + 1; i++)
{
//列处理（去掉边缘几列）
for (j = inputX0; j < m_imgWidth - 
inputW + inputX0 + 1; j++)
{
//指向新DIB第i行第j列的像素的指针
pDst = m_pImgDataOut + lineByte * 
(m_imgHeight - 1 - i) + j;

value=0;
//计算
for (k = 0; k < inputH; k++)
{
for (l = 0; l < inputW; l++)
{
pSrc = m_pImgData + lineByte * 
(m_imgHeight - 1 - i 
+ inputY0 - k) + j - inputX0 + l;
//计算加权平均，保存像素值
value += (*pSrc) * pTemplate[k * inputW + l];
}
}
//乘以系数
value*=f;
//取结果的绝对值
value=(float)fabs(value);
if (value > 255)
{
*pDst = 255; 
}
else
{
*pDst = (unsigned char)(value + 0.5);
}
}
}
} 
}

void CDemoView::OnAveSmooth() 
{
CDemoDoc *pDoc=GetDocument();
ImgCenterDib *pDib=pDoc->GetPDib();

if(pDib->m_nBitCount!=8&
::MessageBox(0,"只处理灰度图像",MB_OK,0);
return ;
}
AvgDialog *dlg=new AvgDialog();
dlg->DoModal();
int TempH, TempW, TempCX, TempCY;
float fCoef;
float Temp[81];
int i;
//设模板为平均模板
for (i = 0; i < 81; i++)
{
Temp[i]=1.0;
}
switch(dlg->m_nIndex)
{
case 0:
//初始化模板窗口
TempH=3;
TempW=3;
TempCX=1;
TempCY=1;
fCoef=(float)(1.0/9.0);
break;
case 1:
//初始化模板窗口
TempH=5;
TempW=5;
TempCX=2;
TempCY=2;
fCoef=(float)(1.0/25.0);
break;
case 2:
//初始化模板窗口
TempH=9;
TempW=9;
TempCX=4;
TempCY=4;
fCoef=(float)(1.0/81.0);
break;
default:
break;
}
CImgEnhance imgnoise(pDib->GetDimensions
(),pDib->m_nBitCount,
pDib->m_lpColorTable, pDib->m_pImgData);
imgnoise.AvgTemplate(TempH,TempW,
TempCX,TempCY,Temp,fCoef);

CMainFrame* pFrame = (CMainFrame *)
(AfxGetApp()->m_pMainWnd);
pFrame->SendMessage(WM_COMMAND, ID_FILE_NEW);

CDemoView* pView=(CDemoView*)pFrame->
MDIGetActive()->GetActiveView();
CDemoDoc* pDocNew=pView->GetDocument();
ImgCenterDib *dibNew=pDocNew->GetPDib();

dibNew->ReplaceDib(imgnoise.GetDimensions
(),imgnoise.m_nBitCountOut,imgnoise.
m_lpColorTable, imgnoise.m_pImgDataOut);
pDocNew->SetModifiedFlag(TRUE);
pDocNew->UpdateAllViews(pView);

Invalidate();
return;
}