地理情報システム

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

1
ArcMapの[高度な編集]ツールバーの[スケール]ツールのようなポリゴンサイズ変更用のArcPyツールはありますか?
ArcGIS 10.3用のPythonスクリプトを書いています。Scale toolArcGISインターフェイスについては知っていますが、そのようなアークピーなコマンドは見つかりません。それが存在します? あなたが写真で見ることができるように、Scale tool異なる作品Buffer tool-それは元の多角形の形を変えます。質問は次のとおりです。 Scale toolarcpy を使用して(ArcGISインターフェイスから利用可能)使用できますか?

3
PostGISが適切にフォーマットされた住所をジオコーディングするのにどれくらいの速度が必要ですか?
PostGISが適切にフォーマットされた住所をジオコーディングするのにどれくらいの速度が必要ですか? PostgreSQL 9.3.7とPostGIS 2.1.7をインストールし、国のデータとすべての州のデータをロードしましたが、ジオコーディングが予想よりはるかに遅いことがわかりました。期待値を高く設定しすぎていませんか?1秒あたり平均3つの個別のジオコードを取得しています。約500万回行う必要がありますが、これを3週間待ちたくありません。 これは巨大なR行列を処理するための仮想マシンであり、このデータベースを横にインストールしたため、構成が少し間抜けに見えるかもしれません。VMの構成の大きな変更が役立つ場合は、構成を変更できます。 ハードウェア仕様 メモリー:65GBプロセッサー:6 lscpuはこれを私に与えます: # lscpu Architecture: x86_64 CPU op-mode(s): 32-bit, 64-bit Byte Order: Little Endian CPU(s): 6 On-line CPU(s) list: 0-5 Thread(s) per core: 1 Core(s) per socket: 1 Socket(s): 6 NUMA node(s): 1 Vendor ID: GenuineIntel CPU family: 6 Model: 58 Stepping: 0 CPU …

3
OSGeo4wシェルスクリプトの実行時のqgis​​.coreのインポートエラー
この投稿とともに、QGISの外部のOSGeo4w Shellでスクリプトを実行しようとしました。しかし、次のエラーが表示されます。 ImportError:qgis.coreという名前のモジュールはありません また、次の投稿を読んで、さまざまなモジュールをインポートしようとしましたが、役に立ちませんでした: QGIS pythonコンソールの外でsextanteアルゴリズムを実行する方法は? QGIS用のシンプルなPythonスクリプトを外部から実行するには(Sublime Textなど)? PyQGISを使用してスタンドアロンPythonスクリプトを作成しますか? 以下は、グリッドを作成し、その上にポリゴンシェープファイルをクリップする簡単なスクリプトです。 注:このスクリプトはテスト済みであり、QGISで実行すると正常に機能します。 ##Test=name import os import glob import sys sys.path.append("C:\Program Files\QGIS Brighton\lib;%OSGEO4W_ROOT:\=/%/apps/qgis;%OSGEO4W_ROOT%\apps\qgis\bin;%OSGEO4W_ROOT%\apps\grass\grass-6.4.3\lib;%PATH%") from PyQt4.QtCore import * from PyQt4.QtGui import * from qgis.core import * QgsApplication.setPrefixPath("C:\Program Files\QGIS Brighton\apps\qgis", True) QgsApplication.initQgis() from os.path import expanduser home = expanduser("~") # Folder path of the …

1
URL経由でフィルターパラメーターをGeoServer WFSに渡しますか?
PHPのCURL関数を使用して、geoServerインスタンスからgeoJSON情報を取得しようとしています。完全なxml GetFeatureリクエストを作成しようとするのではなく、url変数を使用してこれを行っています。 ただし、一部のプロパティフィールドのコンテンツに基づいて結果のサブセットを取得できるようにしたいと思います。 したがって、このURLを使用してすべての結果を取得できます。 http://www.myURL.com/geoserver/namespace/ows?service=WFS&version=1.0.0&request=GetFeature&typeName=layername&outputFormat=json&BBOX=1,2,3,4 結果を「フィールドLIKE値」または「フィールド=値」に制限できますか? Geoserver / ECQLのドキュメントを1時間かけて探しましたが、どれも「この形式でURLに関数を追加する」と明記されていません。その結果、url変数を使用してこれらのタスクを実行できるかどうかは明確ではありませんが、一部の人々はそれを提案しているようです。 GeoServerのWFSとURLパラメーターを使用した、機能するフィルターの明確な例について、誰でも手助けできますか?
17 geoserver  wfs  geojson  php  curl 

