地理情報システム

プロットしたいSpatialPolygonsまたはSpatialPolygonsDataFramesがあります。Rの3つのプロットシステム（ベースグラフィックス、ラティスグラフィックス、およびggplot2）を使用して色を変更するにはどうすればよいですか？ サンプルデータ： library(sp) Srs1 = Polygons(list(Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))), "s1") Srs2 = Polygons(list(Polygon(cbind(c(5,4,2,5),c(2,3,2,2)))), "s2") SpDF <- SpatialPolygonsDataFrame( SpatialPolygons(list(Srs1,Srs2)), data.frame( z=1:2, row.names=c("s1","s2") ) ) spplot(SpDF, zcol="z")

25 polygon  r 

ArcMap 10に読み込まれたフィーチャクラスには約12個のポリゴンがあり、すべて地理WGS 1984にあります。 そのフィーチャクラスの各ポリゴンの各頂点に関連付けられた座標を簡単に取得するにはどうすればよいですか？ 理想的には、スプレッドシート形式でうまく表にしたいです。

25 arcgis-desktop  arcgis-10.0  arcmap  polygon  vertices 

QGISプロフェッショナル認定プログラムについてのうわさを聞いたことがあります。 https://hub.qgis.org/wiki/17/Certificationでドキュメントを見つけましたが、他の記事やWebリンクを見つけることができませんでした。 そのような認証のためのトレーニングとテストを提供する努力がまだあるかどうか誰もが知っていますか？ その場合、どのように/どこでプロセスを開始しますか？ 私は米国にいるので、米国内またはリモートで利用可能なプログラムに関する情報は特に役立ちます。 最も徹底的な回答が提出されてから1年半が経ちました。更新があったのだろうか？

25 qgis  education  training  certification  gis-professional 

この質問は2011年に投稿され、QGISで凡例エントリに複数行が許可されるようになったと誰かが言っていますが、その方法については誰も言及していません。wein 2.8.1を実行していますが、テキストを手動で使用する以外に、複数行で凡例エントリを取得する方法がわかりません。

25 qgis  legend  print-composer 

ユニットが何から返されているフロートののかしらST_Distance。 ドキュメントには次のように書かれています： ...投影された単位の2つのジオメトリ間のデカルト最小距離（空間参照に基づく）。 これらの予測単位は何ですか？ ジオメトリは次のフィールドに保存されますgeometry(Point,4326)。

24 postgis  postgis-2.0 

無効なポリゴンジオメトリから発生する「TopologyException：Input geom 1 is invalid」自己交差エラーは広く議論されています。しかし、R機能だけに依存する便利な解決策をWeb上で見つけられませんでした。 たとえば、私はここでmap("state", ...) Josh O'Brienの素敵な答えを次の出力から「SpatialPolygons」オブジェクトを作成することができました。 library(maps) library(maptools) map_states = map("state", fill = TRUE, plot = FALSE) IDs = sapply(strsplit(map_states$names, ":"), "[[", 1) spydf_states = map2SpatialPolygons(map_states, IDs = IDs, proj4string = CRS("+init=epsg:4326")) plot(spydf_states) この広く適用されているデータセットの問題は、以下のポイントで自己交差が発生することです。 rgeos::gIsValid(spydf_states) [1] FALSE Warning message: In RGEOSUnaryPredFunc(spgeom, byid, "rgeos_isvalid") : Self-intersection at or …

24 r  polygon  rgeos  self-intersection 

