• 应用连续高斯模糊后得到的σ是多少?


    在这里插入图片描述
    对图像应用多个连续的高斯模糊效果与应用单个较大的高斯模糊效果相同,后者的半径是实际应用的模糊半径的平方和的平方根。例如,应用半径为6688的连续高斯模糊与应用半径为1010的单个高斯模糊产生的结果相同,因为62+82=10sqrt{6^2+8^2} =10

    但是我找不到任何证据,为什么会这样呢?

    而且我还发现,在某些代码中,人们会考虑将两个连续的高斯模糊 σ1sigma_1σ2sigma_2用只是一个模糊σ=σ12+σ22sigma=sqrt{sigma_1^2+sigma_2^2}代替。

    我们如何证明这一结论?


    证明如下:


    可以使用一些简单的卷积理论来证明。首先,回想一下卷积:f(gh)=(fg)hf *(g * h)=(f * g) * h
    接下来回想一下高斯模糊图像II就是简单地对高斯核函数GG进行卷积,G(x,yσ)=(2πσ2)1exp(x2+y22σ2)G(x, y | sigma)=left(2 pi sigma^{2} ight)^{-1} exp left(-frac{x^{2}+y^{2}}{2 sigma^{2}} ight)

    所以高斯模糊两次就相当于卷积两次:IB=G1(G2I)=(G1G2)I=GII_{B}=G_{1} *left(G_{2} * I ight)=left(G_{1} * G_{2} ight) * I=G * I,我们知道GG是一个高斯核,因为两个高斯的卷积是一个高斯。

    现在我们只需要证明:G(x,yσ)=G(x,yσ12+σ22)=G1(x,yσ1)G2(x,yσ2)G(x, y | sigma)=G(x, y | sqrt{sigma_{1}^{2}+sigma_{2}^{2}})=G_{1}left(x, y | sigma_{1} ight) * G_{2}left(x, y | sigma_{2} ight)

    一种方法是通过定义:计算
    G(x,yσ)=G1(τ,ξσ1)G2(xτ,yξσ2)dτdξG(x, y | sigma)=iint_{-infty}^{infty} G_{1}left( au, xi | sigma_{1} ight) * G_{2}left(x- au, y-xi | sigma_{2} ight) d au d xi
    这将最终等于期望的结果(参见)。

    但是,有一些使用简单概率论的简便方法。回想一下,两个独立随机变量的总和给出了一个随机变量,其密度等于两个总和随机变量的卷积。如果AN(μA,σA2)A sim mathcal{N}left(mu_{A}, sigma_{A}^{2} ight)BN(μB,σB2)B sim mathcal{N}left(mu_{B}, sigma_{B}^{2} ight)是独立的,那么C=A+BC=A+BAN(μA+μB,σA2+σB2)A sim mathcal{N}left(mu_{A}+mu_{B}, sigma_{A}^{2}+sigma_{B}^{2} ight)多元推广也是如此。

    请注意,如果X1N(0,σ12I2)X_{1} sim mathcal{N}left(0, sigma_{1}^{2} I_{2} ight), X2N(0,σ22I2)X_{2} sim mathcal{N}left(0, sigma_{2}^{2} I_{2} ight),那么它们有密度函数分别是G1G_1G2G_2。因此总和Z=X1+X2Z=X_{1}+X_{2}具有G=G1G2G=G_1*G_2给出的密度函数。

    但是我们知道ZN(0+0,σ12I2+σ22I2)=N(0,(σ12+σ22)I2)Z sim mathcal{N}left(0+0, sigma_{1}^{2} I_{2}+sigma_{2}^{2} I_{2} ight)=mathcal{N}left(0,left(sigma_{1}^{2}+sigma_{2}^{2} ight) I_{2} ight),因此,ZZ的密度函数为:
    pZ(z)=14π2Σexp(12(z0)TΣ1(z0))=12π(σ12+σ22)exp(12zTz[σ12+σ22])=12π(σ12+σ22)exp(x2+y22[σ12+σ22])=G(x,yσ12+σ22)=G(x,yσ)egin{aligned} p_{Z}(z) &=frac{1}{sqrt{4 pi^{2}|Sigma|}} exp left(-frac{1}{2}(z-0)^{T} Sigma^{-1}(z-0) ight) \ &=frac{1}{2 pileft(sigma_{1}^{2}+sigma_{2}^{2} ight)} exp left(-frac{1}{2} frac{z^{T} z}{left[sigma_{1}^{2}+sigma_{2}^{2} ight]} ight) \ &=frac{1}{2 pileft(sigma_{1}^{2}+sigma_{2}^{2} ight)} exp left(-frac{x^{2}+y^{2}}{2left[sigma_{1}^{2}+sigma_{2}^{2} ight]} ight) \ &=G(x, y | sqrt{sigma_{1}^{2}+sigma_{2}^{2}}) \ &=G(x, y | sigma) end{aligned}

    这里z=(x,y)z=(x, y) and Σ=(σ12+σ22)I2=(σ12+σ22)2|Sigma|=left|left(sigma_{1}^{2}+sigma_{2}^{2} ight) I_{2} ight|=left(sigma_{1}^{2}+sigma_{2}^{2} ight)^{2}

    所以:σ=σ12+σ22sigma=sqrt{sigma_1^2+sigma_2^2}

    [参考内容]
    【1】内容1
    【2】内容2
    【3】内容3
    【4】内容4
    【5】内容5
    【6】内容6

  • 相关阅读:
    记录一次线上优化流程
    php ignite 使用问题记录
    invalid contrller specified 错误分析及解决
    koa 2 的 async 和 await 语法
    koa 2 的安装
    vue 自定义组件 v-model双向绑定、 父子组件同步通信的多种写法
    VS2019专业版和企业版激活密钥
    RE:ゼロから始める PKU 生活 episode 2
    CSP-S 2020 游记
    ioi2021集训队作业
  • 原文地址:https://www.cnblogs.com/Jack-Tim-TYJ/p/12831893.html
Copyright © 2020-2023  润新知