numpy で3次元配列を２次元配列にした場合の動作確認。
arr=np.array([[[1,2],[3,4],[5,6]],
[[10,20],[30,40],[50,60]],
[[100,200],[300,400],[500,600]],
[[1000,2000],[3000,4000],[5000,6000]]])
print(arr.shape)
(4, 3, 2)
print(arr.reshape([4,-1]))
[[ 1 2 3 4 5 6]
[ 10 20 30 40 50 60]
[ 100 200 300 400 500 600]
[1000 2000 3000 4000 5000 6000]]
転置すると、i,j,k が k,j,i になります。簡単。
今日は Geoenvironmental Sciences に関する webinar を見ました。
1,2回目はライブで見ていたのですが、十分に理解できるわけもなく、3回目は Youtube にしました。
個人的には久しぶりとなる、地質に関する話でした。休日の穏やかな時間で見ていたこともあり、リラックスして聞けました。仕事では地質にかかわるヒトたちが CIM によるモデリングに迫られつつも見ぬふりをしていますが、ぜひ動画（で紹介されているようなツール）を見てほしいですね。
紹介されていたツールは GemPy。
https://www.gempy.org/
補間の計算には Theano が利用されています。PyMC3と組み合わせることで、確率論的手法（モンテカルロシミュレーション、ベイズ推論など）をサポートしていることが特徴とのことでした（そういえばPyMC3は動かなかったですね）。https://phreeqc.blogspot.com/2020/11/mcmc.html
textベースの入力ですので EVS と同様です。地質の３次元可視化という点では GEORAMA の方が操作性は良いと思われます。が、PyMC3が動けば強力な統合ツールになりそうです。地質って、不確実性もセットで表示すべきですからね。
webinar で、変化の速さを感じたのがココ↓。
https://youtu.be/dfKPxjSge90?t=2430
出典はわかりませんが、概ね理解できますし、利用してきたライブラリ等です。偶然、個人で選択していたツールが他者からも評価・整理されていると、安心します。なんとか、変化にはついていけていたと。一方で、ギリギリだったと複雑な思いも。
今日は webinar を見ただけですが、文献も出ています。
https://gmd.copernicus.org/articles/12/1/2019/
一度、読んでみましょう。
連休中に、スペクトルの平滑化プログラムを書いておこうと思いつきました。
大崎本にサブルーチンが載っていますので、それを利用すれば完成するでしょうと高を括っていました。
が、本を持ち帰っていませんでした。
どうしましょうか、と検索してみると、大崎総合研究所から6月にソースコードが公開されていました。（正誤表もあります。）
http://www.ohsaki.co.jp/activity/download/index.html
昔の有名な図書ですので、プロは既にアレンジしたコードをお持ちです。が、私のようなアマにはありがたい。感謝です。
サブルーチンをダウンロードして、プログラムに組み込み、動かしてみました。が、正しく動きません。どこだ？と探しているうちに、連休が終わりました。
そういえば以前、倍精度に変更したサブルーチンをプロからいただいていました（今回も、連休中に頂きました）。それに組み替えると、改善。さらに、引数を変数に変えると正しく動きました。後者はともかく、前者はあるあるでしたね。感謝。
昨日、お礼がてら雑談していると、周波数によってバンド幅を変える方法もあると教えていただきました。バンド幅を決めて平滑化すると、高周波側と低周波側で見え方が変わってきます。それが改善されるのでしょう。良いですね。
といっても、それを使うのに妥当なスペクトルか判断できない知識レベルですので、今回は通常の平滑化で完成。バグ取りに数日かかりましたが、喰わせるデータや単位のチェックもあわせてできたので良しとしましょう。
https://doi.org/10.1016/j.enggeo.2021.106322
Large-scale field model tests of landslide dam breaching
天然ダムの越流決壊を実験で確認した文献でした。
そういえば、決壊した後の形状はよく見ますが、越流初期からのプログレッシブな形状を追って見たことがありません。再現計算では流量と最終形状で評価しますし、そもそも映像がないので確認できないというのが実情でしょう。
この文献を読んでいて気づいたのが、側壁の侵食形状。水位から上で側岸崩落が発生しています。途中はそうなのかもしれませんね。
側方侵食速度、崩落角度も合わせこみパラメータの一つです。このような実験が増えることで、また一つ、この分野の経験式の積み上げが増えていくのでしょう。
地下水の年代測定において、トリチウムの有用性が低下してから、ずいぶん経過しました。
https://phreeqc.blogspot.com/2017/04/blog-post_10.html
その間、他の同位体を利用していたため、動向のチェックを怠っていました。
久しぶりに実務で利用された方のグラフを見ると、壊変前の値でプロットされておらず、壊変後の値でプロットされえていました。ま、そういう方法もあるよね、と思いつつスルーしていました。が、先日、文献を読んでいて同様のグラフを見かけました。
https://doi.org/10.1016/j.apgeochem.2021.105070
Insight into watershed hydrodynamics using silica, sulfate, and tritium: Source aquifers and water age in a mountain river
出典はUSGS
https://doi.org/10.3133/sir20195090
Tritium as an indicator of modern, mixed, and premodern groundwater age
Scientific Investigations Report 2019-5090
年代を決定できなくても、mixed か premodern かがわかるだけでも価値があります。今後はこのような使い方が主流になるのでしょう。
レイリー減衰でα、βを決定する考え方は、道路橋示方書に書かれています。
考慮すべき固有モードに対して安全側になるよう周波数、減衰比（減衰定数）を定めるような流れです。単純な構造、均一の部材であれば動解まで必要ないのですが、多様な形状・材料で橋が構成されているため、計算が必要になると理解しています。
これ、指定する周波数によっては、α、βの値が大きく変化します。結果にどの程度影響するのでしょう？SoilPlus で単純なモデルを作って計算してみました。
ESRI Japan さんの旧サポートサイトが12月でクローズする模様。
新しいサポートサイト MyESRI へ乗り換えないといけないのですが、紆余曲折を経て、ESRIの営業さんにお願いすることになりました。
その過程で、ライセンスの統合を進めました。驚くことに、2本目以降の年間保守料が半額近くまで下がりました。本体は微々たるものの、エクステンションの値下げ幅が大きい。もっと早く気付けばよかった。
数百の xml から 数項目の情報を抜き出して表にしたい、そんなプログラムを作ってくれるプログラマーがいないか？と他部署の先輩が探されていました。
それなら、先日の CIM スクリプトをダウングレードするだけだなあと思い、手を上げました。
30分ほどで出来上がったのですが、問題はPython（を使えない先輩とその部下）。
なので、重いし遅いのですが、実行ファイルに変換することにしました。
https://phreeqc.blogspot.com/2017/12/python-exe.html
そのまま変換すると 300MB の exe が出来上がりました。
これはダメ、軽くならないかと調べてみると、ありました。
【悲報】PyInstallerさん、300MBのexeファイルを吐き出すようになる
https://lets-hack.tech/programming/languages/python/pyinstaller-issue/
なるほど、全部入りでしたか。それなら仮想環境を毎回作るしかなさそうです（面倒ですが）。
手順
Jupyter の file - save as で .py 保存。
仮想環境確認
conda info -e
仮想環境作成
conda create -n pyinstaller
conda activate pyinstaller
インストール
conda install pyinstaller
今回利用するライブラリのインストール
pip install ●●
コンパイル
pyinstaller 〇〇.py --onefile --clean
変換された exe は 30MB。確かに軽くなりましたが、それでも重い。しかも遅い。
やはり、皆さんに覚えてもらうしかないでしょう。いえ、Python でなくてもよいのですが、何かしら一つくらい言語を習得しておかないと困ると思うのですが。手を貸すからダメなのかな。うーん。
