使用 gnuplot 在标签中显示列值

Display a column value in label with gnuplot

我有一个程序可以生成一个文件,其中包含粒子在盒子中的位置,数据还包含盒子内的压力(比方说每秒)。使用 gnuplot,我在循环中为粒子设置动画,但我想在标签中显示压力。可能吗?

这是文件。压力在第 4 列

我的代码:

# Titre du graphique
set title "Mouvement de particules d'Hélium (R = 3e-2 {/Helvetica=10 5) \n sur une surface de      4 {/Helvetica=10 5^{2}}"
set xlabel "L {/Helvetica=10 (5)}"
set ylabel "L {/Helvetica=10 (5)}"
set style line 2 lc rgb 'blue' pt 7 # cercle
set style fill solid 0.4
unset key
set size square
infile = 'particules.dat'
#stats infile using 1:2 name 'data' nooutput
k=21
do for [i=0:k] {
labels(x,y) = stringcolumn(3)
stats [*:*] infile using 5 name 'pression' nooutput
set label 1 sprintf("Pression : ", pression) at screen 0.2,0.9
plot infile index i using 1:2 with circles ls 2
pause 1.0
}

不太优雅,但你可以使用多图。使用它作为循环的核心:

set multiplot

set size 1,1
set origin 0,0

set border
set xtics; set ytics

set xlabel "L {/Helvetica=10 (5)}"
set ylabel "L {/Helvetica=10 (5)}"
set title "Mouvement de particules d'Hélium (R = 3e-2 {/Helvetica=10 5) \n sur une surface de      4 {/Helvetica=10 5^{2}}"

plot infile index i using 1:2 with circles ls 2

set border 0
unset xtics; unset ytics
unset xlabel; unset ylabel
unset title

set origin -0.35,0.395
plot infile index i u (0):(0):4 w labels

unset multiplot

您可能需要使用当前的 origin 值。

正如@Tom Solid 最初建议的那样,您仍然可以通过绘图 with labels 来完成。 “技巧”是使用 offset。这将给出一个标签,例如向上偏移 1.2 个字符大小。显然,有了这个,标签甚至被放置得比图形区域大。检查 help labeloffset 选项,它也将在绘图样式 with labels 中工作,尽管它没有在文档中明确提及(参见 help labels)。 此外,如果将其限制为一行,即 every ::0::0,则可以避免将许多标签打印在彼此之上,检查 help every.

如果您想要动画 GIF 文件,可以使用带有选项 animateterminal gif,但是不幸的是图像质量不是很好。我建议创建单独的 PNG 框架,例如使用 term pngcairo 并使用其他软件将它们组合成动画 GIF。

代码:

### animation with label from data
reset session

# create some random test data
set print $Data
do for [p=0:10] {
    do for [i=1:50] {
        print sprintf("%.3f %.3f %.2f %.2f",rand(0)*2-1,rand(0)*2-1,0.03,p/10.)
    }
    print "\n"; print "\n"
}
set print

set size ratio -1
set xrange [-1:1]
set yrange [-1:1]
unset key
set style fill transparent solid 0.3
set title "Mouvement de particules d'Hélium (R = 3e-2 Å)\nsur une surface de 4 Å^2"
myLabel(col) = sprintf("Pression: %.2f",column(col))

do for [p=0:10] {
    plot $Data u 1:2:3 index p w circles ls 1, \
         ''    u (-1):(1):(myLabel(4)) index p every ::0::0 w labels offset 0,1.2
    pause 1.0
}

# additionally create an animated GIF file
set term gif size 400,400 animate delay 100 font ",10"
set output "myHeAtoms.gif"
do for [p=0:10] {
    plot $Data u 1:2:3 index p w circles ls 1, \
         ''    u (-1):(1):(myLabel(4)) index p every ::0::0 w labels offset 0,1.2
}
set output
### end of code

结果1:(wxt终端输出单个PNG)

结果 2: GIF 动画(好吧,图像质量不错,我不知道为什么这些点的大小明显不同):