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

それはばかげた質問のように見えますが、それについて明確な答えを見つけることができません：どのユニットがジオパンダ/シェイプを使用してオブジェクト間の距離/面積を計算しますか？

13 python  shapely  geopandas 

私は持っています GeoPandas 0.2.1 py27_0 fiona 1.7.0 np110py27_1 Python 2.7.10 Anaconda 2-4.1.1-Windows-x86環境にインストールされます。GeoDataFrame入力データセットを読み取ってデータを操作して構築できますが、出力データセットを保存しても座標系は保持されません。 import geopandas as gpd world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres')) world.to_file(driver='ESRI Shapefile',filename=r'C:\GIS\Temp\world_out.shp') world_out.shp特徴およびソースシェープファイルのスキーマを含むが、座標ないシステムが定義されていません（<Unknown>）。world_out.prjファイルが0キロバイトで、何が含まれていません。 out = gpd.read_file(r'C:\GIS\Temp\world_out.shp') print out.crs >>> {} 純粋なフィオナを使用してこの操作を実行するようにテストしましたが、正常に機能し、crs情報を保持します。 import fiona datasrc_path = r'C:\Program Files (x86)\ArcGIS\Desktop10.4\TemplateData\TemplateData.gdb' with fiona.drivers(): with fiona.open(datasrc_path,layer='cities') as source: meta = source.meta meta['driver'] = 'ESRI Shapefile' meta['layer'] …

13 python  fiona  geopandas  pandas  anaconda 

GDALをインストールしてコマンドラインで確認したところ、正常に機能していて、コマンドラインでgdal2tiles.pyを実行すると機能しますが、Python IDLEで同じファイルgdal2tiles.pyを実行すると、「osgeoという名前のモジュールはありません。 「」エラー。お手伝いできますか？

13 python  gdal  osgeo  osgeo4w  pyqt4 

最小バウンディングボックスに似たアルゴリズムを作成しようとしています（ただし、最終的には似たようなものにならない場合があります）。この場合、角度はパラメーターとして渡され、角度を指定すると、すべてのポイント/ポリゴンをカバーする最小の長方形が必要になります。これまでの私の考えは、ポイントの中心を見つけることです（セントロイドアルゴリズム）。そこから、パラメーターの角度と同じ角度の2本の平行線と、それらに垂直な2本の線を作成します。次に、反復を使用して、これらの線がすべてのポイントを含むまで外側（反対方向）に移動します。また、厳密な最小バウンディングボックスである必要はなく、おおよその動作です（各反復ステップのサイズに依存すると思います）。 これが私のコードです。すでにすべてのポリゴンを1つに分解しています。次に、凸包を使用して頂点を減らします。次に、すべての頂点をリストに入れます-これがまだ役立つかどうかわかりません... a = layer.getFeatures() for feat in a: geom = feat.geometry() a = geom.convexHull() vertexId = QgsVertexId() vertices = [] b = a.constGet().nextVertex(vertexId) while b[0]: vertices.append(b[1]) b = a.constGet().nextVertex(vertexId) 注：ある時点で、ボックスの角度を渡す必要があります。QGIS 3を使用していますが、これをPythonで作成する必要があります。レイヤー「レイヤー」には、他のすべてのポリゴンのディゾルブされたポリゴンという1つのジオメトリがあります。アクセスするために反復が必要ない場合があります。 詳細/情報をお伝えする必要がある場合はお知らせください。

12 python  pyqgis  geometry 

