c - 单次图像的均值和方差

Question

我正在尝试在 opencv 中使用图像 (hXw) 上的 3X3 窗口计算均值和方差...这是我的代码...这有任何准确性问题吗？或者是否有任何其他有效的方法可以在一次通过。？

int pi,a,b;

for(i=1;i<h-1;i++)
{
    for(j=1;j<w-1;j++)
    {   int sq=0,sum=0;
        double mean=0;
        double var=0;
        for(a=-1;a<=1;a++)
        {
            for(b=-1;b<=1;b++)
            {
                pi=data[(i+a)*step+(j+b)];
                sq=pi*pi;
                sum=sum+sq;
                mean=mean+pi;
            }
        }
        mean=mean/9;
        double soa=mean*mean;//square of average
        double aos=sum/9;//mean of squares
        double var=aos-soa;//variance
    }
}

Answer 1

关于计算效率，我建议在傅立叶域而不是时间(图像)域中使用卷积来执行此操作。请记住，卷积是傅立叶域中的简单乘法。就像在时间序列中，频谱密度函数是方差分解为频率的函数一样，可以将其扩展到图像的二维。应该比嵌套的 for 循环好得多。

我现在没有密码。但这种技术已被用于对象检测或图像配准的“快速模板匹配”等算法中。