本記事はGitHub、nbviewer、Colabでも公開しています。 ※Colabの解説記事はこちら 導入 Lasファイルを処理する場合、PDAL - Point Data Abstraction Libraryを使用することが多いです。PDALはC++で書かれたライブラリですが、Pythonバインディングも提供され…
反響の多かった過去記事：（未解決）測量の高解像度を意識して、浅水流一次元計算でdxをかなり小さくすると上手くいかない件について - 趣味で計算流砂水理 Computational Sediment Hydraulics for Fun Learningについて、スタッガード格子の離散化の問題点…
タンクモデルの話 以下の論文が興味深い www.jstage.jst.go.jp 各種課題 お互い粛々と進めていきましょう。 その他 最近の話題のuvがPythonパッケージ管理で使えそう。またuvxもかなり使えそう。 ブログについて基礎水理の部分を強化したい。
分析的なことを行っていますが、メインの話はWebベースでの図化方法（Python）に関してだったりします。 TL;DR 国土数値情報 平年値メッシュデータのGeoTIFF変換、グラフ化についてPythonコードを詳しく示しています。 2010,2020年の平均年降水量の比較を行…
Delft3D 4をたまに使用しますがライセンスの更新方法を忘れるので備忘録です。 Delft3D 4 (structured)とは 公式ページより Delft3D 4 (structured) は、オランダの研究機関 Deltares によって開発された、水環境の数値モデリングを行うための統合型ソフトウ…
3ヶ月ぶりの勉強会でしたので色々と。 論文紹介 土砂水理に関して、有用な文献のない2000年代以降を中心に整理する。 植生水理、植生消長が主になると思われる。 寄付 順次進める 各種課題 お互い粛々と進めていきましょう。
はじめに 使っている人がほとんどいないことで有名なPython in Excelについて（だって、Pythonが使える人はExcelなんて使わないでしょ）、諸事情で初めて使ったときにつまづいた表題について備忘録としてまとめておきます。 参考サイトも少ないので、誰かの…
computational-sediment-hyd.hatenablog.jp 上記過去記事のリバイスです。サンプルコードも追記したので参考にしていただけると良いかと思います。 潜水橋とは(Copilot) 特徴 利点 欠点 日本での例 佐田沈下橋（四万十川） 流れ橋（上津屋橋）（木津川） 「…
雑談ベースかつ表に出せない内容が多いので議事メモはなしで。 引き続き頑張りましょう。
お疲れ様でした。 事務連絡 Dropboxの運用 寄付の話 浅水流の陰解法（SIMPLE法） 実河川計算では精度とコストの関係からスタッガードスキームの方が有利。コーディングが面倒だけど。冨永先生のモデルが参考になる。 氾濫解析のような大半に流れが無い場合は…
社内環境等でプロキシが設定されており、 vscodeで「ユーザ設定」でプロキシ設定を行っても、スクレイピング等がプロキシを抜けれない場合があります。 その時の最終手段です。 手順1:コマンドプロンプトを開く 手順2: 以下のコマンドを実行する set HTTP_PR…
QGISでレイアウトを作りすぎて管理しきれなくなり、地図帳にまとめるときに作ったPythonスクリプトです。 いろいろ応用が効くので、今後QGISのPythonコンソールを使うときのために上げておきます。 Copilotに大分助けられました。 from qgis.core import Qgs…
本記事はGitHub、nbviewer、Colabでも公開しています。 ※Colabの解説記事はこちら モチベーション 前提条件 川幅～水深～無次元掃流力の関係 式の導出 計算手順 計算結果 川幅、水深、無次元掃流力と掃流砂量の関係 式の導出 計算結果 おまけ：MPM式（1.5乗…
Pythonファイル内で別のPythonファイルをimportする際、 同じフォルダ内やサブフォルダ内にある場合は、以下でimportできます。 # myscript.pyをimportする import myscript # サブフォルダの場合 import folder.myscript また、上記以外のフォルダのPython…
stackoverflow.com nbconvert.readthedocs.io jupyter notebookファイル（.ipynb）をpython scriptファイル（.py）に変換するコマンドラインは以下のとおりです。 jupyter nbconvert --to script [YOUR_NOTEBOOK].ipynb notebook内で実行する場合は!をつける…
概要 インストーラーのダウンロード インストール おまけ：WindowsのVScodeからWSLのPython環境を使用する 概要 WSLのubuntu22.04上でのanaconda, minicondaのインストール方法についての備忘録です。 他のLinuxのディストリビューションでも概ね同様と思い…
お疲れ様でした。 非圧縮性流体のSPH法 論文 A reliable SPH(2) formulation for Darcy–Forchheimer–Brinkman equation using a density-based particle shifting in the ALE description | Computational Particle Mechanics が興味深い。 近年は非圧縮性流…
我々の業界でも地形データの解像度の増加に伴ってGeoTIFFの使用頻度は増えています。 点群密度の増加によってGeoTIFFの容量も大きなるため、データベースの構築やアプリケーションでの使用を考えると、適切な圧縮アルゴリズムを選択する必要があります。 過…
本年もよろしくお願いします。 導入 公式の方法を使用時の不具合 暫定的解決策 GitHub 導入 Pythonユーザーであれば、.pyファイルのようにJupyter Notebooksファイル（.ipynb）をimportしたい場面はよくあります。 その方法は、以下の公式サイトで紹介されて…
勉強会お疲れ様でした。 今回は雑多な話が多かったので、メモ的に残しておきます。 中山間部の河道計画での河床変動計算(土砂流出を含む)の必要性 水位を陰的解くことでコロケート格子の保存則の問題を解決できないか。 三次元解析で本間の式のような経験則…
Click here for the English page TL;DR OneDrive(自動同期時)上のExcelではVBAのThisWorkbook.Path、ActiveWorkbook.Pathなどが上手く動作しません。 解決策としては、本記事に示したfuncOneDrive.bas(https://computational-sediment-hyd.github.io/onedri…
TL;DR オープンデータ（地形データ）とオープンソース（DELFT3D）によって、誰でも実河川の水理解析ができるようなった。 はじめに 兵庫県「全県土の高精度3次元データ」について（Copilot） DELFT3Dとは？（Copilot） 計算条件 計算範囲 標高、河床地形 計…
勉強会お疲れ様でした。 一般座標系平面二次元河床変動計算を書く DELFT-3Dについては概ね特性が理解出来かも。 OpenFOAMについて分析を進める予定。 引き続き作業を進める。 AMR法（Adaptive Mesh Refinement、適合格子細分化法） 進捗なし。引き続き、作業…
久しぶりのOpenFOAMです。 Windows11-WSL2:Ubuntu22.04にOpenFOAM11をインストールしますが、今回は社内環境等を想定してプロキシ環境下での方法をまとめておきます。 他のUbuntuのバージョンでも概ね同様かと思います。 Ubuntuのプロキシの設定 参考サイト …
GISのデータ形式として主流になりつつあるGeoPackageについて、Pythonでの取り扱う方法に関するメモです。 GeoPackageとは（Copilot） 参考資料 PythonでGeoPackageで取り扱うためのライブラリ サンプルデータ GeoPackageへの書き込み GeoPackageからの読み…
勉強会お疲れ様でした。 一般座標系平面二次元河床変動計算を書く DELFT-3Dが参考になりそう。スキームを含めて分析を行う。 渦動粘性係数の設定方法について確認する。実河川での計測値は整理済み。計測方法を調べておく。 ダムブレークによるスキームチェ…
実河川での渦動粘性係数の観測結果に関する論文を備忘録として簡単にまとめておきます。 渦動粘性モデルについて 鉛直方向の渦動粘性係数の分布 横断方向（川幅方向）の渦動粘性係数の分布 流下方向の渦動粘性係数の分布 GitHub 参考文献 渦動粘性モデルにつ…
たまには古典的な水理学の話でも。 概要 本間の越流公式 なぜ未だに多用されるのか 表記上の注意 GitHub 参考文献 関連記事 概要 本間の越流公式とは越流堤上下流の水位と通過流量の関係を示した実験式です。 以下の2報で発表されたのが最初であり、現在でも…
GitHub上の.pyファイルを直接インポートできると便利だなと思ったのでコードを書いてみました。 GitHub上の.pyファイルのインポート方法 参考：GitHubの.pyのソースコードのurlの取得方法 GitHub 参考サイト GitHub上の.pyファイルのインポート方法 基本的に…
勉強会お疲れ様でした。 やる気がでてきました。 一般座標系平面二次元河床変動計算を書く DELFT-3Dが参考になりそう。スキームを含めて分析を行う。 ダムブレークによるスキームチェック 変数定義点、カップリングの影響検証 高速解法が必要、陰解法の適用…
