地理情報システム

地図製作者、地理学者、およびGIS専門家のためのQ&A

4
QGISまたはPostgreSQL、またはその他の無料アプリケーションを使用して空間自己相関を実装しますか?[閉まっている]
閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか?この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 3か月前に閉鎖されました。 GISは初めてです。 PostgreSQLデータベースには2つのレイヤーがあります。 最初のレイヤーには、ヨーロッパでイベントが発生したポイントが含まれています。2番目には、ヨーロッパのシェープファイルが含まれています。 QGISを使用して、ヨーロッパのマップ上のポイントをマップできました。空間オートコレレーション(Moran i)を実装したいと思います。 GeoDAを使用してこれを行うデモをいくつか見ましたが、単一のシェープファイルのみを処理するようです。 QGIS、PostgreSQL、または他の無料アプリケーションを使用してこれを実装する方法について、誰かが正しい方向に私を指すことができますか?

2
ジオリファレンスの結果を改善しますか?[閉まっている]
閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか?この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 2年前に閉店。 背景これは、異なる座標系で他のデータレイヤーと組み合わせて再表示するために、裸のラスターマップをジオリファレンスすることに関する2番目の質問です。前の質問は、メタデータを含まない任意のマップ画像をQGISプロジェクトに変換することです 問題私の目標は、このマップをジオリファレンスすることです。 これはPlate-Carréeではないようです。そこで、QGISで、いくつかの妥当なコントロールポイントを作成しました。完全を期すために、下部に追加しました[参照:1]。プロジェクトファイルEPSG:4326と同じターゲットSRSをQGIS Georeferencerに提供します。Helmertと多項式変換では非常に悪い結果が得られますが、薄いプレートスプラインでは妥当な画像が得られます(結果として得られる地理推定値がコントロールポイントを通過します)。ただし、この結果でさえ、たとえば高緯度では貧弱です(日本の北のロシア沿岸を参照)。これは、Natural Earthの背景を使用したQGIS画面のスクリーンショットです。 代替パスでの結果とコントロールポイントを参照してください。私はMapWarperではるかに簡単に使用できるツールと同様の行使しようとしたhttp://mapwarper.net/maps/758#Preview_Map_tab私はおそらく事実によるもの(貧しい結果を取得する場所を追加したコントロールポイントが少ないこと)。 一言で言えば質問 良いジオリファレンスを取得するのに欠けているトリックはありますか? この投影はすぐに認識可能ですか? で不明な古い図面上の座標系、gdaltransform実際に元のマップを生成するために使用される投影パラメータを暴くことを目標に、いくつかのターゲットSRSに、いくつかの座標点を変換することが示唆されます。私はこのような何かを試してみました:ポイントの私のQGISリストを保存した後、私は経由して、スペース区切りロング/ラッツのリストを取得するには、いくつかの文字列処理をしたcat eurasian-steppe-gcp.points | tail -n+2 | cut -d, -f1-2 | sed 's/,/ /'> tmp.txtとgdaltransformに結果のファイルを入力:gdaltransform -s_srs EPSG:3785 -t_srs EPSG:4326 < tmp.txtスイッチングs_srsとt_srsフラグ(プロジェクトはEPSG:4326を使用します)。幸運なことを望んで、暗闇で撮影していることを知っているので、出力を理解できなくても驚かなかった。誰かがこの方法を使用してソースマップの投影と投影パラメーターの最適な推定値を見つける方法を拡張できますか?この背後にある私の考えは、無数のコントロールポイントを配置して適切なジオリファレンスを作成するのではなく、すべての一般的な座標系をループするだけで、より少ないコントロールポイントでほぼ完璧なジオリファレンスを取得する方が簡単だと思いますか?テスト対象の各CRSについて、他のすべてのポイントに対する各ポイントの相互検証が含まれますか? プロセスを自動化できるように、このアルゴリズムまたはジオリファレンスのいずれかを理解したいと思います---この問題に常に遭遇し、コンテンツ作成者がマップを統合されることのない一回限りの作成として扱うのをやめるまで他のコンテンツで、私は停止する予定はありません。 参照資料 [参照:1] QGIS GCPファイル: mapX,mapY,pixelX,pixelY,enable 142.632649100000009,54.453595900000003,505.941176470588232,-95.220588235293974,1 154.934252200000003,59.559921699999997,536.411764705882206,-52.779411764705742,1 80.080158100000006,9.657192300000000,291.558823529411711,-322.661764705882206,1 10.448442600000000,57.819128900000003,21.676470588235190,-103.926470588235134,1 34.007173000000002,27.761438299999998,101.117647058823422,-244.852941176470466,1 50.950890399999999,11.862196600000001,171.852941176470495,-313.955882352941046,1 29.713217199999999,60.024133200000001,90.779411764705799,-92.499999999999829,1 60.000000000000000,0.000000000000000,208.308823529411683,-362.382352941176350,1 69.867506500000005,66.639146199999999,224.088235294117567,-33.191176470588061,1 27.276107100000001,71.049154799999997,89.147058823529306,-21.764705882352814,1 140.000000000000000,0.000000000000000,536.955882352941217,-362.926470588235190,1 20.000000000000000,0.000000000000000,43.441176470588132,-362.926470588235190,1 …

