我需要使用gnuplot创建一个带线和点的平行图。 不幸的是使用
set style line 1 lc rgb 'blue' lt 1 lw 2 pt 7 ps 0.75 pi 5
或在绘图命令中设置pt ps和pi
Gnuplot不插入任何一点。
我需要插入点来区分具有相同颜色的两条线。 它们必须是相同的颜色,因为它们属于同一个“家族”。因此,我不能简单地改变颜色:(
我该怎么做?
谢谢
答案 0 :(得分:3)
遗憾的是,Gnuplot并不支持这种“开箱即用”的功能,不过如果你需要一些快速的一次性生产情节,那么如何做到这一点有点肮脏。更具体地说,可以获得当前的Gnuplot源代码,以一种或多或少令人满意的方式引入此功能,并使用修改后的Gnuplot版本生成所需的图形。
为了说明这个想法,让我们假设我们想要并行绘制以下数据文件data.dat
仅包含三个“集合”:
1 2 3
2 1 1
3 3 3
受到默认Gnuplot演示parallel.dem
的启发,可能会尝试使用这样的脚本test.gpl
实现此目的:
set title "Parallel Axis Plot" font ",15"
set terminal pdf
set output 'test.pdf'
set border 0
unset key
set xrange [] noextend
unset ytics
# use x-axis tic positions to label the axes
set xtics 1 format "axis %g" scale 0,0
# turn on axis tics for the parallel axes
set for [i=1:3] paxis i tics
plot 'data.dat' using 1:2:3 with parallel
这里,所有的线都是相同的类型和颜色,这是相当混乱的。作为一个小型升级,让我们按如下方式修改输入数据文件:
1 2 3 5
2 1 1 6
3 3 3 7
并使用稍微不同的绘图命令
plot 'data.dat' using 1:2:3:4 with parallel lc var
这会产生
结果数字看起来好一些,但是我们说我们确实坚持要有积分的情节。只需添加pt 1
之类的内容即可。此外,它会产生一条警告消息:"test.gpl", line 17: warning: No pointtype specifier allowed, here
表明parallel
样式根本不喜欢任何点规范。
为了对此进行部分改进,让我们按以下步骤进行:
cd ${HOME}
#prepare a sand box directory
mkdir gpl
cd gpl
#download the latest version
wget -O gnuplot_latest.tgz http://sourceforge.net/projects/gnuplot/files/latest/download?source=files
tar -xzvf gnuplot_latest.tgz
cd gnuplot-5.0.0
./configure --prefix=${HOME}/gpl/local
make && make install
在这里,我们假设我们的系统已经拥有成功构建的所有必要先决条件。例如在Ubuntu上,这个最小的情况应该只需要libpango和libreadline的开发包。
现在,我们需要对${HOME}/gpl/gnuplot-5.0.0/src
中存储的一些源文件进行一些处理。即,我们必须
说服Gnuplot,密谋风格parallel
喜欢积分。为此,只需修改114
中的行gp_types.h
并替换
PARALLELPLOT = 32*PLOT_STYLE_BITS + PLOT_STYLE_HAS_LINE
与
PARALLELPLOT = 32*PLOT_STYLE_BITS + PLOT_STYLE_HAS_LINE + PLOT_STYLE_HAS_POINT
介绍了在parallel
样式中绘制点的能力。为此,必须在文件plot_parallel
中找到函数graphics.c
。默认情况下,它看起来像这样:
static void
plot_parallel(struct curve_points *plot)
{
int i, j;
int x0, y0, x1, y1;
for (i = 0; i < plot->p_count; i++) {
/* rgb variable - color read from data column */
check_for_variable_color(plot, &plot->varcolor[i]);
x0 = map_x(1.0);
y0 = AXIS_MAP(PARALLEL_AXES+0, plot->z_n[0][i]);
for (j = 1; j < plot->n_par_axes; j++) {
x1 = map_x((double)(j+1));
y1 = AXIS_MAP(PARALLEL_AXES+j, plot->z_n[j][i]);
draw_clip_line(x0, y0, x1, y1);
x0 = x1;
y0 = y1;
}
}
}
让我们一点点拉皮条:
static void
plot_parallel(struct curve_points *plot)
{
int i, j;
int x0, y0, x1, y1;
int point_type;
struct termentry *t = term;
for (i = 0; i < plot->p_count; i++) {
/* rgb variable - color read from data column */
check_for_variable_color(plot, &plot->varcolor[i]);
point_type = plot->varcolor?((int)plot->varcolor[i]-1):plot->lp_properties.p_type;
x0 = map_x(1.0);
y0 = AXIS_MAP(PARALLEL_AXES+0, plot->z_n[0][i]);
(*t->pointsize)(plot->lp_properties.p_size);
(*t->point)(x0, y0, point_type);
for (j = 1; j < plot->n_par_axes; j++) {
x1 = map_x((double)(j+1));
y1 = AXIS_MAP(PARALLEL_AXES+j, plot->z_n[j][i]);
draw_clip_line(x0, y0, x1, y1);
(*t->point)(x1, y1, point_type);
x0 = x1;
y0 = y1;
}
}
}
就是这样!现在使用make && make install
再次编译Gnuplot(从目录${HOME}/gpl/gnuplot-5.0.0
调用)。使用${HOME}/gpl/local/bin/gnuplot ${HOME}/gpl/test.gpl
生成的输出(使用绘图命令plot 'data.dat' using 1:2:3:4 with parallel lc var ps 1.5
)应该类似于:
gp_types.h
中的修改确保不会忽略最终绘图中使用的磅值规范ps 3
(并且可以在plot_parallel
函数中访问)plot_parallel
函数中,我们需要引入一个指向全局终端变量的指针,即struct termentry *t = term;
。然后将其用于实际绘制点。在for
循环中迭代输入数据文件中的各行,我们用
`point_type = plot->varcolor?((int)plot->varcolor[i]-1):plot->lp_properties.p_type;`
这确保了如果我们不使用lc var
,则点类型由pt
(或默认值)确定。如果lc var
处于活动状态,则点类型将从输入数据文件中的相应数据列中获取。
(*t->pointsize)(plot->lp_properties.p_size);
设置用ps
(或默认值)指定的所需磅值。(*t->point)(x0, y0, point_type);
的调用用点来增加单个线段。