質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.48%
Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

Q&A

解決済

1回答

451閲覧

データの抽出をループ処理したい

noob4

総合スコア21

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

0グッド

0クリップ

投稿2022/06/04 04:20

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')

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

自己解決

自己解決しました。

python

1import glob 2import os.path 3 4data = np.loadtxt(('conf2.txt'), comments='#', unpack=True) 5xrange1 = data[0] 6xrange2 = data[1] 7 8path = "data/" # use your path 9all_files = glob.glob(path + "*.dat") 10 11path2 = './out' 12if not os.path.exists(path2):#ディレクトリがなかったら 13 os.mkdir(path2)#作成したいフォルダ名を作成 14 15for k in range (len(all_files)): 16 twotheta, intensity = np.loadtxt((all_files[k]), comments='#', unpack=True) 17 file_path = all_files[k] 18 19 x = twotheta 20 y = intensity 21 xcnt = (np.where(x < xrange1)) 22 xcnt = (np.max(xcnt)) 23 x = x[xcnt:] 24 y = y[xcnt:] 25 xcnt2 = (np.where(x > xrange2)) 26 xcnt2 = (np.min(xcnt2)) 27 x = x[:xcnt2] 28 y = y[:xcnt2] 29 30 x_y = np.stack([x,y],axis = 1) 31 32 name = f'{k}' +f"{os.path.basename(file_path)}" 33 np.savetxt(os.path.join(path2, name + ".dat"),x_y,fmt='%.9f')

投稿2022/06/04 04:53

noob4

総合スコア21

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.48%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問