用 GNUPlot 畫 surface plot

這次因為要處理 SVOM 接收站的背景測試，用了相當克難的方式收集到了大約一個星期的無線電背景資料。因為要看各頻率在不同時間的雜訊強度變化，因此要用 dynamical spectrum。之前也否經嚐試著用 IDL 來畫，但是效果並不理想，因此才試著用 GNU plot 畫看看。

GNU plot 畫 surface 的方法有兩種，第一種是用 pm3d，這種情況下，每一行 data 之間需要有空行。需要打的指令如下：

set pm3d at s

set view map

unset surface

unset ztics

set palette rgbformulae 22,13,-31

set xlabel "Days from Jun 1, 2009"

set ylabel "frequency(MHz)"

set term postscript color

set output "dynamical.ps"

splot "dynamical.txt" using 1:2:3

另一種是不用 splot，只需要用 plot 即可。這種情況下，每一筆資料間不需要有空行，不過要先打

set style data image

接著用 plot 就可以了。不過這其中還有一些需要設定的東西，還在摸索中...

2009/08/17 更新：要設定 z 方向的範圍，可用

set zrange [0:40]