如何建立一个合理的Hilbert谱情节与R（计算环境）

Question

我使用的 EMD 包R.该包具有用于显示一个希尔伯特谱（与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函数的源代码，以限定具有用于调色板参数自己的绘图功能。