我有每小时风速数据和每小时表面波高数据,并希望看看这些数据是否相关以及两者之间是否存在时滞。该数据目前在两个独立的熊猫数据框架中。我尝试过使用scipy.correlate和np.xcorr,但它们只返回零或数据大小作为结果。任何帮助将不胜感激。
nsamples = wind.speed.size
wind.speed -= wind.speed.mean(); wind.speed /= wind.speed.std()
wave.HEIGHT -= wave.HEIGHT.mean(); wave.HEIGHT /= wave.HEIGHT.std()
xcorr = correlate(wind.speed, wave.HEIGHT)
dt = np.arange(1-nsamples, nsamples)
recovered_time_shift = xcorr.argmax()
print dt
print "Recovered time shift: %d" % (recovered_time_shift)
非常感谢,