趣味で計算流砂水理 Computational Sediment Hydraulics for Fun Learning

数値計算とか河川工学とかプログラミングのことを書いています

MENU

pythonモジュール:rioxarrayをアップデートしたら結構変わっていたのでよく使うところをまとめました

R

このブログではおなじみのpythonモジュール:rioxarrayを最新版にアップデートすると、主要なメソッドが変わっていたのでまとめておきます。 rioxarrayを使っている主な記事は以下です。 衛星画像処理用スクリプト:備忘録 - 趣味で計算流砂水理 Computation…

数値解析の基礎:10分解説

R

個人的な勉強会で作った資料です。 詳しい方から見たらふざけるな的な内容ですが10分解説なので許して下さい。 本記事はGitHub、nbviewer、Colab、Marp(スライド形式)でも公開しています。 ※Colabの解説記事はこちら ※Marpの解説記事はこちら 数値解析のア…

備忘録:windowsでHEICをJPGに変換する。

R

iphoneでおなじみのHEICファイルですが、他ソフトでの使用のためJPGファイルに変化する機会は多いです。 いつもはGIMPなどを使っていましたが、 Windowsでは、Microsoft Storeの「Image Converter」で一括変換できることを知りました。 参照:【2023年最新】…

備忘録:hvplot 0.8.Xが不具合がでるので0.7.Xを使う(2023/7/21時点)

R

hvplotの最新版0.8.Xで地図系の図化を行うと不具合がでます。 当面はhvplot 0.7.Xを使うしかなさそうです。 condaの場合 conda install -y -c pyviz hvplot==0.7.3 pipの場合 pip install hvplot==0.7.3 こういうのはちゃんとGithubのissueに上げるべきなん…

備忘録:conda自体のアップデートのエラー回避

R

久しぶりにconda自体のアップデートを行った時に引っかかったので。 公式やエラー文には以下のように書いていますが何故か通りません。 conda update -n base conda -c defaults 参考サイトより以下のコマンドで上手くいきました。エラーは結構出ますが問題…

Pythonの基礎:10分解説

R

個人的な勉強会で作った資料です。10分でPythonを解説するという無茶な話です。 本記事はGitHub、nbviewer、Colab、Marp(スライド形式)でも公開しています。 ※Colabの解説記事はこちら ※Marpの解説記事はこちら インデント データ型 データ型:複合的な型 …

Google Colaboratoryの基礎

R

個人的な勉強会で作った資料です。 本記事はGitHub、nbviewer、Colab、Marp(スライド形式)でも公開しています。 ※Marpの解説記事はこちら Google Colaboratoryについて Colabの起動 Colabの実行 90分ルール・12時間ルール 90分ルール 12時間ルール Colabの…

備忘録:marpの参考サイトのリンク

R

最近marpをよく使うので参考サイトをまとめておきます。 marpって開発者は日本人なんですね。 marpとは(chatGPTによる解説) Marp for VS Codeとは(chatGPTによる解説) Marp for VS Codeの使い方 Marpの画像の設定方法 Marpのデザイン エラー:Marp for V…

前○時間平均の降雨データをmatplotlibのbarで正確に書く

R

本記事のデータはGitHub - computational-sediment-hyd/storageFunctionにあります。 github.com 観測された降雨データの多くは前○時間平均で整理されています。 前○時間平均とは、例えば、1時間間隔の観測データで1/1 12:00の観測値が、1/1 11:00~12:00の…

流量ハイドログラフはなぜガンマ関数で近似できるのだろう

R

流量ハイドログラフはなぜガンマ関数で近似できるのでしょうか。 勉強不足でよくわからないです。詳しい方教えて頂きたいです。 アンサーを書きました。 computational-sediment-hyd.hatenablog.jp 参考資料 角屋、流出解析手法 (その2) 単位図法とその進展 …

備忘録:QGIS ポリゴンをラインで分割する

R

年に5回くらい、QGISを使って、ポリゴンをラインで分割することがあるのですが、毎回ググっているのでメモを残しておきます。 QGIS 3.22.5-Białowieżaの場合、 「プロセッシングツールボックス→ベクタオーバーレイ→線で分割」 リンク

備忘録:Ubuntu、プロキシ環境でのssl認証の設定(ssl_verify):(おまけ:intel oneAPIをUbuntu 18.04にインストールする2023年1月時点)

R

