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

Question

我正在使用 R 的 EMD 包。该包有用于显示希尔伯特谱（用 hilbertspec 计算）的谱图函数。然而，输出确实很模糊并且是黑白的。

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

Answer 1

灰度级在频谱图函数 (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

另一个选项是使用源的频谱图函数来定义您自己的绘图函数，该函数具有调色板的参数。