是否有任何看起来(模糊)像这样的过滤器?

信息处理 过滤器
2022-02-09 14:37:11

我正在尝试确定特定非线性过程对我的数据的影响。从模拟中,我可以测量“理想”和“处理”图像的功率谱并获取它们的比率,(我认为)应该对应于某种过滤器或过滤过程。

那么,是否有任何过滤器或过滤过程看起来(模糊地)像所附的图?这两个图是线性线性(左)和对数对数(右)空间中的相同数据。

过滤函数图

更多细节:我有一个获取时间流数据的管道,其中时间流上的每个点都与图像中的一个点相关联。在创建图像之前对时间流执行一些“清理”(PCA 减法)。我的目标是确定时间流处理对结果图像的影响,最终目标是在没有可比时间流数据的图像上复制效果。数据来自测辐射热计阵列,这些是毫米波长的天空图像。

这是模拟图像。左为加工前,右为加工后。边缘周围的尖点是没有时间流数据的地方。

文字描述的图像

另请参阅此 ipython 笔记本http://nbviewer.ipython.org/4393835/进行一些分析。

由于时间流过滤过程完全有可能无法由图像空间中的任何过滤器表示,因此解决“我如何近似这些过滤器函数”的答案是合适的。

1个回答

我将使用有限脉冲设计来近似滤波器,其中 FIR 滤波器的 wrights 是根据滤波器通过的频谱计算的。您可以通过选择所需频率响应的样本(从上图中)并进行样本的离散傅里叶逆变换来做到这一点。结果将是所需滤波器的脉冲响应。然后将滤波器实现为使用脉冲响应样本值作为权重的加权和。所以你的过滤器将是A0+A1Z1+A2Z2... 在哪里A0是脉冲响应的第一个值,A1是第二个等等。 Zn表示 n 的样本延迟(来自滤波器结构的 Z 域表示)。这种方法有时被称为 FIR 滤波器设计的“频率采样”方法。我已经使用这种方法来创建包含特定光谱内容的时间序列。Google FIR 滤波器设计,频率采样。