タグ付けされた質問 「python」

Pythonは、多くのGISプログラムで使用されているオープンソースのインタプリタ型プログラミング言語です。


13
PythonでGDALをWindowsにインストールしますか?
WindowsでPythonを使用してGDAL / OGRをインストールする方法を説明できる人はいますか? Windows Vistaを使用していますが、Webサイトの情報を追跡しようとしましたが、exeファイルがバインドされていないようです。 誰かが必要なファイル/フォルダーへのリンクを含めて、プロセスを説明できますか? 私は今、minGWでgdalセットアップを実行しようとしましたが、これも失敗しました:

11
Arcpyを使用する代わり
私はほとんどすべてのpythonジオプロセシングにESRIのArcpyサイトパッケージを使用しているようです。ESRIの功績として、これらは非常に優れたツールスイートであり、多大な成果を上げることができます。ただし、ESRI Arcpyドメイン外のジオプロセシングスクリプトも作成したいと思います。たとえば、ラスターをポリゴンにクリップする場合、ESRIの次のスクリプトから始めます。 # Import system modules import arcpy from arcpy import env from arcpy.sa import * # Set environment settings env.workspace = "C:/sapyexamples/data" # Set local variables inRaster = "elevation" inMaskData = "mask.shp" # Check out the ArcGIS Spatial Analyst extension license arcpy.CheckOutExtension("Spatial") # Execute ExtractByMask outExtractByMask = ExtractByMask(inRaster, inMaskData) ...

23
Python GDALパッケージには、pip経由でインストールするときにヘッダーファイルがありません
私は尋ねたスタックオーバーフローでこの質問を、私は、このサイトが存在していたが実現する前に。このサイトがより適切であるため、ここでも質問します。 pip install gdal仮想環境(Ubuntu)内のpipからgdalをインストールしようとしています。見つからないため失敗するcpl_port.h extensions/gdal_wrap.cpp:2853:22: fatal error: cpl_port.h: No such file or directory compilation terminated ただし、GDALは正しくインストールされ、ヘッダーファイルはにあります/usr/include/gdal/cpl_port.h。pipがヘッダーファイルを見つけるために設定する必要があるGDALの環境変数はありますか?

3
等高線図のポリゴンを平滑化しますか?
レベルのすべてのポリゴンが使用可能な等高線図を次に示します。 すべての頂点を正確な位置に保持したままポリゴンを滑らかにする方法を尋ねましょう。 実際、輪郭はグリッドデータの上に作成されるため、グリッドデータを平滑化することをお勧めします。したがって、結果の輪郭はより滑らかになります。ガウスフィルターなどの平滑化機能は小さなデータパックを削除し、3番目の変数の範囲、たとえば、アプリケーションで許可されていない高さを変更するため、これは私の希望どおりに機能しないことに注意してください。 実際、私は2Dポリゴン(凸、凹、自己交差などのあらゆるタイプ)の平滑化を合理的に無痛(コードのページを忘れる)で正確に行うことができるコード(好ましくはPythonで)を探しています。 参考までに、ArcGISにはこれを完全に行う機能がありますが、サードパーティの商用アプリケーションを使用することは、この質問に対する私の選択ではありません。 1) Scipy.interpolate: ご覧のとおり、結果のスプライン(赤)は満足のいくものではありません! 2) ここで与えられたコードを使用した結果はここにあります。うまく機能していません! 3) 私にとって最善の解決策は、次の図のようなもので、1つの値のみを変更することで、正方形が徐々に平滑化されることです。あらゆる形のポリゴンを平滑化するための同様のコンセプトを期待しています。 スプラインがポイントを通過する条件を満たします。 4) Pythonでの彼のデータに関する「whuberのアイデア」の行ごとの実装です。結果が良くないので、おそらくいくつかのバグがあります。 K = 2は災害なので、k> = 4の場合です。 5) 問題のある場所の1つのポイントを削除すると、結果のスプラインはwhuberと同じになりました。しかし、なぜこの方法がすべての場合に機能しないのかは疑問です。 6) whuberのデータの適切なスムージングは​​、次のようにすることができます(ベクトルグラフィックソフトウェアで描画) 4): 7) いくつかの象徴的な形については、Pythonバージョンのwhuberのコードの結果を参照してください。 なお、この方法は、ポリラインのために動作しないようだということ。コーナーポリライン(輪郭)については、緑色が必要ですが、赤色になります。等高線マップは常にポリラインであるため、これに対処する必要がありますが、閉じたポリラインは私の例のようにポリゴンとして扱うことができます。また、アップデート4で発生した問題がまだ対処されていないというわけでもありません。 8)[私の最後] これが最終的な解決策です(完璧ではありません!): 星が指す領域について何かしなければならないことを忘れないでください。おそらく私のコードにバグがあるか、提案された方法はすべての状況を考慮し、望ましい出力を提供するためにさらなる開発が必要です。

