趣味で計算流砂水理

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

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

MENU

2021-01-01から1年間の記事一覧

地理情報付きラスタ画像から任意座標のデータを抽出するpythonスクリプト

R

最近よく使うスクリプトなので備忘録的にまとめました。 この記事のポイント ラスタ画像から任意座標のデータの抽出はrioxarrayを使うと超簡単 この手の方法は航空機、ドローン等の高密度測量データから計算データセットを作成する際には必須 environment サ…

最近読んで真似したいと思った河床変動計算の論文

R

研究会お疲れ様でした。 最近読んで真似したいと思った河床変動計算の論文を共有しておきます。 三次元元計算 + 河床変動 www.sciencedirect.com 少し話しましたが流れの三次元元計算と河床変動計算(多分平面二次元)で流路の出現、消失を計算してます。 か…

強調用のcssを導入しました

R

以下のサイトを参考に強調のcssを導入しました。 はてなブログで文章を枠で囲むカスタマイズ - Random Life Blog 【CSS】おしゃれなボックスデザイン(囲み枠)のサンプル30 以下のように記事に書くと、 <div class="box27"> <span class="box-title">ポイント</span> <p>ここに文章</p> </div> こんな感じに表示されます。 …

質問回答:HoloViewsの使い方について

R

以下のような質問を頂いたので回答します。 HoloViewの使い方等はどこで勉強されたのでしょうか?良質な参考書等はございますでしょうか? 質問の回答 それでもHoloViewsは魅力的 こんな人にはHoloviewsが向いてないかも Pythonが書けない オブジェクト指向…

岸・黒木による実験式を修正する

R

