二维卷积c代码

二维卷积c代码

二维信号的卷积原理请参考另外一篇文章：http://blog.csdn.net/carson2005/article/details/43702241

这里直接给出参考代码：

void Conv2(int** filter, int** arr, int** res, int filterW, int filterH, int arrW, int arrH)
{
    int temp;

    for (int i=0; i<filterH+arrH-1; i++)
    {
        for (int j=0; j<filterW+arrW-1; j++)
        {
            temp = 0;
            for (int m=0; m<filterH; m++)
            {
                for (int n=0; n<filterW; n++)
                {
                    if ((i-m)>=0 && (i-m)<arrH && (j-n)>=0 && (j-n)<arrW)
                    {
                        temp += filter[m][n]*arr[i-m][j-n];
                    }
                }
            }
            res[i][j] = temp;
        }
    }
}

void Conv2Test()
{
    int filterW = 3;
    int filterH = 3;
    int arrW = 5;
    int arrH = 5;
    int resW = filterW + arrW - 1;
    int resH = filterH + arrH - 1;

    int** pFilter = new int*[filterH];
    for (int i=0; i<filterH; i++)
    {
        pFilter[i] = new int[filterW];
    }

    int** arr = new int*[arrH];
    for (int i=0; i<arrH; i++)
    {
        arr[i] = new int[arrW];
    }

    int** res = new int*[resH];
    for (int i=0; i<resH; i++)
    {
        res[i] = new int[resW];
    }

    pFilter[0][0] = 1;
    pFilter[0][1] = 3;
    pFilter[0][2] = 1;
    pFilter[1][0] = 0;
    pFilter[1][1] = 5;
    pFilter[1][2] = 0;
    pFilter[2][0] = 2;
    pFilter[2][1] = 1;
    pFilter[2][2] = 2;

    arr[0][0] = 17;
    arr[0][1] = 24;
    arr[0][2] = 1;
    arr[0][3] = 8;
    arr[0][4] = 15;
    arr[1][0] = 23;
    arr[1][1] = 5;
    arr[1][2] = 7;
    arr[1][3] = 14;
    arr[1][4] = 16;
    arr[2][0] = 4;
    arr[2][1] = 6;
    arr[2][2] = 13;
    arr[2][3] = 20;
    arr[2][4] = 22;
    arr[3][0] = 10;
    arr[3][1] = 12;
    arr[3][2] = 19;
    arr[3][3] = 21;
    arr[3][4] = 3;
    arr[4][0] = 11;
    arr[4][1] = 18;
    arr[4][2] = 25;
    arr[4][3] = 2;
    arr[4][4] = 9;

    printf("pFilter: 
");
    for (int i=0; i<filterH; i++)
    {
        for (int j=0; j<filterW; j++)
        {
            printf("%d ", pFilter[i][j]);
        }
        printf("
");
    }
    printf("*********************************************** 
");

    printf("arr: 
");
    for (int i=0; i<arrH; i++)
    {
        for (int j=0; j<arrW; j++)
        {
            printf("%d ", arr[i][j]);
        }
        printf("
");
    }
    printf("*********************************************** 
");


    Conv2(pFilter, arr, res, filterW, filterH, arrW, arrH);

    printf("res: 
");
    for (int i=0; i<resH; i++)
    {
        for (int j=0; j<resW; j++)
        {
            printf("%d ", res[i][j]);
        }
        printf("
");
    }
    printf("*********************************************** 
");
}