4

我正在编写具有可变半径（标准偏差）的高斯模糊，即图像的每个像素都使用不同的内核进行卷积。计算高斯模糊的标准技术在这里不起作用：FFT、轴分离、重复框模糊——它们都假设内核对于整个图像是相同的。

现在，我正在尝试使用以下方案来近似它：

使用由一组 N 组轴对齐矩形 R k和系数 α k定义的分段常数函数 f(x,y) 逼近高斯核 K(x,y) ：

    f(x,y) = ∑<sub>k=1 N α k ·χ R k (x,y)

让 g(x,y) 成为我们的图像，然后

    ∬<sub>ℝ<sup>2 K(x,y)·g(x,y) dxdy ≈ ∬<sub>ℝ<sup>2 f(x,y)·g(x,y) dxdy = ∑< sub>k=1 N α k ·∬<sub>R k g(x,y) dxdy

RHS 上的积分是矩形上的简单积分，因此可以通过预先计算整个图像的部分和来在恒定时间内计算。

生成的算法在 O(W·H·N) 中运行，其中 W 和 H 是图像的尺寸，N 与近似误差成反比 (AFAIK)。

剩下的部分是找到一个好的逼近函数 f(x,y)。当给定矩形数量 N（最小化误差）或给定误差（最小化矩形数量）时，如何找到高斯的最佳近似值？

4

1 回答 1

0

给定矩形的位置和大小，计算出系数应该相当容易，所以真正的问题是计算出矩形的放置位置。

由于您正在逼近高斯，因此将我们的注意力限制在其中心与高斯中心重合的矩形上似乎至少是合理的，因此我们实际上只有一个一维问题 - 计算出一组嵌套矩形的大小如果你的纵横比不是统一的，我认为要么是正方形，要么类似于高斯。

这可以通过动态规划来解决。假设你从外到中工作。在第 N 阶段，您已经制定了一个 nxk 表，该表为您提供来自 1,2...N 个外部像素环的最佳近似误差，用于最多 1,2,..k 个不同的矩形，以及最内层矩形的大小对那个最好的错误负责。要计算阶段 N+1，您需要考虑到目前为止最里面的矩形的所有可能大小，为外部区域贡献 x 个像素环。您计算出该矩形的 alpha，它最适合新环中的像素，并且其外部的环不留给外部矩形。使用已经计算出的表中的值，您知道当您留下多达 k 个外部矩形来覆盖这些区域时，您将得到的最佳可能错误，因此您可以计算出现在 N+1 像素环所贡献的最佳总误差。这允许您为 N+1 个外部像素填写表格条目。当您进入该区域的中间时，您将能够为整个区域制定最佳解决方案。

于 2011-09-25T08:14:53.573 回答