波形データの加工（Z変換・高速フーリエ・ウエーブレット変換など）

fMRIのデータの加工を行った。

というわけで、処理用のPythonスクリプトのメモ

1.Z変換

手書きでもいいけど、それは論文にするときに。

import scipy.stats as SS

x = [1,2,3,4,5]
zx = SS.zs(x)
print zx

2.高速フーリエ変換

実数型の高速フーリエ。

import scipy.fftpack as FF

x = [適当な配列]
fx = FF.rfft(x)
print fx

3.離散型ウェーブレット変換

離散型ウェーブレット変換というもの。

まだ、利用方法は勉強中だけど、一応使い方だけ。

import pywt
"""
pywtがなかったら
sudo apt-get install python-pywt
でインストール 
"""
x = [適当な配列]
(hw, lw) = pywt.dwt(x, "好きなWavelet")
print hw
print lw

フーリエ変換とかZ変換とかは、使い方が分かるからそのまま解析に使ってみて

Waveletはいろいろ実際の使い勝手を勉強しよう。