QGISの重み付きポイントから密度マップを作成する最良の方法は何ですか？

ArcGISカーネル密度ツールの機能を繰り返したいのですが、どうすればよいですか？

簡単にインストールできるフリーでオープンソースのソフトウェアへの回答を制限したいと思います。つまり、QGISはすべてのプラットフォームでGRASSを使用して簡単にインストールできるため、QGIS + GRASSは問題ありませんが、SAGAはそうではありません（残念ながら、理想的なツールである）。

QGISの保護地域全体の野生生物の密度マップを作成しようとしています。次の例の密度マップは、ArcGISで野生生物のポイント観測から作成され、グループ密度（ベクトルファイルのフィールド）の生のカウントを使用して、カーネル密度（SpatialAnalyst）ツールで各グリッドセルの密度を重み付けします。選択した検索半径とグリッドセルサイズ：

密度に関する以前の質問では、GRASS v.kernelツールを使用してArcGISカーネル密度ツールを模倣することが提案されましたが、v.kernelは同じ仕事をしません。マニュアルを見て、（成功して）密度マップを作成した後、v.kernel関数はポイント密度でのみ機能し、各ポイントに重みを付けるために各ポイント（rawカウントなど）に変数を与える機会がないようです。

更新

*.surf.*GRASSには、密度サーフェスの作成に役立つさまざまなツールがあるようです。これらのツールは、重み付け列またはZ値を受け入れるか、ラスタで実行されます。@underdarkはv.surf.rstを提案しました-そして、「zcolumn」は私の重み（カウント）変数になります-しかし、特定のグリッドサイズを作成するか、特定の半径を使用するようにツールに依頼する方法がわかりません。

v.surf.rstまたはその他の方法の使用方法に関する提案

サンプルデータ

x,y,count
431250,8707500,0
418750,8707500,5
413750,8707500,3
411250,8707500,1
408750,8685000,0
411250,8685000,0
416250,8685000,0
416250,8682500,6
411250,8682500,3
418750,8680000,0
433750,8677500,3
421250,8677500,0
423750,8675000,1
431250,8672500,0
428750,8672500,2
426250,8672500,2
423750,8670000,0

そのマニュアルページによると、GRASSコマンドr.resamp.filterがために行いますポイントデータを表すラスタのArcGISが対応致します正確に何ポイントレイヤー：使用filter=box「単純な」ラスターのためのオプションを使いfilter=gauss他のArcGISカーネル用のオプションを選択します。-nフラグを使用して、ヌルの伝搬を回避します。

カーネル密度の推定値（別名「ヒートマップ」）はデータの補間ではないことに注意してください。 位置xでのKDEの値は、xの近くの単位面積あたりの値「Z」の量を推定します。（半径または「帯域幅」は、「近く」の意味を定量化します。）Zの値は、マップ上のすべての可能な場所で定義する必要はありません。たとえば、Zは人などの存在を表す場合があります。その場合、KDEは人口密度を示します。Zの値は、マップ全体で連続的に変化する必要もありません。補間のために、それは想定しているZをはすべての場所で定義され、データは指定されたポイントでのZの値の観測値であることを示します。補間器は、他のすべての点で観測されていないZの値を予測しようとします。とき、これは理にかなってZは、たとえば、温度や圧力であるが、とき、通常は無意味であるZは、データが完全国勢調査のときに何かの存在を記録したり。（後者の場合、ある地域の道路密度マップが意味するものと、道路以外のエリアにまたがる道路を「補間」する意味を理解する方法を考えてください。）

SAGAのモジュール「カーネル密度推定」はあなたが探しているものです。

QGISにSAGAモジュールインターフェースをインストールし（メニュー：[プラグイン]-> [Fetch Pythonプラグイン]）、モジュールを使用します。幸運を！

カーネル密度よりもArcMapのポイント密度に近いGRASS GISを使用した非常に単純な方法：

1. r.in.xyz指定されたラスター解像度でmethod = sum を使用して、ポイントをラスターマップにインポートします（で設定g.region）。

2. 使用r.neighbors方法と地図=（デフォルト）平均を平滑化し、検索半径を設定するためのオプションのサイズを使用します。

（現時点ではGRASSにアクセスできないので、実際にはこれを試していない！）

v.surf.rstの詳細なガイダンスをリクエストしたので、ここに私の入力を示します

最初に、グリッドサイズについて-プラグイン-> GRASS-> Edit Current GRASS regionを使用して、出力解像度を設定します。v.surf.rstからの出力にはその解像度があります。

半径については、「張力」がパラメータのようです。私はこのアルゴリズムの専門家ではありませんが、マニュアルから読むと、これは関連するビットのようです

「...高張力は「ポイント間の距離を広げる」と各ポイントの影響範囲を狭め、低張力は「距離を減らす」とポイントは長距離にわたって相互に影響します。）

したがって、半径パラメーターを使用するのとほぼ同じように、張力パラメーターを使用できます。

サンプルデータから、v.surf.rstの結果は以下のようになり、カウントを補間の重みとして使用していることを考えると妥当なようです。

私はそれをテストしていませんが、QGIS Contributorsリポジトリには「Rを使用したホーム範囲推定」と呼ばれるプラグインがあります。これには、カーネル（密度？）計算が含まれます。うまくいけば、これが最良の方法だと思います。Rは、カーネル密度を計算する実際の統計的方法を実行します。

Rをインストールしている場合は、プラグインをインストールして試してみることができます。

qgis以外でJavaプログラミングを少し行うことに同意する場合は、この密度マップビルダーライブラリを使用できます。

constructorを使用するとHeatMapBuilder(int w, int h, int[][] pts, int[] weights)、必要に応じて各ポイントに重みを付けることができます。出力イメージはgetImage()メソッドで取得し、でディスクに保存できますImageIO.write("mymap.png")。

出力の例を次に示します。

スムージング強度とカラーパレットを変更することができます。