地理情報システム

スナップに関するQGISマニュアルも含め、多くの投稿とチュートリアルを読みました。[設定]> [スナップオプション ] でスナップオプションを設定します。スナップしたい2つのレイヤーを選択し、モードを頂点に設定し、0.00002マップユニットの許容値をチェックして、intを回避します。トポロジ編集を有効にします。 レイヤー1のポリゴンを移動し、近似頂点をレイヤー2の隣接する近似頂点に取得しました。1：500の縮尺では、正しくスナップされているように見えますが、1：1の縮尺にズームインすると、それらは適切に配置されません。私は何かを逃したり、正確にしようとしていますか？ 入力いただきありがとうございます。

13 qgis 

閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか？この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 3年前に閉店しました。 GISの内部とPostGISなどの空間データベースシステムについて、例を使用して背景の前提なしに詳細に説明する良いチュートリアルはありますか？私は特に次のような質問への回答を探しています： GISのジオメトリとは何ですか？ それはどのように表されますか？ 緯度/経度座標が指定されている場合、ST_contain PostGISのようなものを呼び出してポリゴンに緯度/経度座標が含まれているかどうかを確認できる状態にするために、実行する必要がある操作は何ですか？ 私は一般的なリレーショナルデータベースシステムである程度の経験がありますが、GISまたは空間データベースパラダイムのバックグラウンドはありません。

13 postgis  shapefile  geometry  spatial-database  references 

Rの空間オブジェクトが与えられた場合、すべての要素をクリップして境界ボックス内に配置するにはどうすればよいですか？ 私がやりたいことが2つあります（理想的には両方を行う方法を知っていますが、どちらも私の現在の問題に対する許容可能な解決策です-ポリゴンシェープファイルを米国本土に制限します）。 境界ボックス内に完全にはない各要素をドロップします。 これbbox()<-は論理的な方法のようですが、そのような方法は存在しません。 無限小要素（ポリゴン、ラインなど）が境界で切り取られるように、真のクリップ操作を実行します。 sp::bbox割り当て方法を欠いているので、私が作ってみた唯一の方法は、使用するのだろうoverかgContains/ gCrossesSpatialPolygonsと連携して、新たなバウンディングボックスの座標を持つボックスを含むオブジェクト。次に、ポリゴンオブジェクトをクリップするとき、含まれているものとクロスされているものを把握し、それらのポリゴンの座標を変更して、ボックスを超えないようにする必要があります。または何かのようなgIntersection。しかし、確かに簡単な方法がありますか？ バウンディングボックスには多くの問題があり、通常、関心領域を定義するポリゴンへの空間オーバーレイが望ましいことを知っていますが、多くの場合、バウンディングボックスはうまく機能し、よりシンプルです。

13 r  clip  extents 

次のような垂直ラベル（新しい行の各ラベル）を作成する方法： 1998 150 pipe name 水平ではありません： 1998 150 pipe name

13 qgis  labeling 

国勢調査のフローデータ、ジョブ密度分布、およびトランスポートネットワークに基づいて、通勤者が通勤中にたどる可能性のある経路を分析しています。 現時点では、分析にpgRoutingおよびQ​​GISのpgRoutingアドオン（ここで説明）を使用しています。 しかし、私は完全な複製であることを取るすべてのステップが欲しい（これに関する優れた記事を参照して、ここでは）、私は、コマンドラインからの私の分析を行い、将来の参考のためにスクリプトファイルを保存することを好むだろう。 R-geoパッケージは、無料で、軽量で、安定しているため、この要件を十分に満たします。小さなスクリプトファイルをダウンロードすると、私が行ったことをすべてやり直すことができるはずです（Rpubsでのこの無関係な例については、こちらを参照してください））。 これを行うために特定したRパッケージは次のとおりです。 2012年12月に公開されたgdistanceは有望に見えますが、osmトランスポートネットワークとの互換性の兆候はありません e1071はかなり抽象的で、おそらくやりたいことに対してやり過ぎだ spatstatパッケージのpairdist.lpp関数は有望に見えますが、ドキュメントはほとんどありません これらに飛び込む前に、同様のジレンマに直面したかもしれない人々に尋ねる価値があると考えました。 Rでルーティング解析を行うことは可能ですか、それともPostGISを適切に学習してpgRoutingですべて行うのが最善でしょうか？

13 pgrouting  r  routing  osm2pgrouting 

たとえば、次のようにPythonプラグインを使用して追加できQgsVectorLayerましたQGIS。 vl = self.iface.addVectorLayer(uri.uri(), layerName, self.dbConn.getProviderName()) しかし、どうすれば削除できますか？

13 qgis  python  vector  qgis-plugins 

