データを標準出力からそのままグラフ化したい

シェル芸に長けた仕事の仕方をしてる人は往々にしてそういう状況に巡り合いますよね？
私はあります。特に環境が整っていない組込系の仕事をしていると。
(それ以外の人は素直にmatplotlib/seabornでもしてください)

なので、標準入力からデータを受け取ってgnuplotからグラフ化までを行うbash(awk)を作りました。

環境

gnuplot 5.4 patchlevel 0
GNU Awk 5.1.0, API: 3.0 (GNU MPFR 4.1.0, GNU MP 6.2.0)
GNU bash, バージョン 4.4.12(3)-release (x86_64-unknown-cygwin)

グラフ表示するウィンドウはx11を指定していますが、環境によって適当に変更してもらえれば良いと思っています。

使い方

とりあえず表示

$ seq 10 | awk '{print $1,$1**2}'
1 1
2 4
3 9
4 16
5 25
6 36
7 49
8 64
9 81
10 100
$ seq 10 | awk '{print $1,$1**2}' | gplt

optionを使う

一応使えるoptionはhelpで表示されるようにしていますが、メモ書き程度なので、なんとなく使ってみてoption内容を理解するようにしてください。
基本的にはoption指定の後ろに数字なりテキストなりを入れて、どんどん設定入力していくようなイメージになります。

$ gplt --help
[x|y|y2]range,[x|y|y2]r  -X  X  : range setting
[x|y|y2]label,[x|y|y2]l  X      : label text setting
[x|y|y2]tics,[x|y|y2]t  X       : scale setting
[x|y|y2]format,[x|y|y2]f  X     : axis decimal digit num. e.g. X=2->1.00
title,t  X                      : graph title
key,k  X                        : w lp key X. legend label
axis                            : , u X:X axis. y2axis
font  X                         : font size ratio(default 14)
out                             : png file output
file  X                         : png file name
size  X                         : png iamge pix size ratio(default:640,480)

例
pngファイルに出力する

$ seq 10 | gplt out
$ ls
001.png

fileでファイル名を指定しなければ001.pngからの連番で作成されていきます。
scriptと絡めて使うと便利です。

X軸のレンジを指定する

$ seq 10 | gplt xrange 0 5

複数データを表示し、第二軸に設定する

$ seq 10 | awk '{print $1,$1*2,$1^2}' | gplt , u 1:3 axis

よく使うのはこのあたりです。
最終的なoutputとしてグラフデータを使う場合はlabel/key/title設定とか整形したりします。

$ seq 20 | awk '{print $1,$1*2,$1^2,$1/2}' | gplt w lp pt 5 k data1 , u 1:3 axis k y2axis yr 0 50 y2r 0 500 , u 1:4 w p pt 7 ps 4 k data3 xl xaxis yl yaxis t test_graph

コード

内容のほとんど(引数のoption解析)がawkなのでawkでシンタックスハイライトしています。

長いので折りたたみ

gplt

#!/bin/bash
# gnuplot 4.0 upper
# awk gawk 4.0 upper

###############
### setting ###
###############
### X window font size (8,10{11},12{13},14,17,18{20},24{25})
terminal_font_size=14

### default graph size
def_height=640
def_width=480

### deault graph size ratio
graph_size=1.1

### default graph terminal on gnuplot 5.0 over
def_terminal="x11"


# gnuplot version check
gnuplot_ver=`gnuplot -V | awk '{print int($2)}'`
# gnuplot_ver=4  # ver4.0debug 
if [ $gnuplot_ver -lt 5 ]; then
    # 5未満なら組込PC上(Xwindow無し)の使用
    gnuplot4_flag=1
    terminal="png"
    gnuplotcmd="gnuplot"
else
    # 5以上ならローカルPCでXwindow前提の使用
    gnuplot4_flag=0
    terminal=$def_terminal
    gnuplotcmd="gnuplot -p"
fi

# for png font size setting
if [ $terminal = "png" ]; then
    term_font_set="giant"
else
    term_font_set="font\",$terminal_font_size\""
fi

filename='tmp_01.dat'
optionfile='tmp_02.gp'
waferfile='tmp_03.dat'
displayfile='tmp_04.png'
xlinefile='tmp_05.dat'
ylinefile='tmp_06.dat'

