前回の報告 (11/09 2013)以後、 現在の CVS 版に入れられた主な機能について紹介します (ChangeLog, manual の更新部分等より)。
- 一般的な機能
- 新しい描画スタイル with parallelaxes の追加 (12/26)
- 新しい smooth オプション mcsplines の追加 (12/24)
- plot/splot の notitle にタイトル文字列をつけた場合の問題の修正 (12/23)
- set/unset に関するメモリリークの問題の修正 (12/21)
- 文字列値のループ変数に対する問題の修正 (12/21)
- 軸に関するソースの改良 (12/19)
- 「set history」の新設 (12/16)
- "history !N" で N 番目のコマンドを再実行できるように (12/14)
- key のフォントが splot で効いていなかった問題の修正 (12/11)
- unset tics が機能していなかった問題の修正 (11/22)
- plot title の後ろに個々に {no}enhanced をつけられるように (11/14)
- terminal ドライバ
- emf: フォント名の初期化に関する問題の修正 (12/12)
- tek.trm: 新しい出力形式 "sixel" の追加 (12/12)
- wxt: 拡張文字列のセンタリングに関する問題の修正 (12/12)
- qt: 構造体の初期化の問題の修正 (12/11)
- qt: Qt アプリケーションへの埋め込み用の改良 (12/11)
- emf: pattern fill の実装 (12/11)
- emf: 矢軸や point (5, 6) の線幅が太くならない問題の修正 (12/11)
- pnt: linecolor だけが指定されたときは linetype は変更しないように (11/09)
- qt: マウスボタン操作のタイミング処理の再改良 (11/09)
- その他
- volatile データに関する条件コンパイルを削除 (12/23)
- VMS 用の修正 (tek.trm も含む) (12/22)
- automake 1.14 用の修正 (12/21)
- bitmap.c のライセンスの改正 (の許諾) (12/15)
- libgd の configure に関する修正 (12/12)
- UTF-8 文字列を含むソースが SJIS 環境でコンパイルできない問題の修正 (12/11)
- qt terminal を Qt4 か Qt5 のいずれかで作成できるように (11/14)
- defined() をドキュメントからは削除し本体には一応まだ残す (11/14)
- デモ
- parallel.dem: with parallelaxes のデモの追加 (12/26)
また報告が 2 ヶ月ほど開きましたが、今回もいくつかの更新が入っています。
出力形式として、tek.trm に sixel 出力形式というものが追加されました。 これは、DEC の古い端末でサポートされていたグラフィック形式らしいですが、 よくは知りません。
また、emf 出力形式に対して、 「情報やメモ (11/28 2013)」 で報告したいくつかの問題に対する修正が追加されました。 ただし、fill pattern に対しては、 それが作る EMF ファイルの fill pattern は LibreOffice ではうまく表示されない、という問題があるようです (もしかすると LibreOffice 側の問題 ?)。
また、今回は history に関する修正もだいぶ行われて、 シェルの history と同様に、N 番目のコマンドを再実行したり、 「set historysize」が廃止されて、 それを拡張した「set history」が新設されています。
また、plot title に関して、 個別に {no}enhanced をつけれるようにしたこと、 notitle と title が完全に同じオプションを許すようにしたこと、 などの修正もありました。 この後者は、plot ... title "hoge" と書いたものを、 「title」の部分だけ「notitle」としてもちゃんと動くように、 ということのようです。
mcsplines という smooth オプションが追加されましたが、 これは「区分的単調 3 次スプライン」(piecewise monotonic cubic splines) というもので、元のデータ点の単調性、凸性を保持するもののようです。
そして、新しい描画スタイルとして、 「with parallelaxes」というものが追加されました。 これは parallel axes plot、parallel coordinates plot などと呼ばれるもののようで、 日本語では「平行座標描画」などと言われるようです。 統計分野で利用されるもののようですが、 これは x, y の 2 軸 (または x, y, z の 3 軸) を垂直に置くのではなく、 多次元のデータ軸を平面に平行に並べ、 データの組 (データファイルの 1 行のデータ) を各軸に取って それを隣り合う軸毎に線分で結んだものです。 データ行分の独立した折れ線ができます。 サンプルは、例えば以下を参照してください。
- 本家で公開している CVS 版の parallel.dem のサンプル
http://gnuplot.sourceforge.net/demo_cvs/parallel.html
個人的にはあまり見たことがありませんが、 3 次元以上のデータの相関を見るにはよさそうですね。 ただ、なんとなくですが、 データ軸の並べ方などを試行錯誤する必要がありそうな気がします。
(cf. 「情報やメモ (01/29 2014)」, 「情報やメモ (07/22 2014)」, 「情報やメモ (07/22 2014)」, 「情報やメモ (05/03 2016)」)
Gnuplot Q&A 掲示板 (2673 番の記事) で、 emf terminal に関する問題点がいくつか指摘されました。 それらについて報告します。
もともとは、2 ちゃんねるの掲示板「gnuplot を使おう。その 3」 (アドレスは リンクリスト 参照) で報告されていたものをまとめたもののようです。 報告されている問題点はいくつかあり、まとめると、
- pointtype 6 の円 (中を塗らない) の線が lw で太くならない (2 ちゃんねる 200, 203, 204, 269 番、Gnuplot Q&A 掲示板 2673 番)
- 塗り潰し矢先を持つ arrow の軸が lw で太くならない (2 ちゃんねる 269 番、Gnuplot Q&A 掲示板 2673 番)
- set multiplot title と set output を繰替えすと、 2 回目の multiplot title の文字が出ない (2 ちゃんねる 269 番、Gnuplot Q&A 掲示板 2676 番)
- fill pattern がパターンでなく solid の濃さでの表現になっている (2 ちゃんねる 261, 268, 269 番)
のような感じです。 まず、1 つ目の pointtype 6 の円と 2 つ目の arrow の軸ですが、 以下のようなスクリプトで問題が起こります。
set term emf
set output 'temp.emf'
set arrow from -8,14 rto 5,-5 head size 5,20 empty lt 1 lw 6
set arrow from -8, 9 rto 5,-5 head size 5,20 filled lt 1 lw 6
plot x w p lt 1 lw 6 pt 4 ps 2 , x+3 w p lt 1 lw 6 pt 6 ps 2
- バージョン 4.6.3 での上のサンプルの出力
(EMF 画像の gif でのスクリーンショット)
これらは、いくつかの箇所で lw を見ずに 固定した線幅を設定していたことによる問題でした。 pointtype 6 だけでなく 7 の円 (7 は塗り潰しの円) も影響を受けます。 これらは多分以下の差分で修正できます。
- emf-linewidth-1.diff (1275 Byte)
- このパッチを当てた CVS 版の出力 (EMF 画像の gif でのスクリーンショット)
3 つ目の set multiplot title の問題は、 以下のようなスクリプトで問題が起こります。
set term emf
set output 'temp1.emf'
set multiplot title 'hoge'
plot x
unset multiplot
set output 'temp2.emf'
set multiplot title 'hoge'
plot x
unset multiplot
- バージョン 4.6.3 での上のサンプルの 2 つ目の EMF ファイルの出力
(EMF 画像の gif でのスクリーンショット)
- emf-setfont-1.diff (1701 Byte)
- このパッチを当てた CVS 版の出力 (EMF 画像の gif でのスクリーンショット)
4 つ目の pattern fill の問題は、test コマンドで確認できます。
- バージョン 4.6.3 での test コマンドの出力
(EMF 画像の gif でのスクリーンショット)
EMF ではパターンでの塗り潰しが少し面倒くさいこともあり (情報も少ない)、 これまで実装されてこなかったためのものです。 今までは色の濃さでごまかしていたようですが、 他の出力形式のパターンもだいたい同じ形に揃いましたので、 EMF もそれに合わせるべきだろうと思います。 以下の差分でそれを実装してみました。 全部一から作ったわけではなく、 例えばビットパターンなどは win terminal (src/win/wgraph.c) から拝借しています。
- emf-patfill-1.diff (3520 Byte)
- このパッチを当てた CVS 版の出力 (EMF 画像の gif でのスクリーンショット)
なお、現在の CVS 版は、xtics のデフォルトの形式の変更 (%h) に伴い、 × (積) を表す Unicode 文字列が一つ src/util.c に入っていて、 それが Shift_JIS 環境の場合に問題を起こしてしまいますので、 Shift_JIS 環境の VC++ などでコンパイルする場合には、 以下のパッチが必要になると思います。
Gnuplot Q&A 掲示板で情報を紹介してくださった ryo-naka さん、 どうもありがとうございました。 一応いずれのパッチも、すでに本家の方に送ってありますが、 採用されるかどうかはまだわかりません。
(cf. 「情報やメモ (12/27 2013)」, 「情報やメモ (07/18 2014)」, 「情報やメモ (07/22 2014)」, 「情報やメモ (05/03 2016)」)
前回の報告 (09/09 2013)以後、 現在の CVS 版に入れられた主な機能について紹介します (ChangeLog, manual の更新部分等より)。
- 一般的な機能
- save コマンドで軸の fix, extend 等が保存されなかった問題の修正 (11/06)
- デフォルトに戻すための "set xrange [*:*] extend" (11/06)
- 軸名に関する処理の改良 (11/06)
- 3D データの 4 列目が色の場合の Inf, NaN のデータ処理に関する修正 (10/29)
- 等高線の色の割り当てに関するバグの修正 (前回の改良も無効に) (10/27)
- key の箱を拡張文字列も含まれるように拡大 (10/25)
- タブ区切りの CSV ファイルの最初のフィールドが空の時の処理の修正 (10/25)
- 新しいエラーハンドラルーチン (gp_exit) の導入 (10/25)
- 16bit DOS 用のコードを削除 (10/25)
- 出力形式が初期化される前にマウスからイベントを送らないよう (10/23)
- "gnuplot -" 時でエラー後のセッションの継続の問題 (10/23)
- "test palette" コマンドに関する古い問題の修正 (10/21)
- パレット色の割り当てルーチンの改良 (10/21)
- object の属性として "clip/noclip" を新設 (10/18)
- "show object N" の表示が正しくなかった問題の修正 (10/17)
- さまざまな object のクリッピングルーチンの改良 (10/17)
- key title のリセット時のメモリリークの問題の修正 (10/13)
- "set pm3d top/bottom" と " set pm3d interpolate N,M" を同時に使えるように (10/13)
- 格子線に関するデフォルト値の一部変更 (10/13)
- 対数軸で、副目盛/格子線が主目盛/格子線に乗ることの抑制 (10/13)
- "smooth" で不要な列指定があったときに警告を出すように (10/13)
- 軸の刻みのデフォルトの書式用の書式 %h を新設 (10/06)
- 未定義変数をマクロで展開しないように (10/02)
- クリッピングに関するソースの修正 (10/02)
- with labels で文字列と点の記号それぞれで variable color を使えるように (10/01)
- error bar のクリッピング処理に関する問題の修正 (09/30)
- 軸範囲の内部構造の単純化 (09/26)
- 軸ラベルの書式文字列を動的割り当てに変更 (09/26)
- 色の内部構造 (RGB 色か linetype か) を RGB 色に統一 (09/26)
- RGB 色の設定ルーチンの改良 (09/25)
- RGB 色をサポートしてない出力形式でも等高線を set color で描くように (09/23)
- set color をサポートしない出力形式用の空ルーチンを提供 (09/23)
- RGB color をサポートしない出力形式で不透明な塗り潰しをしないよう修正 (09/20, 09/23)
- rtics の初期化に関する問題の修正 (09/20)
- CVS 版で histogram のデフォルトが clustered にならない問題の修正 (09/20)
- 極座標モードの自動縮尺の際の rmax の評価を修正 (09/17)
- 文字列評価の際に致命的エラーが起きる問題の修正 (09/14)
- call の引数評価の際にそれを文字列と見るように変更 (09/14)
- reset session で NaN と pi も再初期化し、 ユーザ定義関数をすべて破棄するように (09/12)
- set style increment と reset コマンドに関する問題の修正 (09/11)
- xticlabels(COL) で目盛が出ないなどの問題に関する修正 (09/09)
- tc variable がデフォルトで不正に色のバーを出していた問題の修正 (09/09)
- set style increment user でない時に lc variable、 with labels tc variable が lt に合わない問題の修正 (09/09)
- terminal ドライバ
- qt: linecolor だけが指定されたときは linetype は変更しないように (11/09)
- qt: マウスボタン操作の応答のタイミング処理の改良 (要 Qt5) (10/27)
- qt: persist 時にホットキー 'i' を処理できるように (10/27)
- qt: データを初期化の際にだけ生成し終了ハンドラで破棄するように (10/26)
- unixpc: 削除 (10/25)
- gd, cairo: ドライバがエラーを標準出力に送らない、exit() を呼ばないよう (10/23)
- wxt: パイプ入力がマウスによって一部消えてしまう問題の修正 (10/11)
- enhanced モードをサポートする出力形式ではそれをデフォルトに (10/06)
- emf: 円の塗り潰し (pt 7) が黒境界で描かれていた問題の修正 (09/17)
- cairo: メモリリークに関する修正 (09/09)
- その他
- aquaterm 用の configure 等の修正 (09/20)
- デモ
少し前回の報告から間が開きましたので、 色々変更が入っていますが、あまり大きなものはないと思います。 ただ、見た目に影響を与える変更もなくはありませんので、 従来と出力が変わった、という場合もあるとは思います。
例えば、今回軸の刻みのデフォルトの書式が、 従来の %g から、新設された %h に変更されました。 %h は、%g が指数表記をする場合には「e+XX」のように表記するのに代わり、 x10^{XX} のように表記するもののようです。 確かにそちらの方が良く使われるような気もしますし、 実際そのように表現されるように自前で設定されていた場合も多いと思います。
クリッピング処理も、ルーチンの改良、ルールの変更などの他に、 オブジェクトに対するクリッピングオプションも新設されています。
その他には、極座標モードの自動縮尺の際の rmax の評価や、 RGB 色に関する変更、lt, lc 回りの修正なども入っています。 enhanced モードをサポートする出力形式では、 それがデフォルトになったのも少し便利かもしれません。
また、今回 16bit DOS 時代の古いコードが削除され、それとともに unixpc という M68010 ベースの AT&T 3b1, AT&T 7300 Unix PC 用の古いドライバも破棄されました。 Wikipedia を見ると、1980 年代位のもののようで、 Microsoft もまだ MS-DOS 3.X 位の時代です。 画面出力もかなり原始的な時代のもののようですが、 コアなファンはかなりいるらしく インターネットにはまだかなり情報があります。 そんなコンピュータで gnuplot がどのようなグラフを出力していたのか 若干興味はありますが、さすがにそういう情報は見あたりません。 ドライバのソースを見ると、 720x300 ビットマップを仮定していて、 線種も 4 種類をサポートしているようですので、 もしかすると意外に綺麗なグラフ (dumb よりはずっとましなもの) を書いていたのかもしれません。
(cf. 「情報やメモ (12/27 2013)」)
gnuplot-4.6.4 がリリースされています。 配布物に含まれる NEWS には、以下のように書かれています。
- 新規: wgnuplot で fit の当てはめを Ctrl-Break で中断するように
- 変更: CVS ファイルの空フィールドを「不正」ではなく「欠損」と扱うように
- 変更: 'using' や 'title' で複数の先頭列 (column header) を参照できるように
- 変更: デフォルトの "make install" では install-info は実行しないように
- 変更: 極座標軸が自動縮尺の場合、なるべく原点を中心に置くように
- 修正: svg と canvas 出力形式のマウス動作が座標軸と反対になっていた
- 修正: システムによって emf のフォントの初期化に失敗していた
- 修正: 日時データ列を column(N), column("HEAD") で参照できるように
- 修正: qt 出力形式の plot with labels の拡張文字列要素のトグル動作
- 修正: columnstacked histograms で生成される key の色/パターン
- 修正: ^C の 2 度打ちで X サーバによる wxt セッションのハングアップを解消
- 修正: win 出力形式でウィンドウのリサイズ後に描画境界の調整の問題
- 修正: ある条件のもとではコマンド入力時にマクロが展開させない問題
- 修正: 数のみからなる文字列を複素数ではなく整数に展開するように
- 修正: 'set grid front' が軸の 0 点の初期配置に失敗していた問題
- 修正: x11 の persist モードでかなり低い精度のマウス座標が表示されていた問題
- 修正: $# ("call" の引数の個数) の解析。これはコメントではない!
- 修正: pngcairo 出力形式で画像の切り出しのメモリリークの問題
- 修正: "lc variable" は linestyle ではなく linetype の色を繰り返すように
- 修正: rtics がたまに 0 の長さで書かれていた問題
マイナーバージョンアップではだいたいそうですが、 見てわかる通り、大きな変更はありません。 CVS 版で既に色々取り入れられている機能はありますが、 それらは次のメジャーバージョンアップ (4.8.0) で採用されるのだろうと思います。
前回の報告 (08/11 2013)以後、 現在の CVS 版に入れられた主な機能について紹介します (ChangeLog, manual の更新部分等より)。
- 一般的な機能
- gnuplot -c script ARG1 ARG2 ... で load でなく call が行えるように (09/08)
- call 回りのコードの修正、$# の処理に関する問題の修正 (09/08)
- gnuplot が定義する変数が削除されないようにグループ化 (09/07)
- パレット関数に関するメモリリークの修正 (09/07)
- reset の新しいオプション session (リセット後初期化ファイルを実行) (09/07)
- reset でデフォルトのマウスの設定を復帰し、軸オプションは受けつけない (09/01)
- object の保存の前に unset object を入れるように修正 (08/28)
- set style histogram, newhistgram でタイトルフォントと文字色を指定できるよう (08/28)
- load, pause <time> などの際にもマウスイベントを処理できるように修正 (08/23, 08/26)
- test コマンドの原点の設定 (set origin) に関する修正 (08/23)
- fit の縮尺されたパラメータ誤差のχ自乗値に関する修正 (08/20)
- r 軸の小目盛 (minor tics) のサポート (08/20)
- r 軸に関する GPVAL_ 値が欠けていた問題の修正 (08/13)
- terminal ドライバ
- x11, xlib: ファイルからの入力時のマウスイベント処理への対応 (08/25)
- x11: -persist モード時のマウス座標の精度を上げる修正 (08/24)
- vgagl, ggi: ^C とマウス動作の更新に関する修正 (08/23)
- qt, wxt, x11: multiplot 等に対するマウスのレイヤー選択に関する修正 (08/19, 08/20)
- x11, win: Layer に関する修正 (08/19)
- win: WIN64 (x86_64) への対応 (08/16)
- その他
- make 時のタイムスタンプの仕組みを改善 (08/12)
- デモ
- games/nibbles.dem, games/xixit.dem: ゲームデモの追加 (08/31)
- histgrams2.dem: タイトル、ラベル位置などの修正 (08/29)
多分一番目につくのはゲームデモでしょう。 いずれも矢印キーなどをホットキーとして使うもので、 nibbles は snake (古くからの Unix ユーザならわかるでしょうか) の変種、 xixit はテトリスの変種のようです。 このようなことが gnuplot でできてしまいますが、 それに合わせてマウス処理に関する修正が行われています。
新設オプションの「reset session」も、 通常の reset よりも便利な場合が多いと思います。
もう一つの大きな変更が gnuplot の新しいコマンドラインオプションである -c です。 gnuplot のコマンドラインでスクリプトとして指定すると、 それは通常 load コマンドで読み込んだ状態になるのですが、 それを load でなく call に変えて、 その後のオプションを gnuplot ではなく call したスクリプトへのオプション とみなすものです。 これにより、今までは間接的にしか呼べなかった call 用のスクリプトを 直接 gnuplot から呼べるようになりますし、 gnuplot スクリプトの先頭に
#!/usr/local/gnuplot -cのようなものを入れて実行権限を与えることで、 オプションを使用する gnuplot スクリプト (call スクリプト) を 直接実行することができるようになります (ただし、私はまだ試していません)。
なお、この「gnuplot -c」の機能は、山本さんの 「gnuplot の精義」(第 2 版) の 6.4 節に書かれていることの示唆から実現したものです。 実は、最近の gnuplot の更新には、 他にもいくつかこの本に関連した修正、追加などがあります。 山本さん、どうもありがとうございました。
(cf. 「情報やメモ (11/09 2013)」)
2 ちゃんねるの掲示板「gnuplot を使おう。その 3」(アドレスは リンクリスト 参照) に書かれた情報について、多少気になるものがありましたので、 少し書いておきます。
-
FreeBSD のパッケージから gnuplot+ がなくなった、 という記事がありました (152)。
確かに最新の ports リストからはなくなっているようです。 私が使っている FreeBSD-8.2 にはまだありますから、 少し前の ports を使って入れる、という手はあるかもしれません。
http://www.freshports.org/math/gnuplot+/ を見ると、今年の 2 月頃になくなっているようです。削除理由は
This package is obsolete. The author does not recommend to use it.
とありますが、 確かに 3.7 系列の gnuplot に戻るよりも 今の 4.6 系列の gnuplot を使う方が賢明だと思います。 マニュアルによれば gnuplot+ ができる機能は以下の通りですが、 現在の gnuplot でも多くが対応できています:- postscript terminal の plus 機能
(文字列部分に LaTeX ライクな命令が使える)
現在は epslatex, cairolatex, lua, mp terminal などで、 グラフと文字部分を分離した出力を生成し、 文字部分は LaTeX に渡すような仕組みがありますから、 ほぼそれで代用できるでしょう。
または psfrag package を使って、Postscript 内の文字を LaTeX の文字に LaTeX 側で置きかえる、という手もあります。
- 日本語対応 (postscript, x11, pbm, png および epson, hp 等のプリンタ)
現在は、postscript, x11, png (gd ライブラリ) では問題なく日本語が使えます。 ただし、pbm や epson, hp のプリンタに関する日本語化は 多分対応していません。
- GNU readline によるコマンド履歴の保存
既に組込まれています (configure 時に選択)
- LANG の値や set locale によるヘルプファイルの切り替え
その機能はありませんが、環境変数 GNUHELP に ヘルプファイルを指定できますので、 gnuplot 実行中に切り替えるのは無理ですが、 起動時にヘルプファイルを切り替えることなら可能です。
- postscript terminal の plus 機能
(文字列部分に LaTeX ライクな命令が使える)
-
qt terminal はどういいのか、という記事がありました (155)。
多分、その前の 153 の記事のサンプルスクリプトで、 対話型 terminal への復帰に set term qt としてあったための質問だと思います。
もちろん、153 の記事のサンプルは、「set term qt」でなくて 「set term wxt」でも「set term win」でも「set term x11」 でもいいわけですが、 qt terminal のいいところは、と言われると微妙ですね。 クロスプラットフォーム、という点では wxt とそう変わりません。 qt terminal が wxt terminal よりも優位な点、というのも思いつきません。 グラフウィンドウのメニューが日本語化されていることくらいでしょうか。
当然ですが、KDE 環境のように Qt がすでに使える環境にあり、 wxWidgets を使うにはインストールしないといけない、 という場合には優位だろうとは思います。
-
最近の gnuplot では set term png で直接 png が出力できると知った、 といった記事がありました (156)。
この記事の意図とは多分違いますが、png terminal について少し書きます。 実は、これまで「png terminal」という形で利用できた (できる) terminal は以下の 3 種類のドライバのものがありました (あります)。
- png.trm (1995 年 11 月 (gnuplot-3.6beta) - 2004 年 8 月 (gnuplot-4.1)): libpng を使用する png terminal。
- gd.trm (2002 年 4 月より png をサポート (gnuplot-3.8) - 現在): gd library (libgd) を使用する png terminal。
- cairo.trm (2007 年 10 月より png をサポート (gnuplot-4.3) - 現在): cairo ライブラリを使用する pngcairo terminal。
最後の cairo 由来のものは「pngcairo terminal」という名前で、 正確には「png terminal」ではないのですが、 gd.trm が有効でなくて cairo.trm が有効な場合、 「set term png」で「pngcario terminal」が起動するようになっています。 だから、png terminal だと思って使っているものが、 実は pngcairo terminal である、という可能性はあります。 それは set term を実行した際に表示される「Terminal type set to ...」 の部分を見れば確認できます。
-
e-100 より小さい値をプロットすると 1 位の大きい値になってしまうのですが、 という質問がありました (159)。
CPU のバグでは、といった返信もありましたが、 これは実際に貼り付けてあったデータを見ると
0.000000000000E+00
となっていて、0 以外のデータの部分の 'E' が抜けていますので、 多分勘違いだと思います。 実際、- の前に E をおぎなうと正しく小さい値が表示されます。 E-201 くらいだとまだ倍精度の制限にはぶつかりませんから 特に問題なく表示されます。
0.000000000000E+00
0.000000000000E+00
0.435473392210-201
0.156895824551-183
.... -
4 列のデータから、中心の x 座標、y 座標、半径、「透過度」を読み取って 円を描きたい (with circles) という質問がありました (162)。
その記事によれば「透過度」というのは fillstyle の solid に指定する値のことで、「solid 0.15」を「solid variable」 に変えて列データから読み取りたい、ということのようです。
163 番の記事でも指摘されていますが、 その「透過度」(というよりも「塗り潰し密度」という方が多分適切) は variable のようにしてデータファイルから読み取ることはできません。 ソースを改良すれば可能かもしれませんが、ちょっと難しそうな気がします。
また、少し試してみたのですが、solid の後ろに 「0.9-0.8」のような数式を書いても正しく 0.1 と認識されるようですが、 「f(x)」のような関数を書くとエラーになるようです。 もし関数が使えるなら、plot for で 1 行ずつ書かせればなんとかなるかと思ったのですが それもだめなようです。
となると、gnuplot 内部からやることはあきらめて、 私なら awk などでやるでしょうか。 例えばデータ file.dat が
1.0 1.0 0.7 0.5
のように、中心の x, y 座標、半径、塗り潰し密度、となっているとして、 以下のような awk スクリプト mk1.awk を作って gnuplot スクリプトを吐かせます (もちろん awk 内部で gnuplot へのパイプを作って、 直接 gnuplot へ命令を渡すことも可能です):
1.0 2.0 0.7 0.8
2.0 2.0 0.7 1.0
2.0 1.0 0.7 0.2
これを使って、BEGIN {
print "set term wxt";
print "set size square";
print "set xrange [0:3]";
print "set yrange [0:3]";
}{
if (NR==1) printf "plot 'file.dat'";
else printf ",\\\n ''";
printf " using 1:2:3 every ::%d::%d", NR-1, NR-1;
printf " notitle with circles lc rgb 'blue'";
printf " fs transparent solid %s", $NF;
} END { print "" }
のようにして gnuplot スクリプト mk1.gp を作成し、% awk -f mk1.awk file1.dat > mk1.gp
のようにするわけです。% gnuplot -persist mk1.gp- 上のサンプルのスクリーンショット
シェルスクリプトのヒアドキュメントを使えば その中で全部の作業を行うことができますし、 上に書いたように awk スクリプト内部だけで行うこともできます。
- 上のサンプルのスクリーンショット
-
gnuplot コマンドを行末に \ を置いて複数行書いた場合、 コマンドヒストリも複数行に分かれて保存されてしまうが、 それを 1 行にできないか、という質問がありました (166)。
確かにそうなっている方がいいかなという気もします。 ただ、かなり長いコマンド行になってしまうと、 編集もやや面倒かもしれません。 簡単にするには、外部でフィルタで変換することでしょう。 多分 perl や ruby などでやると楽なんだと思いますが、 例えば awk でやるなら (ヒストリファイルを .gnuplot_history とします):
のようにすればいいわけです。 長いですが、シェルスクリプトにしておけば楽でしょう。 しかも、いっそそのシェルスクリプトを gnuplot の初期化ファイル (~/.gnuplot) 内で ! を使って gnuplot の起動時に実行しておくようにすれば 意識する必要がなくなります。% awk '{s=$0; while(s ~ /\\$/) { sub(/\\$/, "", s); getline; s=s " " $0 } print s}' .gnuplot_history > tmpf
% mv tmpf .gnuplot_historyただし、これは gnuplot のの再起動時に有効なだけで、 現在実行中の gnuplot では有効にはなりません。 そこまでやろうとすると、やはりソースを改変しないといけないでしょう。 有益な気もしますので、時間があったら少し考えてみたいと思います。
(cf. 「情報やメモ (07/18 2014)」)
前回の報告 (06/29 2013)以後、 現在の CVS 版に入れられた主な機能について紹介します (ChangeLog, manual の更新部分等より)。
- 一般的な機能
- arrow の矢先サイズ調節を行わない "fixed" オプションを新設 (08/09)
- arrow の境界線を描かない "noborder" オプションを新設 (08/08)
- 線分、点の端でのクリッピング処理に関する問題の修正 (08/08)
- 空の columnheader 文字列に関する問題の修正 (08/08)
- grid が off, front の際の zeroaxis の初期化に関する問題の修正 (08/08)
- columnhead(N) を using, title の部分文字列として使えるように (08/06)
- 数字 (と空白) のみの文字列は整数として変換するように (08/06)
- do ブロック内でマクロを利用した場合の問題の修正 (07/22)
- 複素誤差関数をいくつか実装 (libcerf を使用) (07/12, 07/14, 07/15, 07/18)
- Airy 関数の精度を上げる修正 (07/12)
- %s で秒データを読むときの小数以下の秒の処理の修正 (07/01, 07/03)
- 「print ${datablock}」で改行を出力するよう修正 (07/01)
- using でかっこが足りない場合明示的なエラーを表示するよう修正 (07/01)
- set [xy]label rotate parallel が逆向きの軸に起こす問題の修正 (06/29, 07/03)
- 2 以上のレベルの目盛り (tics) の新設 (06/29, 06/30)
- terminal ドライバ
- pslatex: 箱付きテキスト用マクロ定義の再定義に関する問題の修正 (07/25)
- svg: set fontpath の設定を使わないように修正 (07/18)
- win: 箱付きでない文字列のサイズに関するコードの修正 (07/06)
- win: 文字列選択が行の終わりを越えた場合の問題の修正 (07/06)
- win: Ctrl-C でクリップボードにコピーできるように (07/06)
- tgif: arrow を tgif の命令でなく gnuplot のルーチンを使用するように (07/06)
- tgif: 点のサイズを他の出力形式に合わせて半分の大きさに (07/05)
- epslatex, cairolatex: 箱付きラベルのサポート (07/05)
- その他
- wxWidgets に関する configure オプションの修正 (08/09)
- マクロに関する configure オプションの削除 (08/09)
- src/gadgets.c の RCSid が長い間狂っていたことの修正 (08/06)
- pgnuplot の削除 (07/06)
- configure の wxWidget 環境のチェックに関する修正 (07/06)
- emacs lisp ファイルのインストールに関する修正 (07/05)
- Mac 用の条件コード (#ifdef _Macintosh) の削除 (07/02)
- docs/doc2texi.el の terminal セクションのレベルの修正 (07/02)
- xemacs 用の docs/doc2texi.el 等の修正 (07/02)
- ドキュメントの修正 (binary matrix の x と y が実際には逆) (07/02)
- コンパイル時の警告対応のコードの修正 (07/01)
- アセンブラソースの削除 (07/01)
- デモ
- macros.dem: 画面表示部分の修正 (08/09)
- arrowstyle.dem: オプション noborder の追加に対応する修正 (08/08)
- cerf.dem: libcerf を使用した複素誤差関数用のデモ (07/14, 07/15)
今回の重要な変更は以下の通りです。
- arrow の矢先サイズ調節を行わない "fixed" オプションの新設
デフォルトでは、非常に短い arrow は、 矢先も小さくするようになっているそうですが、 size fixed を指定すると小さくしなくなります。
- arrow の境界線を描かない "noborder" オプションの新設
arrowstyle で filled を使用している場合用ですが、 これを指定すると、矢先の周囲の線を描かなくなります。 実は、2 点を結ぶ arrow は、この矢先の周囲の線の幅の分だけ、 矢先が終点から少しずれてしまうそうです。 この矢先の境界線を描かないようにすると、 矢先の先端が丁度終点に一致するようになります。
- columnhead(N) を using, title の部分文字列として使えるように
これにより、一つの plot title 内で、 複数の columnhead(N) を使用するようにできるようになるようです (例えば「plot 'file' using 1:2 title columnhead(1) . ' vs ' . columnhead(2)」のように)
- libcerf による複素誤差関数をいくつか実装
「複素誤差関数」などを実装する libcerf なるものがあるようです (cf. http://apps.jcns.fz-juelich.de/doku/sc/libcerf)。
libcerf は complex.h が利用可能 (C99) な環境用のライブラリですが、 無理矢理やれば少し古い環境 (例えば Solaris 9) でも インストールできなくはありません (cf. 「Unix に関するメモ等 (07/19 2013)」)。 ただし、その場合は gnuplot のコード (src/libcerf.c) も修正が必要です。
- 2 以上のレベルの目盛り (tics) の新設
従来、set xtics の level は 0 (major), 1 (minor) の 2 つしかありませんでしが、 これがより大きな数字レベルも指定できるようになりました。 現在使用可能なレベルは多分 5 くらいまでです。 2 以上のレベルは明示的な指定しかできませんが、 minor tics (level 1) とは違い、ラベルをつけることが可能です。
- tgif: 点のサイズを他の出力形式に合わせて半分の大きさに
tgif terminal の点のサイズは、 他の出力形式に比べてやや大き目だったようです (何で見つけたのか忘れました)。 他の出力形式に合わせてやや小さくしました (ついでにバグも修正してあります)。
- tgif: arrow を tgif の命令でなく gnuplot のルーチンを使用するように
これまでは tgif terminal では、 arrow は tgif の命令に直していましたが、 それだとグラフ上の矢の見た目が他の出力形式と変わってしまいますし、 矢の見た目や矢の形の調整機能は、「tgif の矢」よりも「gnuplot の矢」 の方がずっと上なので、 それを、線分の描画で実現する「gnuplot の矢」に変えました。
その代わり、生成される obj ファイルでは 矢が一つの命令で実現されるのではなく、 複数の折れ線 (多角形) で描かれることになってしまいましたので、 やや善し悪しですが、多分 tgif の矢の方がありがたい、 と思う方は少ないと思います。
また、機能の追加などとは違いますが、 マニュアルの「binary matrix」の部分に一つ修正が入りました。 これは実はかなり長く存在した問題でした。 それについて少し説明します。
これまで binary matrix のマニュアルには 以下のように書かれていました:
===========
よって、非一様な matrix データの構造は以下のようになります:
<N+1> <y0> <y1> <y2> ... <yN>
<x0> <z0,0> <z0,1> <z0,2> ... <z0,N>
<x1> <z1,0> <z1,1> <z1,2> ... <z1,N>
: : : : ... :
これらは以下のような 3 つの数字の組に変換されます:
<x0> <y0> <z0,0>
<x0> <y1> <z0,1>
<x0> <y2> <z0,2>
: : :
<x0> <yN> <z0,N>
<x1> <y0> <z1,0>
<x1> <y1> <z1,1>
: : :
===========
ところが、実際にこのようなデータを描画すると、
この x 座標と y 座標が逆の形で出力されていました。
バグだと思いそれを修正する (x と y を入れかえる)
コードを送ったこともあるのですが、
影響が大きいせいか採用されませんでした。
それがこの度、逆に現実の現象に合った形にドキュメントを修正する形で
ようやく「バグが修正され」ました。
現在は、上のドキュメントの前半は以下のようになっています:
===========
よって、非一様な matrix データの構造は以下のようになります:
<N+1> <x0> <x1> <x2> ... <xN>
<y0> <z0,0> <z0,1> <z0,2> ... <z0,N>
<y1> <z1,0> <z1,1> <z1,2> ... <z1,N>
: : : : ... :
===========
この問題は、山本さんの
「gnuplot の精義」
でも指摘されていました (「3.4.3 マトリックスバイナリー」)。
ようやく修正されたことになります。
(cf. 「情報やメモ (09/09 2013)」)
前回の報告 (06/03 2013)以後、 現在の CVS 版に入れられた主な機能について紹介します (ChangeLog, manual の更新部分等より)。
- 一般的な機能
- load コマンドで macro を 3 段階まで展開するように (06/28)
- 多角形のクリッピングに関する問題の修正 (06/28)
- fit で可能な独立変数の最大数を 5 から 12 に拡張 (06/27)
- 等高線ラベルにも陰線処理を (06/27)
- 陰線処理と等高線色に関する問題の修正 (06/24)
- key title に関する修正 (06/19)
- every で負の値を指定した場合はデフォルトの設定に (06/18)
- using の指定列数のチェックに関する修正 (06/14)
- timefmt に関するテストに関する修正 (06/14)
- show fit で前回の fit コマンドを報告するように (06/14)
- fit 用に仮変数名 set/show/unset dummy を改良 (06/12, 06/14)
- show/unset datafile {commentschars|missing|...} が動作しない問題の修正 (06/10)
- fit の出力ログの変数名のバグの修正 (06/05)
- "文字列"[x:y] の x,y に実数が使えるように (06/05)
- UTF-8 文字列への対応用の修正 (06/05)
- terminal ドライバ
- qt: 文字列の原点に関する修正 (06/28)
- win: console モード版の作成環境に関する修正 (06/24)
- win: バイナリモードでファイルを開くのに合わせた修正 (06/23)
- context: standalone で mono の場合に明示的に color を無効に (06/22)
- x11: XKeycodeToKeysym を XKbKeycodeToKeysym に変更 (06/16)
- win: 折れ線の始点と終点が同じ場合に閉じた多角形にするように (06/11)
- win: レイヤーのサポート (06/11)
- win: 箱付きラベル (文字列の回転もあり) のサポート (06/11)
- win: コンソール版 gnuplot で Ctrl-C で fit を中断するように (06/05)
- qt: Qt5 への対応 (06/05)
- その他
- マニュアルの誤植の修正 (Voigt 関数, %j,%M,%S, kdensity2d など) (06/10)
- Qt4/Qt5 の両環境用のコンパイル環境の整備 (06/05, 06/12)
- デモ
- contours.dem: ラベルの陰線処理に関する修正 (06/27)
- fitmulti.dem: 6 次元 fitting のテストサンプル (06/14, 06/27)
前回同様、fit に関する修正がかなり入っているようです。
また、上の修正にある「"文字列"[x:y] の x,y に実数が」 というのは、もちろんこの x, y は整数しか意味を持たないのですが、 その x, y として列データから取ってきた値を使えるようにするためだそうです。 もちろん、その実数の整数部分を利用するようになっています。
実は最近、gnuplot-4.6 対応になった山本さんの 「gnuplot の精義」 の第 2 版が出版されましたが、 それを読んで気がついたところなどを本家に報告して、 それが修正されたものも上に含まれています。 例えば、
- show/unset datafile {commentschars|missing|...} が動作しない問題の修正
- マニュアルの誤植の修正 (Voigt 関数, %j,%M,%S, kdensity2d など)
- every で負の値を指定した場合はデフォルトの設定に
などがそれです。また、上には書いてありませんが、 マニュアルの every の項目に「: で終わる書き方は不可」 という文言を追加したりもしています。 every は、デフォルトでいい部分は空欄にして 「:3::2」のような指定が可能なのですが、「:」で終わっていると それに続くオプションを解析した上でエラーになるので、 少しわかりにくいエラーメッセージが表示されてしまいます (私は初めて知りました)。 残念ながら、それを許すよう修正するのは原理的に難しそうなので、 そのための説明をマニュアルに追加したわけです。
「every で負の値」というのも上と少し関連した話なのですが、例えば
plot 'data' using 3:2 every pp::ps::pe with points ps 2のようなものが書かれたスクリプトを、gnuplot で
gnuplot -persist -e 'pp=2; ps=3; pe=100' file.gpのように、表示する点のサンプリング間隔、開始点、終了点を 変えながら実行している場合、 それぞれの値をデフォルトの値にしたいときは、 間隔や開始点は 1 や 0 を指定すればデフォルトになるのですが、 終了点のデフォルト、すなわち最後の点の番号は、 データに依存してしまって簡単に設定できません。 そこで、簡単に
gnuplot -persist -e 'pp=-1; ps=-1; pe=-1' file.gpのようにすれば (間隔や開始点も含めて) デフォルトに戻せるように、としてみたわけです。
「show/unset datafile {commentschars...}」の話はバグの修正です。 例えば gnuplot 4.6 では「show datafile」とすると commentschars の設定も表示されるのですが、 「show datafile commentschars」だと出ません (マニュアルにはある)。 また、マニュアルに書かれている「unset commentschars」(これ自身間違い) も、「unset datafile commentschars」も効かず、 他の missing や separator などのオプションも一緒にデフォルトにしてしまう 「unset datafile」しか使えなくなっていました。 現在の CVS 版では「show datafile commentschars」 「unset datafile commentschars」が使えるようになっています。 missing や separator なども同様に show, unset が個別に使えるようになっています。
(cf. 「情報やメモ (08/11 2013)」)
最近「自分で作成した gnuplot スクリプトでオプションをつけて実行したい」 という話をあるところで見ました。
Unix では、シェルスクリプトに実行許可属性を与えて、 1 行目に「#!」から始まる行をつけることで、 それをコマンドのように使えるようになります。 この話は、これを利用したもののようで、 例えば gnuplot の call 用のスクリプト (名前は file.gp とします):
#!/usr/local/bin/gnuplot
plot sin(x)
pause -1 "Hit ENTER key."
if ( "$#" < 2 ) exit
set term $0
set output "$1"
replot
set out
% ./file.gp 'png' 'file.png'のように実行したい、ということのようです。
もちろん、現在の gnuplot にコマンドラインから与えることができるのは 「load」コマンドで実行できるスクリプトだけなので、 このような「call」コマンド用のスクリプトをこのように実行しても うまくいきません。実際には、
% /usr/local/bin/gnuplot ./file.gp 'png' 'file.png'が行われるだけなので、gnuplot は 'png' や 'file.png' も gnuplot スクリプトだと認識してしまいます。
しかし、ソースを改造すれば gnuplot にそのようなことさせることも できなくはありません。 実際そのようなパッチを書いてみたのですが、 ふと気がついて、どうせ Unix ならば他の方法もあるのでは、 と考えてみたら、以下の 2 つの方法を思いつきました。
- シェルスクリプトのヒアドキュメントを利用する
これは、私がいつもよくやる手で、file.gp に代わるシェルスクリプト (ここでは csh script とします) file.csh を以下のように書けば、 上に書いたことが容易に実現できます:
引数の個数チェックは、シェルスクリプト側で行い、 それを viewonly という変数を介して gnuplot で処理しています。 また、シェルスクリプトへの引数 $1, $2 は、 直接ヒアドキュメント内にそう書けばそれがそのまま展開されますので、 普通の文字列として利用できます (シェルスクリプトでは $0, $1 ではなく $1, $2)。#!/bin/csh -f
set viewonly = 0
if ( $#argv < 2 ) set viewonly = 1
gnuplot <<EOF
plot sin(x)
pause mouse "Push mouse RIGHT button."
if ( $viewonly ) exit
set term $1
set output "$2"
replot
set out
EOFなお、gnuplot にヒアドキュメントでスクリプトを流すと、 それが標準入力からの入力となるので、 そのスクリプト内で標準入力を見る「pause -1」が使えません。 よって、上では「pause mouse」に置き換えています。
- 中間スクリプト (プログラム) を利用する
これは、gnuplot スクリプト側では、直接 gnuplot を呼び出すのではなく、 この中間スクリプト (プログラム) を呼び出しておいて、 その中間スクリプト (プログラム) 側で gnuplot を実行し、 与えられたオプション等は gnuplot が認識する形でその gnuplot に渡す、 という方法です。
中間プログラムは、シェルスクリプトでも、perl スクリプトでも、 C のプログラムでも、何でもいいでしょう。 中間プログラムが gnuplot を呼ぶときには、 実際に call を使う形式で呼び出す方法もあるでしょうし、 -e オプションを利用して、 内部変数を設定するコマンドをスクリプトの前に書く、 という方法もあるでしょう。 以下に、call を使って gnuplot を呼び出す中間 csh スクリプト (gpexp.csh) の一例を紹介します:
これが呼び出されるときは、 このスクリプト内部では、 最初の引数 ($1=$argv[1]) は対象となる gnuplot スクリプトで、 それに続く引数は、gnuplot スクリプトに与えたオプションなので、 基本的には引数全体 ($argv) の頭に call をつけた文をほぼ#!/bin/csh -f
set gp = /usr/local/bin/gnuplot
if ( $#argv < 1 ) exit
set s = ""
while ( $#argv > 0 )
set s = "$s '$argv[1]'"
shift
end
$gp -e "call $s"
gnuplot -e "call $argv"
のようにすればいいだけです。 しかし、call スクリプト自身は引用符で囲まなければいけませんし、 また call スクリプトの引数も、ドキュメントに書かれている通り、 引用符で囲まないと動作が保証されませんので、 それらを一々引用符で囲む作業を行っているのが while から end までの部分です。 そして、実際はそうやって作った文字列 s を call に渡す、 という文字列を gnuplot -e で実行しているだけです。なお、このスクリプトに実行許可属性を与えて、 /usr/local/bin/gpexp.csh としたとしても、 これがシェルスクリプトなので、これを使用する gnuplot スクリプト側で
のようにして使うことはできません。 これは、/usr/bin/env を利用して、#!/usr/local/bin/gpexp.csh
...
のようにして使います。 上に書いた file.gp も、その先頭行をこれに変えて、 実行許可属性を与えれば、#!/usr/bin/env /usr/local/bin/gpexp.csh
...% ./file.gp 'png' 'file.png'
のようにして使えるようになります。
前回の報告 (04/30 2013)以後、 現在の CVS 版に入れられた主な機能について紹介します (ChangeLog, manual の更新部分等より)。
- 一般的な機能
- fit のχ二乗和に Kahan-Bubuska アルゴリズムの Neumaier 変種を利用 (06/02)
- splot や contour が無意味なスタイル (histgram 等) で指定された場合に points で描いていたのを、警告をして無視するように変更 (05/31)
- fit の収束、中断後にユーザ変数をすべて更新するように (05/27)
- fit のファイルアクセスと文字出力に関するコードの修正 (05/23)
- columnstack histograms での linetype に関する問題の修正 (05/23)
- 日時データに列のヘッダ文字列指定が効かなかった問題の修正 (05/23)
- 多角形のクリッピングに関する修正 (05/22)
- データに \r が含まれるときの問題に関する修正 (05/22)
- set style textbox の新設 (05/19)
- 空の塗り潰し多角形を描画しないように修正 (05/18)
- 名前付きの列がなかった、部分的に一致した場合の警告に関する修正 (05/18)
- set contour で splot with labels のサポート (05/15, 05/18)
- set fit ... default の読み込みに関する問題の修正 (05/14)
- set fit wrap <column> で長い fit レポートの折り返しができるように (05/14)
- fit 時のユーザエラー変数の設定に関する修正 (05/12)
- set clabel に代わる set cntrlabels の新設 (05/11, 05/15, 05/18)
- set fit limit_abs で収束条件を変更できるように (05/09)
- hypertext point に関するオプションの順序依存がないように (05/08)
- splot の hidden3d での等高線の線種の選択の修正 (05/08)
- データファイルに関する変数の初期化の問題の修正 (05/08)
- fit のχ自乗値の状態変更 (余計な繰り返しの防止) (05/08)
- fit の関数未定義エラーの情報の表示 (05/08)
- fit の 0 割り算エラーが起きる問題の修正 (05/08)
- fit のχ自乗の増加を -1 でなく * で表示 (05/08)
- set fit のいくつかのオプションのデフォルトをキーワード 'default' で (05/08)
- set fit のオプション script の追加 (05/08)
- fit 用のユーザ変数 FIT_P (χ自乗検定の p-値) の追加 (05/08)
- fit に関する文字列処理、NaN データ処理の修正 (05/08)
- ユーザ変数 GPVAL_LAST_FIT の追加 (05/08)
- replot の際の plot/splot のチェックの修正 (05/06)
- 入力データの Inf, NaN の処理の問題の修正 (05/05)
- 多角形のクリッピングルーチンの修正 (05/03)
- set multiplot の自動配置のスキップ、戻りを行う next|previous オプション (05/03)
- terminal ドライバ
- win: ファイルを常にバイナリモードで開くように (06/02)
- win: EMF やクリップボードへの出力はアンチエイリアスを切るように修正 (06/02)
- win: プリンタ (PRN) に出力するとクラッシュする問題の修正 (06/02)
- wxt: ウィンドウを xkill などで閉じた場合に起きる問題の修正 (05/31)
- emf: フォントの初期化に関する問題の修正 (05/20)
- svg: テキスト背景の不透明塗り潰しのサポート (05/19)
- gd/post/x11/xlib/qt/wxt/cairo: set style textbox のサポート (05/19, 05/21)
- qt: center の位置合わせに関する問題の修正 (05/18)
- qt: enhanced 文字列要素を描画要素グループに追加するように (05/18)
- cairo: enhanced モードでの save/restore 後の y 座標に関する問題の修正 (05/17)
- canvas: パターン塗り潰しのサポート (05/08)
- qt: 多角形描画の空の塗り潰しを白から背景色に (05/06)
- canvas: マウスの出入りにともなう不要な再描画の防止 (05/06)
- x11: 画像を受けとる前にグラフを作成する問題の修正 (05/05)
- その他
- with labels に関する図をドキュメントに追加するように (05/31)
- C99 関数 isblank() の削除 (05/08)
- デモ
- all.dem, enhancedtext.dem: UTF-8 でない post terminal 上の 拡張文字列デモが止まらないように (05/24)
- contours.dem, datastrings.dem: 箱付き文字ラベルを追加 (05/19)
- contours.dem: set cntrlabels に関するデモの更新 (05/15)
- transparent.dem: canvas terminal でのパターン塗りのサポート (05/08)
- fit.dem: p 値の出力 (05/08)
- complex_trig.dem: エンコーディングに関する修正 (05/08)
- html/webify_canvas.pl: 出力 HTML にゴミを出さないよう修正 (05/08)
- imageNaN.dem: 入力データの Inf, NaN の処理の修正に伴なう修正 (05/05)
一応前回の報告から 1 ヶ月くらいしか経っていないのですが、 今回はかなり変更などが入りました。 大きな改良は以下のあたりでしょうか。
- set contour で splot with labels が使えるように
- fit に関するオプション、ユーザ変数、表示などのさまざまな改良
- set multiplot の自動配置向けオプション next,previous
- set cntrlabel の追加 (set clabel の改良版)
- textbox スタイルの追加 (set label boxed, with labels boxed に影響)
- canvas terminal のパターン塗り潰しのサポート
Gnuplot Q&A 掲示板 (2653 番の記事) で報告があった、 emf terminal の test 出力のフォントの問題についても修正が採用されています。
いずれも、ほぼ想像がつくもの、 マニュアルを見ればだいたいわかるものでしょうから、 今回は具体例は省略したいと思います (実は面倒くさいだけ ^^;)。
河野＠ロスアラモス研 さんが立ち上げている 「Gnuplot - not so Frequently Asked Questions -」 (http://t16web.lanl.gov/Kawano/gnuplot/) というサイトにお世話になった方も多いのではないかと思いますが、 最近そのサイトがアクセスできなくなっているようです。 以下にその情報があります。
- Gnuplotting
Not so Frequently Asked Questions site down (11/11 2012) - News グループ comp.graphics.apps.gnuplot
Dr Kawano's not-so-faq website not available (11/09 2012)
出版されるのであればそれはそれでうれしいですが、サイトがないのは残念です。 そのため、上にも書いてありますが、現在はミラーがいくつかあるようです。
- ミラーサイト 1:
http://www.ualberta.ca/~xz10/gnuplot/ - ミラーサイト 2:
http://web.archive.org/web/20100801013014/http://t16web.lanl.gov/Kawano/gnuplot/ - ミラーサイト 3 (古めのもの ?):
http://folk.uio.no/hpl/scripting/doc/gnuplot/Kawano/
このサイトは日本語版だけでなく英語版があるせいもありますが、 外国でそういうことが話題になること、 および次々とミラーができていることが、 河野さんのサイトの重要性を物語っているように思います。
前回の報告 (04/02 2013)以後、 現在の CVS 版に入れられた主な機能について紹介します (ChangeLog, manual の更新部分等より)。
- 一般的な機能
- 変数をファイルに保存する際にデータブロックを無視するように (04/27)
- fit のバックアップ処理が上書きせずにエラーを表示するように (04/27)
- fit でパラメータファイルが存在する/しない場合の処理の修正 (04/27)
- DOS の CSV データ (\r\n) も変換なしに読めるように (04/26)
- CSV データの抜けを histogram 用に不正でなく欠けと判断するように (04/26)
- インラインデータの終了処理 (\n, \0) の修正 (04/26)
- 多角形の長方形クリッピング用の Sutherland-Hodgman アルゴリズムの実装 (04/25)
- fit の一行要約機能 (04/24)
- fit のパラメータが 2 つ以上の場合のみ相関行列を表示するように (04/24)
- show variable の余計な改行を抑制 (04/24)
- fit で収束しなかった場合に新たにメッセージを出すように (04/24)
- set fit に出力用オプション quiet, results, brief, verbose を追加 (04/24)
- set key font '[何か]' ; reset でおきるメモリーリークの修正 (04/22)
- set datafile separator で 2 文字以上を指定できるように (04/21)
- set fit に (FIT_* の代用の) オプション maxiter, limit を追加 (04/20)
- 大きな行列に対する stat コマンドのメモリ割り当ての問題の修正 (04/12)
- with labels では set label と異なり at, front, back を受け付けないように (04/09)
- x2,y2 軸の自動縮尺の初期状態に関する問題の修正 (04/07)
- 描画の refresh 時に軸の範囲に余計な小さな因子を追加しないように (04/07)
- 自動縮尺機能の設定のチェックの改良 (04/07)
- set range に noextended オプションの追加、 オプションの順序依存がないように修正 (04/06)
- 時間座標の原点を 2000 年元日から 1970 年元日へ変更 (04/05)
- refresh 時に自動縮尺の "reverse" の元の状態を復帰 (04/04)
- terminal ドライバ
- win: 初期ファイルの設定用の UI の追加 (04/24)
- win: マウスによるグラフの回転のスピードアップ (04/24)
- win: title オプションが適正かのチェック (04/22)
- dumb: グラフのアスペクト比を変更するオプション aspect の新設 (04/24)
- wgnuplot で長い fit を Ctrl-Break で中断できるように (04/22)
- win/wxt: コンソールを上に出すルーチンのコードの修正 (04/22)
- qt: 画像描画を ON/OFF するようなピクスマップを追加 (04/22)
- wxt (cairo): バックスラッシュをフォント名のエスケープ文字に使えるように (04/22)
- wxt: 初期化に失敗した場合に安全に終了するように (04/22)
- x11: 凡例上での左クリックでグラフの ON/OFF が可能に (04/21)
- cairolatex: fillcolor の設定が無視される問題の修正 (04/10)
- x11: リサイズの計算を全体の高さよりグラフ領域の高さを元にするように (04/08)
- デフォルトの出力形式リストから gpic, hp2623a, hp2648 を削除 (04/08)
- wxt: MS-Windows 上でスペースレイズコンソールを実装 (console 版 gnuplot 用) (04/04)
- win: -persist モードでの問題の修正 (04/04)
- post: level2 の画像が 64k バイトを越えた場合に起きる問題の修正 (04/02)
- post: level1 のときの level2 の機能をスキップする部分の改良 (04/02)
- その他
- emacs-lisp ファイル, figure_missing.pdf のインストール環境の修正 (04/25)
- emacs-lisp ファイルに関するインストールディレクトリの問題の修正 (04/10)
- デモ
- fit.dem: pause やグラフなどの削減、タイトルの変更 (04/24)
- fit.dem: set fit のオプション形式の変更に伴なう修正 (04/20, 04/24)
今回は、fit 関係の修正がだいぶ入りました。 あまり fit は使わないので詳しくはよくわかりませんが、 環境変数だったものが set fit などで指定できるようになり、 gnuplot スクリプトで閉じた処理が行いやすいように なっているのではないかと思います。
また、時間の原点が、gnuplot 独自の「2000 年元旦からの秒数」から、 Unix 標準の「1970 年元旦からの秒数」に変更されました。 標準的なものになったので歓迎すべきことだと思いますが、 もしかするとこれまでの gnuplot スクリプトを 修正しないといけない場合もあるかもしれません。
set datafile separator では、これまでは 1 文字指定しかできなかったのが、 "*|" (区切りは '*' または '|') のような指定ができるようになりました。 これもある種のデータに関しては前処理が少し不要になるような 便利な機能なのかもしれません。
(cf. 「情報やメモ (06/03 2013)」)
gnuplot-4.6.3 がリリースされています。 配布物に含まれる NEWS には、以下のように書かれています。
- 新規: Windows 上のコンソール版 gnuplot のスペースレイズコンソール
- 変更: wgnuplot で -persist モードでテキストウィンドウを開かないように
- 修正: Windows で -persist モードが不正だった問題
- 修正: Windows の wxt 出力形式で -persist モードがゾンビプロセスを残した問題
- 修正: "set term ... mono" の後は linetype の色を隠すように
- 修正: latex 出力形式でグラフと文字の色の変更の同期
- 修正: 大きな行列に摘要する stats コマンドに関するメモリ割り当て
見てわかる通り、今回は MS-Windows 用の修正が主です。 4.6.1, 4.6.2 用の MS-Windows 版バイナリは結局リリースされませんでしたが、 今回の 4.6.3 は最初から MS-Windows 版バイナリ (gp463-win32.zip) も、インストーラ形式 (gp463-win32-setup.exe) のものも、 ソースの公開と同時にリリースされています。
ところで、「スペースレイズコンソール」というのは、 コンソール版の gnuplot.exe に関する修正 (実装) のようです。 これをコマンドプロンプトから起動した場合、 グラフをプロットさせるとまずは wxt terminal や win terminal のグラフウィンドウが前面に出てコマンドプロンプト画面を隠しますが、 その状態からコンソール入力に戻るには、 今まではグラフウィンドウをマウスで閉じるか、 コマンドプロンプトをマウスで前面に出す必要がありました。 それが、スペースキー一発でコマンドプロンプトを前面に出し、 コマンド入力ができるようになるようになったのですが、 それが「スペースレイズコンソール」らしいです。
なお、マニュアルは 4.6.2 と 4.6.3 では変更がありませんので、 今回は 4.6.3 用のマニュアルの訳は行いません。 もうしわけありませんが、必要ならば 4.6.2 用のもので代用してください。
報告が遅れましたが、gnuplot-4.6.2 がリリースされています。 配布物に含まれる NEWS には、以下のように書かれています。
- 新規: "bind" コマンドでマウスボタン 1 にユーザコマンドを定義できる
- 新規: hidden3d の pm3d 曲面で (hidden2.dem のように) 他の描画要素を消す
- 新規: 初期設定ファイル ~/.gnuplot をスキップするコマンドラインオプション -d
- 新規: シェル呼び出し中に plot '<&N' でファイルデスクリプタ N から描画
- 変更: "unset term" は元のデフォルト出力形式 (GNUTERM) を復帰する
- 変更: plot コマンドに最後につけた不要なコンマは無視する
- 変更: 一様な binary matrix データの場合に入力が早くなる特別なルーチン
- 変更: 動作が対話型か非対話型かのチェックの変更
- 変更: 軸 (zeroaxis) を grid, border と同じ階層で描画
- 変更: yerrorbars で 2 列のみのデータを受けつける (x 座標は暗黙)
- 変更: using "COLUMNHEADER" で文字列が完全に一致することの義務付け
- 修正: aquaterm での rgbimage 描画スタイルでのレンダリング
- 修正: gd terminal での set_font(",newsize") でのフォントサイズの変更
- 修正: qt terminal のフォントメトリック
- 修正: qt terminal の -persist でのマウス動作、リサイズ動作
- 修正: x11 terminal の線幅のリセットに失敗する問題
- 修正: x11 なしでコンパイルしたときに -persist オプションが失敗する問題
- 修正: -persist で "set view map" グラフで 3D マウス動作を継続する問題
- 修正: 文字列入力解析でのエラー後に解析の状態をリセットする
- 修正: 非常に長い刻みラベル書式に関するバッファオーバーフローの問題
- 修正: 'using 1:2 "<format>"' による書式化された入力処理に失敗する問題
- 修正: カッコ囲みの中から現在のコマンドファイルを exit できるように
- 修正: x 軸の刻みラベルの回転時のスペースの評価
- 修正: volatile データの refresh 時の部分範囲指定 (例: xrange [*:MAX]) の問題
最近の CVS 版に入った機能などもいくつか追加されているようです。 ところで、gnuplot-4.6.0 はリリース後に MS-Windows 版バイナリや、 MS-Windows 用インストーラなどもリリースされたのですが、 gnuplot-4.6.1 のものはついに出ませんでした (もしかしたら気がついていなかっただけかも)。 4.6.2 用のものは、まだ公式ページには「not yet」と書かれていますが、 テスト版が http://www.gnuplot.info/development/binaries/ に置かれているようです。 そのうちに正式にリリースされるだろうと思います。
前回の報告 (01/26 2013)以後、 現在の CVS 版に入れられた主な機能について紹介します (ChangeLog, manual の更新部分等より)。
- 一般的な機能
- 3D の部分範囲指定の初期化の問題の修正 (04/01)
- histogram での "smooth frequency" が スタイルを上書きしないように (04/01)
- monochrome 時の linetype に関する動作の修正 (03/21)
- plot の fillcolor オプションの動作に関する修正 (03/21)
- マウス動作の改訂、グラフの中心でなく現在位置でのズーム (03/17)
- 要求した列ヘッダがない場合に警告を出すように (03/14)
- set fit に [no]prescale オプションが追加 (03/12)
- refresh の際に範囲を正常に延長するように (03/09)
- with yerrorbars で 2 列データをサポート (03/03)
- 片側のみの範囲指定が volatile データの refresh 時に起こす問題の修正 (03/01, 03/09)
- 個別の範囲指定が対数軸で起こす問題の修正 (02/28)
- 標準入力データで volatile の際 replot を refresh と解釈するように (02/28)
- 対話か非対話かの判別の改良 (02/21)
- 対数軸の底を線形軸でも残すように修正 (02/19)
- x 軸の回転したラベルは位置を右合わせするように (02/19)
- zeroaxis を axis や格子線と同じ階層に描くように (02/18)
- 新しいファイル指定: '<&N' (ファイルデスクリプタ番号指定) (02/17)
- binary matrix で大きなデータを読み込む場合の改良 (02/17)
- -e オプション時に stdin からの入力を待たないように変更 (02/17)
- 02/03 の修正で Gnuplot.pm perl module に起きる問題への対処 (02/14)
- set format による空の文字列の初期化に関する問題の修正 (02/08)
- 非 X11 gnuplot で -p を使ったときの問題の修正 (02/03)
- 初期設定ファイルの読み込みをスキップするコマンドラインオプション -d/--default-settings の新設 (02/03)
- cat | gnuplot 時の mouse zoom で起きる問題に関する修正 (02/02)
- terminal ドライバ
- post: 画像の /FlateDecode 圧縮のサポート (04/01)
- svg: 描画タイトルを title 属性にするように (03/29)
- qt: zoom や replot 後のループに関する問題の修正 (03/29)
- x11: ウィンドウマネージャによりリサイズ時に起きる問題の修正 (03/09)
- win: 関数 (ungetch) の再定義に関する問題の修正 (03/02)
- imagen: デフォルトのドライバ群から削除 (02/26)
- aquaterm: rbgimage に関する問題のため 2012 04/24 のパッチを修正 (02/24)
- cgi: 削除 (02/22)
- wxt: 多くの notitle なグラフが title つきのグラフの前に 描かれる場合に起きる問題の修正 (02/19)
- cairo/canvas/svg/wxt: 線分の端 (round/butt) に square を新設 (02/18)
- その他
- graphics.c の整理、分割 (02/27)
- lasergnu スクリプトの削除 (02/26)
- 起動画面から Build System の行を削除 (02/26)
- bf_test に関する Makefile の整理 (02/17)
- デモ
- fit.dem: set fit [no]prescale に関する修正 (03/12)
いくつか補足をします。
- (03/17) マウス動作の改訂は、ズームをグラフの中心を元に行うか、 マウスカーソルの現在位置を元に行うか、のあたりの変更のようです。 ただ、試してみたのですが、よくわかりませんでした。
- (03/12) set fit の prescale オプションというのは、 Marquardt-Levenberg ルーチンの実行前に パラメータの正規化を行うものだそうです。 パラメータの初期値のスケールがバラバラで、 大きさかなり違いがあると Marquardt-Levenberg ルーチンはあまりうまくいかないようで、 以前のマニュアルにも 「パラメータの大きさはあまり違いすぎてはいけません。 絶対値が最も大きいパラメータと最も小さいパラメータの比が大きい程 当てはめの収束は遅くなります。」と書かれていて、 その場合は自前でスケーリングするようヒントが書かれていました。 それを gnuplot 側でやらせるようなオプションのようです。
- (03/03) 2 列データの場合、yerrorbars は x 座標を暗黙 ($0) とし、 1 列目を y、2 列目を yerror とします。
- (02/19) 対数軸の底を線形軸でも残すように、 というのは軸の format 指定の話です。 軸の見出しとして指数表記をする場合、 %L, %l がその指数と仮数になります。この底の値を変えたい場合、 マニュアルでは「現在の対数尺の底での」とあるのですが、 実は logscale 時にしかその底の値は指定できないようになっていて、 さらに線形軸に戻したときに (unset logscale) その底の値は 10 に戻ってしまい、 つまり線形軸では %L 等の底を変更できないようになっていました。 今回、logscale から線形軸に戻したときにも その際の底の値を保持するようになりましたので、 線形軸でも 2 の巾、100 の巾などの指数表記の軸目盛りが 容易に作れるようになっています。
- (02/14) Gnuplot.pm perl モジュールの問題は、 コマンドラインにファイル名が指定されていない場合、 従来は gnuplot は常に対話モードになっていましたが、 その判別ルーチンを改善したようです。 多分 02/21 の改良とも関連すると思いますが、 現在どのようになっているのかは詳しくは調べていません。
- (04/01) Postscript terminal に /FlateDecode 圧縮なるものが あらたにサポートされましたが、これは、set term post のオプションで level3 を指定したときに有効になるもののようです。 ビットマップ画像が PNG 画像として埋め込まれるようになり、 生成される出力のサイズが大幅に改善されるようです。 PNG は、libgd か libcairo のいずれかを使用するようです。
(cf. 「情報やメモ (04/30 2013)」)
先日、「情報やメモ (01/22 2013)」 に、xtics の目盛りラベル文字列をユーザ定義文字列値関数とできれば、 とか書きましたが、良く考えてみたら目盛りラベルも、目盛り位置も、 今の set for の機能を使えばかなり容易に カスタマイズできるんでした。 こんな具合です:
s(x) = sprintf("%.1f=%d+%.1f",x,floor(x),x-floor(x)) # 目盛りラベルの文字列値関数
p(x) = x - 0.3*sin(0.5*pi*x) # 目盛り位置の関数
set xrange [0:5]
unset xtics
set for [i=0:5] xtics ( s(p(i)) p(i) )
plot x
- 上のスクリプトの画像
なお、この手の話は、すでに 4 年以上前にしていました (^^; (「情報やメモ (09/09 2008)」 「情報やメモ (09/11 2008)」 参照)
前回の報告 (12/21 2012)以後、 現在の CVS 版に入れられた主な機能について紹介します (ChangeLog, manual の更新部分等より)。
- 一般的な機能
- bind のコマンド文字列として文字列値変数を使えるように (01/25)
- set timefmt で文字列値変数を受けつけるように (01/24)
- pm3d 曲面の hidden3d 処理を他のオブジェクトも隠すよう改良 (01/21)
- 1 行関数 handle_edge_fragment() のインライン化 (01/08)
- palette が使用されていないときにも線色が適用されるよう (cairolatex 等の出力形式用) (01/08)
- アスペクト比に関する問題の修正 (01/07)
- 不正なコマンドが検出された場合の入力バッファの扱いに関する修正 (01/06)
- readline 用の標準入力に関する修正 (01/06)
- 頂点が 2 つ未満の多角形描画をスキップするように修正 (01/05)
- グラフの描画と multiplot の設定の認識順に関する修正 (01/05)
- set timemft で軸指定ができるように (12/29)
- 文字列を前に出すよう plot title と軸ラベルの描画ルーチンを分離 (12/27)
- terminal ドライバ
- quote_str() が next/openstep 以外では使われないことに関する修正 (01/25)
- ai/vgagl/corel/djsvga: フォントに文字列値変数を使えるよう修正 (01/25)
- x11: 文字列の右、中心合わせのフォントサイズに関する修正 (01/18)
- pslatex: TeX 出力の文字色をレイヤ毎にリセットするように (01/08)
- x11: アスペクト比に関する問題の修正 (01/07)
- qt: GE_fontprops でフォント情報を送ろうとしないように (01/05)
- x11: 点線描画、replot の処理に関するコードの修正 (01/05)
- x11: GE_fontprops に関する入力待ち処理に関する修正 (01/05)
- x11: multiplot 時のマウス処理、replot に関する修正 (01/05)
- emf: 垂直方向のセンタリングがフォントサイズに合うように (12/27)
- wxt: 垂直方向のセンタリングとオーバーサンプリングに関する修正 (12/27)
- x11: 前に要求したフォントリストを保存するように (12/25)
- x11: UTF-8 encoding 時のエラーに関する修正 (12/25)
- その他
- configure 時の disable-volatile-data 指定に関する修正 (01/04)
- gnuplot.texi の作成に関する修正 (12/29)
- クロスコンパイル時の CC, CXX の定義に関する修正 (12/27)
- デモ
- hidden2.dem: pm3d 曲面の隠面処理の改良に合わせた修正 (01/21)
bind で文字列値変数を受けつけるように、というのは、
STRING = "(なんらかの gnuplot コマンド)"とできるようにした、ということのようです。 確かに、複雑なコマンド列を bind したい場合には便利かもしれません。
bind (キー) STRING
(cf. 「情報やメモ (04/02 2013)」)
x11 や wxt terminal で、グラフウィンドウ上で 'h' を打つと、 そのグラフウィンドウ上で利用できるマウス機能、 キーボード機能の一覧が表示されますが、 その一覧 (ヘルプ) は日本語化されていません。 時間がとれたら、日本語表示が可能かどうか検討したいと思いますが、 とりあえず、以下にその訳をあげておきます。
| マウス (<B1> = 左、<B2> = 中、 <B3> = 右) | |
|---|---|
| <B1> 2 回 | `clipboardformat` を用いてクリップボードに座標を出力 (キーの '3', '4' を参照) |
| <B2> | `mouseformat` (キーの '1', '2' を参照) を用いてグラフに付箋付け、 `set mouse labels` が ON ならラベル描画 |
| <Ctrl-B2> | `set mouse labels` が ON ならカーソルに近いラベルの削除 |
| <B3> | ズーム範囲のマーク (2D 描画と maps のみ) |
| <B1-Motion> | 視方向の変更 (回転)。軸のみの回転には <ctrl> を併用。 |
| <B2-Motion> | 視方向の変更 (伸縮)。軸のみの伸縮には <ctrl> を併用。 |
| <Shift-B2-Motion> | 垂直移動 -- xy 平面の変更 |
| <wheel-up> | 上スクロール (Y の増加方向へ) |
| <wheel-down> | 下スクロール |
| <shift-wheel-up> | 左スクロール (X の減少方向へ) |
| <shift-wheel-down> | 右スクロール |
| <control-wheel-up> | グラフの中心方向へズームイン |
| <control-wheel-down> | ズームアウト |
| <shift-control-wheel-up> | X 軸のみズームイン |
| <shift-control-wheel-down> | X 軸のみズームアウト |
| キー | |
| Space | gnuplot 対話ウィンドウを上に |
| q | 描画ウィンドウを閉じる |
| a | 自動縮尺 (set autoscale keepfix; replot) |
| b | 境界描画の ON/OFF |
| e | replot |
| g | 格子線 (grid) の ON/OFF |
| h | help (これ) |
| l | 対数軸の ON/OFF (plot では y 軸、splots では z と cb 軸) |
| L | カーソルが一番近い軸への対数軸の ON/OFF |
| m | マウスモードの ON/OFF |
| r | 十字線 (ruler) の ON/OFF |
| 1 | ひとつ前の mouseformat に |
| 2 | ひとつ次の mouseformat に |
| 3 | ひとつ前の clipboardformat に |
| 4 | ひとつ次の clipboardformat に |
| 5 | 十字線原点からの動径表示 |
| 6 | 対話ウィンドウへの冗長表示 |
| 7 | ratio を巡回的に変更 |
| n | 次のズーム履歴へ移動 |
| p | 前のズーム履歴へ移動 |
| u | ズーム解除 |
| Right | 右回転 (splot のみ)、<shift> で速く回転 |
| Up | 上回転 (splot のみ)、<shift> で速く回転 |
| Left | 左回転 (splot のみ)、<shift> で速く回転 |
| Down | 下回転 (splot のみ)、<shift> で速く回転 |
| Escape | 選択中のズーム領域の解除 |
文章の間違いらしきものも見つけましたので後で報告したいと思いますが、 普段これらの機能はほとんど使っていなかったので、 今回これらを試してみたら結構面白かったです。 特にマウスホイールによるスクロールは、2D だけでなく splot でも有効なので、実行してみるとちょっと不思議な感じがします。 また、キーの '5' (キーの 'r' の後にやる) もこんな機能があるとは知りませんでした。 グラフ上で座標を確認したい人には便利かもしれませんね。
2 ちゃんねるの掲示板「gnuplot を使おう。その 3」(アドレスは リンクリスト 参照) に書かれた情報について、多少気になるものがありましたので、 少し書いておきます。
-
「情報やメモ (08/23 2012)」 にも書いた logscale の切り替えの話題ですが (86,88-94)、 元の質問者が詳細な情報を書いたところ、 bind で if を使ったトグルスイッチ風の割り当てを行うサンプルが紹介され、 それで解決したようです (95-96)。
bind に複雑な式を割り当てるのはなかなか高度な技ですね。 gnuplot の付属のマニュアルにも、 同じようなトグルスイッチ風の少し複雑な bind の例として、 以下のようなものが書いてあります:
bind "ctrl-r" "v=v+1;if(v%2)set term x11 noraise; else set term x11 raise"
対数軸の切り替えについては、「Gnuplot in Action」 (「gnuplot に関する文献等」を参照) の .gnuplot (gnuplot の初期設定ファイル) のサンプルの中 (Chapter 12. Macros, scripting, and batch operations) に、 丁度 95 番の記事の例に似たようなものが載っています:is_log_x = 0;
この本では、変数名を見てわかるように、x 軸、y 軸の対数軸の切り替えを 別々に行うためにこのようなものを書いているようです。 なお、95 番の記事は上のような replot ではなく、 直接 plot 命令が書いてありましたが、 このように replot を使う方が汎用性は高くなると思います。
bind 'x' "if(is_log_x) is_log_x = 0; unset logsc x; replot; \
else is_log_x = 1; set logsc x; replot -
x 座標が 80, 90, 100 などのときに x 軸の目盛りラベルを指数表記にしたとき、 その指数部分を 101 と固定して、 「8.0x101, 9.0x101, 10x101」 となるようにすることはできるか、 set format x "%l2.1lx10^%L" ではうまくいかず 最後のものが 1.0x102 となってしまう、 という質問がありました (110)。
これに対して、10x101 という書き方は 指数表記にする意味はないのでは、 軸名の横に「x 10」のように書いて 「using ($1/10):2」で x 座標を実際には 1/10 にしてラベルをごまかしたら、 という回答がありましたが、 元記事の方はデータではなく、関数描画なのだそうです。
指数表記というよりも 10 で割った数を軸の見出しにしたいなら、 「(x 10)」のように軸に書く方が確かに普通だと思います。 データを 1/10 にするやり方は、もちろん関数描画でも行えます。 「plot [8:10] f(x*10) title "f(x)"」のようにするだけです。 そうすれば x 座標は 8 から 10 まで動きますが、 関数には 80 から 100 までの値が代入されることになります。
なお、昔から考えているのですが、 現在の gnuplot では文字列値関数が定義できますから、 軸の目盛りラベルの文字列もユーザ定義文字列値関数とできれば、 もっと目盛りラベルの自由度が上がるような気がしています。 set format との共存などがちょっと難しそうですが、 明示的な set xtics などがだいぶ楽になるだろうと思いますので、 実現できればなぁと思っています (誰かやりませんか ?)。
また、マニュアルでは %l, %L は「現在の軸の底に基づく」 となっているのですが、 その「現在の軸の底」を変更できるのは対数軸の場合だけで、 線形軸では底は 10 しか使えないようになっています。 これを、線形軸でも任意の底が使えるように修正が可能か (例えば 50x1003 や 1.2x27 などが可能になる)、 後で考えてみたいと思います。
-
0 から 20 までの、10000 個のデータが 1 列 (10000 行) 並んでいるデータがあるが、 それの 5 幅刻みのヒストグラムを書きたいが、 gnuplot でできないか、それとも別な方法で数を数えるしかないか、 という質問がありました (117)。
これに対して、直接ヒストグラムを書くのは無理、 基本的にグラフ描画ツールだから、 といった意見がありましたが、 元記事の方はうまく探したみたいで、 「smooth frequency」(gnuplot 4.0 以降) を利用した ヒストグラムを描くスクリプトを上げています (123)。
gnuplot に付属する smooth.dem にもサンプルがありますが、 smooth という本来データを補間するための機能の中に、 gnuplot-4.0 から smooth frequency や smooth cumulative (4.4.0 以降) などの簡単な統計的なグラフを書くためのオプションが追加されています。 これを利用するとヒストグラムはそれなりに書くことができます。 smooth frequency は、x 座標が同じデータ点を一つに集約して、 それらの y 座標の合計を y 座標とするデータにまとめます。 例えば、
1 2
のようなデータの場合、smooth frequency は、
2 3
1 4
3 5
2 61 6
にしてくれるわけです。
2 9
3 5よって、例えば 0 以上 20 未満の範囲の データに対するヒストグラムを描く場合は、 using を利用して以下のようなデータに変換して読ませます:
- データが 0 以上 5 未満ならば x 座標を 2.5 (0 と 5 の真ん中)、 y 座標を 1 とする
- データが 5 以上 10 未満ならば x 座標を 7.5、y 座標を 1 とする
- データが 10 以上 15 未満ならば x 座標を 12.5、y 座標を 1 とする
- データが 15 以上 20 未満ならば x 座標を 17.5、y 座標を 1 とする
データとこの x 座標の関係は階段関数ですから、 floor() を使えば簡単に実現できます:
using (floor($1/5.0)*5.0+2.5):(1)
y 座標は 1 なので、(1) とかっこが必要です。そしてそのような入力に対して smooth frequency を使えば、 y 座標の 1 の合計、すなわちそのようなデータの個数が 同じ x 座標毎にまとめられるので、結局、
2.5 (0 以上 5 未満のデータの個数)
のようなデータを描画するのと同じことになるわけです。 あとは、これを with boxes で書けばめでたくヒストグラムのできあがりです (s f は smooth frequency):
7.5 (5 以上 10 未満のデータの個数)
12.5 (10 以上 15 未満のデータの個数)
17.5 (15 以上 20 未満のデータの個数)plot 'data' u (floor($1/5.0)*5.0+2.5):(1) s f w boxes fs s 0.5
123 の記事にはより洗練されたスクリプトがのっていますので、 全体のスクリプトはそちらを参照してください。なお、元記事の人は、さらに
using (floor($1/5.0)*5.0+2.5):(1):(0.5) with yerrorbars
として y 誤差線をつけたいが、 これだと高さ 1 のところに点が描かれるだけ、と書いていますが、 これだと smooth frequency がないので 当然 y 座標は高さ 1 になるだけですし、 誤差線も 1 万行だと y 座標のスケールから考えて 0.5 という幅の誤差線は、 点になっても当然です。 しかし、これを例えばusing (floor($1/5.0)*5.0+2.5):(1):(100) s f with yerrorbars
としてうまくいくかというとそうではありません。table 出力を出させてみればわかりますが、 現在の smooth frequency は 2 次元データしか見ておらず、 3 番目以降のデータは y 座標と同じにしてしまうのです。 ただ、これも本来は、y 座標と同じ処理を 3 番目以降のデータにもほどこすのが筋のような気がします (実際、そのようなコードも入っているのですが、 smooth frequency ではその部分を使っておらず、 y 座標で上書きしてしまいます) ので、 時間のあるときに修正可能か考えてみたいと思います。
-
MS-Windows から Xming で Linux の gnuplot を起動して、 その X の画面を Alt+PrintScreen でキャプチャして MS-Windows 上で PowerPoint などへコピーペーストしているのだが、 その X のウィンドウに枠やタイトルバーがついているのをどうにかできないか、 毎回トリミングしないといけないので、 という質問がありました (128)。
それに対して、gnuplot で png 画像にしたら、 という意見がありましたが、 set term png して set output して MS-Windows に転送して、 と結構手間がかかる、という話でした。 それに対しても、作業ディレクトリを Samba にしたら、 という回答がきていました。
一連の作業の手順が同じであれば、 スクリプトを使うとか bind を使うとかで解消できるような気もします。
例えば、上 にあげた 95 番の記事のような bind を利用して、 「set term png ; set out 'hoge.png' ; replot ; set out」 をなんらかのキーに割り当てて、そのキーで png にし、 WinSCP のような Unix と MS-Windows 間の GUI の scp ソフトを使って コピーする (その後ファイル名を適当に変更する)、 などとすれば多少手間は減らせます。
または、Xming はよく知りませんが、 単純に Unix 上であれば、ウィンドウマネージャの機能を使って gnuplot のグラフのフレームを消すことはできなくはありません (ただし使っているウィンドウマネージャに依存します)。 例えば、私は普段 fvwm (2.4.X) を使っていますが、 fvwm の初期化ファイルで 「Style "noframe" NoTitle, NoHandles, BorderWidth 0」 のようにしておいて、 gnuplot 側では、「set term x11 title "noframe"」 としてグラフを描かせると、 gnuplot のグラフウィンドウの枠が完全に消えます。
Unix 上でウィンドウを画像にキャプチャする際は、 xwd とか ImageMagick の import などを使いますが、 import はデフォルトでは枠はつけません (逆に -frame で枠つきで画像化する)。 そういうものを Unix 上で利用する、という手もあるかもしれません。
-
wxt, qt は X11 に比べて何が良くなるのか、 という質問がありました (133)。 Solaris の古い環境で ライブラリなどのインストールが大変そうだから、 使う価値があればやってみようかと思う、という話のようです。
それに対して、自分で試したら、 マウスホイールを使って範囲指定ができるけど 最近の gnuplot なら wxt 以外でもできるのでは、 最近の Linux はデフォルトが wxt で x11 にするとかなり見劣りする、 MS-Windows でまずは試してみたら、 といったような回答がありました。
wxt terminal 上で特別にできること、 というのは実はそんなにあるわけではありません。 むしろ、pango や cairo などのライブラリにより グラフのレンダリング、文字のレンダリングが綺麗になるところが 大きな特徴でしょうか。 フォントも fontconfig 管理になるので x11 に比べれば指定するフォント名などのフォント管理は 多分だいぶ楽でしょうし、 UTF-8 文字列なども x11 terminal に比べれば 適切に処理されやすくなっているだろうと思います。
なお、そこから派生して、いつまで古い Solaris を使っているんだ、 といったような意見もでていたようですが、 実はうちでは Solaris 9 のマシン (Ultra10) が健在で、 gnuplot に関する作業のほとんどをその上でやっていたりします (^^)
その Solaris に wxt, qt も入れてありますが、 Solaris 9 だとデフォルトの gtk とかのライブラリが古いので、 確かに wxt のインストールはそれなりに苦労しました ( 「Unix に関するメモ等 (10/10 2006)」 参照)。 ただ、特に wxt でないと困る、ということはないので、 普段は私は x11 terminal を使用しています。 -
[0:2*pi] で sin(x), [-2*pi:0] で cos(x), [-2*pi:2*pi] で tan(x) を描画するにはどうしたら、 という質問がありました (140)。
FAQ の一つのような気がしますが、 3 項演算子を使って不要な部分を 1/0 で切り捨てる、 という標準的な回答がついていましたが、 「情報やメモ (12/21 2012)」 で紹介した通り、 最近の CVS 版ならばグラフに個別に範囲指定をすることでも 行うことができます。
set yrange [-2:2]
set grid
plot [-2*pi:2*pi] tan(x) lt 3, [0:2*pi] sin(x) lt 1, [-2*pi:0] cos(x) lt 2- 上のスクリプトの画像
- 上のスクリプトの画像
-
pm3d map 上に x=1 のような直線を描きたいが、 nohead の arrow だとカラーマップの下に入って見えなくなるし、 媒介変数表示だと pm3d map では表示されない、 という質問がありました (144)。
それに対して、描画の順番を入れかえたら、という意見がありましたが、 set arrow なので splot より前に指定せざるをえない、 と返ってきていました。 そして、直線は外部ファイルでデータを作って、 それを描画させる、という強引な解決策を取ったそうです。
まず、媒介変数表示で 「set parametric ; set pm3d map ; splot 1,u,0」のようにしても 表示されないのは、多分 set pm3d map で塗る面がないからでしょう。
また set arrow だと上書きされる、というのには解決法があります。 gnuplot-4.0 以降では set arrow に front, back のオプションがあり、 front を指定するとグラフの最上位に表示されますから消されなくなります:
set pm3d map
set arrow 1 from 1,0,0 to 1,10,0 nohead front lt -1
splot sin(x+y)- 上のスクリプトの画像
- 上のスクリプトの画像