vbaのプログラムを、とりあえずはっときます。でもこれあってるのかなー？ 続きは、また今度 Sub べき乗分布サンプリング() 'y = 1/x の場合 Randomize x = 1 '初期値 For i = 1 To 6000 newx = x + WorksheetFunction.NormSInv(Rnd()) '新しい候補 If 1 <= …

「google earthで表示」を右クリックしてリンクアドレスをコピーして、そのURLの中の「output=nl」を「output=kml」に変更してアクセスするとkmlファイルがダウンロードされます。↓のプログラムをkml2csv.jsという名前で保存して、それに↑のKMLファイルをド…

自然は放っておくと、自然に平衡状態（自己組織化臨界状態と同義語？）になるらしいです。その状態では、大小さまざまな出来事が、べき乗則に従って予測不可能なタイミングで起きます。例えば、地震の発生や、砂山の雪崩現象や、自然林の林冠ギャップの分布…

10mメッシュ標高DEMなどの高解像度のラスタデータを広域に計算しようとすると、とても時間がかかって待ってられません。 そこで、マルチコアCPU(Core2Duoとか)をフル活用して、短時間でラスタ演算できるGDALとOpenMPを使った並列プログラミングの方法を紹介…