8
Pythonを使用して投影座標を緯度/経度に変換しますか?
これは基本的な質問だと思いますが、解決策を見つけたり認識したりできないようです。 このサイトは戻ります Point: X: -11705274.6374 Y: 4826473.6922 例として最初のキー値000090で検索する場合。これは空間参照であ​​り、それが何であるかがわかります。 Pythonを使用してこれを緯度と経度に変換する方法の手順または例を探しています。

17
一般的なGISの目標を念頭に置いてPythonプログラミングを学習しますか?
あなたの意見では、GISの目標を念頭に置いてPythonを学ぶのに最適な本/サイトは何ですか? 「最高」とは、次のことを意味します。 それほど長くない(本) わかりやすい(本/サイト) 実用的な例(書籍/サイト) Python用ArcPyサイトパッケージを使用してArcGIS for Desktopをカスタマイズする方法を学習するための具体的な回答については、別のQ&A:ArcPyを学習するためのリソースを教えてください。

1
NumPyなしでOGRポイントの下でGDALラスターのピクセル値を取得しますか?
私は、風景全体の野生の花粉媒介者の豊富さの計算モデルに取り組んでいます。モデル自体は完成しており、私は現在、後処理ステップに苦労しています。 次のようなGDAL受粉器供給ラスターがあります(明るい色は、受粉器への訪問回数が多いことを意味します)。 そして、ランドスケープ上のサンプル位置を表すポイントのOGRシェープファイルがあります。 私はこれらのポイントの下のピクセルでいくつかの分析を実行しようとしていますが、そうするには、ポイントの下のピクセルの値を抽出できる必要があります。 PythonでOGRとGDALのみを使用して、ポイントの下のピクセルの値を抽出することは可能ですか?ReadAsArray()出力ラスターは非常に大きい(メモリに収まるには大きすぎる)ため、ラスター全体をを介してメモリに読み込むことは避けたいと思います。 私はこの投稿に気づきましたが、これは似ていますが、コマンドライン呼び出しが必要です。
45 python  gdal  ogr  pixel 

3
ArcMapライセンスなしでPythonを使用してShapefileのコンテンツを表示することはできますか?
ArcMapライセンスがなくてもPythonを使用してシェープファイルの内容を見ることができるかどうか疑問に思いました。状況は、ESRIソフトウェアだけでなく、多くの異なるアプリケーションからシェープファイルを作成できることです。空間参照、フィーチャタイプ、属性名と定義、およびシェープファイルのフィールドの内容をチェックし、それらを許容値のセットと比較するPythonスクリプトを作成したいと思います。組織にESRIライセンスがない場合でも、このスクリプトが機能するようにします。このようなことを行うには、ArcPyを使用する必要がありますか、それともArcPyを使用せずにシェープファイルを掘り下げることができますか?
40 python  shapefile 

5
Pythonコードからgdal_translateを呼び出す方法は?
いくつかのgdal APIを使用gdal_translateしてPythonコードから呼び出すことは可能ですか?ファイルシステムからgdal_translate.exeを実行するだけではなく、コード内でなんらかの方法で呼び出すので、gdal_translate実行可能ファイルが存在する正確なディレクトリを知る必要はありません。
40 python  gdal 