pip（pandas、scipy、simplekml）を使用していくつかのpythonライブラリをインストールしましたが、QGISをロードしようとするたびにこのエラーが発生します。 Warning: loading of qgis translation failed [/usr/share/qgis/i18n//qgis_en_US] Warning: loading of qt translation failed [/usr/share/qt4/translations/qt_en_US] Warning: QCss::Parser - Failed to load file "/style.qss" Warning: QVariantMap DBusMenuExporterDBus::getProperties(int, const QStringList&) const: Condition failed: action Warning: QVariantMap DBusMenuExporterDBus::getProperties(int, const QStringList&) const: Condition failed: action QH6248 qh_lib_check: Incorrect qhull library called. Caller uses reentrant …

12 qgis  python  installation  ubuntu 

現在、rasterioを使用してSentinel-2 AWS .jp2イメージファイルを開くと「アクセス拒否」エラーが発生します（現在はSinergiseによってホストされており、こちらにあります）。 私のコードは次のようになります。 import rasterio access_key = '*****************' secret_access_key = '***********************' region_name = 'eu-central-1' Session = rasterio.env.Env(aws_access_key_id=access_key, aws_secret_access_key=secret_access_key, region_name=region_name) url = 's3://sentinel-pds/tiles/10/S/DG/2015/12/7/0/B01.jp2' with Session: with rasterio.open(url) as dataset: print dataset 応答は次のとおりです。 --------------------------------------------------------------------------- CPLE_AWSAccessDenied Traceback (most recent call last) <ipython-input-1-4f6da7274649> in <module>() 12 13 with Session: ---> 14 with rasterio.open(url) …

12 python  remote-sensing  sentinel-2  rasterio  amazon-web-services 

以下のコードを使用して、数百万のGPSポイントの国（および場合によっては州）を見つけています。現在、コードは1ポイントあたり約1秒かかりますが、これは非常に遅いです。シェープファイルは6 MBです。 ジオパンダは空間結合にrtreeを使用しており、信じられないほど効率的であると読みましたが、これはここでは機能しないようです。私は何を間違えていますか？私は毎秒1000ポイントかそこらを望んでいました。 シェープファイルとcsvはここからダウンロードできます（5MB）：https ://www.dropbox.com/s/gdkxtpqupj0sidm/SpatialJoin.zip ? dl=0 import pandas as pd import geopandas as gpd from geopandas import GeoDataFrame, read_file from geopandas.tools import sjoin from shapely.geometry import Point, mapping,shape import time #parameters shapefile="K:/.../Shapefiles/Used/World.shp" df=pd.read_csv("K:/.../output2.csv",index_col=None,nrows=20)# Limit to 20 rows for testing if __name__=="__main__": start=time.time() df['geometry'] = df.apply(lambda z: Point(z.Longitude, z.Latitude), axis=1) …

12 python  spatial-join  geopandas  rtree 

経度と緯度の行を使用して、Pandas Data Frameからシェープファイルを構築したいと思います。 私はcsvファイルを持っているので、パンダで処理して、処理しやすいデータフレームを作成します 行ごとにループを作成せずにそれを行うことは可能ですか？

12 python  shapefile  pandas 

非常に大きなデータセットを持つQGISの一連の住所ポイントの最も近い道路からの距離を計算する最良の方法は何ですか？（QGIS 2.0.1、Win XP） 以前の同様の質問に対する回答を読みました。 QGISでの線と点の間の距離の測定 PyQGISで垂直線を描く？ QGISでポイントとポリゴン間の最小距離を計算しますか？ /programming/16448402/distance-of-point-feature-to-nearest-polygon-in-r それにもかかわらず： MMQGISプラグインは、ラインまでの最短（垂直）距離ではなく、ラインネットワークの最も近いノードまでの距離を検出します。 RとPythonが提案するスクリプトは、小さなデータセットで正しく機能し、モデラーに実装できますが、私の仕事では、530ポイントと約135'000のフィーチャを持つ地域のストリートレイヤーがあり、すべてのフィーチャをサイクリングしてから、すべてのポイントには膨大な時間がかかります QGIS 2.0.1でGRASSを動作させるのに問題があります（QGISのGRASSが動作しません（Windows XP））。

12 qgis  python  r  distance 

バッファーレイヤー（緑のポリゴン）があり、バリア（青い線）を横切るたびに2つのポリゴンに分割します。私は「splitGeometry」メソッドを使用しようとしましたが、動作させることができません。これまでの私のコードはこれです： while ldbuffprovider.nextFeature(feat): while barprovider.nextFeature(feat2): if feat.geometry().intersects(feat2.geometry()): intersection = feat.geometry().intersection(feat2.geometry()) result, newGeometries, topoTestPoints=feat.geometry().splitGeometry(intersection.asPolyline(),True) result（エラー）に対して1を返し、newGeometriesに対して空のリストを返します。どんな助けも大歓迎です。

12 qgis  python  pyqgis  splitting 

gdal-devメーリングリストからの移行： 2013年9月2日月曜日午後7時9分、David Sheanは次のように書いています。 こんにちはリスト、私は配布用の単一のNetCDFファイルとして同一の投影/範囲/解像度を持つGTiffラスターの時系列をパッケージ化しようとしています。過去1時間、オンラインドキュメントを参照し、gdal_translate、gdalbuildvrt、gdalwarpで遊んでみましたが、成功しませんでした。 既存のgdalコマンドラインユーティリティを使用してこれを行う簡単な方法はありますか？NetCDF Python APIを使用したカスタムソリューションに頼る前に、私は尋ねたいと思いました。 ありがとう。-デビッド 2013年9月3日火曜日、午前10時15分、エティエンヌトゥーリニーは次のように書いています。 あなたが望むのはおそらくgdalの範囲外です。gdal_translateがそれらを単一のファイルに入れるためには、巧妙なメタデータ管理が必要になります... gdal_translateを使用してそれらをすべてnetcdfに変換してから、python-netcdf4（numpy / scipyのものではない）を使用して時間ディメンションにスタックすることをお勧めします。 2013年9月3日火曜日午前7時55分、「Signell、Richard」は次のように書いています。 David、質問をGIS stackexchangeグループ/gis// に投稿する場合、 参考になるサンプルコードを提供します。 -リッチ ===================== 更新日9/3/13 17:04 PDT 入力データセットの1つに対するgdalinfoの出力を次に示します。 gdalinfo 20120901T2024_align_x+22.19_y+3.68_z+14.97_warp.tif Driver: GTiff/GeoTIFF Files: 20120901T2024_align_x+22.19_y+3.68_z+14.97_warp.tif Size is 10666, 13387 Coordinate System is: PROJCS["unnamed", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]], PROJECTION["Polar_Stereographic"], PARAMETER["latitude_of_origin",70], …

12 python  gdal  netcdf 

私はExcelのX、Yポイントをシェープファイルに変換するPythonコードを書いています。このプロセスでは、shhet1からxyポイントを読み取る必要があります。Excelワークブック（97-2003）のsheet1からデータを読み取るプロセスに含める必要があるステップは何ですか。 次のようにコードを書きました。 import arcpy from arcpy import env import Xlrd env.workspace="E:\" input_table="123.xls\Sheet1" しかし、機能していません。

12 python  excel 

私たちの組織は、ジオプロセシングワークフローをPostGISに移行することを検討しています。現在、ArcGISを使用しており、ModelBuilderで使用される多数のカスタムPythonツールを使用しています。ほとんどのデータをPostGISに移動して、さまざまなアプリで使用できるようになりましたが、現在、そこでデータ処理を実行することも理にかなっていますか。 ソフトウェアと互換性があるようにデータを処理します。顧客がソフトウェアを購入し、データを提供すると、ソフトウェアでの使用に最適化されるように処理します。これには、さまざまな品質の入力データを処理するためのさまざまなツールを構築する必要があります。特定の形式またはスキーマでデータを受け取ることは期待できないため、入力フィールドを出力フィールドにマップし、単一のフィールドを複数のフィールドに解析し、複数のデータセットをマージするなどのツールを構築します。フィールドの連結、および他の多くの一般的な操作。PostGISは、すべての処理ニーズを完全に実行できるようです。 PostGISを使用してデータ処理を行う場合、組織、使用するツールなどについて何かアドバイスはありますか？ QGIS python処理と組み合わせて使用​​しますか？ 非Web処理にPython ORMを使用している人はいますか？GeoGISにはPostGIS用のPython ORMがあるため、GeoDjangoの使用に傾倒しています。PostGISを使用してデータを処理する最初のテストでは、Pythonコードに多くの大きなSQLテキストブロックがあり、GeoDjango ORMが管理しやすく読みやすいコードの作成に役立つと考えています。また、PostGISと同様に対話するGeoAlchemy ORM があり、DjangoほどWeb固有ではないようです。 QGISやArcGISを使用している人と同じくらい、PostGISを使用してジオプロセシングを行っている人のことを聞いたことがありません。

12 python  postgis  geoprocessing  geodjango  geoalchemy 

マネージャーから、スクリプト/プログラミングを理解している従業員のための何らかのコンピテンシーメジャーを考案するように依頼されました。私は独学のプログラマーであるため、スクリプトの作成に関しては、あるレベルの理解を測定する方法について正式に理解していません。 ArcGISとQGISでPythonを使用しています。MS Office製品でVB.Netを使用し、WebマッピングのためにJavaに手を出し始めました。 スクリプトを作成する際に、従業員のスキルレベルを判断するための何らかのマトリックス/メソッドを使用しているGISマネージャーがいるのではないかと思っています。

12 python 

QGISの開発がここ数か月でどのような大きな一歩を踏み出したのか、本当に驚きました。新しいアイコンセット、改善された印刷コンポーザ、およびQGIS 2.0で見られる多くの変更は本当に素晴らしいです！ 他の多くのクールな改良点の中には、この気の利いた青いタスクバーもありますが、これは非常に魅力的です。 スクリプトやプラグインのためにPythonでこのバーに対処することは可能ですか（たとえば、計算終了後に小さなメッセージを表示するために）？多くのQGIS開発者もこのQ＆Aサイトを訪れていると思うので、ここでこの質問をしています。 PS：バーはQGIS開発でのみ表示され、一部の操作でのみ表示されます。

12 qgis  python  pyqgis