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

PostGISサーバーからの地理データを処理するプログラムの作成を開始し、GDALを使用します。 あなたの経験によると、それを使ってプログラムするのに最適な言語はありますか？PerlとPythonでプログラミングできますが、再利用するコード、サンプル、ライブラリなどがより多く含まれている言語を知りたいです。

20 python  postgis  gdal 

Pythonで開発し、OSGEOのGDALを使用して、ラスターとシェープファイルを操作および操作しています。 ポイントフィーチャを含むシェープファイルを取得して、サーフェスラスタに補間したいです。現在、私は「RasterizeLayer」メソッドを使用しています。このメソッドは、ポイントフィーチャからラスター（すべてのnodata値で設定されている）に値を焼き付けますが、すべてのタッチされていないピクセルを「nodata」値として残します。したがって、チェッカーボードタイプのラスターが残っています。 RasterizeLayerを使用した後の内容： 最終製品に必要なもの： 私が探している関数は、arcgisscripting importの 'Spline_sa（）'として知られていると思います。 GDALには同様の機能がありますか、または希望する出力を得るための別の方法がありますか？

20 python  raster  gdal  interpolation 

Pythonを使用して、大きな円距離を作成する必要があります。これは、数値であり、できればクライアント側のマップを描画するために使用できる何らかの「曲線」です。曲線の形式（WKTであろうと座標のペアであろうと）については気にしませんが、データを取得したいだけです。 そこにはどんなツールがありますか？何を使うべきですか？

20 python  distance  polyline-creation  great-circle 