QGISポイントベクターレイヤーのポイントフィーチャ用にカスタムラベルを作成しようとしています。それを行うためにPythonプラグインを作成したいと思います。ラベルの後ろに色付きの長方形を表示する必要があります（または！）。ラベルは長い文字列であり、ラベル付けエンジンの式と条件を使用して作成されます。 最も簡単な方法は、ユーザーがアクティブレイヤーのラベルを好きなように設定し、プラグインがレイヤーのフィーチャからラベルテキストを取得し、テキスト文字列を使用してボックスを作成できるようにすることだと思いました。私の問題は、QGIS APIでこれらのテキスト文字列にアクセスする方法が見つからないことです。ラベル付けエンジンによって生成されたラベル文字列はQGISにどのように保存され、どのようにして取得できますか？

13 qgis  pyqgis  labeling 

ユーザーがテーブルから表形式の統計を取得できるQGISツールはありますか？テーブルを開いて、興味のあるフィールドを選択し、フィールドヘッダーを右クリックして、これに似たもの（ArcGISから取得したスクリーンショット）を取得できると期待していました。 具体的には、ArcGIS Frequency and Summary Statisticsツールに代わるQGISを探しています。これは、表形式のデータについて取得できる他の多くの統計情報を提供します。 質問した直後に、この機能を見つける方法を説明するこのQGISチュートリアル（セクション4.2を参照）を見つけました。Vector > Analysis Tools> Basic statistics。また、そのすぐ横に[ベクター]> [分析ツール]> [一意の値の一覧表示]があります。これにより、一意の値のリストが取得されますが、各一意のアイテムの頻度（カウント）の追加の列も確認したいと思います。この同様のツールは「すぐに使用可能」ですか？

13 qgis  statistics 

私はインド太平洋に生息する種（〜20〜-65度の経度）の凸包を構築しています。私が遭遇した問題は、180度の経度線の両側にある種の場合、再構築された船体が、レイヤーの両端を構成する非常に小さなポリゴンにポイントを含めるのではなく、レイヤー全体に広がることです。この動作の原因は、この線がレイヤーの端を表しており、QGISがベクトルマップをラップしていないことです（下の画像を参照してください...私の「評判」は今ではそれを含めるのに十分です、ありがとう）。 いくつかの考えとグーグルの後、これらは私が遭遇した問題への可能なアプローチの一部です。ただし、そのメリットや実装についてはあまりよくわかりません。 中央子午線を0（グリニッジ）から150度に変更します。これをどのように行うことができますか？これは単にマップの表示を変更するだけですか（それ自体が有用です）、それとも正しい凸包の構築に貢献しますか？ 座標系を+/- 180度形式から0〜360度形式に変更します。どうやら、これらの線に沿った何かは、PostGISの「ST_Shift_Longitude」関数を使用して実現できます。これは、負の経度に360度を追加します。これを行う簡単な方法はありますか？また、この関数は180度中心のマップを生成します。これは世界の終わりではありませんが、150度中心のマップの方が優れています（〜インド太平洋の中心）。 お時間を割いてお世話になり、長い投稿をありがとうございます。

13 qgis  convex-hull 

このマップをロビンソン図法で投影したい： library(ggmap) world <- map_data("world") ggplot() + geom_path(data = world, aes(long, lat, group = group)) そして、投影を「ロビンソン」に変更したいと思います（前の質問への回答からのアドバイスに従ってください：ウィキペディアからの地球規模の気候地域マップはどの投影を使用していますか？ このプロジェクションのデフォルトの実装を見つけるのに苦労しました。proj4ライブラリを使用するために次のことを考えました。 library(proj4) robinson <- project(cbind(world$long, world$lat), proj = "+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs") 私は以下を含む多くのアプローチを試みました： # using ggmap::get.map() get_map("world", projection = mapprojection(robinson)) # using ggplot2::coord_map coord_map(projection = robinson) # and sp::coordinates: library(sp) …

13 coordinate-system  r  proj  ggmap 

OpenLayers 2.12を使用して、HTTP基本認証が有効になっているサーバーからWMSレイヤーを表示しようとしています。 JavaScriptコードのURLパラメーターにユーザー名とパスワードを入力して、認証を処理しようとしました。レイヤー作成の例： myLayer = new OpenLayers.Layer.WMS('background', 'https://username:password@ws.nls.fi/rasteriaineistot/image?', { layers: 'background', bbox: '-380188,6249943,1347312,8226943' }, { displayInLayerSwitcher: true, isBaseLayer: false, projection: 'EPSG:3067', visibility: true }); もちろん、資格情報はJavaScriptコードに保存され、すべてのブラウザーで機能するわけではないため、これは安全ではありません。Internet Explorer 8では、OpenLayers.jsを指すセキュリティエラーが発生し、マップの表示は一切拒否されます。Firefox 13は、キャンセルできる認証ダイアログをポップアップ表示します（その後、マップが正しく表示されます）。Chrome 23では、認証は問題なく動作するようです。 HTTPベーシック認証をURLでエンコードし、例のようにOpenLayersに渡すことで、クロスブラウザー方式でHTTP基本認証を処理できないことを確認できますか？ HTTP基本認証を処理してユーザーに透過的に機能する（認証ポップアップが表示されない）別の方法を提案できますか？この問題を回避するには、おそらく何らかの種類のプロキシサーバーを使用してください。

