r - R 中的图像平滑

Question

如何在 R 中平滑这张图片，以便只保留两个峰值？

如果这是一维数据，我会对其进行运行均值或拟合回归函数。但我没有找到有关在二维矩阵上应用这些方法的非常具体的信息。例如，我尝试使用 stats 包中的 filter() 。

我也考虑过克里金法，但这更多的是插值，是吗？

Answer 1

spatstat 包包含一个应用高斯模糊的函数 blur()。这在某种程度上使图像变得平滑，大部分噪声消失并且两个主峰清晰可辨。

效果如下图所示，相当显着，尤其是在3d图中。

生成图片的代码是:

library(jpeg)
library(fields)
library(spatstat)

picture <- readJPEG("~/Downloads/spectrogram.png.jpeg")
picture2 <- as.matrix(blur(as.im(picture), sigma=6))

layout(matrix(c(1:4), nrow=2))
image.plot(picture, col=gray.colors(50), main="original image", asp=1)
image.plot(picture2, col=gray.colors(50), main="blurred with sigma = 6", asp=1)
drape.plot(1:nrow(picture), 1:ncol(picture), picture, border=NA, theta=0, phi=45, main="original spectrogram")
drape.plot(1:nrow(picture), 1:ncol(picture), picture2, border=NA, theta=0, phi=45, main="blurred with sigma = 6")