Gnuplot在平行图中插入点

时间:2015-10-06 07:09:58

标签: gnuplot

我需要使用gnuplot创建一个带线和点的平行图。 不幸的是使用

set style line 1 lc rgb 'blue' lt 1 lw 2 pt 7 ps 0.75 pi 5

或在绘图命令中设置pt ps和pi

Gnuplot不插入任何一点。

我需要插入点来区分具有相同颜色的两条线。 它们必须是相同的颜色,因为它们属于同一个“家族”。因此,我不能简单地改变颜色:(

我该怎么做?

谢谢

1 个答案:

答案 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

然而,这会产生可能类似的情节: enter image description here

这里,所有的线都是相同的类型和颜色,这是相当混乱的。作为一个小型升级,让我们按如下方式修改输入数据文件:

1 2 3 5
2 1 1 6
3 3 3 7

并使用稍微不同的绘图命令

plot 'data.dat' using 1:2:3:4 with parallel lc var

这会产生

enter image description here

结果数字看起来好一些,但是我们说我们确实坚持要有积分的情节。只需添加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中存储的一些源文件进行一些处理。即,我们必须

  1. 说服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
    
  2. 介绍了在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;
        }
    
        }
    }
    
  3. 就是这样!现在使用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)应该类似于:

  4. enter image description here

    说明

    1. gp_types.h中的修改确保不会忽略最终绘图中使用的磅值规范ps 3(并且可以在plot_parallel函数中访问)
    2. plot_parallel函数中,我们需要引入一个指向全局终端变量的指针,即struct termentry *t = term;。然后将其用于实际绘制点。
    3. for循环中迭代输入数据文件中的各行,我们用

      确定相应的点类型
      `point_type = plot->varcolor?((int)plot->varcolor[i]-1):plot->lp_properties.p_type;`
      

      这确保了如果我们不使用lc var,则点类型由pt(或默认值)确定。如果lc var处于活动状态,则点类型将从输入数据文件中的相应数据列中获取。

    4. call (*t->pointsize)(plot->lp_properties.p_size);设置用ps(或默认值)指定的所需磅值。
    5. 最后,类型(*t->point)(x0, y0, point_type);的调用用点来增加单个线段。