地理情報システム

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

1
PyQGISファイルロックを解放しますか?
pyQGISでファイルロックのリリースをトリガーするものは何ですか? を呼び出して(一時的に使用される)いくつかのデータソースを削除しようとしていますQgsVectorFileWriter.deleteShapeFileが、それを行う前にQGISを終了する必要があります。ソースをQgsVectorLayerオブジェクトにロードしました。ソースを削除する前に、これらのオブジェクトとそれらへの参照をすべてガベージコレクションする必要がありますか?これを強制する方法はありますか? 失敗する最小限のコードサンプルを作成できました。実行する前に、一時ディレクトリが空であることを確認してください。 from qgis.core import * import processing, os, gc project_temp_dir = "C:/Path/To/My/Dir/" layer1_path = project_temp_dir + "layer1.shp" layer2_path = project_temp_dir + "layer2.shp" input_layer = QgsMapLayerRegistry.instance().mapLayersByName('in_layer')[0] if not input_layer.isValid(): raise Exception("Failed to grab input layer") # Create layer 1 err = QgsVectorFileWriter.writeAsVectorFormat(input_layer, layer1_path, "utf-8", input_layer.crs()) if err != QgsVectorFileWriter.NoError: …
16 pyqgis  lock 

4
QGISでOD道路距離マトリックスを作成するには?
出発地と目的地(OD)のコストマトリックスを作成するために必要な手順を探しています。 私が探しているのは、道路距離のマトリックスです。距離行列ツールは、ユークリッド距離のみを生成します。QGIS Roadgraphプラグインはから最短道路距離生成1人の Nの宛先に起源を。必要なのは、すべての出発地からN個の最も近い目的地までの距離を含むファイルです。私には多くの起源があるので、ロードグラフで分析を繰り返すことは選択肢ではありません。再度、感謝します。 私は、ゼロから始めて段階的な説明を探しています(つまり、3つのシェープファイル:出発地、目的地、道路網)。

5
バッファなどのジオプロセシング用のPythonライブラリ(ArcPy以外)を探していますか?[閉まっている]
閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか?この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 8か月前に閉鎖されました。 ArcPyを除き、シェープファイルを使用してバッファー/交差などのジオプロセシングを実行できるPythonライブラリはありますか?

4
Pythonを使用して属性に基づいて機能をフィルターしますか?
Pythonを使用してQgisで属性(arcobjectsのIqueryfilterに似ています)で機能を取得する方法は?すべての機能を取得して手動でフィルタリングする代わりに、where句を使用してフィルタリングするオプションはありますか? 例:「Counties」というフィールド名があります。5万を超える機能があります。つまり、時間がかかるため、すべての機能を取得してフィルタリングすることはできません。だから私はアークオブジェクトでiqueryfilter.whereclause = 'Counties = Norwich'を使用してそれをクエリすることができます。PyQgisでも同様のことが必要です。
16 qgis  pyqgis 

2
Postgisでポリゴンのポイントを数える
単純な問題があります。ポリゴンのセット内のポイントの数をカウントしたいのです。 私はすでにSQLを持っていますが、実際にポイントを含むポリゴンのgidのみを返します。 私のテーブル:19.000行のポリゴンレイヤーと450行のポイントレイヤー。 次のSQL select grid.gid, count(*) AS totale FROM grid, kioskdhd3 WHERE st_contains(grid.geom,kioskdhd3.geom) GROUP BY grid.gid; 実際にポイントを含む320個のポリゴンのみを返します。しかし、ポイントの数が0であっても、すべてのポリゴンが返されるようにします。 もちろん、それは私のWHERE節に関係しています。st_containsのどこに置く必要がありますか? ありがとう、カーステン

1
NumPy配列を使用してビッグデータジオプロセシングを最適化するにはどうすればよいですか?
NumPy配列を利用してジオプロセシングを最適化する方法を学ぶことに興味があります。私の仕事の多くは「ビッグデータ」に関係しており、ジオプロセシングでは特定のタスクを完了するのに数日かかることがよくあります。言うまでもなく、これらのルーチンを最適化することに非常に興味があります。ArcGIS 10.1には、arcpyを介してアクセスできる次のような多くのNumPy関数があります。 NumPyArrayToFeatureClass(arcpy.da) RasterToNumPyArray(arcpy) TableToNumPyArray(arcpy.da) 例として、NumPy配列を利用した次の処理集中型ワークフローを最適化するとします。 ここでの一般的な考え方は、ベクトルとラスタベースの操作の両方を移動する膨大な数のベクトルベースのポイントがあり、その結果バイナリ整数ラスタデータセットが得られるということです。 このタイプのワークフローを最適化するためにNumPyアレイをどのように組み込むことができますか?

7
2つのマップを素敵なWebインターフェイスと比較しますか?
LeafletまたはOpenLayersでは、2つのマップをバーで比較して、2つのマップをスライドさせることはできますか? 例は次のとおりです。 http://www.esri.com/services/disaster-response/hurricanes/hurricane-sandy-the-aftermap http://www.arcgis.com/apps/StorytellingSwipe/index.html?appid=c62112f6ba834227978d8f1cf2b6f5af http://www.abc.net.au/news/specials/qld-floods/(地図ではありません) 私はそれが可能であることを知っています、しかし、私はそれをする方法を知りません。誰でもアイデアはありますか?例?

