google-site-verification: google3bd66dd162ef54c7.html
FC2ブログ

台風19号の通過に伴う気圧変化

大型の台風19号が近くを通過しそうだったので、気圧の変化の記録を取ってみました。

やり方は以前と同じで、LPS-25Hを使ったアナログ表示の気圧計から10分間隔でシリアルに流しているデーターを、TeraTeram を使ってタイムスタンプ付きでパソコンに記録しました。データーがPCに記録出来ればこっちのもの、エクセルで簡単にグラフにすることが出来ます。

▼気圧の値をPCに記録中
記録開始

▼台風19号の気圧変化
2019年台風19号の気圧変化@稲城市
規模が大きい台風だったので、気圧が下がってから元に戻るまで1日半くらいかかっています。

記録されていた最低の気圧は965hPaで、私がこれまでに観察した中で一番低い値でした。なお、気圧の値は海面気圧に換算した値です。ちなみに、台風の経路図を見るとほぼ中心、つまり目が通過していて、目の通過中は風が収まっていました。

▼気圧計の針がスケールオーバー
最低気圧付近
この気圧計の最低側は980hPaまでしか目盛りが無いのですが、それを超えて(下回って)針が振れました。まあメカニカルなストッパーがあるわけではないので、機械的には問題ありません。

◆まとめ
台風の中心が通過するなんて、めったにないこと、しかもかなり大型の台風です。その気圧変化をきっちりと記録出来てよかったです。そうは言っても、被害にあった方がいらっしゃるので、よかったと書くのは不謹慎ではあります。

気圧変化のグラフは滑らかな変化では無く、変曲点があったりコブのような不連続な変化があるようです。以前記録したデーターでもそんな感じでした。これって、ひょっとしたら台風の気圧分布には螺旋形の溝みたいな構造があるのかも知れません。銀河の渦には腕のような構造がありますが、それと同じようなことがもし台風の渦にもあったら、興味深い話です。

日の出日の入り時刻で描く太陽のアナレンマ (RStudioでグラフ化)

久しぶりにアナレンマの話です。この観察を始めてからもうすぐ3年が経過するのですが、今でもしぶとく続けています。というか、何もしなくても勝手にデーターが ambient に蓄積されるので手間がかかるわけではありません。

なお、アナレンマの観測の話はアナレンマのカテゴリにまとめてあるので興味のある方はご覧ください。

▼アナレンマの観察ユニット
アナレンマ観察装置
太陽光発電パネルの電圧をこのユニットで測定し、日の出日の入り時刻を割り出しています。得られた情報はWiFi経由で ambient に送っているだけなので、とても簡単な仕掛けになっています。

さて、突然アナレンマの話が出てきた理由ですが、最後にグラフを公開してからかなり時間が経ったということが第一にあります。あと、最近パソコンの入れ替えをやったのでいろいろアプリを入れ直しているのですが、その中に RStudio があり、久しぶりにこれを使ってみたかった、というのがもう一つの動機です。

▼最後に作成したアナレンマのグラフ
以前のアナレンマプロット

この図を作ったのは2017年の12月
で、エクセルの散布図でプロットしています。このグラフでは月別にプロットの色を変えることで時間の経過を表現していますが、本当はもっと滑らかに色を変化させたかったです。しかし、エクセルではこれくらいの表現が限界だと思います。そういうことはVBAやPythonなどでプログラムを書けば確実に出来るはずですが、環境の立ち上げが面倒なのでそのままになっていました。

前述したように、今回新しいPCにもR (RStudio) をインストールしたので、これを使って最新の観察データーでアナレンマの作図を行ってみました。R のグラフの表現力は優れているので、欲しかったグラフが簡単に作れることを期待しました。

以下おさらいを兼ねて、データーのダウンロードから順に作業を振り返ります。

▼Ambient で公開中のグラフ
ambientのアナレンマのページ
開始からもうじき3年。天候の影響でノイズが載っていますが、太陽の南中時刻の変化がはっきりと捉えられています。

データーは ambient のサーバーに自動アップロードしているので、そこからダウンロードします。なお、このデーターは、ambient の中に私が作った公開チャンネル 日の出・日の入り時刻でアナレンマを描く (チャネルID: 200)に入っているので、グラフを見るだけならて誰でも閲覧出来ます。

昔は一括でデーターのダウンロードができたのですが、さすがにサーバーへの負荷が重過ぎるのでしょう。現在は一度にダウンロードできるのは最大35日に制限されています。ローカルに過去分のデーターを保存しておき、新しく追加された分だけダウンロードして全データーを準備します。

現時点で約1300点のデーターがありますが、これをエクセルで前加工してCSVのファイルに変換し RStudioに読み込みます。

▼RStudio
Rstudioでグラフ化

RStudioから散布図作成コマンドの plot で作図しました。そのコマンドはたったこれだけで済みました。(不等号>は全角に修正しています)

> data = read.csv("c:/Myfiles/Rdata/anaremaSorted.csv")
> par(bg=rgb(0.8, 0.8, 0.8))
> plot(data$MidTime, data$DayLength, col = hsv(data$dayN/365, 1, 1), pch=20, tck=1)

最初の行でデーターをファイルから読み取り、2行目でグラフの背景を灰色に指定、3行目で散布図を書いています。ここですばらしいのは、要素の塗分けを hsv( , , )パラメーターで色相を直接指定出来ることです。色相は0から1の範囲の値で指定するので、年間の日付連番を365で割るだけで、1年周期の色環でプロット点を塗り分けることが出来ました。

▼アナレンマ (2019年10月10日時点)
最新のアナレンマのプロット

以前作ったものに比べ、データー数が増えたと共に、色の変化が滑らかになって格好良くなりました。

◆まとめ
久しぶりにRでグラフを作ってみたのですが、思いのほかうまくいきました。世界中で多くの人が使っているので、コマンドのオプションが豊富で、欲しいと思った機能はほとんどそろっている感じです。VBAでやらなかったのは正解でした。
カレンダー
09 | 2019/10 | 11
- - 1 2 3 4 5
6 7 8 9 10 11 12
13 14 15 16 17 18 19
20 21 22 23 24 25 26
27 28 29 30 31 - -
プロフィール

ラジオペンチ

Author:ラジオペンチ
電子工作を中心としたブログです。たまに近所(東京都稲城市)の話題など。60過ぎて視力や器用さの衰えを感じつつ日々挑戦!
コメントを入れる時にメールアドレスの記入は不要です。なお、非公開コメントは受け付けていません。

記事が気に入ったらクリックを!
最新記事
カテゴリ
最新コメント
リンク
FC2カウンター
検索フォーム
月別アーカイブ
RSSリンクの表示
QRコード
QRコード