このところのコロナ感染の急拡大が気になって、変異型のコロナの感染力が原因かもしれないが、以前行ったデータからの予測を試みた。大阪府の場合、3月の初めから徐々に増え、ここへ来て大幅増になっている。前週との比は1.3,　1.5,　2.3,　1.7の4週にわたって増え、平均すると1.7近い。そこで、この1.7の増加で、グラフを作成してみた。一方、前週比が0.8で抑えられた場合のグラフも同時に作ってみた。スタートの時点を800人とすると2週間後に2200人、520人となっている。
4週間後の予測に変えてみたら、青色6336人、赤色335人、
ちょっと恐ろしい数字になっている。このままだといわゆる

大渕の夏目さんから2日ほど前に桜情報をもらっていたので、中新井池の御衣黄桜を見てきました。天気が曇りだったので、ちょっと残念でしたが、緑色の花は、高貴で、新鮮でした。ほぼ満開を迎えているようでした。
あまり見る機会が無いので、印象に残りました。
一枚だけスマホ写真を載せてみます。

地元、夏目さんの案内で、大渕地区の歴史を訪ねてグループウォーキングに参加しました。2回目です。
今回のテーマは「百万遍と軽便鉄道」
コースは大渕地区農村改善センター→銭瓶堂→六所神社→軽便鉄道・野賀駅跡→野賀熊野神社→野賀天神社→中新井池→元の場所
でおおよそ8000歩歩きました。
軽便鉄道跡はとても歩き易くよく整備されていました。百万遍の堂はまだ念仏行事が続いていると言うことや、神社もそれぞれ地元の子どもさん達できれいに掃除されていることを聞き、伝統が息づいていることを感じつつ、楽しく会話しながら回りました。
写真は、地蔵菩薩、銭瓶堂

コロナがまだ思うようには収束せず、家にいることの多い毎日ですが、空を眺める楽しみはあります。しかし、最近はあまり見上げていません。

3月の惑星は火星以外あまり見やすくないようです。

スコープタウンからのメール情報です。
図版つきの案内が以下のサイトで見ることができるようです。

https://scopetown.jp/2021/02/5173/

おすすめの天文現象は以下の通りです。
以下はメールからの引用です。
-----------
・火星とM45プレアデス星団（すばる）
３月上旬　宵〜22時過ぎ、西の空で火星とM45プレアデス星団（すばる）が近づきます。一番近づくのは３月４日頃で、その前後数日が

ちょっと子どもの疑問に答えてみようかなと思った。
昔住んでいた居た古屋の雨戸に節穴が空いていて、寝ていると障子に見事に景色が映っていた。まさにピンホールカメラ。
だから窓でなくてもピンホールで外の景色は全部見える（暗いけど）。それは光がその穴を直進してくるからだろう。
象は穴の点を結んで、点対称になる。だから景色は上下左右が対称になっている。(逆さま）

では、窓の外の景色はどうか。基本は同じだろう。今度は窓が大きいので、障子には明るくまどがボーと映るが、景色は映らない。
それは、通過する射線は窓を通って、いくつもの線となる。だから

以前、気象庁のデータを元に気温や降水量のグラフを描いてみたが、気になって自分で降雨量を計算できないものか考えてみた。
多分、いろいろと間違いを含んでいると思います。そのまま信用しないで、ご自分でお確かめ下さい。

降水量の計算をする。水の蒸気圧を使う。Kはケルビン(絶対温度で、0℃＝273Kである、理科年表を見る。hPaはヘクトパスカル）
288Kで16.886hPa、298Kで31.403hPa,290Kでは19.192hPaと出ている（もちろん湿度100％の蒸気量）。15℃から17℃になると13.6％ほど水蒸気量が増すことになる。
状態方程式を使う、空気1.0ｍ3中の水蒸気量を計算してみる。式の中のｗの水分量でkg

東京都の感染者予測を一度グラフにしたので、静岡県もできるのかなと考えたが、静岡県の場合、何だか難しそうである。
理由は、クラスターがどのように起きるのかで数字が大きく振れる。地域にもよると言うので、本当に大雑把ながら、ＮＨＫのまとめグラフの動きを見て考えてみた。一か月単位で見ると、増えたのは11月後半50人前後、それが12月後半は35人ぐらいそして1月後半は60人越えくらい。
これは1週間の平均感染者データからだ。だから、このところは一か月で2倍の増え方、11月と12月では0.5倍というわけで、
一か月単位で2倍と0.5倍の増え方。当然2月の後半は120人前後と30人

今日は、強風でしたが、地元夏目さんの案内で安倍晴明塚などグループで歩きました。
赤石を積み上げた、晴明塚初めて知りました。晴明のすすぎ井戸跡、晴明硯井戸跡、大須賀歴史民俗資料館、中新井池（御衣こう、オニバスで有名）を巡るウォーキングでした。
民俗資料館は昔生活の場で活躍していたかまど、脱穀機、むしろ織り機など、、、

東京都の感染者がどのくらい増えるのか、心配になって、単純なモデルで予測してみようと考えた。緊急事態宣言がでるようだけど、もちろんこれによって数値は変わるものと思う。
専門家の使う実効再生産数と言う意味が正確には分からないので、ここでは、11月から12月一杯までの東京の一週間の平均を出して、その増加割合を計算してみた。その結果おおよそ1.2倍になっていたので、2割増と考えて、一週間単位の周期としてこれを複利計算して、増加を予測してみた。1割減になったときも参考にグラフ化してみた。つまりy=1.2のｘ乗で現在の一日の感染者が1000人としてスタート。70日後

今回は厚生省のデータでなく、NHKのデータを拾ってみた。こちらは各都道府県ごとの一連のデータ。直接NHKのURLにアクセスしてweb上ではダウンロードできるが、直接pythonプログラムからcsvを読み込もうとすると、エラーになる。この原因はURLの長さに問題があるようだ。
仕方が無いので練習のため、自分のHPのサイトにアップロードした。
これならデータ取得できる。NHKのデータは都道県別の一連データだ、全ての県が入っているcsvなので、取り出すのにやっかい。もっとうまい方法があるかもしれないけど、今回は7435から7788までで何とか静岡県のデーターは取り出せる。NHKのデータは毎日