物理学者(ポスドク)による日々の研究生活のメモ書きです ( python/emacs/html/Japascript/シェルスクリプト/TeX/Mac/C言語/Linux/git/tmux/R/ポケモンGO)
こんな感じのテキストデータがあって、これを読み込む
■ 過去記事 : Python を用いて実験データを読み込み、可視化する
これで1つのテキストファイルをndarrayに収納はできた
あとは、複数のファイルを読み込んで1つにまとめる
listにしといて、appendを使ったほうが早いのかもしれんけど、とりあえず動いたので良し
(2022/08/7 追記)
このhstackで2つのarrayを結合している
https://deepage.net/features/numpy-stack.html#基本的な結合の仕方
(2023/02/01 追記)
上の例だと、ファイルの通し番号順に上からarrayを作ってくれない
1つ目に開いたファイルの中身がarrayの一番最初に来て欲しい時は・・・
1 2
3 4
5 6
みたいな2列のデータを読み込みたい
(2023/06/09 追記)
ある任意の列だけ読み込みたいときは、
np.loadtxt(fname, dtype=str, unpack=Truc, usecols=[0, 2, 3])
のように書く
これで、0, 2, 3列目のみ取り出せる
loadtxtでテキストファイルを読み込むときは、すべての列数が同じことが必要みたいで
通常は5個しか要素がないのに、ある行だけ6個要素があるとエラーになる
必要な列だけ書いて読み込むことで、そのエラーを回避することができる
ツイート
1 21
2 22
3 23
4 24
5 25
2列目のみ解析に使う、2列目のみ使うので1列目は捨てる2 22
3 23
4 24
5 25
import numpy as np
fname="hoge.txt"
a = (np.loadtxt(fname, unpack=True))[1]
unpackのオプションがなんで必要かはいかの記事を参照fname="hoge.txt"
a = (np.loadtxt(fname, unpack=True))[1]
■ 過去記事 : Python を用いて実験データを読み込み、可視化する
これで1つのテキストファイルをndarrayに収納はできた
あとは、複数のファイルを読み込んで1つにまとめる
import numpy as np
for i in range(100):
fname="hoge%d.txt" % (i)
a = (np.loadtxt(fname, unpack=True))[1]
if i == 0:
b = a
else:
b = np.hstack((b, a))
これでbにすべての結果が入るfor i in range(100):
fname="hoge%d.txt" % (i)
a = (np.loadtxt(fname, unpack=True))[1]
if i == 0:
b = a
else:
b = np.hstack((b, a))
listにしといて、appendを使ったほうが早いのかもしれんけど、とりあえず動いたので良し
(2022/08/7 追記)
このhstackで2つのarrayを結合している
https://deepage.net/features/numpy-stack.html#基本的な結合の仕方
(2023/02/01 追記)
上の例だと、ファイルの通し番号順に上からarrayを作ってくれない
1つ目に開いたファイルの中身がarrayの一番最初に来て欲しい時は・・・
1 2
3 4
5 6
みたいな2列のデータを読み込みたい
for i in range(7, 21):
fname = "hoge%d.txt" % i
# ファイルサイズが0以上かチェック
if os.path.getsize(fname) > 0:
aa, bb = np.loadtxt(fname, unpack=True)
if i == 7:
a = aa
b = bb
else:
a = np.hstack((a, aa))
b = np.hstack((b, bb))
else:
print("error at %s" % fname)
fname = "hoge%d.txt" % i
# ファイルサイズが0以上かチェック
if os.path.getsize(fname) > 0:
aa, bb = np.loadtxt(fname, unpack=True)
if i == 7:
a = aa
b = bb
else:
a = np.hstack((a, aa))
b = np.hstack((b, bb))
else:
print("error at %s" % fname)
(2023/06/09 追記)
ある任意の列だけ読み込みたいときは、
np.loadtxt(fname, dtype=str, unpack=Truc, usecols=[0, 2, 3])
のように書く
これで、0, 2, 3列目のみ取り出せる
loadtxtでテキストファイルを読み込むときは、すべての列数が同じことが必要みたいで
通常は5個しか要素がないのに、ある行だけ6個要素があるとエラーになる
必要な列だけ書いて読み込むことで、そのエラーを回避することができる
PR
この記事にコメントする
プロフィール
HN:coffee
職業:物理屋(自称)
趣味:映画鑑賞、登山
出身:大阪府の南の田舎
自己紹介:
import MyProfile
import coffee_pote from TWITTER
import amazonのほしい物リスト from WISH_LIST
print "先月子供が産まれました!"
# 最終更新 2022/10/25
職業:物理屋
趣味:映画鑑賞、登山
出身:大阪府の南の田舎
自己紹介:
import MyProfile
import coffee_pote from TWITTER
import amazonのほしい物リスト from WISH_LIST
print "先月子供が産まれました!"
# 最終更新 2022/10/25
カテゴリー
ブログ内検索
リンク
相互リンク募集中です
(Twitterにてお知らせください)
Demo scripts for gnuplot version 5 (gnuplotのさまざまなデモ画像と作り方がまとめられている、眺めているだけでできるようになった気分になれる)
gnuplotスクリプトの解説 (米澤進吾さんの個人ページ、gnuplotと言えばこのかた)
gnuplot のページ (Takeno Lab、うちのブログがリンクされていたのでリンク返し)
Twitterから映画の評価が分かる & 映画の鑑賞記録が残せる coco (映画の感想をまとめられるサイト、いつもお世話になっています)
Astronomy Picture of the Day Archive (天文や宇宙関連の最新の話題について画像とともにNASAが説明しているページ)
今日のほしぞら (任意の時刻の空で見える星を表示してくれる、国立天文台が管理している)
GNUPLOTとアニメーション (応用の項目の「見せてあげよう!ラピュタの雷を!!」あたりからすごすぎる)
読書メーター (読んだ本をリストできる便利なサイト)
flickr難民の写真置き場 (20XX年、flickrは有料化の炎に包まれた。あらゆるflickr無料ユーザーは絶滅したかに見えた。 しかし、tumblr移住民は死に絶えてはいなかった。)
教授でもできるMac OS X へのLaTeX, X11, gccのインストレーションと環境設定 (阪大の山中卓さんのwebページ、タイトルにセンスが溢れている、内容は超充実してる、特にTeX関連、学振DCとかPDの申請書類作成時にはお世話になっております)
英語論文執筆用の例文検索サービス (とんでもないものを見つけてしまった・・・・ arXivに収録されている 811,761報の 英語論文から,例文を検索するための検索エンジン)
Demo scripts for gnuplot version 5 (gnuplotのさまざまなデモ画像と作り方がまとめられている、眺めているだけでできるようになった気分になれる)
gnuplotスクリプトの解説 (米澤進吾さんの個人ページ、gnuplotと言えばこのかた)
gnuplot のページ (Takeno Lab、うちのブログがリンクされていたのでリンク返し)
Twitterから映画の評価が分かる & 映画の鑑賞記録が残せる coco (映画の感想をまとめられるサイト、いつもお世話になっています)
Astronomy Picture of the Day Archive (天文や宇宙関連の最新の話題について画像とともにNASAが説明しているページ)
今日のほしぞら (任意の時刻の空で見える星を表示してくれる、国立天文台が管理している)
GNUPLOTとアニメーション (応用の項目の「見せてあげよう!ラピュタの雷を!!」あたりからすごすぎる)
読書メーター (読んだ本をリストできる便利なサイト)
flickr難民の写真置き場 (20XX年、flickrは有料化の炎に包まれた。あらゆるflickr無料ユーザーは絶滅したかに見えた。 しかし、tumblr移住民は死に絶えてはいなかった。)
教授でもできるMac OS X へのLaTeX, X11, gccのインストレーションと環境設定 (阪大の山中卓さんのwebページ、タイトルにセンスが溢れている、内容は超充実してる、特にTeX関連、学振DCとかPDの申請書類作成時にはお世話になっております)
英語論文執筆用の例文検索サービス (とんでもないものを見つけてしまった・・・・ arXivに収録されている 811,761報の 英語論文から,例文を検索するための検索エンジン)
最新記事
(11/20)
(05/17)
(05/17)
(05/09)
(03/05)
(02/29)
(02/21)
(02/21)
(02/21)
(02/21)