グライダーパイロットからは、一定の間隔でGPSを修正するという形で多くのフライトデータを取得しています。飛行経路を分析し、グライダーのパイロットがサーマルを見つけたときに行う「旋回」の開始と終了を検出したいと思います。 理想的には、アルゴリズムは線上の開始点と終了点を与え、1つの「円」を定義します。これらのポイントは、GPS修正の1つと等しくなる可能性があり、補間する必要はありません。 飛行経路に沿って歩き、回転数を確認し、グライダーが旋回しているかどうかを判断するための基準を設定することができました。 PostGIS拡張機能を備えたPostgreSQLを使用しているため、この問題に対するより良いアプローチがあるかどうか興味がありました。私はすでに2つの線分の角度を計算する手順を持っています： CREATE OR REPLACE FUNCTION angle_between( _p1 GEOMETRY(PointZ,4326), _p2 GEOMETRY(PointZ,4326), _p3 GEOMETRY(PointZ,4326) ) RETURNS DECIMAL AS $$ DECLARE az1 FLOAT; az3 FLOAT; BEGIN az1 = st_azimuth(_p2,_p1); az3 = st_azimuth(_p2,_p3); IF az3 > az1 THEN RETURN ( degrees(az3 - az1)::decimal - 180 ); ELSE RETURN ( degrees(az3 - …

24 postgis  postgresql  algorithm  trajectory 

FWToolsのogr2ogrを使用してESRI Shapefileに変換するGeoJsonファイルがあります。 使ってみた ogr2ogr -F "ESRI Shapefile" filename.shp geojsonfile.json OGRGeoJSON 出力を得た エラー6：シェープファイルレイヤーにStringList型のフィールドを作成できません。エラー1：シェープファイルタイプに非ポリゴン（POINT）ジオメトリを書き込もうとしました。エラー1：レイヤーOGRGeoJSONの翻訳に失敗した後、翻訳を途中で終了する 提供するGeoJSON入力は、標準のgeojson形式です。 geojsonのコンテンツは{共有するために、features配列の一部のみを取りました} {"id":1223, "properties":{"name":"3","level":2}, "map_version":11,"type":"FeatureCollection","x_id":1234,"obj_type":"LevelGeometry", "features":[ { "my_area":[131.40093253528894,33.58879853190332,6.515867411008615,3.490549810987222,-0.6237271428108215], "id":1984791, "location":{"type":"Point","coordinates":[131.4009325,33.5887985]}, "properties":{"display_name":"xyz"}, "type":"Feature", "geometry":{"type":"Polygon","coordinates":[[[131.4009508,33.5888314],[131.4008895,33.5887946],[131.4009148,33.5887653],[131.4009761,33.5888020],[131.4009508,33.5888314]]]}, "obj_type":"Geometry" } ] }

24 geojson  ogr2ogr  fwtools 

QGIS 2.18を使用しています。地図上でナビゲーション用のセクターライトを作成する必要があります。 ライトセクターデータは、becon_id、開始度、終了度、および色を与えるフィールドとして、500以上のブイ、ビーコン、灯台がマップに表示される必要があるシェープファイルにあります。各ビーコンには、多くの行があり、それぞれが1つの明るいセクター（たとえば、白いセクター）を記述します。 最終的な結果は次のようになります。適切な色のライトセクター、カラーフィールド（RGW）の文字としてマークされた色、およびブイ/ビーコン/灯台から100mから1000mの点線。 これはほとんどの場合、ルールベースのシンボルとして作成されるはずですが、Pythonが必要ですか？ 以下は、114〜154度の緑のセクター、154〜168度の白のセクター、168〜237度の赤のセクター、緑セクター237〜314度、白色セクター314〜320度、赤色セクター320〜337度（何らかの理由で、0は北ではなく南）：

24 qgis  python  labeling  symbology  style 

PostGISで新しいデータベースを作成したいので、現在のデータベースが使用されている間にそれをロードできます。ドキュメントによると PostGISのパッケージ化されたディストリビューション（特にPostGIS 1.1.5以降のWin32インストーラー）は、PostGIS関数をtemplate_postgisと呼ばれるテンプレートデータベースにロードします。PostgreSQLインストールにtemplate_postgisデータベースが存在する場合、ユーザーやアプリケーションは、単一のコマンドを使用して空間的に有効化されたデータベースを作成できます。 私の場合、これはそうではないようです： $ createdb -T template_postgis my_spatial_db createdb: database creation failed: ERROR: template database "template_postgis" does not exist 過去には、プライマリgisデータベースをコピーしてから、すべてのテーブルの内容を削除することをいじりました。より良い方法がなければなりません。誤ってドロップした場合はどうしますか？

24 postgis  postgresql  database  create 

地理的特徴（ライン、ポリゴン、およびそれらの同等のマルチパート）がantimeridian（±180°経度）にまたがり、GeoJSONとしてクライアントWebアプリケーションに送信および受信する必要がある場合、それらを保存するためのベストプラクティスは何ですか？ 私は、Postgres / PostGISデータベースのサポートを受けて、サーバー側のWeb APIで作業を開始し、熱帯低気圧の履歴と風の半径を予測します。太平洋の多くのサイクロンは、時には寿命の間に複数回、反時計回りを通過する不幸な傾向があります。 ニュージーランドのアンティリディアンの近くに住んでいる私は、地域のデータでこの問題に何度も対処し、対処戦略を立てましたが、実際にベストプラクティスと考えられるものを見つけたいと思います。残念ながら、antimeridianでタグ付けされた既存の質問はないため、関連する質問を検索するのは困難です。私がこの問題に苦労しているのを見た質問はすべて、非常にアプリケーション固有のアドバイスを求めているようです。この質問では、境界のない地球に広がるGeoJSONポリゴンの場合の反時計回りについて簡単に説明します。この質問は、私が尋ねていることにかなり近いものです。 空間的データベースに歴史的サイクロンと予測サイクロンを保存する必要がありますが、反時計回りに問題があると予想しています。たとえば、緯度/経度(0,179)で始まり、で終わる線(0,-179)は、その方向に関してあいまいです。それは、反時計回りを横切る短い経路を取るか、惑星全体を「包む」かです。このようなパスを空間データベースにどのように保存する必要がありますか（具体的には私はPostGISを使用していますが、ソリューションが十分に汎用的であることを望みます）。私が持っているいくつかのアイデア： フィーチャジオメトリに変更を加えず、クライアントアプリケーションへの曖昧さをシフトします。 antimeridianを横切るフィーチャを、antimeridianでブレークするマルチパートジオメトリに分割します。（GeoJSON仕様は名前付きCRSをサポートしています。） 仕事以外のグローバル予測などの不連続を持っていない別のサイクロン流域や海洋地域のために サイクロンが惑星全体を巡回することが一度も観察されていないという事実を利用して(90,-90) 、360°の位相でオフセットされた緯度範囲から始まるサイクロンの座標を保存します（他の-180-180°を維持） サイクロンがアフリカの南端の南にある可能性が非常に低いという事実を利用して、経度30°でブレークを使用します（上記のマップのように）。 座標がEPSG 4326の有効範囲を超えて拡張できるようにします。たとえば、反時計回りを通過するすべてのフィーチャに対して、180度以上-180度未満です。 TopoJSONと同様に、デルタエンコーディング（例：で始まり(0,-179)、次の座標は-3緯度西です）。PostGISにデータを保存するときにこれを実装するかどうか、または実装する方法がわかりませんが、これはクライアントアプリケーションにデータを送信するための優れたソリューションです。 何らかの形式のベクトル表記または極座標。（かなり難しくて珍しいようです。） これらのうち、アイデア2〜5は一般的なものではないので好きではありませんが、私の特定のアプリケーションにとって意味があるので好きです。太平洋のデータのみを処理するアプリケーションの場合、それらは非常に理にかなっている可能性があるため、オプションとして完全に割引したくはありません。 アイデア6および7は、Tom MacWrightのブログから取り上げられました。このブログは読む価値はありますが、反時計回りに関しては決定的なものではありません。 アイデア4はGeographicaGS 'GeodesicLinesToGISPythonによって使用されfiona.transform.transform_geom、360°の反時計回りオフセットで使用されています。次に、FionaはOGRを使用してい-wrapdatelineます。これは非常に堅実な先例であり、実際にはかなり一般的だと思います。 ストレージの問題に関連して、そのような機能をクライアントアプリケーションに送信する方法、およびアプリケーションがそこにポストバックされるデータを考慮する方法を検討する必要があります（たとえば、太平洋のサイクロンの予測トラックを変更する人間予報士）。交換形式はおそらくGeoJSONになりますが、そうである必要はありません。 残念ながら、GeoJSON仕様は、反時計回りの問題について明示的ではありません。ウィキペディアからこれ： 多くの地理的ソフトウェアライブラリまたはデータ形式は、世界を長方形に投影します。多くの場合、この長方形は180番目の子午線で正確に分割されます。これにより、180番目の子午線上で単純なタスク（エリアやラインを表すなど）を実行できなくなることがよくあります。いくつかの例： GeoJSON仕様では、仕様で180番目の子午線の処理に言及していないため、180番目の子午線を横切る線の表現は、世界中を回っていると解釈できます。 OpenStreetMapでは、領域（ロシアの境界など）は180番目の子午線で分割されます。 私の読書では、GeoJSONには反時計回りに広がる特徴の特定の標準表現はなく、意図的に曖昧なままになっています（マルチパートジオメトリによって問題が解決される可能性があります）。同様に、OpenStreetMapにはantimeridianにジオメトリ分割がありますが、そのような分割されたフィーチャがマルチパートであるか、実際に個別のレコードであるかはわかりません。 これは、特にバウンディングボックスまたはこの行にまたがる他の空間要求を作成するという観点からだけでなく、入力の解析とサニタイズおよびフィーチャジオメトリの更新においても、かなり問題があります。これが、私が従うことができるベストプラクティスを決定しようとしている理由です。

24 postgis  geojson  api  antimeridian 

ArcGISのカバレッジ、シェープファイル、ジオデータベースで使用される空間データストレージ方法の違いについて疑問に思っていました。カバレッジが最初の形式で、その後にシェープファイル、そしてジオデータベースが続きました。では、これらの新しい形式のシェープファイルとジオデータベースの改善点は何ですか？ 誰かがそれを例で説明してくれたら素晴らしいと思います。

24 arcgis-desktop  shapefile  esri-geodatabase  esri-coverage-format 

MXDには、デフォルトツールのレイヤー、シンボル、構成に関する情報が含まれていることを知っています。50 MBを超えるものがいくつかあり、読み込みに時間がかかり始めています。 MXDのファイルサイズを増やす主なものは何ですか？サイズを劇的に増加させる具体的なものはありますか？ 私がまったく証明していない私のアイデアのいくつかは次のとおりです。 複雑なレイヤーシンボル レイヤーごとに複数のラベルクラス 定義クエリ ソースからレイヤーをプルすると、ファイルパスが長くなります（C：/MyFolder/Data/Data2/SubContainer/YouGetThePoint.shp） レイヤーのスケール範囲 レイヤーをグループに配置する これらのアイデアのいずれかに妥当性はありますか？

24 arcmap  arcgis-desktop  mxd  file-size 

既存のGeocoding API（Google、Mapquestなど）とその機能を比較するための適切な更新されたソースはありますか？ 彼らの言語サポートにもとても興味があります。たとえば、Googleには、JavaScript、Java、Pythonなどのジオコーディングサービス用のライブラリがあります。

24 google-maps  java  geocoding  api  mapquest 

私は現在、デフォルトのシンプルシンボル機能またはポイントオブジェクトの絵文字フォントを使用しています。QGIS用のSVGシンボルライブラリは、含まれているものより少し拡張されていますか？ 私はこのサイトhttp://mapicons.nicolasmollet.comを見つけました。それはいくつかのアイコンを持っていますが、探しているものではありません。これらは、実際にはGoogleマップにより適しています。

24 qgis  symbology  svg