前提・実現したいこと
現在、csvファイルから読み込んだデータをフーリエ変換し、Excelに出力するプログラムを作成しています。
csvファイルの読み込みの際にはエラーが発生していませんでしたが、Excelファイルに結果を出力するコードを書き加えた際にエラーが発生しました。
出力するExcelファイルのA列に周波数fq[:int(N/2)+1]の値を、B列からは、各列にData[:int(N/2)+1]の値を入れたいと思っています。
発生している問題・エラーメッセージ
UnicodeDecodeError: 'utf-8' codec can't decode byte 0x83 in position 0: invalid start byte
該当のソースコード
Python
1import numpy as np 2import pandas as pd 3import os 4import matplotlib.pyplot as plt 5import openpyxl 6 7file_name = input("Please input filename: ") 8 9Sampling_frequency = int(input("Please input Sampling_frequency: ")) 10 11Number_of_iterations = int(input("Please input Number_of_iterations: ")) 12 13color_list=["black","gray","lightcoral","brown","darkred","red","tomato","sienna","sandybrown","darkorange","tan","gold","darkkhaki" 14 ,"yellowgreen","olivedrab","darkturquoise","lightblue" 15 ,"steelblue","cornflowerblue","navy","blue","indigo","darkviolet","purple","magenta","crimson","lightpink"] 16 17data=pd.read_csv(os.path.join(r"Z:/private/結果/4月/",file_name),header=None,usecols=list(range(1,n+1)),names=[f"{n}回目" for n in range(1,n+1)]) 18 19book=openpyxl.load_workbook(r"Z:/private/結果/4月/検証.xlsx",encoding="shift_jis") 20sheet=book["0%"] 21#データ数 22N=len(data) 23#サンプリング周波数 24fs=Sampling_frequency 25#窓関数 26window=np.hamming(N) 27#窓関数の補正値 28acf=(1/(sum(window)/N)) 29#時間軸のデータ作成 30t=np.arange(0,N*(1/fs),1/fs) 31#周波数軸のデータ作成 32fq=np.linspace(0,fs,N) #周波数軸 linspace(開始,終了,分割数) 33 34fig = plt.figure(figsize=(6,4)) 35 36for n in range(1,Number_of_iterations+1): 37 38 #2次元配列から1次元配列に変更 39 data2=data[str(n)+"回目"].values 40 #窓関数後の信号 41 data_w=window*data2 42 #読み込んだデータをフーリエ変換 43 data_w_FFT=np.fft.fft(data_w) 44 #FFTの複素数結果を絶対値に変換 45 data_w_FFT_abs=abs(data_w_FFT) 46 #窓補正 47 Data=(data_w_FFT_abs)*acf 48 plt.plot(fq[:int(N/2)+1],Data[:int(N/2)+1],color=color_list[n],label=str(n)+"回目") 49 df=pd.DataFrame(Data[:int(N/2)+1],fq[:int(N/2)+1],columns=[str(n)+"回目"]) 50 print(df) 51 52 sheet["A1"]=fq[:int(N/2)+1] 53 sheet["B1"]=Data[:int(N/2)+1] 54plt.xlabel("frequency[Hz]",fontsize=14) 55plt.ylabel("signal intensity",fontsize=14) 56plt.xlim(0,2000) 57ax = fig.add_subplot(111) 58ax.legend(prop={"family":"MS Gothic"}) 59plt.show() 60book.save(r"Z:/private/結果/4月/検証.xlsx",encoding="shift_jis")
試したこと
他の文字コードを試してみましたが、同様のエラーが見られました。
補足情報(FW/ツールのバージョンなど)
openpyxl 3.0.0
となっていました。
エラー箇所は、
data=pd.read_csv(os.path.join(r"Z:/private/結果/4月/",file_name),header=None,usecols=list(range(1,n+1)),names=[f"{n}回目" for n in range(1,n+1)])
book=openpyxl.load_workbook(r"Z:/private/結果/4月/検証.xlsx",encoding="shift_jis")
この2文どちらかが考えられます。
回答1件
あなたの回答
tips
プレビュー