私はなぜかSAGA（System for Automated Geoscientific Analyses）が好きで、QGISで処理をするときは、ほとんどSAGAプラグインを使います。先日PCがとんだので、QGISをインストールし直したら、SAGAが全然動きませんでした。で、そのときの対応方法です。

QGIS SAGAのエラー

https://gis.stackexchange.com/questions/209902/problem-with-saga-in-qgis-2-14-6

によると、SAGA2.2.3までしか対応していないとのことです。私がインストールしたQGIS2.18 64bitでは、SAGA2.3.2でした。

修正方法1 ⇒ 駄目でした。

上のサイトにあるように、SAGAをプラグインを2.2.3に変更すれば良いと考えて、2.2.3をダウンロードしてきて、プラグインのディレクトリを差し替えてみました。

でも、全く動きませんでした。ちなみに、QGIS2.14では、この対応により解決しました。

修正方法2

https://gis.stackexchange.com/questions/217414/saga-algorithms-not-working-on-qgis-2-18-on-mac-osx

にあるように、saga コマンドの変更により、"INTERPOL"が"RESAMPLING"が変更になったようです。それが、プラグインに反映されていないことが原因です。

なので、上記サイトを参考に修正しました。

対象フォルダに移動

私の環境では以下でした。

c:\OSGeo4W64\apps\qgis\python\plugins\processing\algs\saga\

bashで以下を実行

mv description description_old
mkdir description
cd description_old
for i in *.txt; do sed "s/|INTERPOL|/|RESAMPLING|/g" $i > ../description/$i ; done

こんなときは、bashが便利ですね。

以上でSAGAが使えるようになりました。

2017-10-10

seabornで日本語表示

何かとバタバタしていて全然書けないです。今月は書きます。

seabornで日本語表示

windowsが変わるごとにフォントが変わるのでもう大変です。多分この方法が一番ラクかなと思います。

ここを参考にしました。

ここからIPAexフォントをダウンロードします。
解凍し、C:\Anaconda3\Lib\site-packages\matplotlib\mpl-data\fonts\ttfにコピー

これだけで、下準備はO.K.です。

後はpython上で

import seaborn as sns
sns.set(font='IPAexGothic', context='notebook', font_scale=1.5, style="whitegrid", palette="deep")

として、普通にmatplotlibでグラフを書けばO.K.です。

windows10と同じ、Yu Gothic UIのttf場所が全然わかりません。

2017-09-28

Upwind Flux Schemeの問題点と修正

2d Upwind Flux Schemeの計算で厳しい条件の場合、全然回らなくて苦労していたのですが、1dでいろいろテストして、ようやく見つけました。
Upwind Flux Schemeに限ったことではないですが、水面の不連続面が出てきた場合、厄介なことになります。コロケート格子特有かもしれませんが、下のgifのように水面はなんとか解けますが、どんどん水が増えてしまいます。

f:id:SedimentHydraulics:20170928205415g:plain

そこで、水面と河床高の関係より、不連続面を判断し、そこの運動量フラックスをキャンセルすることにより、連続性は確保されるようになりました。また、二次元でもテストしてみます。

f:id:SedimentHydraulics:20170928205504g:plain

一次元だとここまで急変する断面はないですが、平面二次元だと普通にでてくるので、注意が必要です。とゆうか、こうゆう話ってどこにも書いてないですよね。あえて挙げると以下を参考にしました。

http://www.impact-project.net/publications/wp3_pub7.pdf

さくっと書いてますが、見つけるまでに相当苦労しました。最初はドライベッドの問題だと思い、そこばかりを気にしていたのに全然違いました。

gist4faa1ea740856745c1892422d06de8fd

2017-09-12

続：改めて不等間隔格子の河川流一次元計算

冷静に考えたら、座標変換ですね。

基礎式です。

$\frac{\partial A}{\partial t} + \frac{\partial \xi}{\partial x}\frac{\partial Q}{\partial \xi} = 0$

$\frac{\partial Q}{\partial t} + \frac{\partial \xi}{\partial x}\frac{\partial }{\partial \xi} \left( \frac{Q^2}{A} \right) = -gA\frac{\partial \xi}{\partial x}\frac{\partial H}{\partial \xi} -gAI_e$

$\frac{\partial \xi}{\partial x} = \frac{1}{\Delta x_{i+1} + \Delta x_{i}}\left(\frac{ \Delta x_{i}}{ \Delta x_{i+1} }(\xi_{i+1} - \xi_{i} ) + \frac{ \Delta x_{i+1}}{ \Delta x_{i} }(\xi_{i} - \xi_{i-1} )\right)$