4
ArcGIS DesktopまたはQGISを使用してモザイク処理中に黒いエッジを削除しますか?[閉まっている]
閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか?この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 3か月前に閉鎖されました。 3つの画像をモザイク化します。各画像には黒い背景があります。背景は、画像をモザイク化した後も残っているため、画像間のギャップのようなものです。 ArcGISまたはQGISを使用してこの黒い色を削除するにはどうすればよいですか。

2
ggmapでデンマークの地図を作成する方法
次の2つのパッケージをインストールしました。 library(ggmap) library(maps) 世界地図を作成するためのコードをオンラインで見つけました(以下を参照)。国レベルでプロットをズームするにはどうすればよいですか?たとえば、デンマーク。 ggplot()+ borders("world", colour="gray50", fill="gray50")
13 r  ggplot2  ggmap 

1
GRASS処理はQGISで動作しませんか?
GRASSアルゴリズムを処理ツールボックスから機能させるには、いくつかの問題があります。 QGIS 3.4(MacOS High Sierra)を新規インストールすると、最初にpythonが見つかりません-カスタムPATH変数を作成して解決しました。次に、タブとスペースの一貫性のない使用について文句を言いますが、これもQGISアプリ内のgrass74ファイルを編集することで解決しました。今、私は私が使用するGRASSアルゴリズムにかかわらずこのエラーを受け取っています: WARNING: Default locale settings are missing. GRASS running with C locale. Default locale not found, using UTF-8 Traceback (most recent call last): File "/Applications/QGIS3.app/Contents/MacOS/grass/grass74", line 2029, in <module> main() File "/Applications/QGIS3.app/Contents/MacOS/grass/grass74", line 1880, in main set_paths(grass_config_dir=grass_config_dir) File "/Applications/QGIS3.app/Contents/MacOS/grass/grass74", line 622, in set_paths os.environ['MANPATH'] = sys_man_path …
13 qgis  grass  installation  mac 

2
QGIS 3.0のRスクリプト?
QGIS 3.0の処理オプションにはRスクリプトが含まれていないようです(SAGA、GRASS、GDALとともに利用可能なものとしてリストされていません)。これはこのバージョンのQGISの変更ですか、インストールに問題がありますか?
13 qgis  r  qgis-3 

6
QGIS 3.0でのジオリファレンス?
QGIS 3.0(OS XとWindows 10)にアップグレードしたばかりで、ジオリファレンスの方法がわかりません。ソフトウェアから、GDALジオリファレンスはまだインストールされており、削除できなかったため、「ラスター」カテゴリにあるべきだと言われました。しかし、メニューの他の部分にも見当たりません。 新しいプロファイルが作成されましたが、同じ問題があります...:

1
MacOSからQGISをアンインストールする方法
多くの人が知っているように、QGISのMacOSインストーラーはWindowsバージョンほどターンキーではありません。基本的に3〜4セットのパッケージをインストールする必要があり、そのうち1つだけがMacアプリディレクトリ(GDALなど)に格納されます。QGISを完全にアンインストールする場合、GDALなどのパッケージをアンインストールする方法を教えてください。
13 qgis  mac 

