opencv下面滤波程序最后两行没看懂什么意思,望各路前辈解答

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;

http://blog.csdn.net/WINCOL/article/details/6103884
高斯滤波器就是抑制服从正太分布的噪声么