in the light
gnuplot
PyXPlot은 gnuplot의 확장 interface로 그림을 그릴 수 있게 해주는 프로그램입니다. 내부적으로는 Python과 PyX package를 이용하죠. 그림은 eps파일로 만들어 gv로 보여줍니다. 위의 그림과 아래의 그림은 아래의 명령으로 각각 Gnuplot과 PyXPlot에서 그린 결과를 보여줍니다.
set xlabel ‘Time (s)’
set ylabel ‘Amplitude’
set xrange [0:10]
p sin(x)
pyxplot
그림 그려주는 속도는 Gnuplot이 훨씬 빠릅니다. 하지만, Gnuplot 결과물은 presentation이나 paper에 쓰기에 부적합하다는 생각이 듭니다. 경험상 간편한 presentation 준비를 위해서는 PyXPlot을 쓰는 것도 괜찮더군요. Paper를 위해서는 PyX와 Postscript를 다 쓰는 것이 좋고요^^.
Gnuplot은 리눅스에서 텍스트파일에 저장된 값을 빠르게 그림으로 그려주는 프로그램입니다.
다양한 기능을 가지고 있지만, 제 경우에는 주로 수치해석 후 결과 확인용으로 씁니다.
gnuplot으로 그림을 그릴 때에는 command line 상에서 ‘gnuplot’이라고 치고 들어가서 gnuplot 명령어들을 이용하여 그림을 그리고 ‘q’를 입력하여 빠져나옵니다.
그런데 간단히 결과를 확인해보기 위해서 gnuplot에 들어가서
p 'file1' w l,'file1' u 1:3 w l,'file1' u 1:4 w p
또는
set grid
set xrange[:10]
set log y
p 'file1' w l,'file2' w l,'file3' w l
과 같이 매번 치려니 귀찮다는 생각이 들었습니다.
그래서 이 script를 만들게 되었습니다. 처음에는 ruby로 만들어 사용하다가 안 쓰는 기능은 빼고 python script로 바꿨습니다.
이 script를 사용하면 command line 상에서 바로 gnuplot 명령어를 사용하여 그림을 그릴 수 있습니다. 사용 방법은 아래와 같습니다.
./gnuqp.py [options] filename1 [u 1:2] [w l], filename2 [u 1:2] [w l], filename3 ...
실행파일 이름, 몇 가지 setting 관련 옵션들, 이후에는 gnuplot의 plot 명령어를 입력합니다.
filename1
Empty filename[2,3,...] will be replaced by the filename1
두 번째 위치부터는 파일명을 생략하면 첫 번째 파일명으로 대체합니다. 하나의 파일에서 여러 column들을 그릴 때 편리합니다.
u 1:2 : columns you want to plot
w [lp..]: line style- line, point, dot or impulse ..etc (default: w l)
plot 명령어의 옵션들 중에는 using (columns)과 with (line style)만 지원합니다. 그 외의 명령은 잘 안 쓰더라고요^^.
위의 옵션을 주지 않았을 때 기본적으로 ”with line” 옵션으로 그립니다.
-p : do not run gnuplot. just print the gnuplot command
-c : no comma seperation - the arguments are filenames seperated with a blank- use with glob pattern
-l : set logscale y
-g : set grid
-x[:10] : set xrange [:10]
-y[1:5] : set yrange [1:5]
위의 옵션들은 gnuplot의 setting을 간편하게 하기 위해 만들었습니다.
-p 옵션을 붙이면 gnuplot의 명령어만 출력하고 그림은 안 그립니다.
-c 옵션을 붙이면 파일들을 기본 옵션(with line)으로 그립니다. 이 때 파일명들 사이의 “,”를 생략하고 파일명만 씁니다. command line상에서 glob pattern을 이용하여 여러 그림을 그릴 수 있도록 하기 위한 옵션입니다.
나머지 gnuplot setting들은 위의 설명으로 충분할 것이라 생각합니다.
./gnuqp.py -p "file1", u 1:3 w p
-> p 'file1' u 1:2 w l, 'file1' u 1:3 w p
./gnuqp.py -p -c file.00*
-> p 'file.0010' w l,'file.0020' w l,'file.0030' w l
Download: gnuqp.py
이 글은 스프링노트에서 작성되었습니다.
