pythonでwavファイルをフィルタ処理した後に書き出しができない.

Question

前提・実現したいこと

wavファイルを読み込んで、scipy signal decimate()を用いたLPF込みの間引きを行い、wavファイルを書き出したいと考えているのですが、以下のようなエラーが出ており書き出すことができません。

発生している問題・エラーメッセージ

 File "bandpass.py", line 109, in <module>
    write(wav_filename_out,rate_d,data_d)
  File "/Applications/anaconda3/lib/python3.8/site-packages/scipy/io/wavfile.py", line 788, in write
    fmt_chunk_data = struct.pack('<HHIIHH', format_tag, channels, fs,
struct.error: required argument is not an integer

該当のソースコード

import numpy as np
import matplotlib.pyplot as plt
import sys
import scipy.io.wavfile
from scipy.io.wavfile import write
from scipy.io.wavfile import read
from scipy.fftpack import rfft,irfft,fftfreq
from scipy import signal

##音声ファイル読み込み##
args = sys.argv
wav_filename = args[1]
rate, data = scipy.io.wavfile.read(wav_filename)
time = np.arange(0,data.shape[0]/rate,1/rate) #時間軸作成

print(len(time))
print(len(data))
print("data_MAX:",np.max(data))
print("data_min:",np.min(data))


#dataを高速フーリエ変換
fft_data = np.fft.fft(data)
freqList = np.fft.fftfreq(data.shape[0], 1.0/rate) #周波数取得

fft_data = fft_data/22050 #サンプル数/2

######################      scipy.signal.decimate      #############################
data_d = signal.decimate(data, 10)  # 間引き
rate_d = rate/10 #間引き後のサンプリング周波数
time_d = np.arange(0,data_d.shape[0]/rate_d,1/rate_d) #間引き後の時間軸

print("data_d:",len(data_d))
print("time_d:",len(time_d))
print("data_d_MAX:",np.max(data_d))
print("data_d_min:",np.min(data_d))

#dataを高速フーリエ変換
fft_data_d = np.fft.fft(data_d)
freqList_d = np.fft.fftfreq(data_d.shape[0], 1.0/rate_d) #周波数取得

fft_data_d = fft_data_d/22050 #サンプル数/2

######################################################################################



##音声ファイル書き込み##
data_d = data_d.real.astype(np.int16)
args = sys.argv
wav_filename_out = args[2]
data_d.flags.writeable = True
write(wav_filename_out,rate_d,data_d)

試したこと

今回使用しているwavファイルは16bitのサンプリングレートが44.1kHzのものを使用しています。

補足情報（FW/ツールのバージョンなど）

python3

実行時：python3 file名　入力ファイル　出力ファイル

Comments

[melian]

最終行の write(wav_filename_out,rate_d,data_d) の第２引数である rate_d は整数でないと上記のエラーが発生することになります。実際の rate_d の値はどうなっていますか？

[iface]

回答ありがとうございます。
data_d_MAX: 21400.095955183147
data_d_min: -21494.69185424528
でしたので整数にすると無音ではありますが書き出すことができました。

Answer 1

rate_d = rate/10 #間引き後のサンプリング周波数

を以下に変更しましょう。

rate_d = rate//10 #間引き後のサンプリング周波数

Comments

[iface]

変換したところ書き出すことができました。ありがとうございます。
しかし、なぜか無音の状態で書き出されるみたいです。

[iface]

コード上の書き出しのところで
data_d = data_d.real.astype(np.int32)
になっていたことが原因でした。