サンプリング周波数を調整する方法

Question

こんにちは。
bitalinoをpythonから動かしてリアルタイム処理することを目標にいろいろやっているのですが、
出力されたcsvファイルを見ていると、サンプリング周波数が1000になっているのにも関わらずデータ数が1000個しか出力されてきません。
改善策を教えていただきたいです。
また、このプログラムでとれた心拍からピークを取って出力するというところまでをリアルタイム処理でしていきたいと考えています。
乱文になりましたが、ご教授お願いいたします。

import matplotlib.pyplot as plt

import time

from bitalino import BITalino

# The macAddress variable on Windows can be "XX:XX:XX:XX:XX:XX" or "COMX"
# while on Mac OS can be "/dev/tty.BITalino-XX-XX-DevB" for devices ending with the last 4 digits of the MAC address or "/dev/tty.BITalino-DevB" for the remaining
macAddress = "98:D3:A1:FD:5B:E8"

#This example will collect data for 5 sec.
running_time = 10
    
batteryThreshold = 30
acqChannels = [0, 1, 2, 3, 4, 5]
samplingRate = 1000
nSamples = 10
digitalOutput = [1,1]

# Connect to BITalino
device = BITalino(macAddress)

# Set battery threshold
device.battery(batteryThreshold)


# Read BITalino version
print(device.version())
    
# Start Acquisition
device.start(samplingRate, acqChannels)

start = time.time()
end = time.time()
f = open('hogehoge2.csv','w',)

while (end - start) < running_time:
    # Read samples
    data = device.read(nSamples)
    s = data[:, 7]

    for i in range(0,9):
        print(s[i])

    f.write(str(s[i]))
    f.write('\n')
    #f.close()
    end = time.time()
    plt.plot(data[:, 7])
f.close()

# Turn BITalino led on
device.trigger(digitalOutput)
    
# Stop acquisition
device.stop()
    
# Close connection
device.close()

Comments

[y_waiwai]

問題は何でそれをどうしたいんでしょうか
問題が見えません

[N.R]

筋電を計測したいのでサンプリング周波数を1,000に設定したいのです。
また、ECGまで計測したいからです。

[y_waiwai]

で、それがどうなってるんですか？
問題がなにかわかりません

[N.R]

わかりにくくてすいません。
csvファイルで筋電の値が出力されてくるのですが、サンプリング周波数を1000にして10秒間計測してもデータ数が1000個しか出力されてこないのです。samplingRateを100にしようが1000にしようが、100のまま固定されている。というのが問題です。

[thkana]

そのセンサーは秒あたりどれくらいのデータを取得できるのですか?
秒10個のデータしか取れないセンサーなら、1000個/秒とろうとしても無理、ということになりますが。

[N.R]

返信ありがとうございます。

Bitalino向けのOpenSignalsというソフトウェアではサンプリング周波数を10〜1000まで調整することができているので、センサー自体は大丈夫だと考えています。
よろしくお願いします。

[ozwk]

nSamples を100とか1000にするとか
plt.plot(data[:, 7])をとりあえず消してみるとどうなります?

あとcsvの出力は1000個ピッタリですか?

samplingRateを10にしたらどうなります?

[N.R]

お返事遅れて申し訳ありません。
plt.plot(data[:, 7])を消しました。

nSampleを100・samplingRateを1000にするとcsvの出力は100個でした。
nSampleを1・samplingRateを1000にすると
print(s[i])
IndexError: index 1 is out of bounds for axis 0 with size 1
というエラーがでてきます。
csvファイルは何回か計測してみましたが、997~998個出力されます。

samplingRateを10にするとcsvファイルは10個出力されてきました。
よろしくお願いします。

[YT0014]

f.write(str(s[i]))
f.write('\n')
上記のファイル書込みですが、ご提示いただいたコードでは、forのループの範囲外となっているようですが、実コードも同様でしょうか？