4
Pythonを使用したシェープファイルの簡単なテーマ別マッピング?
ArcGIS / ArcPyを使用せずにPythonで地理データを視覚化し、地図を作成したいと思います。 インターネットで、Pythonを使用して主題図を作成する方法を見つけました。 私が試したコードは次のとおりです。 import shpUtils import matplotlib.pyplot as plt shpRecords = shpUtils.loadShapefile('C:\\Users\\shapefile.shp') for i in range(0,len(shpRecords)): x = [] y = [] for j in range(0,len(shpRecords[i]['shp_data']['parts'][0]['points'])): tempx = float(shpRecords[i]['shp_data']['parts'][0]['points'][j]['x']) tempy = float(shpRecords[i]['shp_data']['parts'][0]['points'][j]['y']) x.append(tempx) y.append(tempy) plt.fill(x,y) plt.axis('equal') plt.title("Testing") plt.show() ただし、これを実行すると、ランダムな色が表示されます。 シェープファイルの特定の列を視覚化する場合、同様のコードを使用してこれを実装するにはどうすればよいですか? これは、彼が色の使用法についてのみ議論している上記のリンクでは非常に不明確です... デカルトやPySALなど、これを実現するためにおそらく追加のモジュールが必要でしょうか?

7
ポイントレイヤーとラインレイヤーの最近傍 [閉まっている]
#qgisと#postgisの間のstackoverflowとircでこの質問を数回しました。また、実際の答えなしにそれをpostgisにコーディングまたは実装しようとしました。 プログラミング(できればpython)を使用して、ポイントレイヤーから、ラインレイヤーまたはポリゴンレイヤーの最も近いラインに投影するまで、ラインを描画します。 現在のところ、私のデータのほとんどはESRIの形状とpostgis形式です。ただし、私は主にshp + qgisユーザーなので、postgisソリューションには近づかない方がいいでしょう。 理想的な解決策は、Pythonまたは同様のライブラリでGDAL / OGRを実装することです GDAL / OGRライブラリを使用して、どこから始めるべきですか?ソリューションプランを提供することは可能でしょうか? NetworkXを使用して最近傍分析を実行できますか? これは実際に可能ですか? 簡単な場合、ポイントは投影ポイントではなくセグメントのエンドポイントに接続できます

4
arcpyを使用せずにシェープファイルをジオプロセシングするためのPythonスクリプトの例
arcpyに基づいていないPythonスクリプトを使用して、属性によるシェープファイルのクエリ、選択からの新しいレイヤーの作成、ポリゴンの面積の計算、ポリゴンのポイントへの変換などを行います。 他のPythonモジュールまたはライブラリを使用するコード例はありますか?私はarcpyを使用してこれを簡単に行うことができますが、他のオプションを調べたいと思いました。

3
ArcGIS“ in_memory”ワークスペースの適切な構文と使用法は何ですか?
in_memoryArcGISでワークスペースを最適に使用する方法がわかりません。このテーマで見つけた最良のリソースは、ESRIの「インメモリワークスペースの使用」です。残念ながら、ヘルプセクションはあまり有益ではありません。また、このサイトには、こことここのテーマについて簡単に触れた投稿があります。 私の特定の質問: in_memoryArcGIS / arcpyベースのスクリプトでワークスペースを使用するための適切な使用法と構文は何ですか? あるin_memoryワークスペースは、例えば、使用してレイヤーを作成し、同じarcpy.MakeFeatureLayer_management()? in_memoryスクリプトの最後にワークスペースを削除するなどの標準はありますか?

10
ベクターデータセットのすべてのフィーチャを簡単にシフトする方法は?
シェープファイルを作成し、すべてのフィーチャの頂点が一定量シフトしているとしましょう。任意のシフトですべてのフィーチャをシフトする最も簡単な方法は何ですか(したがって、頂点の(x、y)位置)?この修正を適用するファイルがたくさんあるので、Bash / OGRの回答が望ましいでしょう:) 最後に、いい機能があるので、私はこのためにSpatialiteを使用することになりましたShiftCoords。しかし、スレッドは非常に有益でした!皆さんありがとう!

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