if [ $# -eq 0 ]; then

sed 's/,/ /g' | awk 'substr($1,1,1)!~/#|%/' > $filename

$gnuplotcmd <<EOF
set term $terminal size $def_height*$graph_size,$def_width*$graph_size $term_font_set
set out '$displayfile'
set nokey
set grid
plot '$filename' w lp pt 6
EOF

else

# help disp
if [ $1 = "--help" ] || [ $1 = "--h" ] || [ $1 = "-h" ] || [ $1 = "-help" ] || [ $1 = "help" ]; then
    echo "####################"
    echo "### graph option ###"
    echo "####################"
    echo "[x|y|y2|cb][range|r] (min) (max)    : range setting. e.g. xr 0 5 -> 0to5"
    echo "                     {abs (absmax)} : range abs setting. e.g. xr abs 100 -> -100to100"
    echo "[x|y|y2|cb][label|l]  X             : label text setting. e.g. xl \"test\""
    echo "[x|y|y2|cb][tics|t]  X              : scale setting. e.g. xt 2->2 4 6..."
    echo "[x|y|y2|cb][format|f]  X            : axis decimal digit num. e.g. xf 2->1.00"
    echo "[title|t]  X                        : graph title. If no file name is set and the file is to be output, the title will be used as the file name."
    echo "[key|k]  X                          : w lp key X. legend label"
    echo "[keypos|kp] [left|l][right|r]"
    echo "            [top|t][bottom|b]       : legend label position e.g. kp r b"
    echo "axis                                : using y2axis e.g. , u X:X axis"
    echo "log [x|y|y2|cb]                     : axis logscale"
    echo "font [(3)|(giant|g|2)|(large|l|1)|(medium|m|0)|"
    echo "      (small|s|-1)|(tiny|t|-2)|(-3)]: font size e.g. font l. e.g. font 2"
    echo "size  X                             : png iamge pix size ratio(default:640,480)"
    echo "all {no \"keystr\"}                   : col all data plot. X is key title number. XX is key title str."
    echo "    {[h|s][s|split]}                : vsplit arranges charts vertically. hsplit arrange charts horizon."
    echo "[x|y]line X                         : line plot"
    echo "wafer [2|4|6|8|12]                  : wafer circle [inch]"
    echo "bar {rot} {label}                   : bar graph. rot is xtics dir rot. label is data label."
    echo "bubble {scale X} {refoff}           : bubble chart. scale is ref chart scale size."
    echo "map {gray|viridis}                  : color map graph u 1:2:3."
    echo ""
    echo "#####################"
    echo "### output option ###"
    echo "#####################"
    echo "out                                 : png file output"
    echo "file  X                             : png file name"
    echo "svg                                 : for markdown"
    echo "sixel|term|terminal                 : terminal display"
    echo ""
    echo "####################"
    echo "### other option ###"
    echo "####################"
    echo "fit {poly X} {txt|text|dat}         : poly fitting. degree X.  txt is fitting log output \"fit.log\""
    echo "    {gauss {height} {peak} {sigm} {base}}"
    echo "                                    : gauss fitting. "
    echo "headder                             : "
    exit;
fi

sed 's/,/ /g' | awk 'substr($1,1,1)!~/#|%/' > $filename

# record number get
nf=`awk 'NR>1{print NF;exit}' $filename`

echo $@ |\
awk '\
    function adb(s, n1, n2, n3, n4, n5){return out = sprintf(s, n1, n2, n3, n4, n5) out}
    function ad(s, n1, n2, n3, n4, n5){return out = out sprintf(s, n1, n2, n3, n4, n5)}
    function adop(s, n1, n2, n3){return w = w sprintf(s, n1, n2, n3)}
    function number_judge(x){return x~/^[-+]*[0-9]+\.*[0-9]*$/}
    function natural_no_judge(x){return x~/^[+]*[1-9]+[0-9]*$/}
    function range_set(s,  min_tmp, max_tmp){
        min_tmp = max_tmp = ""
        if($(i+1)=="abs" && number_judge($(i+2)) && $(i+2)>0){
            # arg1 abs -> range
            min_tmp = -$(i+2)
            max_tmp = $(i+2)
        }else{
            if(number_judge($(i+1))){
                min_tmp = $(i+1)
                next_flag++
            }
            if(number_judge($(i+2))){
                max_tmp = $(i+2)
                next_flag++
            }
        }
        gsub("range", "", s)
        gsub("r", "", s)
        ad("set %srange [%s:%s];", s, min_tmp, max_tmp)
    }
    function tics_set(s){
        gsub("tics", "", s);
        gsub("t", "", s)
        ad("set %stics %s;", s, $(i+1))
    }
    function label_set(s){
        gsub("label", "", s)
        gsub("l", "", s)
        ad("set %slabel \"%s\"", s, $(i+1));
        if(!gnuplot4_flag)ad(" noenhanced")
        ad(";")
    }
    function format_set(s){
        gsub("format", "", s)
        gsub("f", "", s)
        return format[s] = $(i+1)
    }
    function maxminget(  rc, rcr, k, n, i){
        if(maxminget_flag)return
        i = 0
        while(getline rc < filename > 0){
            n = split(rc, rcr);
            if(!maxminget_flag)for(k=1; k<=nf; k++)mint[k] = maxt[k] = rcr[k]
            i++
            for(k=1; k<=nf; k++){
                data[i, k] = rcr[k];
                if(rcr[k]>maxt[k])maxt[k]=rcr[k]
                if(rcr[k]<mint[k])mint[k]=rcr[k]
            }
            maxminget_flag++
        }
        nr = maxminget_flag
        close(filename)
    }
    function wafer_setting(){
        if(out!~/xrange/ && out!~/yrange/)ad("set size square;")
        if(out!~/xrange/)ad("set xrange [-%d*1.1:%d*1.1];", wafer_size, wafer_size)
        if(out!~/yrange/)ad("set yrange [-%d*1.1:%d*1.1];", wafer_size, wafer_size)
        if(out!~/xlabel/)ad("set xlabel \"x[mm]\";")
        if(out!~/ylabel/)ad("set ylabel \"y[mm]\";")
        ad("wafer_size=%s;", wafer_size)
    }
    BEGIN{
        axis_list["x"]++
        axis_list["y"]++
        axis_list["y2"]++
        axis_list["cb"]++
        # TICS_LIST
        tics_list["t"]++
        tics_list["tics"]++
        for(i in axis_list)for(k in tics_list)TICS_LIST[i k]++
        # RANGE_LIST
        range_list["r"]++
        range_list["range"]++
        for(i in axis_list)for(k in range_list)RANGE_LIST[i k]++
        # LABEL_LIST
        label_list["l"]++
        label_list["label"]++
        for(i in axis_list)for(k in label_list)LABEL_LIST[i k]++
        # FORMAT_LIST
        fromat_list["f"]++
        fromat_list["format"]++
        for(i in axis_list)for(k in fromat_list)FORMAT_LIST[i k]++
    }
    {
        fname = "001"
        font_size = 1
        wafer_size = 150
        def_scale = 10
        ad("set nokey;")
        ad("set ytics nomirror;")
        for(i=1; i<=NF; i++){
            if(next_flag && next_flag--)continue
            if($i == "title" || $i == "t"){
                ad("set title \"%s\"",$(i+1));
                if(!gnuplot4_flag)ad(" noenhanced")
                ad(";")
                title_str = $(i+1)
                next_flag++
                continue
            }
            if($i == "file"){
                fname = $(i+1)
                fname_flag++
                next_flag++
                continue
            }
            if($i == "size" && number_judge($(i+1))){
                graph_size *= $(i+1)
                next_flag++
                continue
            }
            if($i in TICS_LIST){
                tics_set($i)
                next_flag++
                continue
            }
            if($i in RANGE_LIST){
                range_set($i)
                continue
            }
            if($i in LABEL_LIST){
                label_set($i)
                next_flag++
                continue
            }
            if($i == "keypos" || $i == "kp"){
                gsub("set nokey;", "", out);
                trans_dir["r"] = trans_dir["right"] = "right"
                trans_dir["l"] = trans_dir["left"] = "left"
                trans_dir["t"] = trans_dir["top"] = "top"
                trans_dir["b"] = trans_dir["bottom"] = "bottom"
                i1 = $(i+1)
                i2 = $(i+2)
                if(i1 in trans_dir){
                    i1 = trans_dir[i1]
                    next_flag++
                    ad("set key %s", i1)
                    if(i2 in trans_dir){
                        i2 = trans_dir[i2]
                        ad(" %s", i2)
                        next_flag++
                    }
                    ad(";")
                }else if(i1 == "off"){
                    ad("unset key;")
                    next_flag++
                }
                continue
            }
            if($i == "key" || $i == "k"){
                gsub("set nokey;", "", out);
                adop(" title \"%s\" ", $(i+1))
                next_flag++
                continue
            }
            if($i == "font"){
                font_size = $(i+1)
                if(terminal != "png" &&
                   term_font_set = font_size == 3  ? "font\",24\"" :
                                   font_size == -3 ? "font\",8\"" : "" &&
                    ++next_flag)continue
                font_size_list["giant"] = font_size_list["g"] = font_size_list[2] = "giant"
                font_size_list["large"] = font_size_list["l"] = font_size_list[1] = "large"
                font_size_list["medium"] = font_size_list["m"] = font_size_list[0] = "medium"
                font_size_list["small"] = font_size_list["s"] = font_size_list[-1] = "small"
                font_size_list["tiny"] = font_size_list["t"] = font_size_list[-2] = "tiny"
                font_size in font_size_list &&
                (term_font_set = font_size_list[font_size]) &&
                terminal != "png" &&
                term_font_set = term_font_set == "giant" ? "font\",18\"" :
                                term_font_set == "large" ? "font\",17\"" :
                                term_font_set == "medium" ? "font\",14\"" :
                                term_font_set == "small" ? "font\",12\"" :
                                term_font_set == "tiny" ? "font\",10\"" : term_font_set;
            }
            if($i in FORMAT_LIST && format_set($i) && ++next_flag)continue
            if($i == ","){
                gsub("set nokey;", "", out);
                N = split(w, T)
                t_flag = w_flag = 0;
                for(k=N; k>=1; k--){
                    if(T[k] == "title")t_flag++
                    if(T[k] == "w")w_flag++
                    if(T[k] == ",")break
                }
                keyn++
                if(!t_flag)adop(" title \"%d\" ", keyn);
                if(!w_flag)adop(" w lp pt 6 ")
                adop(" , \"\" ");
                continue
            }
            if($i == "axis" && adop(" axis x1y2 "))continue
            if($i == "log" && $(i+1) in axis_list && ad("set logscale %s;", $(i+1)))continue
            if($i == "u" && adop(" u %s ", $(i+1)) && ++next_flag)continue
            if($i == "w"){
                i1 = $(i+1)
                if     (i1 == "lp")adop(" w lp pt 6 ");
                else if(i1 == "p") adop(" w p pt 5 ");
                else if(i1 == "l") adop(" w l ");
                continue
            }
            if($i == "pt"){
                N = split(w, T)
                flag = 0
                for(k=N; k>1; k--){
                    if(T[k] == "pt" && natural_no_judge($(i+1)) && (T[k+1] = $(i+1)) && ++flag)break
                    if(T[k] == ",")break
                }
                if(flag){
                    w = ""
                    for(k=1; k<=N; k++)adop(sprintf(" %s",T[k]));
                }
                continue
            }
            if($i == "ps" || $i == "lt" || $i == "lw"){
                N = split(w, T)
                flag = 0
                for(k=N; k>1; k--){
                    if((T[k] == "p" || T[k] == "lp" || T[k] == "l") && ++flag)break
                    if(T[k] == ",")break
                }
                if(flag)adop(" %s %d ", $i, natural_no_judge($(i+1)) ? $(i+1) : 1)
                continue
            }
            # if gnuplot version > 5.0
            if(($i == "hasen" || $i == "hasenn" || $i == "dot") && !gnuplot4_flag){
                N = split(w,T);
                flag = 0
                for(k=N; k>1; k--){
                    if((T[k] == "lp" || T[k] == "l") && ++flag)break
                    if(T[k] == ",")break
                }
                if(!flag)adop(" w lp pt 6 ")
                adop(" dt \"-\" ")
                continue
            }
            if($i == "out" && ++out_flag)continue
            if($i == "all"){
                all_flag++
                gsub("set nokey;", "", out);
                k = i+1
                while(k < NF){
                    if(!natural_no_judge($k))break
                    set_all_key_t[$k] = $(k+1)
                    k+=2
                }
                continue
            }
            if(all_flag){
                if(($i == "vsplit" || $i == "vs") && !hsplit_flag && ++vsplit_flag)continue
                if(($i == "hsplit" || $i == "hs") && !vsplit_flag && ++hsplit_flag)continue
            }
            # 複数列あってもいいかもしれないが一応allじゃないときに有効にする
            if($i == "fit" && !all_flag && nf >= 2){
                fit_flag++
                # とりあえず20個くらいパラメータを初期化 ⇒ gauss以外はいらないかも?
                for(k=1; k<=20; k++)ad("p%d=1;", k);
                if($(i+1) == "gauss"){ # 一次元gauss
                    # 1:Height, 2:PeakPos, 3:Sigm, 4:Base ←収束しないこともあるのでこれらのparam_initは個別指定可能とする.初期値はある程度求めたあげた方が良いな…
                    for(k=1; k<=4; k++){
                        if(!number_judge($(i+k+1)))continue
                        ad("p%d=%s;", k, $(i+k+1))
                        next_flag++
                    }
                    ad("f(x) = p1 * exp(-((x - p2) / p3)**2) + p4;")
                    ad("fit f(x) \"%s\" via p1, p2, p3, p4;", filename)
                    fx = "ttt*exp(-((x - ttt) / ttt)^2) + ttt"
                    param_no = 4
                    next_flag++
                }else if($(i+1) == "gauss2" && nf >= 3){ # 二次元gauss 表示は重ねるんじゃなくてsplitだな
                    # 1:Height, 2:PeakPos1, 3:Sigm1, 4:PeakPos2, 5:Sigm2, 6:Base
                    for(k=1; k<=6; k++){
                        if(!number_judge($(i+k+1)))continue
                        ad("p%d=%s;", k, $(i+k+1))
                        next_flag++
                    }
                    ad("f(x) = p1 * exp(-((x - p2) / p3)**2 - ((y - p4) / p5)**2) + p6;")
                    ad("fit f(x) \"%s\" via p1, p2, p3, p4, p5, p6;", filename)
                    fx = "ttt*exp(-((x - ttt) / ttt)^2 - ((y - ttt) / ttt)^2) + ttt"
                    param_no = 6
                    next_flag++
                }else if($(i+1) == "poly" && natural_no_judge($(i+2))){  # 多項式一変数
                    param_no = $(i+2) + 1
                    st1 = "f(x) = "
                    st2 = "fit f(x) \"%s\" via "
                    for(k=1; k<=param_no; k++){
                        st1 = st1 sprintf("p%d*x**%d", k, param_no-k)
                        st2 = st2 sprintf("p%d", k)
                        fx = fx sprintf("tttx^%d", param_no-k)
                        if(k != param_no){
                            fx = fx " + "
                            st1 = st1 " + "
                            st2 = st2 ", "
                        }
                    }
                    st1 = st1 ";"
                    st2 = st2 ";"
                    ad(st1)
                    ad(st2, filename)
                    next_flag+=2
                }else if($(i+1) == "poly2" && nf >= 3 && natural_no_judge($(i+2))){ # 多項式二変数 ⇒ これもアルゴ化したい
                    if($(i+2) == 2){ # 二次
                        ad("f(x) = p1*x**2 + p2*x*y + p3*y**2 + p4*x + p5*y + p6;")
                        ad("fit f(x) \"%s\" via p1, p2, p3, p4, p5, p6;", filename)
                        fx = "tttx^2 + tttxy + ttty^2 + tttx + ttty + ttt"
                        param_no = 6
                        next_flag+=2
                    }else{          # 一次(平面)
                        ad("f(x) = p1*x + p2*y + p3;")
                        ad("fit f(x) \"%s\" via p1, p2, p3;", filename)
                        fx = "tttx + ttty + ttt"
                        param_no = 3
                        next_flag++
                    }
                }else{
                    # 引数何もない時は一次fit固定
                    ad("f(x) = p1*x + p2;")
                    ad("fit f(x) \"%s\" via p1, p2;", filename)
                    fx = "tttx + ttt"
                    param_no = 2
                }
                # fitting parameterをグラフ上に表示する. 4.0はsprintfが使えないので5.0以上限定
                # 適当に4個ずつくらいで改行を入れてみる？ ⇒ 改行文字が上手く入れられないからlabelを追加していくしかなさそう
                if(!gnuplot4_flag && param_no){
                    st = "set label %d right at graph 0.95,0.05 sprintf(\"f(x) = %s\"" 
                    for(k=1; k<=param_no; k++){
                        st = st sprintf(", p%d", k)
                    }
                    st = st ");"
                    ad(st, ++label_cnt, fx)
                    gsub("ttt", "%g", out)
                }
                continue
            }
            # fit result text output
            # ファイルの有無でfit.logを残すか残さないかshで判断する
            if(fit_flag && ($i == "txt" || $i == "text" || $i == "dat"))print 1 > "fit_log_get"
            if($i == "svg" && ++svg_flag)continue;
            if(($i == "sixel" || $i == "term" || $i == "terminal") && ++sixel_flag)continue;
            if($i == "map" && ++cmap_flag)continue;
            if(cmap_flag){
                if     ($i == "gray")   gray_color_flag++
                else if($i == "viridis")viridis_color_flag++
            }
            if($i == "xline"){
                maxminget()
                printf"%s ",$(i+1) > xlinefile
                for(k=1; k<=nf; k++)printf"%s ", maxt[k] > xlinefile
                printf"\n%s ",$(i+1) > xlinefile
                for(k=1; k<=nf; k++)printf"%s ", mint[k] > xlinefile
                printf"\n\n" > xlinefile
                next_flag++
                xline_flag++
                continue
            }
            if($i == "bubble" && ++bubble_flag)continue
            if(bubble_flag){
                if($i == "scale" && number_judge($(i+1)) && scale = $(i+1) && ++next_flag)continue;
                if($i == "refoff" && ++refoff_flag)continue
            }
            if($i == "yline"){
                maxminget()
                printf"%s %s\n%s %s\n\n", nf==1?0:mint[1], $(i+1), nf==1?nr:maxt[1], $(i+1) > ylinefile
                next_flag++
                yline_flag++
                continue
            }
            # inch
            if($i == "wafer" && natural_no_judge($i+1)){
                wafer_size = $(i+1)
                if(!(wafer_size%2) && wafer_size<=12 && wafer_size>=2)wafer_flag++
                if(!wafer_flag)continue
                wafer_size *= 25.4 / 2
                wafer_plot_point = 1000
                pi = atan2(0, -0)
                for(j=1; j<=wafer_plot_point; j++){
                    p = j * 2 * pi / wafer_plot_point;
                    print wafer_size*cos(p), wafer_size*sin(p), j/wafer_plot_point > waferfile
                }
                continue;
            }
            if($i == "bar" && nf > 1){
                maxminget();
                bar_flag++;
                continue;
            }
            if(bar_flag){
                if($i == "rot")  rot_flag++;
                if($i == "label")label_flag++;
            }
            if($i == "hedder"){
                hedder_flag++
                maxminget()
                # 一行目はヘッダーとして読み込むので#でコメントアウト化
                cmd = sprintf("sed -i -e \"1 s/^/#/g\" %s", filename);
                system(cmd)
                close(cmd)
                for(k=1; k<=nf; k++)hedder_list[k] = data[1, k]
                for(r=2; r<=nr; r++)for(k=1; k<=nf; k++)data[r-1, k] = data[r, k]
                nr--
                continue
            }
        }

        for(i in format)adb("set format %s \"%%.%df\";", i, format[i])

        # terminal type
        adb("set grid;")
        set_width = int(def_width * graph_size)
        set_hight = int(def_height * graph_size)
        if(out_flag){
            # title & fname flag check
            if(!fname_flag && title_str)fname = title_str
            # exist file check
            exst = fname ".png"
            if(getline < exst > 0){
                close(exst);
                exst = sprintf("%s_t.png", fname)
                if(fname~/^[0-9]+$/){
                    while(getline < (exst = sprintf("%03d.png", ++fname)) > 0)close(exst)
                    close(exst)
                }
            }
            close(exst)
            adb("set out \"%s\";", exst)
            adb("set term png size %s, %s %s;", set_hight, set_width, term_font_set)
        }else if(svg_flag){
            adb("set term svg font \", %d\";", set_font)
        }else if(sixel_flag){
            adb("set term sixel size %s, %s font \", %d\";",set_hight, set_width, set_font)
        }else{
            adb("set out \"%s\";", displayfile)
            adb("set term %s size %s, %s %s %s;", terminal, set_hight, set_width, term_font_set, gnuplot4_flag?"":"noenhanced")
        }

        if(w!~/ w /){
            if(cmap_flag)adop(" w p pt 5 ps 1 ");
            else         adop(" w lp pt 6 ");
        }else{
            S = w;
            nw = gsub(" w ", "", S);
            S = w;
            nc = gsub(" , ", "", S);
            if(nw < nc + 1)adop(" w lp pt 6 ")
        }
        if(w~/axis/){
            ad("set y2tics;")
        }
        if(w~/,/){
            N = split(w, T)
            flag = 0;
            for(k=N; k>=1; k--){
                if(T[k] == "title" && ++flag)break
                if(T[k] == ",")break
            }
            if(!flag)adop(" title \"%d\" ", ++keyn);
        }

        # graph type select
        if(cmap_flag){
            ad("set palette defined ")
            if     (gray_color_flag)   ad("( 0 \"#000000\", 1 \"#252525\", 2 \"#525252\", 3 \"#737373\", 4 \"#969696\", 5 \"#BDBDBD\", 6 \"#D9D9D9\", 7 \"#F0F0F0\", 8 \"#FFFFFF\");")
            else if(viridis_color_flag)ad("( 0 \"#440154\", 1 \"#472c7a\", 2 \"#3b518b\", 3 \"#2c718e\", 4 \"#21908d\", 5 \"#27ad81\", 6 \"#5cc863\", 7 \"#aadc32\", 8 \"#fde725\");")
            else                       ad("( 0 \"#000090\", 1 \"#000fff\", 2 \"#0090ff\", 3 \"#0fffee\", 4 \"#90ff70\", 5 \"#ffee00\", 6 \"#ff7000\", 7 \"#ee0000\", 8 \"#7f0000\");")
            ad("set view map;")
            if(wafer_flag){
                wafer_setting()
                # どうせカラーはnmであることが多いのでdefaultはそうする
                if(out!~/cblabel/)ad("set cblabel [nm];")
            }
            ad("splot \"%s\" u 1:2:3 %s palette", filename, w)
            if(wafer_flag)ad(", \"%s\" w l lt -1 lw 2;", waferfile)
            ad(";")
        }else if(bubble_flag){
            # ここでdefaultのスケールを設定
            if(scale == "")scale = def_scale
            tmp_scale = scale
            scalestr = scale
            scale = 15 / scale
            ad("set multiplot;")
            # ver4.0はcirclesに対応していないので、psで対応するため、レンジを固定する
            if(wafer_flag || gnuplot4_flag){
                wafer_setting()
                ad("ref_pos=%s;", wafer_size*0.93)
                if(!refoff_flag && !gnuplot4_flag){
                    ad("ref_scale=%s;", 15)
                    ad("scalestr=\"%s\";", scalestr)
                    ad("plot \"%s\" u (1>0?ref_pos:0):(1>0?ref_pos:0):(1>0?ref_scale:0) every 1 w circles lc \"gray\" fs solid border -1, \"\" u (1>0?ref_pos:0):(1>0?ref_pos:0):(1>0?scalestr:0) every 1 w labels offset 0,0;", filename)
                }
                ad("set parametric;")
                ad("plot [0:2*pi] wafer_size*cos(t),wafer_size*sin(t) w l lt -1 lw 2;")
            }else{
                # 横縦軸のレンジを固定する
                maxminget();
                ad("set xrange [%s:%s];", mint[1]-(maxt[1]-mint[1])/10, maxt[1]+(maxt[1]-mint[1])/10);
                ad("set yrange [%s:%s];", mint[2]-(maxt[2]-mint[2])/10, maxt[2]+(maxt[2]-mint[2])/10);
            }
            # gnuplot ver4 correspondesce
            # 別にこっちに統一してもいいけど5.0verも残しておく
            if(gnuplot4_flag){
                maxminget();
                # gnuplotバージョンによって〇のポイントタイプが異なる。6が中抜き円、Ver4.0は7が中抜き円
                circle_type = 7
                scale /= 4
                # refの表示
                ad("set label %d \"%s\" center at first ref_pos,ref_pos;", ++label_cnt, scalestr)
                scale_plot = scale * tmp_scale
                ad("b=%s;", scale_plot<0 ? -scale_plot : scale_plot);
                ad("plot \"%s\" u (1>0?ref_pos:0):(1>0?ref_pos:0) every ::1::1 w lp pt %d ps b*2 lt %d lw 2 notitle;", filename, circle_type, scale_plot<0?1:2)
                for(i=1; i<=nr; i++){
                    ad("set label %d \"%s\" center at first %s,%s;", ++label_cnt, data[i,3], data[i,1], data[i,2])
                    scale_plot = data[i,3] * scale
                    ad("i=%d;", i-1);
                    ad("b=%s;", scale_plot<0 ? -scale_plot : scale_plot);
                    ad("plot \"%s\" u 1:2 every ::i::i w lp pt %d ps b*2 lt %d lw 2 notitle;", filename, circle_type, scale_plot<0?1:2) # ここで正負の色決めてます
                }
            }else{
                ad("plot \"%s\" u 1:2:($3<0?-$3*%f:0) w circles lc 4 fs solid border -1, \"\" u 1:2:3 w labels offset 0,0;", filename, scale)
                ad("plot \"%s\" u 1:2:($3>0?$3*%f:0) w circles lc \"light-blue\" fs solid border -1, \"\" u 1:2:3 w labels offset 0,0;", filename, scale)
            }
            ad("unset multiplot;")
        }else if(bar_flag){
            # ver4.0はxlabelに行数しか指定できない(xtic()が使えない) set xformat ""してから個別にラベルを配置するか…
            # なので横軸が1列目のデータではなく行数になる…⇒ 1列目を文字列変換すればなんとかなるか？
            ad("set style fill solid border;")
            ad("set boxwidth 0.925;")
            if(out!~/yrange/){
                maxall = maxt[2]
                if(all_flag)for(k=2; k<=nf; k++)if(maxt[k] > maxall)maxall = maxt[k]
                ad("set yrange [:%s];", maxall*1.1)
            }
            if(rot_flag)ad("set xtics rotate by -90;")

            if(all_flag){
                if(label_flag && gnuplot4_flag){
                    maxminget();
                    axis_cnt = 0
                    for(i=1; i<=nr; i++){
                        for(k=2; k<=nf; k++){
                            ad("set label %d \"%s\" center at first %s,%s;", ++label_cnt, data[i,k], axis_cnt++, data[i,k]+maxall*0.02)
                        }
                        axis_cnt++
                    }
                }
                ad("plot \"%s\" u ($0*%d):2", filename, nf)
                if(nf == 3 && !gnuplot4_flag)ad(":xtic(1)")
                ad(" t")
                ad(" \"%s\"", set_all_key_t[1]?set_all_key_t[1]:1)
                ad(" w boxes lw 1 lt 3")
                if(label_flag && !gnuplot4_flag)ad(", \"\" u ($0*%d):2:2 t \"\" w labels offset 0, 1", nf)
                for(i=2; i<nf; i++){
                    ad(", \"\" u ($0*%d+%d):%d", nf, i-1, i+1)
                    if(i == int(nf/2) && !gnuplot4_flag)ad(":xtic(1)")
                    ad(" t")
                    ad(" \"%s\"", set_all_key_t[i]?set_all_key_t[i]:i)
                    ad(" w boxes lw 1 lt %d", i+3)
                    if(label_flag && !gnuplot4_flag)ad(", \"\" u ($0*%d+%d):%d:%d t \"\" w labels offset 0, 1", nf, i-1, i+1, i+1)
                }
            }else{
                if(label_flag && gnuplot4_flag){
                    maxminget();
                    for(i=1; i<=nr; i++)ad("set label %d \"%s\" center at first %s,%s;", ++label_cnt, data[i,2], i-1, data[i,2]+maxt[2]*0.02)
                }
                ad("plot \"%s\" u 0:2", filename)
                if(!gnuplot4_flag)ad(":xtic(1)")
                ad(" w boxes lw 1 lt 3")
                if(label_flag && !gnuplot4_flag)ad(", \"\" u 0:2:2 w labels offset 0, 1")
            }
            if(yline_flag)ad(", \"%s\" u ($1-1):2 t \"\" w l lt -1 lw 3", ylinefile)
            ad(";")
        }else if(all_flag){
            if(vsplit_flag){
                # ad("set multiplot layout %d,1;", nf-1)  # ver4.0ではlayoutは存在しない ⇒ 存在する場合はそのまま使えば？
                ad("set multiplot;")
                ad("set format x \"\";")
                # ver4.0ではscreenは使えないのでベタ書きで設定する。
                # ad("set lmargin screen 0.1;")
                ad("set lmargin 10;") # 左だけ軸の数字がちゃんと入るようにこれくらい空ける
                # ad("set rmargin screen 0.9;") # y2軸使用しない場合は設定不要なはず
                graph_no = nf - 1
                ratio_set = (gnuplot4_flag?1.1:1.15) / graph_size # 係数は経験則
                all_size = 30
                for(i=1; i<nf; i++){
                    # ad("set tmargin screen %s;", 0.95-0.90/(nf-1)*(i-1))
                    # ad("set bmargin screen %s;", 0.95-0.90/(nf-1)*i+0.025)
                    ad("set tmargin %s;", (2.5 + all_size / graph_no * (i-1)) / ratio_set)  # 係数は経験則
                    ad("set bmargin %s;", (3 + all_size - all_size / graph_no * i) / ratio_set)
                    if(i == nf - 1){
                        ad("set format x \"ttt\";")
                        gsub("ttt", "%g", out);
                    }
                    ad("plot \"%s\" u 1:%d t \"%s\" %s", filename, i+1, hedder_list[i+1]?hedder_list[i+1]:set_all_key_t[i]?set_all_key_t[i]:i, w)
                    if(xline_flag)ad(", \"%s\" u 1:%d t \"\" w l lt -1 lw 3", xlinefile, i+2)
                    ad(";")
                }
            }else if(hsplit_flag){
                ad("set multiplot;")
                graph_no = nf - 1
                ratio_set = (gnuplot4_flag?1.1:1) / graph_size
                all_size = 77
                for(i=1; i<nf; i++){
                    # if(i>1)ad("set format y \"\";")
                    ad("set lmargin %s;", (5 + all_size / graph_no * (i-1)) / ratio_set)
                    ad("set rmargin %s;", (5 + all_size - all_size / graph_no * i) / ratio_set)
                    ad("plot \"%s\" u 1:%d t \"%s\" %s", filename, i+1, hedder_list[i+1]?hedder_list[i+1]:set_all_key_t[i]?set_all_key_t[i]:i, w)
                    if(yline_flag)ad(", \"%s\" u ($1-1):2 t \"\" w l lt -1 lw 3", ylinefile)
                    ad(";")
                }
            }else{
                ad("plot \"%s\" u 1:2 t", filename)
                ad(" \"%s\"", hedder_list[1]?hedder_list[1]:set_all_key_t[1]?set_all_key_t[1]:1)
                ad(" %s", w)
                if(xline_flag)ad(", \"%s\" u 1:3 t \"\" w l lt -1 lw 3", xlinefile)
                for(i=2; i<nf; i++){
                    ad(", \"%s\" u 1:%d t", filename, i+1)
                    ad(" \"%s\"", hedder_list[i]?hedder_list[i]:set_all_key_t[i]?set_all_key_t[i]:i)
                    ad(" %s", w)
                    if(xline_flag)ad(", \"%s\" u 1:%d t \"\" w l lt -1 lw 3", xlinefile, i+2)
                    if(yline_flag)ad(", \"%s\" t \"\" w l lt -1 lw 3", ylinefile)
                }
                ad(";")
            }
        }else{
            if(wafer_flag){
                wafer_setting()
                ad("set parametric;")
                ad("set multiplot;")
                ad("plot [0:2*pi] wafer_size*cos(t),wafer_size*sin(t) w l lt -1 lw 2;")
            }
            ad("plot \"%s\" %s", filename, w)
            if(fit_flag)ad(", f(x) w l")
            if(xline_flag){
                flag = 0
                for(i=1; i<=nf; i++){
                    if(!index(out, sprintf(":%d", i)))continue
                    ad(", \"%s\" u 1:%d t \"\" w l lt -1 lw 3", xlinefile, i+1)
                    flag++
                }
                if(!flag){
                    ad(", \"%s\"", xlinefile)
                    if(nf != 1)ad(" u 1:3")
                    ad(" t \"\" w l lt -1 lw 3")
                }
            }
            if(yline_flag)ad(", \"%s\" t \"\" w l lt -1 lw 3", ylinefile)
            ad(";")
        }
    }
    END{
        print out > optfile;
    }' optfile=$optionfile filename=$filename nf=$nf waferfile=$waferfile\
       terminal=$terminal displayfile=$displayfile xlinefile=$xlinefile\
       ylinefile=$ylinefile term_font_set=$term_font_set gnuplot4_flag=$gnuplot4_flag\
       graph_size=$graph_size def_height=$def_height def_width=$def_width

$gnuplotcmd <<EOF
load '$optionfile'
EOF

if [ $gnuplot4_flag -eq 0 ] && [ `echo $@ grep out | wc -l` -gt 0 ] && [ `whereis mogrify | awk 'NF>1' | wc -l` -gt 0 ]; then
    # mogrify -fuzz 10% -trim +repage $out_file
    mogrify -trim +repage $out_file
fi

fi

# png出力する場合にimagemagickで表示するようにする
if [ -f $displayfile ]; then
    if [ $terminal = "png" ] && [ `whereis display | awk 'NF>1' | wc -l` -gt 0 ]; then
        display $displayfile
    fi
    rm $displayfile
fi

# temp file delete
if [ -f $filename ]; then   rm $filename; fi
if [ -f $optionfile ]; then rm $optionfile; fi
if [ -f $waferfile ]; then  rm $waferfile; fi
if [ -f $xlinefile ]; then  rm $xlinefile; fi
if [ -f $ylinefile ]; then  rm $ylinefile; fi
if [ -f fit_log_get ]; then
    rm fit_log_get
elif [ -f fit.log ]; then
    rm fit.log
fi

おわりに

gnuplotとawkしかないような環境でログデータ分析を行う場合に、これひとつつっこんでおけば分析が捗るような気がします。