1
長時間の作業でも目に優しい色
この質問がここのルールに合わない場合はおApび申し上げます。 私は大量のリモートの歩道をデジタル化し、パン、ズームをたくさん行い、これまで2日間行ってきましたが、常にひどい頭痛がします。私が使用している色ですか? 機能を追加する既存の歩道/道路レイヤーがあります。新しいフィーチャを作成する歩道のポリゴン、この歩道の周囲のバッファ、完了した領域をカバーするために使用するプレーンポリゴンレイヤー。 誰かが目の疲れや頭痛を最小限に抑えるために使用するいくつかの代替色を推奨できますか?
17 color 

1
西に面するすべてのビーチのOSMを照会しますか?-QGISが望ましい
この惑星(またはその特定のサブセット)で西向き(+/- 20-30度)にあるすべてのビーチを特定したい。 私の考えは、OSMをデータソースとして使用することです。タグ「natural = beach」を持つすべてのポリゴンを抽出すると、ビーチとしてフラグが付けられているすべてのエリアのデータセット/シェープファイルが得られます。ここまでは順調ですね。 また、海岸線のデータ(OSMやNatural Earthなど)を使用して、ビーチの向きを特定できることも知っています。 小さな不一致を修正するためにビーチの形状をバッファリングし、これらを海岸線と交差させると、形状のどちら側が水に面しているかに関するデータが得られます。交差するすべてのポイントの最大値/最小値を使用すると、方位がわかります。(90度+/- 20-30度なら問題ありません...) そして、それは私が失われた場所です... 私の好みのツールはQGISですが、PostGISも同様です。 このパズルを解く方法を知っている人はいますか?

12
経度と緯度の意味をどのように記憶しますか?[閉まっている]
閉じた。この質問は意見に基づいています。現在、回答を受け付けていません。 この質問を改善したいですか?この投稿を編集して事実と引用で答えられるように質問を更新してください。 5年前に閉鎖されました。 経度と緯度は、私が通常覚えておく必要のない概念です。しかし、私は時々それらを使用しているため、しばらくの間突然意味を学ぶ必要があります。 私の日々の終わりまでその意味を最終的に記憶するための補佐官のメモ(記憶-援助)はありますか? これら2つの用語の語源は何ですか?

2
GDALを使用したPythonでのWGS 1984 Webメルカトル(EPSG:3857)の再投影
GDALを使用してPythonでラスターを再投影しています。後でOpenStreetMapおよび多分Googleマップと一緒にOpenlayersで使用するために、地理的なWGS 84座標からWGS 1984 Webメルカトル(補助球)にいくつかのTIFFを投影する必要があります。私はここからPython 2.7.5とGDAL 1.10.1を使用し、ここからのアドバイスを使用して座標を変換しています(私のコードは下にあります)。つまり、osgeo.osrをインポートし、ImportFromEPSG(code)とCoordinateTransformation(from、to)を使用しました。 私は最初にUTMゾーン29であるEPSG(32629)を試し、この投影されたラスターを取得しました(多かれ少なかれ細かい)ので、コードは正しいようです: それから私はこれとこの質問を読んで見つけたのでEPSG(3857)を使用しましたそれが正しい最近の有効なコードであること。ただし、ラスタは空間参照なしで作成されます。WGS 84データフレームでは遠く離れています(ただし、データフレームをWebメルカトルに切り替えても問題ありません)。 EPSG(900913)出力が地理参照されているが、北に3つのラスタセルについてシフト: ArcGIS(WGS_1984_Web_Mercator_Auxiliary_Sphereにエクスポート)を使用してラスターを再投影すると、結果はほぼ問題ありません。 そして、古いコード102113(41001,54004)を使用すると、結果は完璧です。 すべてのコードを使用した私のテストの要約: 3857: far away up (missing georeference) 3785: far away up (like 3857) 3587: far away right 900913: slightly jumped up 102100: python error 102113: perfect 41001: perfect 54004: perfect ArcGIS (web merc. aux.): good だから私の質問は: 正しいEPSGコードが間違った結果を与えるのはなぜですか? そして、なぜ古いコードがうまく機能するのでしょうか? …