実験式修正シリーズの第2段として、前回の岩垣の式(岩垣による限界掃流力の式の不連続部を解消する - 趣味で計算流砂水理)に続いて、河床波の遷移と有効掃流力と無次元有効掃流力の関係を示した岸・黒木による関係式(水理委員会移動床流れの抵抗と河床形状研…

オープンデータのALB測量よりjupyter上で陰影図を作成する

R

こんなのを作りました。 モチベーション 使用データ グリットテキストデータよりgeotiffを作成 陰影図を作成 図化 github モチベーション グリーンレーザありのALB(Airborne Laser Bathymetry)に河川測量成果について、 九頭竜川(code:8606070001) 、那賀川(…

pythonの高速化の基礎:二次元配列の処理を例に

R

結論:何も考えずnumbaを使え。 モチベーション テストケース 計算結果一覧 基本ケース リスト内包表記 スライス スライス ケース2 numba numba + リスト内包表記 numba + スライス おまけ おまけ1:配列へのアクセスは結構時間がかかる おまけ2:ループの順…

WSL2で外部ネットワークにつながらない理由がSymantec Endpoint Protectionだった話

R

気付くのほぼ1日かかってしまいしました。 pingが通るのにwgetやupdateがつながらなくて、いろんな方法を試したのですが、その原因がまさかのSymantec Endpoint Protectionでした。 WSL2のUbutnu18.04 LTSでapt updateするとconnection timed out Could not …

河川横断測量成果のフラグの説明

R

河川横断測量成果では、下のように横断座標値の前に整数型のフラグがついています。 -0.200,-200.000,3.279,1.374,-9999.000,-9999.000,52,-9999.000,-9999.000,-9999.000,0.000,0,20121105,8101010001,天塩川水系,天塩川 0,-18.000,1.980 0,-9.800,2.040 0,…

欠損したグリッド標高テキストデータからgeotiffを作成するpythonスクリプト

R

とてつもなくユーザーが少ない話ですが。。。 目的 ファイル名から図郭北西座標を取得 元データをラスタ用の二次元配列に変更 データをgeotiff化 おまけ:VRTファイルの作成 Environment 目的 航空測量成果に納められているグリッド標高テキストデータは、 …

備忘録:numba.jitclassのロード方法が変わった?

R

以前に以下の記事でも書きましたが、numbaによるクラスの高速化に使用するjitclassのロード方法が変わったようです。 computational-sediment-hyd.hatenablog.jp 以前のnumba 0.46.0では、 from numba import jitclass でロードできましたが、最新のnumba 0.…

備忘録:地理院タイルをQGISに読み込む

R

すぐに忘れるので備忘録として残しておきます。 地理院タイルをQGISに読み込む場合は、githubで公開されている公式のwmtsをロードしておくと便利です。 github.com WMTSメタデータ: https://gsi-cyberjapan.github.io/experimental_wmts/gsitiles_wmts.xml …

河川断面データ処理ツールの開発:途中段階

R

GUIの作成はつらいな。。。アホほど時間がかかる。 参考サイト Param subobjects — Panel 0.11.3 documentation Param — Panel 0.11.3 documentation Parameters and Widgets — PyViz 0.10.1a3 documentation Dynamicmap — HoloViews 1.14.4 documentation C…

備忘録:パソコンの画面をgif化

R

ストアからインストールできるようになっていました。 forest.watch.impress.co.jp

備忘録:アルファシェイプ

R

2年ぶりくらいに使おうとしたら見つからなかったので。 アルファシェイプとか 凹包(concave hull)(良いサイトがなかった) と言われる幾何関係のアルゴリズムです。 凸包(convex hull)と比べると格段に複雑になります。 パラメータαへの依存度が高くあまり…

備忘録:web上のcsvファイルをpythonで直接読み込む

R

web上のcsvファイル(サンプル)は、 下のようなローカルファイルと同様の方法では読み込むことが出来ません。 fpath = '***.csv' with open(fpath, "r", encoding='shiftjis') as f: l = next(f) そのため、urllib.requestモジュールを使って以下のように扱…

備忘録:新規PCのセットアップ

R

新規PC(ノートPCのスペック : VAIO Z - 趣味で計算流砂水理) を導入したのでセットアップです。 備忘録として書いておきます。 まずは ソフトウェア関係 Firefox Dropboxデスクトップアプリ google日本語入力 LINEデスクトップアプリ スマホ同期 QTTabBar 7-…

備忘録:TeXをms-wordに似せる2:参照図表番号のゴシック体

R

過去記事備忘録:TeXをms-wordに似せる - 趣味で計算流砂水理の続きでちょっと上級編です。 文中で図表を参照するときにゴシック体にしますが、TeXだと結構面倒だったりします。 まずms-wordで本文をMS明朝で図表番号をMSゴシックにしたものが以下です。 こ…

ノートPCのスペック : VAIO Z

R

新しいPCです。 CPU : intel core i7-11370H @ 3.30GHz RAM : 32 GB OS : Windows 10 Pro 64bit (function(b,c,f,g,a,d,e){b.MoshimoAffiliateObject=a; b[a]=b[a]||function(){arguments.currentScript=c.currentScript ||c.scripts[c.scripts.length-2];(b…

岩垣による限界掃流力の式の不連続部を解消する

R

河床変動の数値計算では、よく使用する岩垣による限界掃流力の式について不具合があったので若干修正しました。なんかこういうのって土木っぽい。。。 限界掃流力に関する基礎的研究 モチベーション 岩垣の式は実験式のため、図化すると関数の接続部で不連続…

備忘録:TeXの表中の箇条書き(改行あり)

R

TeXで下のような表中の箇条書き(改行あり)の作り方です。minipageを使った強引なやり方です。もっといい方法があれば教えて欲しいです。 \begin{table}[t] \begin{center} \caption{表中の箇条書き(改行あり)minipage} \footnotesize \begin{tabular}{|p…

備忘録:TeXのローマ数字

R

なんで準備されていないのか不思議です. 以下の方法が最もきれいな気がします.文中で使うと若干詰まりますが. LaTeXで(文中に)ローマ数字を使う方法 - Qiita ここにも書いてますが,ググるとアルファベットで打って文字間隔を調整する方法が引っかかりま…

備忘録:bibtexでエスケープが必要な記号

R

久しぶりにbibtexを使ったのでアンダーバーに気付かずに1時間程エラーを探してました. beiznotes.org

仮説が間違っていたかも:農業用ダムの下流河道の流砂環境へ与える影響

R

一般的な治水機能を持つダムの場合, ダムができると応答の遅い礫河床でも運用後30年も経過すれば,河床材料の粗粒化やみお筋の固定化や植生域の拡大等がみられるます. これは,出水のピークカットとダムによる流砂の遮断によるものです. なので,出水とそ…

備忘録:TeXをms-wordに似せる

R

論文の執筆期間に入りましたが、私が書く系の論文はフォーマットがms-wordしかないものが多く、TeXユーザーとしてはつらいところです。 さらに、私はlualatexを使うので、日本語latexはplatexユーザーがメインなので、フォーマットは一から作ることが多いで…

備忘録:河川コード

R

今年は実河川ものを多くアップロードしていきたいのですが、河川コードを書くことを習慣付けようかと。 (ほとんど誰も使わないのですが。。。。) 鉄板ネタとして「荒川」は一級水系に2つありますが、日本の全河川名でみると32河川あります。 (チェック用…

備忘録:python-selenium-googlechromeで画面コピー

R

グラフをhtmlで作ることが多いのですが、いざ紙ベース資料(??)を作るときに手動でやるのは大変なので、seleniumで画面コピー(スクリーンショット)をやる話です。 以下の参考サイトそのままですが。 Python + Selenium で Chrome の自動操作を一通り - …

備忘録:python-selenium-googlechromeのchromedriverのインストール方法

R

以下のサイトのようにいろいろな方法がありますが、 Python + Selenium で Chrome の自動操作を一通り - Qiita 私はcondaでやるほうが楽だと思います。 ただし、バージョンには注意が必要です。google chromeのバージョンと合わせて下さい。 バージョンの確…