物理学者(ポスドク)による日々の研究生活のメモ書きです ( python/emacs/html/Japascript/シェルスクリプト/TeX/Mac/C言語/Linux/git/tmux/R/ポケモンGO)
なんとなくやってみたくなって調べていたら、良さそうな方法が見つかったので・・・
こんなをgnuplotを使って作ります
まずはデータの準備。
そちらはいろんなところに良い記事がありますので、そちらを参考にどうぞ(なんて不親切な記事なんだ・・・)
リサジュー図形
データさえ生成できれば、あとは以下のようなスクリプトで1枚目の20枚のリザージュ曲線をまとめたようなグラフがパパっと描けます
次はgifアニメの作り方
以前にも記事は書きましたが、今回新たに別の方法を知りました・・・
・以前の記事
【Gnuplot】のみを用いてgifアニメを作る
【gnuplot】gifアニメを作ってみた
・今回見つけた方法
gnuplotで1つのデータファイルからgifアニメを
要するに、グラフを書くためのデータに例えば100行毎に改行があれば
それを index 1 のようにしてその部分だけをグラフ化することができると・・・
(using の親戚らしい)
で、データをそのように作ってあげると・・・
(完全に自分用のメモだなこれは・・・)
ツイート
こんなをgnuplotを使って作ります
まずはデータの準備。
そちらはいろんなところに良い記事がありますので、そちらを参考にどうぞ(なんて不親切な記事なんだ・・・)
リサジュー図形
データさえ生成できれば、あとは以下のようなスクリプトで1枚目の20枚のリザージュ曲線をまとめたようなグラフがパパっと描けます
#!/bin/sh
for((phase=0; phase<=180; phase+=45))
do
f1=10; f2=10;
gene_data ${f1} ${f2} ${phase} > ${f1}_${f2}_${phase}.txt
f1=10; f2=20;
gene_data ${f1} ${f2} ${phase} > ${f1}_${f2}_${phase}.txt
f1=10; f2=30;
gene_data ${f1} ${f2} ${phase} > ${f1}_${f2}_${phase}.txt
f1=20; f2=30;
gene_data ${f1} ${f2} ${phase} > ${f1}_${f2}_${phase}.txt
#echo ${phase}
done
TERMINAL="set terminal png"
#XLABEL="set xlabel\"x\""
#YLABEL="set ylabel\"y\""
OUTPUT=reserge.png
OPTION="w l lw 2 notitle"
gnuplot <<EOF
set terminal png size 1080, 720 fontscale 1.0
${XLABEL}
${YLABEL}
set size square
set output "${OUTPUT}"
set multiplot layout 4, 5
plot "10_10_0.txt" ${OPTION}
plot "10_10_45.txt" ${OPTION}
plot "10_10_90.txt" ${OPTION}
plot "10_10_135.txt" ${OPTION}
plot "10_10_180.txt" ${OPTION}
plot "10_20_0.txt" ${OPTION}
plot "10_20_45.txt" ${OPTION}
plot "10_20_90.txt" ${OPTION}
plot "10_20_135.txt" ${OPTION}
plot "10_20_180.txt" ${OPTION}
plot "10_30_0.txt" ${OPTION}
plot "10_30_45.txt" ${OPTION}
plot "10_30_90.txt" ${OPTION}
plot "10_30_135.txt" ${OPTION}
plot "10_30_180.txt" ${OPTION}
plot "20_30_0.txt" ${OPTION}
plot "20_30_45.txt" ${OPTION}
plot "20_30_90.txt" ${OPTION}
plot "20_30_135.txt" ${OPTION}
plot "20_30_180.txt" ${OPTION}
unset multiplot
EOF
rm ./*.txt
for((phase=0; phase<=180; phase+=45))
do
f1=10; f2=10;
gene_data ${f1} ${f2} ${phase} > ${f1}_${f2}_${phase}.txt
f1=10; f2=20;
gene_data ${f1} ${f2} ${phase} > ${f1}_${f2}_${phase}.txt
f1=10; f2=30;
gene_data ${f1} ${f2} ${phase} > ${f1}_${f2}_${phase}.txt
f1=20; f2=30;
gene_data ${f1} ${f2} ${phase} > ${f1}_${f2}_${phase}.txt
#echo ${phase}
done
TERMINAL="set terminal png"
#XLABEL="set xlabel\"x\""
#YLABEL="set ylabel\"y\""
OUTPUT=reserge.png
OPTION="w l lw 2 notitle"
gnuplot <<EOF
set terminal png size 1080, 720 fontscale 1.0
${XLABEL}
${YLABEL}
set size square
set output "${OUTPUT}"
set multiplot layout 4, 5
plot "10_10_0.txt" ${OPTION}
plot "10_10_45.txt" ${OPTION}
plot "10_10_90.txt" ${OPTION}
plot "10_10_135.txt" ${OPTION}
plot "10_10_180.txt" ${OPTION}
plot "10_20_0.txt" ${OPTION}
plot "10_20_45.txt" ${OPTION}
plot "10_20_90.txt" ${OPTION}
plot "10_20_135.txt" ${OPTION}
plot "10_20_180.txt" ${OPTION}
plot "10_30_0.txt" ${OPTION}
plot "10_30_45.txt" ${OPTION}
plot "10_30_90.txt" ${OPTION}
plot "10_30_135.txt" ${OPTION}
plot "10_30_180.txt" ${OPTION}
plot "20_30_0.txt" ${OPTION}
plot "20_30_45.txt" ${OPTION}
plot "20_30_90.txt" ${OPTION}
plot "20_30_135.txt" ${OPTION}
plot "20_30_180.txt" ${OPTION}
unset multiplot
EOF
rm ./*.txt
次はgifアニメの作り方
以前にも記事は書きましたが、今回新たに別の方法を知りました・・・
・以前の記事
【Gnuplot】のみを用いてgifアニメを作る
【gnuplot】gifアニメを作ってみた
・今回見つけた方法
gnuplotで1つのデータファイルからgifアニメを
要するに、グラフを書くためのデータに例えば100行毎に改行があれば
それを index 1 のようにしてその部分だけをグラフ化することができると・・・
(using の親戚らしい)
で、データをそのように作ってあげると・・・
set terminal gif animate nocrop font "arial,12" fontscale 2.0 animate delay 15 loop 0 size 1080,720
#set terminal gif animate
set output "hoge.gif"
nkmax=10
nk=0
while (nk < nkmax) {
plot [-1:1][-1:1] "${OUTPUT}" index nk w p pt 5 ps 1
nk=nk+1
}
のようにしてgifアニメを作ることができます#set terminal gif animate
set output "hoge.gif"
nkmax=10
nk=0
while (nk < nkmax) {
plot [-1:1][-1:1] "${OUTPUT}" index nk w p pt 5 ps 1
nk=nk+1
}
(完全に自分用のメモだなこれは・・・)
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)
(03/05)
(02/29)
(02/21)
(02/21)
(02/21)
(02/21)
(01/13)
(01/05)
(01/05)