Python：10msごとのサンプリングデータを1msごとのサンプリングデータに作り変えるには？

20msごとのサンプリングデータを1msごとのサンプリングデータに作り変えたい

疑問

20msごとのデータを1msごとのデータに変形する方法がわからない

なぜそれが必要なのか？

もともと1ms毎に取得したデータと10msごとに取得したデータで相互相関をとりたい。
そのために10ms毎のデータを1msに変形させる必要性があると思ったからであります

行動規範の内容に同意します

回答4件

ベストアンサー

どない頑張っても相互相関には1msの解像度はないので、1msのデータを(移動平均などで)10ms間隔にダウンサンプリングしてなまらせるほうがいいような。

投稿2018/09/10 00:48

tachikoma

総合スコア3601

Kalman_taku

2018/09/10 00:58

tacikomaさん、ありがとうございます！無知で大変申し訳ないのですが、相互相関を使うには、何msまでが限度なのでしょうか？

tachikoma

2018/09/10 02:47 編集

粗いほうのサンプリング間隔が限度です。例えば1msと10msの組み合わせの場合だと10msです。これは、10msのほうには1msの変動は含まれません（サンプリングできません）が、一方で1msのほうには10msの粗い周期の信号が含まれているためです。

Kalman_taku

2018/09/10 01:07

なるほどですね。つまり、10msごとのデータをアップサンプリングする方法はあるが、それをしてしまうと、そのアップサンプリング時に追加した値は1ms毎のデータには含まれていないため、10msを1msにアップサンプリングするより、 1msを10msにダウンサンプリングしたほうが良いということでしょうか？

tachikoma

2018/09/10 01:17

そうです。心配なのは、アップサンプリングによるアーチファクトです。それから、ダウンサンプリングするときは、適当なフィルタでスムージングしてからダウンサンプリングしないと、今度はエイリアシングが出る可能性があります。

行動規範の内容に同意します

時系列データなら、ダウンサンプリング前のフィルタ設計が結構大変（移動平均だと周波数特性が平坦ではありません）だと思いますし、フィルタによる位相遅れも発生しますので、位相遅れの無いフィルタ込みのscipy.signal.decimate()を使うのがおすすめです。

https://docs.scipy.org/doc/scipy/reference/generated/scipy.signal.decimate.html

ただ、間引き後のデータにおいても高周波帯域の成分はフィルタである程度減衰しますので、ダウンサンプリングレートはよく検討してください。確認するには、例えば、白色雑音を生成してそれをdecimate()で処理した信号をFFT解析するなどの方法があります。元も白色雑音は平坦なスペクトルを有するはずですが、処理後は高周波帯域が減衰していることが観察できると思います。

投稿2018/09/11 04:46

編集2018/09/11 04:48