2つのジオデータフレームがあります。 import geopandas as gpd from shapely.geometry import Point gpd1 = gpd.GeoDataFrame([['John',1,Point(1,1)],['Smith',1,Point(2,2)],['Soap',1,Point(0,2)]],columns=['Name','ID','geometry']) gpd2 = gpd.GeoDataFrame([['Work',Point(0,1.1)],['Shops',Point(2.5,2)],['Home',Point(1,1.1)]],columns=['Place','geometry']) そして、gpd1の各行について、gpd2の最も近いポイントの名前を見つけたいです。 desired_output = Name ID geometry Nearest 0 John 1 POINT (1 1) Home 1 Smith 1 POINT (2 2) Shops 2 Soap 1 POINT (0 2) Work 私はラムダ関数を使用してこれを機能させようとしています： gpd1['Nearest'] = gpd1.apply(lambda row: min_dist(row.geometry,gpd2)['Place'] , …

20 python  distance  geopandas  pandas 

閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか？この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 閉じた3年前。 Webマッピングアプリケーションの開発に興味があります。私はPythonに比較的精通しているので、Pythonフレームワークを探しています。私のオプションは何ですか？各オプションの長所と短所は何ですか？ 地図上の地理情報を表示し、ユーザーが「情報」を使用して詳細を取得できる単純なオプションを探しています。この時点で、クエリや動的機能は必要ありません。私はよく文書化されたオプションを好むでしょう。

19 python  web-mapping  software-recommendations 

Pythonで2つのポリゴンが交差するかどうかを判断するのに役立つアルゴリズム、高レベルのソリューション、またはライブラリを探しています。 2つの異なる配列に2つのポリゴン（これらは穴のない単一パーツポリゴンです）の頂点があります。ポリゴンは2Dです（XおよびY座標のみ） これらの2つのポリゴンが交差するかどうかを示すブール値を返す関数を作成します。 私はarcpy、またはarcgisこのコンポーネントを使用できないことに注意してください。 これを行うためのアルゴリズムまたはライブラリを提案できますか？

19 python  algorithm  intersection 

WindowsのQGISプラグインでサードパーティライブラリを使用するにはどうすればよいですか？私は顧客のためにrasterioとnumpyを使用するプラグインを開発しましたが、rasterioとnumpyのインストールに問題があります。 実際には、rasterioとnumpyはメインシステムPython（C：\ Python27）にインストールされていますが、それを認識するためにQGIS Pythonが必要です。 どうすればできますか？

18 python  pyqgis  installation  windows 

私はqgis.coreと他のさまざまなqgisライブラリを使用するWindows 7でスタンドアロンのPythonスクリプトを作成しようとしています。 このテーマに関する古い質問をここで見つけました：PyQGISを使用してスタンドアロンPythonスクリプトを作成しますか？こちらのガイドを使用して、アドバイスに従っていました：http : //docs.qgis.org/testing/en/docs/pyqgis_developer_cookbook/intro.html#python-applications そこで提案されているコードは次のとおりです。 set PYTHONPATH=c:\qgispath\python に続く： set PATH=C:\qgispath;%PATH% IDLEは「最初の設定」を好まないことがわかったため、コードは次のようになります。 PYTHONPATH="C:\OSGeo4W64\apps\qgis\python" PATH="C:\OSGeo4W64\apps\qgis\python\qgis\core;%PATH%" import qgis.core そして私は得る： Traceback (most recent call last): File "<pyshell#15>", line 1, in <module> import qgis.core ImportError: No module named qgis.core 最後に;％PATH％ "を削除することを考えると考えられるパスのすべてのバリエーション（少し余分に見える）を試しましたが、それ以上の情報を見つけることができませんでした。間違っている？ 編集： ルイージ・ピレリのコメントを見た後、私は今試しました： >>> PYTHONPATH="C:/OSGEO4~1/apps/qgis/./python\qgis" >>> import qgis.core Traceback (most recent call last): File …

18 qgis  python  pyqgis  windows 

指定した基準の近傍内で、ラスターの各セルの焦点統計を計算しようとしています。 背景-3つのバイナリラスタがあり、それぞれが関心のある単一の植生タイプを表します。調査地域内のセル（近隣のセルの合計/合計）から20 km ^ 2以内の各植生タイプの被覆率を計算したいと思います。問題は、各セルの周囲に単純な円または正方形の近傍を使用できないことです。これを行うと、合計の計算に使用される検索領域に研究領域外の領域が組み込まれるためです。この例外は重要です。統計は生息地モデルの入力として使用され、研究エリア外の地域は生息地の可能性があると見なすことができないため、都市化されているためです。それらを含めると、誤った統計が得られます。だから、私はnは、基準を満たす、希望する近隣のサイズに等しいエリアをカバーするために必要なセルの数によって決定されます。基準は、都市化された地域に該当しないことです。何らかの形のセルオートマトンを使用すべきだと考えています。私はCAで働いたことはありませんが。 私が望むのは、スターターコード、または正しい方向のポイントのようなものだと思います。 以下のコメントに返信： 私の研究サイトの境界にあるセルのこの統計を計算しているとしましょう。調査エリア外のすべてのエリアをゼロに割り当てる（またはNoDataを無視する）と、関心のあるエリアカバレッジの約半分を表す統計が得られます。したがって、約10 km ^ 2のエリアのカバレッジ率、20 km ^ 2エリアの代わりに。私は家の範囲のサイズを勉強しているので、これは重要です。それは動物が風景を見る/使用する方法であるため、近所は形を変えなければなりません。20 km ^ 2が必要な場合は、形状またはホームテリトリーを変更します。[NoDataを無視する]チェックボックスをオンにしない場合、セル出力はNoDataになり、NoDataは役に立ちません。 2014年10月24日現在の「プログラム」 ShapelyとFionaを使用してこれまでに作成したコードは次のとおりです。 import numpy as np import pprint import shapely from shapely.geometry import* import fiona from fiona import collection import math traps = fiona.open('C:/Users/Curtis/Documents/ArcGIS/GIS_Data/occurrence/ss_occ.shp', 'r') study_area = fiona.open('C:/Users/Curtis/Documents/ArcGIS/GIS_Data/Study_Area.shp', 'r') for i in study_area: …

18 python  raster-calculator  spatial-statistics  shapely  fiona 

Pythonでは、GDALを使用して、次のようにラスターの投影をWKT文字列として抽出しました。 wkt = dataset.GetProjection() # wkt is 'PROJCS["GDA_1994_Transverse_Mercator",GEOGCS["GDA_1994",DATUM["GDA_1994",SPHEROID["GRS_1980",6378137,298.2572221010002],TOWGS84[0,0,0,0,0,0,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",117],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",10000000],UNIT["Meter",1]]' WKT文字列を使用して、次のようにSpatialReferenceインスタンスを作成できます。 src = osr.SpatialReference() src.ImportFromWkt(wkt) これは簡単です。srcUTMゾーンなどから投影のいくつかのパラメーターを簡単に抽出できます。しかし、投影の名前のようなものを抽出する方法がわかりません"GDA_1994_Transverse_Mercator"。これは確かに可能になりましたが、Python APIのドキュメントは、すべての用途において存在しない可能性があります。 投影の名前と地理座標系を抽出するにはどうすればよいですか？

18 python  gdal  coordinate-system  python-2.7 

無人航空機（UAV、無人機とも呼ばれます）があり、それを使用して、より小さなエリア（10 km2）でデジタル標高モデルとオルソモザイクを作成したいと思います。私のプラットフォームでは、高度が約120メートル（400フィート）の画像を重ねて表示できます。 いくつかの商用製品は完全なソリューションを提供しているようです。例えば。Pix4DおよびDronemapper。 しかし、私は主にプロセスをより細かく制御でき、場合によっては独自のバッチ処理も設定できるソリューションを探しています。Pythonやその他のオープンソースツールは興味深いものです。 私の選択肢は何ですか？

18 python  dem  unmanned-aerial-vehicle 

別のGISサーバー（GeoServerなど）からのWMS / WFSの呼び出しをサポートし、応答データ（WMS BasemapおよびWFSレイヤー）を写真として保存するためのAPIを備えたPythonオープンソースGISライブラリがあるかどうか疑問に思っています。 推奨事項はありますか？ 入力をありがとう！ 更新： 私がしようとしているのは、OpenLayersをフロントエンドとして、Djangoをサーバーとして使用する地図印刷サービスです。クライアントユーザーがエクステントとレイヤーを設定し、印刷要求（パラメーター、つまり、マップエクステント、レイヤーの名前を参照）をサーバーに送信し、サーバーがこの要求を引き継ぎ、要求パラメーターを使用してWMS / WFSを再度呼び出し、保存します応答をPDFとして、このPDFリンクをクライアントにエクスポートします。 困難な部分は、サーバーがWMS / WFSを呼び出してこれらの応答を結合/オーバーレイする方法です（つまり、WMSは通常ベースマップであり、WFSはフィーチャレイヤーを指すため、これらのマップ/レイヤーを結合する）。画像として。 現在の回答では、urllibは良いもののようですが、これらの応答（WMS、WFS）をどのように組み合わせるかはわかりません。OWSLibも別の良いオプションのように見えますが、それはクライアントプログラミングツールであることを示しています。 その他の入力 感謝する！

18 python  open-source-gis  wms  wfs 

この投稿はコミュニティWikiです。この投稿を改善するには、既存の回答を編集してください。現在、新しい回答は受け付けていません。 Pythonでいくつかのコードを開発していますが、GDAL / OGR pythonバインディングを使用したいと思います。GDALおよびOGR pythonバインディングを使用するためにインストールする必要のある最低限のものは何ですか？ Pythonからのみ使用したい場合、GDAL / OGRをインストールする最も簡単な方法は何ですか？ MacとWindowsの両方に当てはまる答えを探しています。 GDALはPythonパッケージインデックス（PyPi）にリストされていますが、その依存関係（libgdalおよびlibgdalとnumpyの両方のヘッダーファイル）をインストールする方法がわかりません。たぶん私がする必要があるのは、numpyをインストールしてからどこかでlibgdalを取得することです。Windowsのインストールには、libgdalとそのヘッダーファイルに加えて Windowsバイナリが必要ですか？ 使用しているバージョンのPython（2.5、2.6、2.7、3）に適切なPythonバインディングを確実に取得するにはどうすればよいですか？ PythonからGDALとOGRを使用するには、PATH環境変数を変更する必要がありますか、またはサイトパッケージフォルダーと一般的なimportステートメントからすべてにアクセスできますか？ これらの質問はこの質問に関連していますが、答えはしません 。MacOS Xでpython virtualenvにGEOS、PROJ、GDAL / OGR をインストールするWindowsでPythonを使用してGDALをインストールする方法は？ FWTools、OSGeo4W、およびkyngchaosはpythonバインディングでogrとgdalを提供し、優れていることに注意する必要がありますが、FWToolsとOSGeo4Wは既存のpythonインストールのライブラリとして機能するのではなく、独自のpythonとともに出荷され、kyng chaosフレームワークはOS Xのみ。

18 python  gdal  ogr 

Python OGR APIに少し問題があります。私がやろうとしているのは、ポリゴンの外側のリングの各頂点のすべての座標を取得することです。 これは私がこれまでに持っているものです： import osgeo.ogr import glob path = "/home/woo/maps/" out = path + 'output.txt' file = open(out,'w') for filename in glob.glob(path + "*.shp"): ds = osgeo.ogr.Open(filename) layer1 = ds.GetLayer(0) print layer1.GetExtent() for feat in layer1: geom = feat.GetGeometryRef() ring = geom.GetGeometryRef(0) points = ring.GetPointCount() #Not sure what to …

18 python  open-source-gis  ogr 

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コードが間違った結果を与えるのはなぜですか？ そして、なぜ古いコードがうまく機能するのでしょうか？ …

17 python  gdal  coordinate-system  epsg  web-mercator