2
地図上で「ターミナル」都市を見つける
終わりの村/都市を見つけたい。つまり、この村はX kmの距離で最大1つの村にリンクされている必要があります。 私はこのアルゴリズムがどのように呼ばれているのか疑問に思っています(誰かが以前に似たようなものを作成したと思います)。 このためのツールはありますか?この概念はどのように呼ばれますか? 例えば City ------------ Village 1 ------- Village 2 ------ Village 3 - - - [ Mountain] \ \ \ Village 4 | Road | \ in mountains Village 5----------Village6---------Village7 - - - - - 限り山の中の道は長いXkmよりも、Village 3そしてVillage 7私が考えられるであろう、端末彼らは(それぞれV2、V6)1つの村にリンクされているため、。 そうでなければ、私たちが持っているオープンデータを使って構築したいと思います。 私が試したもの: 私の国(ルーマニア)のOSMデータをダウンロードし、村や都市をデータベースにインポートします。データベースのジオロケーション機能を使用して、R kmの半径内にX個以上の村がない村を見つけることができます。 しかし、これは私の場合の解決策ではありません。私の場合、上の例のように村は山の反対側にあるかもしれませんが、良い方法はありません(またはまったくありません)。

1
放射照度分析-GRASSとSAGAの大きな差異
プロットの放射照度値を計算して視覚化したいと思いました。理由はわかりませんが、QGIS 2.18.5のコピーでは、「地形分析->稲妻」で適切なSAGAモジュールが欠落しているため、GRASSの「r.sun」アルゴリズムを選択しました。 結果は驚くべきものでした。解析が行われたラスタが適切にジオロケーションされているにもかかわらず、プロットはポーランド東部ではなく金星に配置する必要があるようです。ここでは、6月21日にほぼ5 kWh / sqメートルを受信することは不可能です。 数値を再確認するために、SAGA 5.0のスタンドアロンコピーを見つけて分析を再実行しました(「Potential Incoming Solar Radiation」アルゴリズム)。今回の結果はより信頼性が高かった(比較のためにQGISにインポートされたスクリーンショットのラスター)。 これらの2つのアルゴリズムは大きく異なりますか? 誰も同じ問題に直面していますか? それでもこの機能のテストのみ。 QGISバージョン: 2.18.5 GRASSバージョン: 7 SAGAバージョン: 5.0.0。 入力:ラスターの標高、勾配、およびアスペクトデータ(3つの個別)。SAGAは、高さラスターのみで実行されました。GRASSはすべて3を使用しました。

3
最大の幅/高さの比率で「斜めの境界ボックス」を作成しますか?
以下の画像に示すように、入力フィーチャのセットから最大の幅/高さの比率(「斜めの境界ボックス」)で囲む長方形を作成したいと思います。 たとえば、PostGIS関数ST_EnvelopeおよびST_Box2DFMEトランスフォーマーBoundingBoxReplacerをBoundingBoxAccumulator使用して、軸平行バウンディングボックスを生成します。 PostGIS、QGIS、またはFMEを使用したアプローチを使用した提案は大歓迎です!
13 qgis  postgis  fme  extents 

3
最長辺を向いたポリゴン内にグリッドを作成する方法は?
行と列の属性が割り当てられたポリゴン内にグリッドを作成するために、「Create Polygon Graticule」というプロセスを使用しています。QGISで実行しているSAGAツールです。セルがNからSおよびEからWである入力ポリゴンの範囲に基づいてグリッドを生成します。グリッドの方向がポリゴンの最長エッジに平行である場合、分析に適しています。参考のために、現在の自動出力と、グリッドをポリゴンに合わせるために手動で変更した出力の画像を添付しました。グリッドは元々ポリゴンの範囲全体をカバーしていましたが、内側にあるセル部分のみを表示するためにクリップされています。 上記の画像は、ツールがデフォルトで出力するものです。 これは、必要なパラメーターに合うように手動で編集されたグリッドの画像です。 私の質問は、自動化された方法で2番目の画像のグリッドのように整列したグリッドを作成する方法はありますか?数千のポリゴンのグリッドを生成する必要があります。

