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

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

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

而且我还发现,在某些代码中,人们会考虑将两个连续的高斯模糊 σ1\sigma_1σ2\sigma_2用只是一个模糊σ=σ12+σ22\sigma=\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}\right)^{-1} \exp \left(-\frac{x^{2}+y^{2}}{2 \sigma^{2}}\right)

所以高斯模糊两次就相当于卷积两次:IB=G1(G2I)=(G1G2)I=GII_{B}=G_{1} *\left(G_{2} * I\right)=\left(G_{1} * G_{2}\right) * 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}\right) * G_{2}\left(x, y | \sigma_{2}\right)

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

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

请注意,如果X1N(0,σ12I2)X_{1} \sim \mathcal{N}\left(0, \sigma_{1}^{2} I_{2}\right), X2N(0,σ22I2)X_{2} \sim \mathcal{N}\left(0, \sigma_{2}^{2} I_{2}\right),那么它们有密度函数分别是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}\right)=\mathcal{N}\left(0,\left(\sigma_{1}^{2}+\sigma_{2}^{2}\right) I_{2}\right),因此,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σ)\begin{aligned} p_{Z}(z) &=\frac{1}{\sqrt{4 \pi^{2}|\Sigma|}} \exp \left(-\frac{1}{2}(z-0)^{T} \Sigma^{-1}(z-0)\right) \\ &=\frac{1}{2 \pi\left(\sigma_{1}^{2}+\sigma_{2}^{2}\right)} \exp \left(-\frac{1}{2} \frac{z^{T} z}{\left[\sigma_{1}^{2}+\sigma_{2}^{2}\right]}\right) \\ &=\frac{1}{2 \pi\left(\sigma_{1}^{2}+\sigma_{2}^{2}\right)} \exp \left(-\frac{x^{2}+y^{2}}{2\left[\sigma_{1}^{2}+\sigma_{2}^{2}\right]}\right) \\ &=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}\right) I_{2}\right|=\left(\sigma_{1}^{2}+\sigma_{2}^{2}\right)^{2}

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

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

posted @ 2020-05-02 22:07  田纳尔多  阅读(114)  评论(0编辑  收藏  举报