地理情報システム

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

2
マップマッチングGPSは道路網をポイントします！
私はPostgresとPostGISを初めて使います。 GPSポイントの単純なマップマッチング（または、必要に応じてパス推論）を試みています。私はpostgres（postgis）を使用しています。 データベースにGPSポイントを含むテーブルがあります： gpspoints（緯度、経度、時間、速度、方位角、ジオメトリ、...） 道路ネットワークのシェープファイルもデータベースにインポートしました： ルート（gid、idrte、version、nomrte、norte、clsrte、geomなど）-リンクの方位角がありません。 postgis関数ST_azimuthを使用してどのように計算できるかわかりません。私が持っているのは、上記の列を含む道路網のシェープファイルだけです。 リンクの方向がGPSポイントの方位角（+または-15度）と一致する場合にのみ、各GPSポイントを（ポイントの周りの20メートルのバッファー内の）最も近いリンクに関連付けて、投影された位置を取得します。それ以外の場合、許容可能な方位角を持つ20メートルのバッファー内で、NEXTの最も近いリンクを検索します。（写真のように！） 投影されたGPSポイントの新しい座標を、テーブル「gpspoints」に「projectedLat」および「projectedLong」として追加する必要があります。 （下の図では、実際のポイントは方向を使用して示されていますが、投影されたポイントには方向マークがありません）
3
Python、GDAL、およびScikit-Imageを使用した画像処理
私は処理に苦労しており、できればここで解決できることを願っています。 林業に適用されるリモートセンシング、特にLiDARデータを使用して作業しています。アイデアは、ツリートップの検出にScikit-imageを使用することです。私はPythonに慣れていないので、次のことをするために大きな個人的な勝利を考えました。 CHMをインポートします（matplotlibを使用）。 ガウスフィルターを実行します（scikit-imageパッケージを使用）。 最大フィルターを実行します（scikit-imageパッケージを使用）。 peak_local_maxを実行します（scikit-imageパッケージを使用）。 ローカル最大値（matplotlibを使用）でCHMを表示します。 今私の問題。matplotでインポートすると、画像の地理座標が失われます。したがって、私が持っている座標は、基本的な画像座標（つまり、250,312）です。必要なのは、画像内の極大ドット（画像内の赤いドット）の下にあるピクセルの値を取得することです。ここで、フォーラムで同じこと（NumPyを使用せずにOGRポイントの下でGDALラスターのピクセル値を取得しますか？）私の場合、ポイントはscikit-imageで計算されました（各ツリーの頂点の座標を持つ配列です）。そのため、シェープファイルはありません。 結論として、私が最後に望むのは、地理座標での各極大点の座標を含むtxtファイルです。例えば： 525412 62980123 1150 ...
2
QGISで二変量カラーランプを作成しますか？
最近見た多くの地図には、連続した推定値間の色の相関に使用できるこれらのきちんとした双方向の色のグラデーションがあります。 以下は最近の出版物からの例です。 クリスチャン・レバー、ピーター・J・バーカーク、ダニエル・ミュラー、ピーター・H・バーブルク、ヴァン・ブティック、ペドロ・J・レイタン、マーカス・リンドナー、トビアス・キューマーレ -ヨーロッパの森林収穫強度パターンのドライバー。エルゼビア。2014年3月1日（アクセス：http : //dx.doi.org/10.1016/j.foreco.2013.12.030）。 これらのグラデーションはで生成されたArcGIS 10.xと思いますが、無料のGISツール（または少なくともLinuxで使用可能なもの）でこのようなカラーグラデーションをどのように再現できるのでしょうか？ 提案やコードスニペットを喜んで聞きます。
1
ベクトルアフィン変換を使用してレイヤー/すべてのフィーチャを移動するQGIS
QGISを使用して、異なるソフトウェアを使用して作成されたいくつかのフライトトラックを分析しています。これで、それらがQGISにあるとき、すべての基準点は飛行場基準点である0,0になります。 しかし、現実の世界では、この飛行場の基準点は、EastingとNorthingの価値です。つまり、Easting = 286499.025 Northing = 6230965.204であるため、他のすべて（滑走路が終了し、トラックポイントは実際の世界とまったく同じになるはずです） いくつかの調査を行ったところ、Vector Affine Transformationプラグインが見つかりました。これを使用して、Easting = 286499.025 Northing = 6230965.204 Easting = XおよびNorthing = Yの基準点にすべてを移動する方法を説明します。 プラグインのマトリックスを理解できません。
1
複数のシェープファイルをPostgisにバッチロードする
shp2pgsqlヘルプによると、オプションdrop、append、createおよびprepareは相互に排他的です。そのため、シェイプからテーブルを作成し、他の複数のシェープファイルを追加する場合、次のような操作を行い、作成モードか追加モードかを示すカウンターを保持します。 cnt=0 for shp in $(ls *.shp); do if [ $cnt -eq 0 ] ; then shp2pgsql -s 27700 -c $shp schema.table_name | psql -h localhost db else shp2pgsql -s 27700 -a $shp schema.table_name | psql -h localhost db fi ((cnt++)) done これは予想どおりに機能しますが、もっと簡単な方法があるのではないかとよく疑問に思います。
2
GeoserverのSLDファイルで単純な数学を使用していますか？
GeoserverのSLDファイルで単純な数学を使用することは可能ですか？データベースからテキストサイズプロパティを読み取り、その値を乗算したいとします。どうすればよいですか？ <sld:TextSymbolizer> ... <sld:CssParameter name="font-size"> <ogc:PropertyName>text_size</ogc:PropertyName> * 2.1 </sld:CssParameter> ... </sld:TextSymbolizer> この単純なアプローチは明らかに機能しないため、これを達成するための何らかの方法があるに違いないと思った。
11 geoserver  sld 
2
QGISで複数の行でポリゴンを分割する方法はありますか？
ポリゴン（.shp）レイヤーとライン（.shp）レイヤーがあります。ポリゴンに重なる線に沿って、ポリゴンをいくつかの部分に分割します。その後、各ポリゴンを個別に編集します（属性テーブル）。 回避策を試してみました：短い距離でラインをバッファリングしてから、ジオプロセシングツールで動作しました（相違点-「Unterschied」）。その後、「マルチパートからシングルパートへ」。それはうまくいきましたが、私はむしろ線でそれをやりたいです。 QGISの拡張機能やタスクを実行する別の方法はありますか？ESRI ArcINFOツール「ライン付きポリゴンの分割」に似ていますか？しかし、私が持っているのはQGIS 2.6だけです。
11 qgis  vector  line 
1
トラックの自動検出
からのオープンエリアにくぼみがあるエリアを与えられます。すべての地形車両が柔らかい表面を走行していると言います。車両は、深さが約10〜20センチメートル、幅が約15〜30センチメートルの平行な窪み（パス）を引き起こし、長さは表面の堅牢性によって異なります。 後の分析に関連するリモートセンシングプラットフォームはどれですか？Quickbird、小型ドローン、LIDAR、航空写真？ パスの文書化に使用できるツール（FME / QGIS / ESRI /その他）で利用可能な手順はありますか？ この質問を単純化するために、エリア内に他のパスがないこと、またはそれらが画像から除外されていることを明確に知っていると仮定します。 完全な自動化は必要ではなく、おそらく不可能です。 これは、トラックがどのように見えるかの例です。
3
ラスターの面積を計算するより正確な方法
私の日常業務では、30秒の解像度での地理的投影でグローバルラスターデータセットの面積を計算するように常に求められます。これらのデータセットは通常、Combine操作の結果です（典型的な例は、国のレイヤーと組み合わせた植生クラスです）。これを行うために、ユニットは、地理投影の各ピクセルの面積が30アーク秒のラスタデータセットを作成しました。このエリアグリッドを使用して、各クラスのエリアを合計するためにzonalstatが実行されます。この面積グリッドがどのように作成されたのかわからないので、このアプローチが等面積投影でラスターを再投影するだけの方が正確であるかどうかを常に疑問に思っていました（簡単なテストから、2つの方法の結果は似ています）。誰か同じような状況を経験しましたか？
3
多角形から多様性の連続測定まで？
まず第一に; 私は成功せずに、同様の質問を見つけようとしました。たぶん、私はGISに慣れていないので、探しているものが正確にわからないからでしょう。誰かが同様の問題を指摘してくれたら、この投稿を削除させていただきます。 特定の国の人口の多様性の「連続」またはラスター（小さなグリッドセル）変数を作成する必要があります。ポリゴンでの民族グループの広がりを示すシェープファイルがあり（図1）、私が探している結果は、各行政単位（AUの場合、この場合、 360ナイジェリアの選挙区）。 図1.ナイジェリアの人口グループポリゴン 私が思いついた解決策は、各AUの各ポリゴンの面積割合を取得し、それから不均一性インデックスを計算することでした。しかし、問題は、行政単位の分布のために、かなりの量の情報を残してしまうことです。図に示すように。2、正方形「a」、「b」、および「c」は同じ「分離指数」を持ちますが、「ホットスポット」に対して同じ位置にないことは明らかです。 図2。 ですから、別のソリューションでグリッドマップを作成し、最も近い境界までの距離を計算できますが、1つの境界のみを共有することは、複数のグループが一緒に住んでいるマップの中央部分にあることとは異なります。 この質問を見つけた後、重心を使用してポリゴンをポイントに変換し、同じ方法を適用できると思います。しかし、真実は、私がこれに慣れていないということです、そして、その質問は本当にはっきりと答えられません。どうしてそんなものを作ることができますか？ 別の例を使用して、次のようなもの（このWebサイトの画像）を作成します。 さまざまな質的特性を持ついくつかのポイントの分布を考えると、各行政単位の「平均不均一性」を推定できる場所から多様性の尺度を取得します。 どうすればできますか？私はRとQGISを使用しているため、どのプラットフォームがソリューションベースであるかは気にしません。
11 qgis  raster  polygon  r 
1
OpenLayers v3で機能の場所を取得する方法
描画インタラクションまたは変更インタラクションを使用した後、変更または追加された機能を返すリスナーがあります。 例 draw.on('drawend', function (event) { // get the feature var feature = event.element; // ...listen for changes on it logStatus(feature.getId()); }); IDを取得する方法は知っていますが、DBに保存する必要があるため、機能の場所（緯度経度）が必要です。どうすればできますか？APIで見つかりませんでした。
1
ジオサーバーの管理者パスワードの変更
ジオサーバー管理者のデフォルトのパスワードを別のパスワードに変更するにはどうすればよいgeoserverですか？ 2つの異なる方法を試しましたが、どれも機能しませんでした。オンラインマニュアルに従って最初に試したのは、data_dirにアクセスすることです。 /usr/local/lib/geoserver-2.5/data_dir/security そして私は見つけることができました <user enabled="true" name="admin" password="digest1:XXXX..."> でusergroup/default/users.xml。しかし、ダイジェストを編集することはできません。 2番目に試したのは、左側のパネルの[パスワード]メニューのWeb管理インターフェイスからパスワードを変更することです。パスワードは変更できますが、マスターパスワードと呼ばれるものを変更しているようです。のパスワードadminは同じままです。 osgeolive8のgeoserver2.5を使用しています。
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.