最近会社のセキュリティがさらに厳しくなってきて、SSLサーバ証明書がない一部のサイトについて、 SSL認証を無効化する設定も必要になったのでまとめておきます。 (以前、pythonのcondaで引っ掛かったの同様です。(参照:備忘録:プロキシ環境下でanaconda…

メモ:d4PDFの5kmダウンスケーリングモデルSI-CAT DDS5TKについて

R

本業(河川屋)で使ったのでまとめておきます。 SI-CAT DDS5TKとは 解析の対象期間 解析の対象エリア データの座標系、地図投映法 河川計画への適用について 公開場所 SI-CAT DDS5TKとは 気候変動適応技術社会実装プログラム SI-CAT(Social Implementation …

雑談:メインPC:VAIO Zがとんだので、VAIO Zを批評します。⇒ 一般の方にはあまりおすすめしません。

R

年明け早々最悪です。厄年の洗礼かな。 VAIO Zが壊れた経緯 VAIO Zを購入した理由 壊れるまでの使用感 修理依頼から完了まで まとめVAIO Zの批評 VAIO Zが壊れた経緯 私のメインPC:VAIO Zが壊れました。(スペックはノートPCのスペック : VAIO Z - 趣味で計…

ダム建設が河道に与える影響:十勝川を例に

R

十勝川について調べていた時に気になったのでまとめておきました。 一応専門としている土砂水理の話なので少しマニアックです。 ダムによる河道への影響 ダム建設後の十勝川の変化 十勝川 札内川 当該分野の研究動向 おわりに 参考図書 十勝川についてのみを…

備忘録:「河道計画検討の手引き」の間違い-潜水橋による堰上げの影響

R

さらに詳しい記事も書きましたのでそちらもご覧下さい。 computational-sediment-hyd.hatenablog.jp 自明な内容ですが5年に1回くらい聞かれるので。 直リンク:http://www.jice.or.jp/cms/kokudo/pdf/tech/material/kadoukeikaku.pdf#page=93 リンク 上下流…

備忘録:プロキシ環境下でanacondaが急にネットワークに接続できなくなった要因

R

プロキシ環境下でanacondaが急にネットワークに接続できなくなりました。 色々調べると、プロキシ環境の変化によりssl_verifyがTrueだと接続できなくなっていたようです。Falseにすると通りました。 conda config --set ssl_verify False 確認すると、 conda…

備忘録:Wordで図とキャプションを分離させない

R

Wordを使い初めて15年くらい経ちますが初めて知りました。 図とキャプションを分離させないようにするには ・キャプションを後ろに置く場合 : 図を配直する段落用のスタイルを定義し、スタイル定義ダイアログの[段落]にある[改ベージと改行]タプで[次の…

備忘録:windowsでwifiの候補リストが表示されなくなったときの対処法

R

windowsのアップデート等(私の場合はネットワークのリセット)でwifiの候補リストが表示されなくなったときの対処法です。 考えられる理由としては、以下のserviseが停止しているためと思われます。 WLAN AutoConfig piyopiyocs.blog115.fc2.com Network Li…

pythonで地域メッシュごとのポリゴンデータからGeoTIFFを作成する:国土数値情報平年メッシュデータを例に

R

TL;DR 地域メッシュポリゴンデータからGeoTIFFを作る手順を解説 pythonのrioxarrayを使うと簡単 テストケースで作成した平均年間雨量のデータが興味深い。 結果だけ参照されたい方は「 国土数値情報:平均年間雨量メッシュデータ」まで読み飛ばして下さい。 …

備忘録:QGISのarea()と$areaの違い

R

たまに使うと忘れているので。 以下のサイトの一部をDeepLしただけです。 gis.stackexchange.com 楕円形と平面、どちらの面積も計算できます。 投影CRSで定義されたジオメトリの平面面積を知りたい場合は、area($geometry) 式を使用します。返される値は、そ…

議事メモ:2022/11/5

R

非構造格子の計算方法 よくわからないので勉強したい 。資料収集から。 貯留関数のパラメータ同定 スパース推定を使えないか。 浅水流モデルの高速化 ステンシル数が少なく並列化では限界がある。 dtを大きくとる、つまり、陰解法を用いるしかない。その場合…

画像処理のど素人がまとめたpythonによるTIFF圧縮形式の比較

R

画像処理なんて全くわからないのですが、必要になったのでまとめておきました。 TL;DR TIFF画像作成時には多くの圧縮形式が準備されている 今回の事例では、deflate、lzma、lzw、zstdによる圧縮でファイルサイズが30~50分の1になった。 はじめに TIFFの圧縮…

コーディングを意識したかなり丁寧なSIMPLE法:スタッガード格子版の離散化の解説

R

以前にSIMPLE法のプログラミングに関する以下の記事を書いたので離散化の記事を書いておきます。 SIMPLE法でもpythonならシンプルに書けるのか ⇒ かなりシンプルになるがやっぱり...(ソースコードあり) - 趣味で計算流砂水理 Computational Sediment Hydraul…

備忘録:python-openpyxlを使うときの注意-旧形式ファイル.xlsは使えません

R

pythonでexcelファイルを処理するときの最もポピュラーなライブラリopenpyxlを使用する場合、旧形式のexcelファイル(拡張子.xls)は使えないので、.xlsxまたは.xlsm形式に変換する必要があります。 大量のxlsファイルを変換する場合は、例えば以下のサイト…

備忘録:numpyでインデックス用の連番を振る

R

数値計算ではインプットファイルなどでインデックス用の連番を振ることよくあります。 テクニックのメモです。 配列を繰り返す 配列0~4を3回繰り返したインデックスの作成方法です。 import numpy as np x = range(5) print(list(x)) # [0, 1, 2, 3, 4] i =…

備忘録:geopandasのGeoDataFrameをKMLで出力する

R

参考サイト:Field names mixed up writing via "KML" driver のとおり、fionaを使って出力します。 ただし、GeoDataFrame内の日本語column名が上手く出力されないようです。 また、KMLはNameFieldが必須のため、設定する必要があります。これも日本語名は上…

備忘録:windows環境でintel fortran無償版の動作環境を作る:2022-9-8

R

windows環境でintel fortranを使う必要があり、仕方なく環境を構築しました。(Linux版はこちら) 基本的には以下の2記事をそのままトレースしているだけです。 Intel oneAPI Base Toolkit + HPC Toolkit (Intel Fortran)のインストール(Windows版) - Qiit…

備忘録:tensorflowのインストール:2022-9-8

R

急にtensorflowを使い始めました。 インストールは以下を参考にしました。 https://www5.dent.niigata-u.ac.jp/~nisiyama/grad/python/anaconda-install-j.pdf#page=11 condaで仮想環境を作ります。今回はGPU無とします。 conda create --name AI python=3.9…

備忘録:海外の洗掘(scour)深予測の簡易式について

R

洗掘深予測式は、護岸の根入れ等を河床下へ埋設する構造物を設計する際に必要となります。 洗掘は大きく分類すると、 一般的な洗掘:General Scour 構造物周り等の局所的な洗掘:Local Scour になります。 前者は現在の技術では河床変動計算による予測が概ね…