地理情報システム

QGISジオメトリツール->ポリゴンの重心がポイント座標をどのように正確に計算しているかを知りたいのですが。たとえば、ポリゴンの境界を小さなポイントに分割し、最大距離を持つ2つのポイントのペア（長さと幅）を2つ取り、これらの2つの線の交点を重心と見なしますか？

10 qgis  polygon  geometry  centroids 

ArcGIS Desktop v10.2.2を使用して、美術館をポイントとして表示するマップを作成しています。Museumsレイヤーにはシンボルがなく、ラベルは[ID]フィールドから取得した、IDが内側にある色で塗りつぶされた円でできています。 [LABEL]手動で作成せずに、博物館のリスト全体をレイアウトマップの凡例に追加し、フィールドのIDと説明（それらの名前）を含むシンボルのリストを作成します。 手動による解決策として考えられるのは、レイヤーシンボルを「一意の値」に設定し、すべての値（100以上）を追加し、各ポイントに同じラベルシンボルを1回設定することです。しかし、もっと良い解決策があることを願っています！ 以下の私の地図を見てください。 レイヤーは「Musei」（ミュージアム）で、左側の凡例にあるように、シンボルはありません。右側には、説明された凡例の記号、[ID]フィールドから取得された内部に赤いテキストが入った暗い灰色の円が表示されます。レイアウトに、IDのあるすべての円と対応するラベルを含む列を追加したいと思います。 -更新- これは、@ FelixIPによって提案されたソリューションに従って使用したフィールド計算機の設定です。 しかし、ポイントは期待どおりに垂直に配置されていません。代わりに、それらはあちこちに少し移動されます。また、10、100、1000、10000のさまざまなステップ値を使用してみました。

10 arcgis-desktop  arcmap  arcgis-10.2  symbology  legend 

