第一节--滤波操作的VC++实现

//-----------------------------------------【滤波操作的VC++实现】--------------------------------------
//1--下面为CImgProcess类添加用于--线性滤波的模板操作，函数Template类似于Matlab中的imfilter方法，它支持任意
// 大小的模板,它支持任意大小的模板(nTempH*nTempW)并且提供了更大的模板设定的自由度
//2--模板的中心--可以由参数int nTempMX和nTempMY任意指定
//3--参数FLOAT *pfArray为指向模板数组的指针，而参数FLOAT fCoef为模板的系数-即例子5.1中的1/9
//4--尽管模板本身是一个二维数组，但是程序实现中，我们将它按行存储作为一个一位数组对待，引用起来同样十分方便，
// 如pfArray[k*nTempW+L]表示模板的第K行的第L个元素
//5--Template函数在一开始就将目标图像初始化为黑色，并且,并且没有填充图像的边界，才有5.2.2小结中，提到的第一种
// 收缩边界的策略，因此，经过函数Template滤波的输出图像的周围都有很细的一圈黑边，如图5.6(f)所示。对于3*3的
// 模板，黑边的宽度为1个像素
//------------------------------------------------------------------------------------------------------
/*******************************************************************************************************
*函数原型--void CImgProcess::Template(CImgProcess* pTo,int nTempH,int nTempW,int nTempMY,int nTempMX,
* FLOAT* pfArray,FLOAT fCoef)
*函数功能:模板操作--该函数用于指定模板的大小(任意大小)来对图像进行操作
*函数参数:1--CImgProcess* pTo---输出图像的CImgProcess指针
* 2--int nTempH---模板的高度
* 3--int nTempW---模板的宽度
* 4--int nTempMY---模板中心元素的Y坐标(<=iTempH-1)
* 5--int nTempMX---模板中心元素的X坐标(<=iTempW-1)
* 6--FLOAT* pfArray--指向模板数组的指针
* 7--FLOAT fCoef-----模板系数
*函数返回值:无
********************************************************************************************************/
//1--GetHeight()--返回CImg实例中的图像每列的像素数目，即纵向分辨率或高度
//2--GetGray()----返回指定坐标位置像素的灰度值
//3--SetPixel()设定指定坐标位置像素的颜色值
void CImgProcess::Template(CImgProcess *pTo,int nTempH,int nTempW,int nTempMY,int nTempMX,
FLOAT *pfArray,FLOAT fCoef)
{
//目标图像的初始化
pTo->InitPixels(0);
//循环变量
int i,j;
//扫描图像进行模板操作
for(i=nTempMY;i<GetHeight()-(nTempH-nTempMY)+1;i++)//高度--相当于先行扫描行，一行一行的进行扫描
{
for(j=nTempMX;j<GetWidthPixel()-(nTempW-nTempMX);j++)//宽度--相当于列扫描
{
//(j,i为中心点)
float fResult=0;
for(int k=0;k<nTempH;k++)
{
for(int l=0;l<nTempW;l++)
{
//计算加权和
fResult+=GetGray(j+1-nTempMX,i+k-nTempMX)*pfArray[k*nTempW+1];//此块计算出来的相当于----梯度幅值
}
}
//乘以系数
fResult*=fCoef;
//取正
//fabs()---取一个数的绝对值---#include<math,h>
fResult=(float)fabs(fResult);//锐化时，可能出现负值
BYTE byte;
if(fResult>255)
byte=255;
else
byte=fResult+0.5;//四舍五入
pTo->SetPixel(j,i,RGB(byte,byte,byte));//因此，最后输出的图像是一个---灰度值范围为0~255的灰度图
}//for j
}//for i
}