LoginGNUPLOTで光度曲線を描く−AIP4Winのlogファイルからグラフ化する−AIP4Winを使って測光しlogファイルに記録した場合の、データファイル化と光度曲線のグラフ化をWindowsPCで行う方法のメモです。1つは、Excelを使いテキストファイルを読み込んで処理する方法で、Excelに慣れた人には方法は簡単ですが、同じ処理を何度も繰り返すといい加減嫌になります。グラフも大きさが一定にならず、グラフタイトルや目盛を毎回毎回付けるのうんざりしています。(VBAマクロを組めばよいのでしょうが、それに慣れるまで大仕事です)そこで、AWK(テキスト処理)とGNUPLOT(グラフ処理)の出番です。この方法はコマンドプロンプトからコマンドを打ち込むのですが、テキストファイルに書いておいたものをコピー・ペーストすれば良いので手間は要りません。PCにインストールします。その具体的な方法は、ここでは省略しますが、検索して見てください、世の中には日本語の丁寧な解説が溢れています。データファイルの整理あるディレクトリ(=フォルダー)例えば、"20130118"にAIP4WinのLogファイルが"LT_GemR.txt"というファイル名で置いてあるとします。長いですが実物を示します。1890.116st8.00000254.RED.FIT<以下省略>実際には、1晩観測すると測光データが延々と続きます。何行あっても以下の処理は、行数に無関係に全く同じように処理出来ます。このディレクトリ内にテキストエディタで、次のAWK用の1行スクリプトを保存して置き、コピーしてコマンドプロンプトにマウス右クリックでペーストします。(最近は私はテキストエディタには、サクラエディター(オープンソース)を使っています。>LT_GemR.datこの1行スクリプト(1行野郎とかワンライナーと呼ぶようです)の意味は次のとおりです。最初の38行分と空行を無視し、次のブロックを抜き出し出力する。第4フィールド(Julian_Day)から56200.0を引いたもの、第7から第10フィールド(V-C1、sigma、C2-C1、sigma)。これらのデータを"LT_GemR.txt"というファイルから抜き出し(計算し)、"LT_GemR.dat"というファイルに書き出す。コマンドプロンプトに、マウス右クリックでペーストすると次のようになり、Enterキーを押します。>LT_GemR.datするとそのディレクトリに次の様な内容のテキストファイル"LT_GemR.dat"が書きだされます。0.0047<以下省略>以上で、測光データファイルが整理できました。この夜の光度曲線を描く上記の作業で、この夜はIc、Rc、V、Bの4バンドで測光し、次のファイル名で保存してあるとします。"LT_GemIc.dat","LT_GemR.dat","LT_GemV.dat","LT_GemB.dat"GNUPLOTのスクリプトファイルを準備するデータファイルが置いてあるのと同じディレクトリに、GNUPLOTの作図のためのpltファイルを作ります。次のような内容のテキストファイルを、例えば"LT_GemBVRI.plt"というファイル名で保存しておきます。点の形、点の大きさ、点の色を表す。GNUPLOTを動かすコマンドプロンプトから次のように打ち込みます。まとめ光度曲線を描くのに行った作業は次の3種類だけです。"LT_Gem20130107IRVBC1.png"何日もの観測データを整理する場合、Excelで行う作業量に比べると圧倒的に少ない手数で、しかも統一のとれたフォーマットのグラフが得られます。