gnuplot でこのような濃淡グラフを書けないでしょうか？

質問

（x座標 y座標 値）という形式のデータで次のようなものがあります

1 1 32.4
1 2 23.4
..........
..........
1 50 65.2
2 1 23.1
2 2 91.9
..........
..........
2 50 73.2
..........
..........
50 49 26.8
50 50 76.8

これを、二次元(x,y)のセル上で
50×50の濃淡プロットを作るには
どうすれば良いでしょうか？

gnuplotでなくてもOKです。
最初はexcelでやろうとしたのですが、思いの外大変でした。

回答No.3

imoriimori · ベストアンサー

＃２の者です。
「Densityplotというのがwebで見つかったのですが、関数のプロットのよう（？）です。
質問のデータを濃淡グラフでプロットするにはどうしたらよいのでしょうか？」

DensityPlotはおっしゃるとおり関数のプロットです。リストデータの濃淡プロットにはListDensityPlotというコマンドを使います。

ここでちょっとがっかりさせてしまうかもしれませんが、ListDensityPlotが扱うデータは {ｘ座標、y座標、値｝のデータ列ではなく、値がずらずらと並んだリストです。
もし手近にmathematicaを試せる環境ならば次のようなことをやってみればわかります。
data = {{1, 2, 3}, {4, 5, 6}, {7, 8, 9}, {10, 11, 12}};
ListDensityPlot[data, Mesh -> False];
１から１２までの濃度値が横３×縦４の升目に反映された図になります。

ですから、ご質問の例では、方法は二つです。
その１：{ｘ座標、y座標、値｝で並んだリストデータを値だけのリストデータに形を変えてやらねばなりません。
値だけのリストデータをmathematicaに書き込めばもちろんいいんですが、こつこつタイプインするのは５０×５０ではちょっと無理ですね。
リストデータを外部のテキストファイルに書き込んであるのなら、それをまずImportコマンドで読み込んで、そのあとリストデータの構造を変えて値だけのリストデータにします（リストデータの編集法は長々しい説明になってしまいますのでご勘弁ください）。
外部ファイルを作るときに値だけの羅列データにして、ずらずらと一列のデータとしてImportし、そのあとPartitionコマンドで50に区切ってやるのも良いです。

その２：{ｘ座標、y座標、値｝のデータ列からInterpolationコマンドで補間関数を作り、その補間関数をDensityPlotで描く。

正直、mathematicaはなぜ{ｘ座標、y座標、値｝のデータ列をスムーズに扱ってくれないのかちょっとうらめしく思うこともあります。ご質問の例ではきれいな直交格子にデータが並んでいるのでまだ良いのですが、座標がばらばらの（つまり直交格子から見れば歯抜けのあるような）データだとややこしいことになってしまいます。

お礼

ありがとうございます！
周りに mathematica に詳しい人がいないので助かりました。
Partitionというので何とかなりそうな気がします。