重複したエントリや重複したエントリがないか、長期間にわたる鳥の観察を確認する必要があります。 異なるポイント（A、B、C）の観測者が観測を行い、紙の地図にマークを付けました。種、観測点、およびそれらが見られた時間間隔の追加データを含むラインフィーチャに取り込まれたライン。 通常、観察者は観察中に電話で連絡を取り合いますが、時には忘れてしまうため、重複した行を取得します。 私はすでにデータを円に接する線に減らしたので、空間分析を行う必要はありませんが、各種の時間間隔を比較するだけで、比較によって見つかったのは同じ個体であることを確信できます。 私は今、Rで以下のエントリを識別する方法を探しています。 重複した間隔で同じ日に行われます そしてそれが同じ種である場合 異なる観測点（AまたはBまたはCまたは...）から作成されたもの） この例では、重複する可能性のある同じ個人のエントリを手動で見つけました。観測点が異なり（A <-> B）、種は同じ（Sst）で、開始時間と終了時間の間隔が重なっています。 ここで、data.frameに「duplicate」という新しいフィールドを作成し、両方の行にエクスポートできるように共通のIDを与え、後で何をするかを決定します。 すでに利用可能な解決策を探してたくさん検索しましたが、種のプロセスをサブセット化する必要があり（できればループなし）、2 + x観測点の行を比較する必要があるという事実については何も見つかりませんでした。 試してみるデータ： testdata <- structure(list(bird_id = c("20150712_0810_1410_A_1", "20150712_0810_1410_A_2", "20150712_0810_1410_A_4", "20150712_0810_1410_A_7", "20150727_1115_1430_C_1", "20150727_1120_1430_B_1", "20150727_1120_1430_B_2", "20150727_1120_1430_B_3", "20150727_1120_1430_B_4", "20150727_1120_1430_B_5", "20150727_1130_1430_A_2", "20150727_1130_1430_A_4", "20150727_1130_1430_A_5", "20150812_0900_1225_B_3", "20150812_0900_1225_B_6", "20150812_0900_1225_B_7", "20150812_0907_1208_A_2", "20150812_0907_1208_A_3", "20150812_0907_1208_A_5", "20150812_0907_1208_A_6" ), obsPoint = c("A", "A", "A", "A", "C", "B", "B", …

10 r  time  comparison  duplication 

Rで2つのラスターグリッドを配置しようとしています。配置したら、それらを一緒に追加できるようにしたいと考えています。 私はを作ることstackがうまくいくかどうかを確認しようとしました： grid_snap <- stack(habi_sdw, Pop_sdw) そして、私は次のエラーを受け取ります： compareRaster（x）のエラー：範囲が異なります ラスターグリッドには次のプロパティがあります。 show(habi_sdw) # class : RasterLayer # dimensions : 9187, 9717, 89270079 (nrow, ncol, ncell) # resolution : 0.00892857, 0.00892857 (x, y) # extent : -28.83706, 57.92186, -36.02464, 46.00214 (xmin, xmax, ymin, ymax) # coord. ref. : +proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs # …

10 raster  r  clip  snapping  resolution 

経由qgis --code myscript.pyでpythonスクリプトを実行し、すぐに終了したいと思います。私はを使用iface.actionExit().trigger()しています。これは、Pythonコンソールから実行したときにQGISを強制終了しますが、渡したスクリプトに入れたときは強制終了しません--code。 すぐに終了する正しい方法は何ですか？QGIS 2.0.1を実行しています アップデート：私も試しましたsys.exit()。QGISはそれをキャッチし、次のウィンドウをポップアップします。 An error occured during execution of following code: execfile('myscript.py') Traceback (most recent call last): File "", line 1, in File "myscript.py", line 14, in sys.exit() SystemExit 更新：os.kill(os.getpid(), 9)動作しますが、それは汚いハックであり、私はもっと良いものを探しています。

10 qgis  python  pyqgis 

多くの地理データセットをPostGISにインポートしていますが、それらのセットは異なりSRIDます。（あるものEPSG:3857、あるものEPSG:4326、あるもの）。 追加を作成したいgeometry columnです。the_geom_mercatorを使用しSRID EPSG:3857、元のgeom列を元の列のSRIDままにします。 PostGIS関数でこれを行うにはどうすればよいですか？

10 postgis  coordinate-system  srid  mercator 

属性を持つシェープファイルをクライアント用の.kmlに変換しました。次に、.kmlと属性がGoogle Earth内で編集されました。更新された.kmlを受け取ったので、.shpに戻す必要があります。問題は、すべての属性がHTML形式になっているため、Google Earthにポップアップ情報ボックスが表示されることです。.shpに変換しようとすると、属性が（Googleで作成された）「説明」フィールドにまとめられます。そして、私の属性と情報はすべて1つのセルにまとめられます。「説明」フィールドからデータを抽出する方法はありますか？私がそれをしたとしても、htmlフォーマットを取り除き、それを使用可能にする簡単な方法はありますか？ これは、私がいくつかのフォーラムで見つけた「Fusion Table」方法を使用して得られるものです。左側にhtmlが表示されています。表示されているよりもさらに進んでいます。 ArcGIS Desktop 10.2を使用しています。

10 arcgis-desktop  shapefile  kml  convert  google-earth 

ラインストリングをポイントでカットする最良の方法をチェックしていました。 シナリオは次のとおりです。たくさんの道​​路があり、交差点でカットされたセグメントが必要です。 わかった 線ストリング（ポイントでカットされていない完全な）テーブル st_intersectionポイントテーブル 交差点テーブルでカットされた独立したラインストリングセグメントを取得する必要があります。 私はPostGIS関数を使用していて、いくつかのアプローチを見つけましたが、それらのすべてが何らかの問題を引き起こしています。 これは私がすでにテストしたものです： 1 ラインテーブル：1行、1200行のst_memunionポイントテーブル：1700行（ポイント） 悪い点：本当に多くの時間とメモリのフラッシュが必要です。同時に複数のテーブルを作成できないため、メモリはそれを処理できません。そしてその結果は汚くて厄介です。私に正しい行番号を与える代わりに、後でそれをクリーンアップする必要があります（ここで説明されている交差点でのラインの分割） CREATE TABLE lines_with_messy_result AS ( SELECT ((ST_DUMP(ST_SPLIT(a.geom,b.ix))).geom) as geom FROM st_union_lines a INNER JOIN lots_of_points b ON ST_INTERSECTS(a.geom, b.ix) ); --then need to clean this up create table lines_segments_cleaned as ( SELECT DISTINCT ON (ST_AsBinary(geom)) geom FROM lines_with_messy_result ); …

10 postgis  postgresql  openstreetmap  splitting  streets 

入力グラフの折れ線グラフの計算に興味があります。後で他のGRASSネットワーク分析ツールを使用するので、できればGRASSを使用します。GRASSドキュメントに折れ線グラフへの参照が見つかりませんでした。検索で使用する別の用語や、そのようなグラフを作成する非表示の関数はありますか？ コンテキスト： 無向グラフGの線グラフは、Gのエッジ間の隣接関係を表す別のグラフL（G）です。 出典：http : //en.wikipedia.org/wiki/Line_graph

10 grass  network  graph 

私は現在、等時性と基礎となるアルゴリズムの分野で働いています。現在問題となっているのは、等時線自体の場合の計算ではなく、結果の視覚化です。 私の等時線アルゴリズムの結果は、ポイントとエッジです。実際、私には有効なソリューションがありますが、3873エッジと1529ノードの場合、時間がかかるようです（2015 Core i7 CPUとかなり高速なSSDを含むLenovo T440sラップトップでは約2.0秒）。秒の代わりに、msec :-)のようなものがもっと欲しいです。 たぶん誰かが、到達可能な領域を視覚化するポリゴンを構築するのに必要な計算時間を短縮するのを手伝ってくれるでしょう。 しかし、待ってください...まず最初に！ これは、私の等時線の計算結果であるエッジの視覚化です。 これらのエッジはPostGISデータベーステーブルに格納され、単純な線ストリングです。 ユーザーに見せたいのは次のようなもの です。画像の最南端と最東端にある切断された領域に注意してください。これらは個別の領域として描画する必要があります（したがって、ここではマージできません:-)） 現在、私はこのクエリを使用しています： SELECT ST_AsGeoJson(St_Transform(ST_Multi(ST_Collect(polygons)), 4326)) AS coverage FROM ( SELECT ST_MakePolygon(ST_ExteriorRing(ST_GeometryN(segments, generate_series(1, ST_NumGeometries(segments))))) AS polygons FROM ( SELECT ST_Union(ST_Buffer("GEOMETRY", 20, 'quad_segs=2')) AS segments FROM my_edges AS a ) AS b ) AS c 私はすでにいくつかの実験を行っており、多くのドキュメントも読みましたが、より良い解決策を見つけることができません。 私の目には、大きな問題はST_Unionの使用です（ドキュメントで述べられているように、この関数は遅くなる可能性があります）。非常に興味深いのは、それをST_Collectで置き換えると、ST_Bufferの計算が遅くなるように見えるため、次のクエリのオールインオールがさらに長くかかりますが、エッジ間の領域は埋められません（ラインの周りにバッファーが作成されるだけです） ）： SELECT ST_AsGeoJson(St_Transform(ST_Multi(ST_Collect(polygons)), …

10 postgis  postgis-2.0  polygonize 

時々、興味のあるポイントを表示するためにマップブックを作成する必要があります。通常のメッシュを使用してページを作成する最初のステップ： a）単一のポイント（25ページなど）が端にあるページがいくつかあり、b）ページが多すぎるため、このソリューションは好きではありません。 最初の問題はコードを使用して簡単に修正できます。ページ範囲の長方形を関連するポイント範囲の中心に移動します。 私はまだそれが好きではありません、ページ数は同じままなので非常に込み合っています。それらはすべて、レポートの複数のコピーで実際のA3用紙ページになることに注意してください。 そこで、ページ数を減らすコードを作成しました。この例では、45から34です。 これが達成できる最良の結果かどうかはわかりませんが、 すべてのポイントをキャプチャするために指定されたサイズの長方形の数を最小限に抑えるためにポイントをシャッフルするための最良の戦略（疑似コード、パブリケーション、Pythonライブラリ）は何ですか？確かに、誰かがゲーム理論、軍事芸術、漁業でそれを発見しました これは元の質問の更新です： これは、必要な実際の範囲とページサイズを示しています。 164ページのうち10ページを表示する拡大ズーム： サンプルポイントフィーチャクラス 長方形のサイズは、制限内に収まるとすぐに変更できます。

10 arcgis-desktop  arcpy  algorithm  data-driven-pages  maximum-coverage-problem 

QGISの「ラベルの移動」ツールに加えて、PostGISビューを使用して動的な引出線を作成しようとしています。 CREATE VIEW leader_line AS SELECT gid, ST_MakeLine(geom, ST_SetSRID(ST_MakePoint(xcord_label, ycord_label), SRID))::geometry(linestring, SRID) AS geom FROM point WHERE xcord_label IS NOT NULL; これはすべてのラベルで問題WHERE ST_X(geom) < xcord_labelなく機能しますが、ラベルの見栄えが悪くなりWHERE ST_X(geom) > xcord_labelます。 誰かがラベルの引き出し線を正しく配置する方法を知っていますWHERE ST_X(geom) > xcord_labelか？ラベルのxmax座標を参照する方法はありますか？

10 qgis  postgis  labeling  coordinates  dynamic-layer 

リーフレットオーバーレイペインを使用してSVGを「上に」配置したリーフレットマップがあります。D3ブラシをオーバーレイペインに追加して、別のグラフを更新します。D3.brush（）とリーフレットの例の同時例を探しましたが、見つからないようです。 私が見つけた最も近いものだった。このhttps://github.com/mbostock/d3/issues/1321と、このhttps://www.mapbox.com/mapbox.js/example/v1.0.0/timeline-scaled-markers/（これはマップではなくチャートをブラッシングしています）。 私は疑問に思っていました： 誰もがどんな例でも知っていた、または これが可能であるか、合理的に簡単であるか誰でも知っていた、または 可能であれば、私が遭遇する可能性のある種類の問題について、誰かがアドバイスをしました。

10 leaflet  d3 

QGISはGDALを使用しているため、PDFとしてエクスポートされたマップは地理空間PDFになりますか？QGISから携帯電話のAvenza PDF Mapアプリにマップをインポートしたいと思っています。以前の投稿では、地理空間PDFおよびGeoTiffとしてエクスポートすることについて言及していましたが、正確な答えやプロセスはありませんでした。

10 qgis  geotiff-tiff  pdf  geospatial-pdf  avenza 

昨年5月にArcGIS Proベータプログラムの発表というタイトルのEsriブログ投稿へのコメントとして、mahenderrajという名前の誰かが尋ねました： ArcGIS ProからArcGISサーバーへのマップサービスとしてマップを共有することはできないと思います。そのためには、通常のArcMapデスクトップ製品を使用する必要があります。 私はこれに対する答えを見たことがないようで、ArcGIS Pro 1.0インターフェースでは、ArcGIS OnlineまたはPortal for ArcGISにマップを公開するのは簡単に思えますが、ArcGIS Serverは言及されていません。 ArcGIS ProからArcGIS Serverにマップサービスとしてマップを公開することはできますか？

10 arcgis-server  map-service  arcgis-pro