python-如何使用scipy.signal.welch和年度数据生成功率谱

Question

我想使用scipy.signal.welch生成功率谱密度图。但是，我的数据是光谱辐射率与天/年的关系，如果我的x轴以天为单位，我不确定如何选择采样率。

我的光谱辐射度图看起来像这样，我试图通过进行快速傅立叶变换（蓝色图）来显示趋势（不是很成功）：

以下代码的构想来自这个问题here

在我的情况下，采样率是多少？由于数据点是在16天的固定时间间隔内获取的，请问这是我的采样率吗？此外，如果选择了一个窗口，即使文档指出以下内容，仍然有必要将nperseg定义为256：

  

nperseg：int，每个段的可选长度。默认为无，但是   如果window是str或tuple，则设置为256，如果window是array_like，   设置为窗口的长度。

from scipy.signal import welch

sig = data_plot_multi_year1["y"]

f, psd = welch(sig,
fs= ????  ,           # sample rate
window='hanning',     # apply a Hanning window before taking the DFT
nperseg= 256,         # compute periodograms
detrend='constant')   # detrend x by subtracting the mean