xyデータを任意のxの範囲で抽出し、新たなxyデータとして保存するコードを作成しました。
このコードをループ処理して、フォルダ内のすべてのxyデータに対して同様の処理を行い、
5つのファイルを読み込んだ場合、新たに5つのファイルを出力するようなコードを書きたいです。
この場合ループ内はどのようなコードを書けばよいのでしょうか。
python勉強中のため、説明だけでなく
サンプルコードを書いていただければ幸いです。
python
1import numpy as np 2import glob 3import os.path 4 5#####config ファイル読み込み 6data = np.loadtxt(('conf2.txt'), comments='#', unpack=True) 7xrange1 = data[0] 8xrange2 = data[1] 9 10path = "data/" #dataフォルダを指定 11all_files = glob.glob(path + "*.dat") #data読み込み 12 13path2 = './out'#書き出しフォルダの指定 14if not os.path.exists(path2):#ディレクトリがなかったら 15 os.mkdir(path2)#作成したいフォルダ名を作成 16 17twotheta, intensity = np.loadtxt((all_files[0]), comments='#', unpack=True) 18 19x = twotheta 20y = intensity 21 22xcnt = (np.where(x < xrange1)) 23xcnt = (np.max(xcnt)) 24x = x[xcnt:] 25y = y[xcnt:] 26xcnt2 = (np.where(x > xrange2)) 27xcnt2 = (np.min(xcnt2)) 28x = x[:xcnt2] 29y = y[:xcnt2] 30 31x_y = np.stack([x,y],axis = 1) 32 33np.savetxt(os.path.join(path2, "x_y" + ".dat"),x_y,fmt='%.9f')
回答1件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。