以下用于图像调整大小的 cpp Bucubic 插值代码有什么问题

Question

我正在尝试使用双三次插值对图像进行上采样，我需要与 opencv 的 cvResize() 函数匹配的准确值，但是以下代码的结果与 cvResize() 的结果不匹配，您可以看一下并帮助我修复错误。

Image *Image::resize_using_Bicubic(int w, int h) {
    float dx,dy;
    float x,y;
    float tx, ty;
    float i,j,m,n;


Image *result=new Image(w, h);
tx = (float)this->m_width/(float)w;
ty = (float)this->m_height/(float)h;
for(i=0; i< w; i++)
{
    for(j=0; j< h; j++)
    {
        x = i*tx;
        y = j*ty;
        dx = float(i-x)-(int)(i-x);
        dy = float(j-y)-(int)(j-y);
        float temp=0.0;
        for(m=-1;m<=2;m++)
        {
            for(n=-1;n<=2;n++)
            { 

                int HIndex,WIndex;
                HIndex=(y+n);
                WIndex=(x+m);

                if (HIndex<0) {
                    HIndex=0;
                }
                else if(HIndex>this->getHeight())
                {
                    HIndex=this->getHeight()-1;

                }
                if (WIndex<0) {
                    WIndex=0;
                }
                else if(WIndex>this->getWidth())
                {
                    WIndex=this->getWidth()-1;

                }


                temp+=this->getPixel(HIndex,WIndex)*R(m-dx)*R(dy-n);
            }
        }

        result->setPixel(j, i, temp);
    }
}
return result;

}

Answer 1

你还没有说结果有多么不同。如果它们非常接近，例如每个 RGB 通道中的 1 或 2 以内，则可以简单地通过舍入差异来解释。

双三次插值的算法不止一种。 Don Mitchell 和 Arun Netravali 进行了分析，并提出了一个公式来描述其中的一些：http:// /www.mentallandscape.com/Papers_siggraph88.pdf

编辑：还有一件事，各个滤波器系数应该相加并用于除以最终值最后，对值进行标准化。我不确定为什么一个是 m-dx ，另一个是 dy-n ，它们不应该是相同的符号吗？

r=R(m-dx)*R(dy-n);
r_sum+=r;
temp+=this->getPixel(HIndex,WIndex)*r;
. . .
result->setPixel(j, i, temp/r_sum);

Answer 2

将：更改

            else if(HIndex>this->getHeight())

为：

            else if(HIndex >= this->getHeight())

并将：更改

            else if(WIndex>this->getWidth())

为：

            else if(WIndex >= this->getWidth())

编辑

还将：更改

    for(m=-1;m<=2;m++)
    {
        for(n=-1;n<=2;n++)

为：

    for(m = -1; m <= 1; m++)
    {
        for(n = -1; n <= 1; n++)