1
分岐した多角形の形状を検出しますか?
何百万ものポリゴンが連続したカバレッジを持つベクターレイヤーがあります。形状に応じて分類する必要があります。私はすでにコンパクト(のような景観生態学からいくつかの形状インデックスを使用しています4piA / P ^ 2)、平均幅(2A / P)、形状番号(P / SQRT(A) )、私も見ました、この答えをするの丸み/コンパクトさを計算しますポリゴン? 私の問題は、これらすべてのメトリックが面積と境界のある比率のみを使用していることです。フラクタル次元インデックスでさえ、面積と周囲長(2ln(0.25P)/ ln(A))のみを使用しています。しかし、面積と周囲長が同じで形状がまったく異なる2つのポリゴンをどのように区別できますか?この分岐ポリゴンAのように: 曲線ストリップBと同じ面積と周囲で描画しようとしました。既知のインデックスはすべて同じです。しかし、私にとっては、単純なストリップ(新月のような曲線を含む)と複雑な分岐形状を区別することが非常に重要です。 多角形Bを意図的に直線ストリップではなく曲線ストリップとして表示します。これは、直線の細長い形状を検出する関連円インデックスを知っているためですが、私の多角形も同じ外接円を持っている可能性があります。凸包を構築し、面積Apolygon / Aconvexの比率を計算しても、ここでは非常に似ているかもしれません。 それでは、ベクターデータで分岐ポリゴンAとポリゴンBを自動的に明確に区別するにはどうすればよいですか?(ラスターに変換するには、非常に小さいセルサイズ、膨大なデータセット、メモリ不足が必要になるため、不可能です)。他のパラメーターを含む他の形状インデックスはありますか?理想的には、この方法は明確に分岐したポリゴンだけでなく、CとDを区別するでしょう: 私の唯一のアイデアは、凸包を構築し、その凸包からポリゴンを消去し、残された(大きな)ピースの数をカウントすることです(レイヤー全体ではなくポリゴンごとにポリゴンを消去します)。これは、ボーダーの複雑さを示している可能性があります。 後でPythonで実装する数学的な解決策/アルゴリズムを歓迎します。

1
QGIS Measure Toolから距離を取得していますか?
入力マスクを使用してSpatialiteデータベースにフィードするプラグインを作成しました。オブジェクトの最小/最大次元を表す2つのフィールドがあります。 iface.actionMeasure().trigger()2つのプッシュボタン「メジャー」でメジャーツール()を呼び出した瞬間に、寸法をクリックし、キーボードでフィールドに値を入力します(QLineEdit)。 測定ツールで右クリック(測定の終了)を行った後、値が自動的にフィールド(最小または最大、QLineEdit)にコミットされる可能性はありますか? 測定ツールの呼び出し(動作します) 左クリック:測定を開始 4.一度に:右クリックして、測定ツールの合計値を最小フィールド(QLineEdit)にコミットします。 更新1: この質問をしたある日、メジャーツール(ST_Length(geom))として使用できるラインストリングレイヤーを作成するというアイデアを思いつきました。利点は、そこからデータを操作および取得することを完全に制御できることです。短所は、プロジェクトに追加の層があることです。 ワークフローは次のとおりです。 測定をクリック:編集を切り替え、「機能を追加」を有効にします オブジェクトの最小寸法を表す線ストリングを作成します 変更の確認、ラインストリングの長さの計算(ST_Length(geom))、QLineEditへの長さのコミット ダイアログを閉じてラインストリングを削除(DELETE FROM) しかし、私の実際の質問は残っています。測定ツールから直接距離値を取得する方法はありますか? アップデート2: Jakobの助けを借りて、Measure Tool を閉じるまで、Measure ToolのTotal値でQLineEditを更新できます。次に、値は「0,000 m」に設定されます。私はこれを使用してこの問題を解決しました: def text_changed(str_measure): if str_measure == '0,000 m': pass else: self.min.setText(str_measure) 動作しますが、これが適切な解決策であるかどうかはわかりません。

2
奇妙な形のポリゴンの重心座標
QGISで特定のポリゴンの重心を計算する際に問題があります(ジオメトリツールを使用):特定のポリゴン-奇妙でかなり長い形状の場合-計算された重心はポリゴンから外れています。たとえば、下の図の赤でマークされたポリゴンを見てください。 その重心(青い矢印で示されています)は、実際には対象のポリゴンの外側にあります。これは、少数の(まだ厄介な)他のポリゴンで発生します。 私はQGISが比較的新しいので、この問題がシェープファイルとポリゴン(WGS 84)の投影によるものかどうかは考えていません。とにかく、ベルギー(私が見ている国)に対応するWGS84 / UTMゾーン31で地図を再投影してみました。 私は過去の議論で同様の質問に対する答えを見つけていません。誤って見落とした場合は申し訳ありません。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.