double D0=30;
void ILPF(CvMat* src, const double D0)
{
int i, j;
int state = -1;
double tempD;
long width, height;
width = src->width;
height = src->height;
long x, y;
x = width / 2;
y = height / 2;
CvMat* H_mat;
H_mat = cvCreateMat(src->height,src->width, CV_64FC2);
for(i = 0; i < height; i++)
{
for(j = 0; j < width; j++)
{
if(i > y && j > x)
{
state = 3;
}
else if(i > y)
{
state = 1;
}
else if(j > x)
{
state = 2;
}
else
{
state = 0;
}
switch(state)
{
case 0:
tempD = (double) (i * i + j * j);tempD = sqrt(tempD);break;
case 1:
tempD = (double) ((height - i) * (height - i) + j * j);tempD = sqrt(tempD);break;
case 2:
tempD = (double) (i * i + (width - j) * (width - j));tempD = sqrt(tempD);break;
case 3:
tempD = (double) ((height - i) * (height - i) + (width - j) * (width - j));tempD = sqrt(tempD);break;
default:
break;
}
//二维高斯低通滤波器传递函数
tempD = exp(-0.5 * pow(tempD / D0, 2));
((double*)(H_mat->data.ptr + H_mat->step * i))[j * 2] = tempD;
((double*)(H_mat->data.ptr + H_mat->step * i))[j * 2 + 1] = 0.0;