4
カスタムリーフレットコントロール
リーフレットのカスタムコントロールを作成するのにいくつかの助けがありますが、リーフレットの最新バージョンで使用されているCSSスタイルへのアップグレードにまだ問題があります。基本的に、必要なのは、ドロップダウンシャドウが周囲にある白いボックス内のアイコンだけです。 ここに私がこれまで取り組んできたものがあります。見て、気軽にフォークしてください: http://codepen.io/DrYSG/pen/zJsiG
16 leaflet 


1
「ノードなし交差点」とは何ですか?
「ノードなし交差点」とは何かを正確に定義できますか?このテキストに誤りがあるさまざまなGISプラットフォームとツールキットがあります。このためのGoogle検索結果は、何が問題なのかを実際に説明することなく「精度を落とすとすべての問題がなくなる」というカテゴリに分類される傾向があります。2つの有効なジオメトリがある場合、なぜこのエラーが発生するのでしょうか? ドキュメントへの参照が最も役立ちます。

4
RESTful APIを介したArcGIS Serverのセキュリティ保護されたレイヤーへのプログラム認証
インターネット上のセキュリティで保護されたマップサービスを公開するArcGIS 10.1サーバーインスタンスがあります。私が必要としているのは、ユーザーがこれらのセキュリティで保護されたWebサービスを表示できるようにするクライアントアプリケーション(現在、ArcGIS Javascript APIの3.3バージョンを使用して構築中)をコーディングすることです。 このオンラインESRIの例は良いスタートだと思います。 私はすべてのマップサービスが彼女/彼に属し、同じユーザー名とパスワードで彼女/彼によってアクセス可能であることをすでに知っているので、ユーザーに各マップサービス認証を要求することはありません。私の考えでは、資格情報のプロンプトは1回だけ表示されるはずです。したがって、JSコードは、何らかのRESTfulログイン呼び出しを介して資格情報を各mapserviceにフィードする必要があります。私には、ArcGIS Server REST APIがそのような呼び出しを提供しているようには見えません...多分間違っています。 したがって、セキュリティで保護されたマップサービスにログインするこの「RESTful」な方法は、ArcGIS Serverで実行可能ですか?ある場合、これを説明するWebリソースへの例またはリンクを提供できますか?


6
ArcGIS Desktopで使用するCSVファイルを準備していますか?
ArcGIS Desktopで使用するCSVファイルを準備するにはどうすればよいですか。 ArcGISは間違ったフィールドタイプを列に割り当て、áやêなどの特殊文字を誤って解釈するため、CSVファイルの使用に問題があるため、お願いします。 Esriフォーラムで、「Col22 = V002 Text」などのフィールドタイプを何らかの形で定義するいわゆるschema.iniファイルがあることを読みました。http://forums.esri.com/Thread.asp?c = 93&fを参照してください= 1149&t = 64464 ディスク上でこれらの.iniファイルを頻繁に見たことがありますが、実際には何が良いのか疑問に思ったことがないので、それはちょっとおかしいです。Rのような他のプログラムはそうしないので、Excelがそのようなメタデータを追加のファイルに保存するのはちょっと変です。 「文字列」型などの適用方法がわからなかったため、すでにこの.iniファイルを操作しようとしてほとんど成功しませんでした。MSサイトにはいくつかの情報があります。http://msdn.microsoft.com/en-us/library/windows/desktop/ms709353%28v=vs.85%29.aspxを参照してください。しかし、解決策が見つかりませんでした。 。 また、この.iniファイルで作業するというアイデアはあまり好きではありませんでした。50列とすると、すべてのフィールド名を定義して入力するのはかなりの作業です。そして、これらの.iniファイルは失われる可能性があります。

3
ラスターデータベースのクエリを高速化する方法は?
私はこれらの列を持つpostgresql / postgisにラスターデータベースを持っています: (ID、rast、data_of_data)。 「ラスト」は、WKT形式のラスターファイルがある列です。WGS84システム(30.424、-1.66)および2002-01-09のポイントのDN値を検索するクエリの例は次のとおりです。 SELECT st_value(rast,(st_GeomFromText('POINT(30.424 -1.66)', 4326))) as val FROM my_table WHERE date_of_data='2002-01-09' これらの種類のクエリを高速化する方法(空間インデックスなど)はありますか?

3
PythonでGeoJSONを生成する
シェープファイルのポリゴンを使用してGeoJSONファイルをプログラムで作成しますが、独自のアプリケーションから属性を追加します。 これは、シェープファイルに対して簡単に実行できます。 def create_data_dayer(self,varlist, data): """ Creates a new shape to contain data about nodes. varlist is the list of fields names associated with the nodes. data is a list of lists whose first element is the geocode and the remaining elements are values of the fields, in the same …
16 python  geojson  ogr  fiona 

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