#はじめに
nim関係で2度目の投稿です。
今回はファイルを読み込んで得られた数値をgnuplotとgnuplotインターフェース(https://github.com/dvolk/gnuplot.nim)を用いてグラフにしていきたいと思います。

nimでgnuplotを使うに至るまではこちらを参考にしてください(なるだろうか?)
https://qiita.com/YagiEL39/items/be050672b2a49757d6de

LD_curve.txt

	鉄筋４本梁				鉄筋２本梁	
	荷重(kN)	たわみ(mm)			荷重(kN)	たわみ(mm)
1	0	0		1	0	0
2	0.283	0.137		2	0.129	0.034
3	0.512	0.219		3	0.422	0.084
4	1.01	0.357		4	0.408	0.158
5	1.561	0.484		5	0.616	0.204
6	1.994	0.57		6	0.877	0.253
7	2.574	0.699		7	1.171	0.299
8	3.09	0.843		8	1.425	0.393
9	3.598	0.988		9	1.582	0.446
10	3.673	1.018		10	1.912	0.594
11	3.831	1.091		11	1.858	0.65
12	3.934	1.215		12	1.844	0.834
13	3.809	1.308		13	1.862	1.025
14	3.974	1.41		14	1.88	1.306
15	4.028	1.515		15	1.872	1.516
16	4.149	2.017		16	1.908	2.001
17	4.217	3.018		17	1.908	2.525
18	4.21	4.012		18	1.894	3.027
19	4.221	5.011		19	1.854	3.496
20	3.963	6.026		20	1.726	4.054
21	3.501	6.62

大学で鉄筋コンクリート試験体破壊(小さいやつだけど)実験をやったときの荷重とたわみのデータです。
今回はこれを使おうかなと思います。pythonをいじっていた時はmatplotlibを利用してこんな感じでやってました。

lord_deflection.py

import matplotlib.pyplot as plt

file=open("LD_curve.txt",'r',encoding="utf-8")
data=file.readlines()
file.close()
kajyu1,tawami1,kajyu2,tawami2=[],[],[],[]
for j,i in enumerate(data):
    if 1<j and j<22:
        i=i.strip("\n").split("\t")
        kajyu1.append(float(i[1]))
        tawami1.append(float(i[2]))
        kajyu2.append(float(i[5]))
        tawami2.append(float(i[6]))
    elif j==22:
        i=i.strip("\n").split("\t")
        kajyu1.append(float(i[1]))
        tawami1.append(float(i[2]))
fig=plt.figure()
ax=fig.add_subplot(1,1,1)
ax.set_title("Lord-Deflection Curve")
ax.set_xlabel("Deflection[mm]",size=14,weight="light")
ax.set_ylabel("Lord[kN]",size=14,weight="light")
ax.set_xticks([i/2 for i in range(14)])
ax.set_yticks([i/2 for i in range(10)])
ax.tick_params(direction = "inout", length = 5, colors = "blue")
ax.grid(which = "major", axis = "x", color = "blue", alpha = 0.8,
        linestyle = "--", linewidth = 1)
ax.grid(which = "major", axis = "y", color = "green", alpha = 0.8,
        linestyle = "--", linewidth = 1)
ax.plot(tawami1,kajyu1,label="4 Rebars")
ax.plot(tawami2,kajyu2,label="2 Rebars")
ax.legend(loc="best")
plt.savefig("result.jpg",dpi=300)
plt.show()

pythonだとこんな感じになりますね。

これを今回はnimでやってみようかなと思います。

#目次
1.　ファイル読み書きの練習
2.　最低限のグラフ描画まで
3.　gnuplotの練習

#1.ファイル読み書きの練習
まずtxtファイルの作成&書き込みをやってみましょう。

file_write.nim

block:
    var f:File=open("test.txt",Filemode.fmWrite)
    defer:
        close(f)
        echo "closed"
    f.writeLine "nimnimしてきたな。"
    f.writeLine "先生に訊いてみよう"
    f.writeLine "先生、nimしてもいいですか?"
    f.writeLine "ダメです。"
    f.writeLine "あああああああああああああ"

で、その結果。

test.txt

nimnimしてきたな。
先生に訊いてみよう
先生、nimしてもいいですか?
ダメです。
あああああああああああああ

次に読み込んでみましょう

file_read.nim

import os
block:
    var f:File=open("test.txt",Filemode.fmRead)
    defer:
        close(f)
        echo "closed"
    for i in 1..4:
        echo f.readLine()

出力はこんな感じになります。
"nimnimしてきたな。
先生に訊いてみよう
先生、nimしてもいいですか?
ダメです。
あああああああああああああ
closed"

ってなります。
じゃあこれを利用してLD_curve.txtの読み込み→そのまま表示までやってみましょう。

ld_curve.nim

block:
    var f:File=open("LD_curve.txt",Filemode.fmRead)
    defer:
        f.close()
        echo "closed"
    echo f.readAll()

でokですね(^ω^)b

そうそう、これは余談ですが、nimはgoみたいにdeferを利用することができます。block内でbreakとかしても最後に必ずdefer内が処理されます。

途中で抜けようが抜けまいが必ずやってほしい処理を記述できます。

deferについては私の理解がまだこの程度なので後で詳しく書いてみたいですね。

#2.最低限のグラフ描画まで
"必要な数値を取り出してseqに入れる→プロットする"の手順で最低限のグラフ描画を行います。

gnu_ld.nim

import strutils,gnuplot
var
    x1,y1:seq[float64]= @[]#鉄筋4本用
    x2,y2:seq[float64]= @[]#鉄筋2本用
block:
    var f:File=open("LD_curve.txt",Filemode.fmRead)
    defer:
        f.close()
        echo "closed"
    while f.endOfFile!=true:#Fileに残りの行がなくなるとtrueになるので注意
        var
            temp1:seq[string]= @[]
        temp1=f.readLine().split("\t")
        try:
            y1.add(temp1[1].parseFloat)
            x1.add(temp1[2].parseFloat)
            try:
                y2.add(temp1[5].parseFloat)
                x2.add(temp1[6].parseFloat)
            except:
                echo "failed"#pythonのpass的なのを探してます。
        except:
            echo "failed"

plot x1,y1,"4Rebars"#最後のstringは凡例
plot x2,y2,"2Rebars"

discard readChar stdin

で、こいつをターミナルで"nim c -r gnu_ld.nim"してやると...
gnuplotの方がこんな感じになります。

#3.gnuplotの練習
gnuplotは触り始めたばっかりなので私はド素人ですが、グラフを描くにあたっていくらか調整していきたいと思います。

gnuplotインターフェースにはcmdプロシージャが用意されていて、ちょっとスマートではないですが、cmd("gnuplotコマンド")でgnuplotコマンドを利用することができます。というかgnuplotインターフェースにはplotとそのスタイル(線とか点とか)をいじるためのプロシージャしか用意されていないので、グラフの詳細設定にはcmdを使うしかないようです。

###3-1.グラフタイトル

cmd "set title 'グラフタイトル'"

でグラフタイトルを付けられます。
###3-2.各軸ラベル

cmd "set xlabel 'x軸ラベル'" 
cmd "set ylabel 'y軸ラベル'"

###3-3.格子線

cmd "set grid"

###3-4.最終的なグラフ
さてさて、色々いじくった後のグラフを表示してみましょう。

gnu_ld2.nim

import strutils,gnuplot
var
    x1,y1:seq[float64]= @[]#鉄筋4本用
    x2,y2:seq[float64]= @[]#鉄筋2本用
block:
    var f:File=open("LD_curve.txt",Filemode.fmRead)
    defer:
        f.close()
        echo "closed"
    while f.endOfFile!=true:#Fileに残りの行がなくなるとtrueになるので注意
        var
            temp1:seq[string]= @[]
        temp1=f.readLine().split("\t")
        try:
            y1.add(temp1[1].parseFloat)
            x1.add(temp1[2].parseFloat)
            try:
                y2.add(temp1[5].parseFloat)
                x2.add(temp1[6].parseFloat)
            except:
                echo "failed"
        except:
            echo "failed"

cmd "set title 'Lord-Deflection Curve'"#グラフタイトル
cmd "set grid"#格子線
cmd "set xlabel 'Deflection[mm]'" # 各軸ラベル
cmd "set ylabel 'Lord[kN]'"
plot x1,y1,"4Rebars"
plot x2,y2,"2Rebars"

discard readChar stdin

で、結果のグラフがこちら!!

ちゃんと出てきましたねー

以上、gnuplotファイル読み込み→グラフ描画してみたでしたー