#1/#2/#3

꽃가루의 비산 상황을 재료로, 시계열 데이터의 연습을.

전제

5.2 시스템 확인 (-c 또는 사용)

–persist(-p) 옵션을 사용합니다.

-p 는 OS X 라면 그 나름대로 설정하지 않으면 사용할 수 없습니다.
-p가 설정되어 있지 않은 환경에서는 ~/.gnuplot 파일을 설정하고, 바꾸고 나서 읽으면 좋을지도( #1 의 마지막으로 접하고 있습니다).

.gnuplot

~/.gnuplot 에 이하의 내용을 기재.

set xdata time
set timefmt "%Y%m%d%H"
set format x "%Y/%m/%d"
set xtics rotate by 90  right

load 파일에 시간축 설정을 기재할 때의 주의

$ cat ./.gnuplot
set xdata time
set timefmt "%Y%m%d%H"
set format x "%Y/%m/%d"
set xtics rotate by 90  right
$ gnuplot -c <( echo 'load ".gnuplot" ; ...') ...

등으로 호출하는 방법은 CentOS에서는 잘 갔지만 내 OS X (brew)에서는 잘 가지 않았다 1

단, ~/.gnuplot 에 같은 기재를 해, load 커멘드를 실행하지 않는 경우에는 잘 간다.

데이터 준비

에서 꽃가루 비산수 데이터를 적당히 다운로드. 2019/02/01부터 03/10까지의 데이터였다.

나는 이바라키 현민이므로 이바라키 3 관측점 데이터를.

format.csv 읽으면 체재는 설명된다.

사용하기 쉽도록 측정국을 1파일씩 분할.

% perl -F, -lane 'open my $out, ">>", $F[0] . q{.txt} ; print $out join " ", $F[2] . $F[3], $F[10]' Data.csv 

측정국 코드는

50810200 쓰쿠바
50810100 미토
50820100 히타치

그렇다고.

우선, 플롯

$ gnuplot -p -c <( echo 'plot for [ i = 1:ARGC ] ARGV[i] u 1:2 w l' ) 50810200.txt 50810100.txt 50820100.txt

하나코 씨

흠. 겹쳐서 맛있어.

축척 같이 그대로, 어긋나고, 조금 체재를 정돈한다

$ gnuplot -p -c <( echo 'array Stat[3]=["tsukuba","mito","hitachi"]; \
set yrange [-0.5:3] ; \
set ytics ( Stat[1] 0, Stat[2] 1, Stat[3] 2 ); \
plot for [ i = 1:ARGC ] ARGV[i] \
    u 1:($2/9000 + 1 * ( i - 1 ) ) w l title Stat[i]\
' ) 50810200.txt 50810100.txt 50820100.txt

2/#1 로부터 새로운 일이라고 하면,

배열(Stat[3])

인덱스가 1부터 시작한다.

ytics (Y 축의 눈금)는 이름 Y 값의 쉼표로 구분 된 집합으로 덮어 쓸 수 있습니다. 물론 xtics도.

다만 mytics등이 효과가 없어져 버린 것 같다.

using 인수 내에서의 사칙 연산은 기본

이 경우는, 9000으로 나누어, 묘화를 어긋나기 위해서, index 분족하고 있다.

정도?

#2



ARG1은 아무런 문제가 없지만 ARGV [1]에 "ariel"이 들어갑니다. 다른 제어 변수 인 w.

'\' 개행을 넣어 시인하기 쉽게 하고 있다. ↩