1
shapelyを使用して、ポイントに最も近いラインセグメントを検索しますか?
バックグラウンド 既知のポイントから、図に示されているように、MultiLineStringのテーブルに対して最も近い周囲の「見える境界」を確立する必要があります。 私はこのサイトをいくつかの用語(例:最小エッジ、最小周囲、最近傍、クリップ、ポリゴンを含む、可視性、スナップ、カットノード、レイトレース、フラッドフィル、内部境界、ルーティング、凹型ハル)で検索しました。このシナリオに一致すると思われる以前の質問を見つけることができません。 図 緑の円は既知のポイントです。 黒い線は、既知のMultiLineStringです。 灰色の線は、既知のポイントからの放射状の掃引を示しています。 赤い点は、放射状スイープとMultiLineStringsの最も近い交点です。 パラメーター PointはMultiLineStringと交差しません。 Pointは常にMultiLineStrings内の名目上中央に配置されます。 MultiLineStringはPointを完全に囲むことはないため、境界線はMultiLineStringになります。 約1,000個のMultiLineStringを含むテーブルがあります(通常は約100ポイントの単一行を含みます)。 考慮される方法論 既知のポイントから一連の線を構築することにより、放射状スイープを実行します(たとえば、1度の増分で)。 MultiLineStringを使用して、各放射状スイープラインの最も近い交点を確立します。 放射状スイープラインの1つがMultiLineStringのいずれとも交差しない場合、これは、MultiLineStringの周囲構造に収まる周囲のギャップを示します。 概要 この手法は最も近い交差点を検出しますが、放射状スイープの解像度に応じて、必ずしもすべての最も近い境界ノードポイントを検出するわけではありません。誰もがすべての周辺ポイントを確立するための代替方法を推奨したり、何らかの形のバッファリング、セクター化またはオフセットでラジアルスイープテクニックを補完できますか? ソフトウェア 私は、SpatiaLiteやShapelyをソリューションに使用することを好みますが、オープンソースソフトウェアを使用して実装できる提案を歓迎します。 編集:実用的なソリューション(@geneの回答に基づく) from shapely.geometry import Point, LineString, mapping, shape from shapely.ops import cascaded_union from shapely import affinity import fiona sweep_res = 10 # sweep resolution (degrees) focal_pt = Point(0, 0) …

1
GDAL Pythonでデータセットを閉じるのはなぜですか?
Python GDALコードには、スクリプトの最後でデータセットを閉じる人がいることがよくあります。Python GDALでデータセットを閉じるのが理にかなっているのはなぜですか?実行しないと、結果はありますか? import gdal # open dataset ds = gdal.Open('test.tif') # close dataset ds = None
17 python  gdal 

8
ArcGIS for Desktopでフィールド名を変更しますか?
ArcGIS 10.2でフィールド名(エイリアスではなく)を変更するにはどうすればよいですか? ArcGISサポートでこのページを既に見つけました: " シェープファイルフィールドの名前変更(ディフェンスマッピング) "が、ArcCatalogのカスタマイズモードカテゴリに" ディフェンスマッピング " がありません! 防衛マッピングにアクセスする方法、またはフィールド名を変更する方法を教えてください。

1
ogr2ogrを使用して属性でフィーチャを選択しますか?
何らかの処理の後にシェープファイルを作成しています。ただし、最後の手順では、シェープファイルの属性による選択を行う必要があります。 コマンドが動作するときにPythonスクリプトで使用することを目的として、シェルでコマンドを使用しています。 ogr2ogr -f "ESRI Shapefile" -select * where ID="1" outfile.shp infile.shp エラーメッセージが表示されます。 FAILURE: Unable to open datasource `Downloads' with the following drivers. 何が間違っているのでしょうか?

6
QGISを使用して、選択したポリゴンの境界ボックス座標を計算します
郡ポリゴンの境界ボックス座標を手動で(そして多少不正確に確信しています)取得しました。郡ポリゴンは、郡形状のより大きなレイヤーからの選択です。QGISプラグインまたは他の関数を使用して、これを手作業ではなく迅速に計算できるようにしたいと思います。これはほぼ確実に可能ですが、Google検索で何も見つかりませんでした
17 qgis  extents 

2
Esriの.ascファイルを理解していますか?
ArcGISヘルプを読む:サポートされているラスターデータセットファイル形式で、16ビットの符号付き整数(離散データ用)と32ビットの浮動小数点ファイル(連続データ用)の両方について、ASCIIグリッドファイルタイプが単一のファイル拡張子.ascでアドレス指定されていることを読みました)。 Esri Grid形式のセクションで次のように述べています。 グリッドは、Esri固有のラスターデータストレージ形式です。グリッドには、整数と浮動小数点の2つのタイプがあります。 次に、両方の形式について説明します。 私の混乱は、次のヘルプセクションEsri ASCIIラスターフォーマットで発生します。そこでは、浮動小数点形式のみに言及しています。だから、私は知りたい: 離散データの16ビット符号付き整数形式も、グリッド以外にラスター形式と見なされますか? 「ESRIグリッド形式-グリッドデータストレージ」セクションでは、.bnd、.hrdなどの他の拡張機能をいくつか挙げています。16ビットの符号付き整数の.ascファイルは、Esriの離散データグリッド形式のASCIIバージョンですか? この質問の動機は、拡張子が.ascのファイルを説明/理解することです。

2
リーフレット:コンテナ座標とレイヤー座標?
マップコンテナ座標とレイヤー座標の違いは何ですか? 私が集められたところでは、地図のコンテナの座標は常にマップビューポートに対して相対的である、としながら、層の座標が常に同じ値でスタートコンテナcoordsのマップがパンされると、彼らが変更されます。 座標/ポイントを参照するAPI: LatLngからレイヤーポイント コンテナポイントへのLatLng

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