13 openlayers-2  javascript  authentication 

データ自体とともに地理空間メタデータを非ESRI空間データベースに格納する（したがって、ダンプ時に移動できる）既存の標準的な普遍的なアプローチ（および、できれば管理をサポートする既存のツール）があります。 テーブルとリレーションに単純に依存するため、PostGIS、Spatialite、Oracle、SQL Serverなどのデータベースで使用できるアプローチを特定したいと考えています。ここで、メタデータとは、データに関する説明情報（つまり、US FGDCまたはISO 19139地理空間メタデータタイプ情報）-BBOXおよび内部データではありません。 ESRIユーザーは、ファイル（Shapefiles）であろうとジオデータベースであろうと、データを普遍的に記述および付随できるいくつかのXML形式を持っています。ただし、ESRIソフトウェアを使用しない場合、既存のオプションは何ですか？はい、もちろん、独自のテーブル、データ構造などを設計できます。しかし、なぜ確実に存在するはずのホイールを再発明するのでしょうか。 更新： Geonetwork（または必然的にサーバーに関係するもの）のような複雑なアーキテクチャコンポーネントは、まさに避ける必要があるものです。また、メタデータは、個別のデータベースとしてではなく、データとともに存在します。要件は以下のとおりであり、最初に述べておかなければなりませんでした。 システム要件：1.アーキテクチャは、QGISとSpatialiteデータベースのみを必要とします-一部は、組織がサーバー上で何かを実行するほど洗練されておらず、何かを購入したり、ビルド/デプロイしたりするお金がないためです。 機能要件：1.データは多くの人々に容易に配布されなければならず、ドキュメントはデータから簡単に分離されてはなりません。つまり、データが何であり、なぜ作成されたかを常に把握できるように、それらは一緒に配布されるべきです。など-データがある場合、ドキュメントがあります。2.データ自体と同様に、メタデータのドキュメントは、直感的なデスクトップツールを使用して、非技術スタッフが簡単に編集および保守できる必要があります。 ユースケース：1.ボビー学生スチューデント（およびGISの学習）は、調査の一環として監視サイトのデータを作成します。2.ボビーは、使用した入力、処理ステップの説明、および他の人がデータの系統を理解するのに役立つその他の情報を記録します。3.ボビーは実際の仕事を得て退職し、彼のデータはCD-ROMにバックアップしたままにします。4. 2年後、誰かがデータを見つけて、データ内のドキュメントを読むことができるため、非常に有用であると判断します。 あなたが洗練された組織から来たなら、あなたは言うだろう。しかし、関連するシナリオは実際、私の世界では非常に一般的です。

13 spatial-database  metadata  standards 

閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか？この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 4年前に閉鎖されました。 ArcGISのさまざまなシナリオで長さと面積を計算する方法を理解しようとしています。フィーチャクラスのフィールドで答えが見つからない理由はわかりませんが、何かを理解しておらず、いくつかの歴史があることを知っていない限り、正確な答えを見つけることができません。疑問符の記入を手伝ってもらえますか？または、なぜこれについてすべて間違っているのか教えてください。） GCS =地理座標系PCS =投影座標系 すべてのリンクは10.1ヘルプドキュメントにあります- フィーチャクラスのshape_lengthおよびshape_areaフィールド a。GCS-？ b。PCS-単純な平面の使用 c。シェープファイルを除き、常に自動更新されますか？はい ArcMap測定ツール a。GCS-デフォルトの測地線、代替案はLoxodromeとGreat Ellipticですが、平面ではありません。面積計算は利用できません！ b。PCS-デフォルトの平面、代替の測地線、ロクソドローム、グレートエリプティック http://resources.arcgis.com/en/help/main/10.1/index.html#//00s500000022000000 属性テーブル計算機 a。GCS-利用不可 b。PCS-プラナー http://resources.arcgis.com/en/help/main/10.1/index.html#//005s00000027000000 フィールド計算ツール（データ管理ツールボックス） a。GCS-測地線線形、利用可能な面積は疑わしい b。PCS-プラナー http://resources.arcgis.com/en/help/main/10.1/index.html#//00170000004m000000 バッファリングツール（およびその他のツール） a。GCS-測地線 b。PCS-平面またはGCS出力を指定 http://resources.arcgis.com/en/help/main/10.1/index.html#//000800000019000000 Javascript APIクライアントサイド a。GCS-測地線の面積と長さの関数 b。PCS-Webメルカトルから地理への変換（またはジオメトリサービスの使用） http://help.arcgis.com/en/webapi/javascript/arcgis/help/jsapi/namespace_geometry.htm Flex APIクライアントサイド a。GCS-測地線の面積と長さの関数、「長さ[または面積]は、カスタムの円柱等面積投影法を使用して計算されます」。これはjavascript apiには記載されていません!! b。PCS-Webメルカトルから地理http://resources.arcgis.com/en/help/flex-api/apiref/com/esri/ags/utils/GeometryUtil.htmlに変換できます ArcGIS Server REST API-ジオメトリサービス a。GCS-測地線 b。PCS-プラナー http://help.arcgis.com/en/webapi/javascript/arcgis/help/jsapi/geometryservice.htm 別の質問、測地線測定とは正確には何ですか？スフェロイドの3Dトリガー式を意味すると思いました（ハバーシン？）。また、面積の計算に使用するのは遅すぎるので、等しい面積の投影が使用されるのはなぜですか？ 長さと面積を決定するときの別の質問は、同じデータ、回転楕円体を使用した測地線計算よりも等面積投影法の方が正確ですか？そして、なぜ簡単に？

