r - 如何使用 R（计算环境）创建合理的希尔伯特谱图

Question

我正在使用 R 的EMD包。这个包有一个频谱图功能，用于显示希尔伯特频谱（用 hilbertspec 计算）。然而，输出实际上是模糊的和黑白的。

此功能似乎没有输出彩色图像的选项。我怎样才能清楚地显示光谱，如果可能的话以彩色显示。

Answer 1

灰度级在spectrogram函数 ( gray()) 中进行了硬编码，要覆盖此设置，您可以使用以下内容：

# define a color palette
colors <- colorRampPalette(c("#007FFF", "blue", "#000077"))  
gray <- function(x) colors(255*x)  # redefine gray palette
spectrogram(test1$amplitude[,1], test1$instantfreq[,1])
gray <- grDevices::gray  # reset gray palette function

另一种选择是使用spectrogram函数的源来定义您自己的绘图函数，该函数具有调色板的参数。