使用 gnuplot 在标签中显示列值

问题描述

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

这是文件。压力在第 4 列

我的代码：

# Titre du graphique
set title "Mouvement de particules d'Hélium (R = 3e-2 {/Helvetica=10 \305) \n sur une surface de      4 {/Helvetica=10 \305^{2}}"
set xlabel "L {/Helvetica=10 (\305)}"
set ylabel "L {/Helvetica=10 (\305)}"
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
}

解决方法

正如@Tom Solid 最初建议的那样，您仍然可以通过绘制 with labels 来实现。 “技巧”是使用offset。这将给出一个标签，例如向上偏移 1.2 个字符大小。显然，标签的放置甚至超出了图形区域。检查 help label 和选项 offset，它也适用于绘图样式 with labels，尽管文档中没有明确提到它（请参阅 help labels）。 此外，如果将标签限制为一行，则可以避免将多个标签打印在一起，即 every ::0::0，勾选 help every。

如果您想要动画 GIF 文件，则可以使用带有选项 terminal gif 的 animate，但不幸的是，图像质量不是那么好。我建议创建单独的 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,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,1.2
}
set output
### end of code

结果 1：（wxt 终端的单个 PNG 输出）

结果 2： GIF 动画（不错，图像质量不错，我不知道为什么点的大小略有不同）：

,

不太优雅，但您可以使用多图。将此作为循环的核心：

set multiplot

set size 1,1
set origin 0,0

set border
set xtics; set ytics

set xlabel "L {/Helvetica=10 (\305)}"
set ylabel "L {/Helvetica=10 (\305)}"
set title "Mouvement de particules d'Hélium (R = 3e-2 {/Helvetica=10 \305) \n sur une surface de      4 {/Helvetica=10 \305^{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 值。