13 arcgis-10.1  arcgis-desktop  coordinate-system 

編集済み： 質問を説明したい。「ポイントA」にいて、「ポイントB」に行きたいとします。このポイントは「at_2po_4pgr」テーブルにはないため、ソース/ターゲットノードではありません。次に、ポイントAおよびB（緑色のポイント）の最も近いノードを検索します。この後、グリーンポイントIDを使用してshortest_pathコールを実行し、「オレンジ」パスを取得します。しかし、最初のケースで実際のパスコスト（距離）を取得するには、「offsetA」を減算し、2番目のケースでde「offset B」を追加する必要があります。赤点と緑点の間の距離を計算するには、次のクエリを実行します。 SELECT * FROM st_distance( ST_GeomFromText('POINT(-3.6963314 42.3498066)',4326), ST_GeomFromText('POINT(-3.6954276 42.3479634)',4326))。 オフセットを追加または減算するタイミングを知るにはどうすればよいですか？ 私の英語でごめんなさい！

13 postgis  openstreetmap  pgrouting  osm2po 

ベクターの背景マップとして使用する必要があるエリアのオルト写真があります。生のTIFF形式で、対応するTFWワードファイル-180 GBのデータを含む2500ファイル（各71,5 MB）として取得しました。座標系はローカルであり、ベクターと一致します（EPSGコードはありませんが、「32805」と命名し、適切な定義を付けました）。 MapInfoでデスクトップで使用する場合は、それらをECW（MapInfoに付属のツールを使用）に変換し、2500個のTIFFファイルを開くのはやり過ぎなので、71個のファイルのみを保持するように大きくします。49個のTIFF（7x7）を1つのECW（35000x35000pixels-最大は約200MB）に統合しました。MapInfoでは非常に高速で非常に高速です。 今、私はかなり混乱しています-GeoServerでそれらを提供する方法は？ 比較のために、1つのTIFFと1つのECWを公開しています。ECWはブラウザーのプレビューではるかに高速です（ECWサーバーのライセンスの問題は認識していますが、問題ではないはずです）。「GeoServer on steroids」というプレゼンテーションを見つけ、ImageMosaic、ImagePyramid、リタイリング、概要の追加などについて読んでみましたが、非常に有益でしたが、どうすればいいのかわかりません。 私の質問は：どうすればいいですか？モザイクまたはピラミッド、そして肯定的な答えの1つである場合、私はあなたのアドバイスまたはいくつかのヒントが必要です。ディスク容量のためにECWになりたいので、サーバーに180GBのTIFFを保持する必要はありません。 データは、ピーク時に最大20人のユーザーが接続されたLANを介して提供されます。SQLサーバーにはそれほど大量のデータはありません。他の情報を見逃してすみませんが、必要に応じて送信します。 Geoserver 2.1.4、Windows 7 32ビット、2GBシステムメモリ、（1.7.0_09（Java HotSpot（TM）Server VM）、ネイティブJAI +ネイティブJAI ImageIO = true Original TIFF gdalinfo D:\75720-47970.tif Driver: GTiff/GeoTIFF Files: D:\75720-47970.tif D:\75720-47970.tfw Size is 5000, 5000 Coordinate System is `' Origin = (7572000.000000000000000,4797500.000000000000000) Pixel Size = (0.100000000000000,-0.100000000000000) Metadata: TIFFTAG_SOFTWARE=Adobe Photoshop 7.0 TIFFTAG_DATETIME=2006:10:09 13:02:57 …

13 geoserver  geotiff-tiff  image-pyramids  ecw  tiling