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

Shapelyを使用して、以下に示すすべての重複しないポリゴン（与えられたポリゴンA、B、C）をキャプチャしようとしています。さらに、反復や交差のテストなどを行わずにそうすることを望んでいます。この質問に対する受け入れられた答えは、PostGISメソッドを表していますが、「ユニオン」は人によって異なることを意味するようです。

15 python  polygon-creation  shapely 

バッファ、ユニオン、インターセクトなどのジオメトリ操作を実行するための純粋なPython（CまたはJava依存関係なし）ライブラリを知っている人はいますか？Google AppEngineと互換性のあるものを探しています。

15 python  open-source-gis 

シェープファイルをインポートするためにPythonでGDALを使用します。ファイルのフィールド名を知りたいのですが、現在の方法は次のとおりです。 fields = [] for i in range(1, layer.GetFeature(0).GetFieldCount()): field = layer.GetFeature(0).GetDefnRef().GetFieldDefn(i).GetName() fields.append(field) しかし、この方法で、最初のレイヤーの機能を取得しています。レイヤーごとに機能が異なる可能性があるということですか？ そうでない場合、これほど深くに入る代わりに、フィールド名を一度に取得することは可能ですか？はいの場合、フィールド名を取得する簡単な方法はありますか？

15 python  shapefile  gdal 

PythonでLASデータを操作する方法を学び始め、他のユーザーがLASファイルをどのように処理するかを知りたいと思いました。ポイントを読み取り（numpy配列を使用しています）、クラス1および2（未分類およびグラウンド）を別の配列にフィルター処理します。私は次のコードを持っていますが、ポイントをフィルタリングすることはできません。 # Import modules from liblas import file import numpy as np if __name__=="__main__": '''Read LAS file and create an array to hold X, Y, Z values''' # Get file las_file = r"E:\Testing\ground_filtered.las" # Read file f = file.File(las_file, mode='r') # Get number of points from header num_points = int(f.__len__()) …

15 python  lidar  numpy  array  liblas 

Windows XPでPython 2.7.3のgdalをインポートしようとしています： >>> import gdal Traceback (most recent call last): File "<stdin>", line 1, in <module> File "C:\Python27\lib\site-packages\gdal.py", line 2, in <module> from osgeo.gdal import deprecation_warn File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 21, in <mo _gdal = swig_import_helper() File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 17, in swit_helper _mod = imp.load_module('_gdal', fp, pathname, description) ImportError: DLL …

15 python  gdal  windows-xp 

Pythonで機能ごとにシェープファイルを分割することは可能ですか？（最良の解決策は、結果のベクトルオブジェクトを一時的にディスクではなくメモリに保存できるソリューションです）。 理由：シェープファイルのいくつかの異なるサブセットでgdal rasterizeLayer関数を使用したい。この関数には、osgeo.ogr.Layerオブジェクトが必要です。 まあ、私は少し試してみましたが、次のように機能します。次のように、フィーチャごとにgdalレイヤーオブジェクトのジオメトリを取得できます。 # Load shape into gdal shapefile=str(vectorPath) layer_source = ogr.Open(shapefile) lyr = layer_source.GetLayer(0) for i in range(0,lyr.GetFeatureCount()): feat = lyr.GetFeature(i) ge = feat.geometry() このジオメトリに基づいてosgeo.ogr.layerオブジェクトを作成する方法を知る必要があります。 明確にするために。プレーンogr / gdalコードの関数が必要です！これは他の人にも興味があるようであり、セカンダリモジュールなしのソリューションが引き続き必要です（ただし、ここから来るソリューションは無料のqgis​​プラグインで使用されます）。

15 python  gdal  rasterization  osgeo  split-by-attribute 

経度、緯度、およびこのポイントの3番目のプロパティ値を持つデータポイントを指定します。プロパティ値に基づいてポイントをグループ（地理的サブリージョン）にクラスター化するにはどうすればよいですか？私はグーグルで検索し、この問題は「空間的制約付きクラスタリング」または「地域化」と呼ばれているようだと考えました。しかし、私は地理データの処理に精通しておらず、どの種類のアルゴリズムが適しているか、このタスクにどのpython / Rパッケージが適しているかについてはわかりません。 私が欲しいものについてより直感的なアイデアを与えるために、私のデータ散布図は次のようだとしましょう： したがって、各ドットは点、xは経度、yは緯度、カラーマップは値が大きいか小さいかを示します。値の位置と類似性に基づいて、これらのポイントをサブリージョン/グループ/クラスターに分割したいと思います。次のように（直感的なアイデアを示すためだけではありません）。 どうすればこれを達成できますか？

15 python  r  clustering  grouping 

私は、QGISが機能を「溶解」するのと基本的に同じことを行う関数を作成しようとしています。私はそれは非常に簡単だと思ったが、明らかにそうではない。だから、私が集めたものから、フィオナとshapelyの使用がここでの最良の選択肢であるはずです。ベクターファイルをいじり始めたばかりなので、この世界は私にとってもPythonにとっても新しいものです。 これらの例では、http： //tinyurl.com/odfbanuでここに設立された郡のシェープファイルを使用して います ：今の私の最善の方法は、に基づいて、以下の通りであるhttps://sgillies.net/2009/01/27/a-more-perfect-union-continued.html。それはうまく機能し、52の状態のリストをShapelyジオメトリとして取得します。この部分を行うためのより簡単な方法がある場合は、お気軽にコメントしてください。 from osgeo import ogr from shapely.wkb import loads from numpy import asarray from shapely.ops import cascaded_union ds = ogr.Open('counties.shp') layer = ds.GetLayer(0) #create a list of unique states identifier to be able #to loop through them later STATEFP_list = [] for i in range(0 , …

15 python  shapely  fiona  geopandas 

データフレームのリストがあります。すべて同じ構造で、すべてポリゴンです。 ジオパンダを使用してそれらをすべて一緒に追加するにはどうすればよいですか？

15 python  geopandas 

rasterioを使用してラスター内のポイントで単一のピクセル値を取得するには、https：//github.com/mapbox/rasterio/pull/275に例があります。 しかし、ラスタ内の単一のポイントで値を抽出するために使用できる、ラスタリオ（CLIではなく）内に直接APIがありますか？ -編集 with rasterio.drivers(): # Read raster bands directly to Numpy arrays. # with rasterio.open('C:\\Users\\rit\\38ERP.tif') as src: x = (src.bounds.left + src.bounds.right) / 2.0 y = (src.bounds.bottom + src.bounds.top) / 2.0 vals = src.sample((x, y)) for val in vals: print list(val)

14 python  numpy  rasterio 

フロントエンドにリーフレット付きのdjangoアプリがあり、座標を渡すことができるようにgeojsonオブジェクトのbboxを計算する必要があります map.fitBounds() 私はいくつかのライブラリを試し、getBoundsリーフレットの特定のオブジェクトのメソッドを試しましたが、などFeatureGroup、境界が定義されていないことを訴えました。 誰も私にこれの簡単な解決策を教えてくれますか？

14 python  leaflet  geojson  django 

Rを使用してtwitterをデータマイニングすることに関するいくつかの研究を行ってきましたが、私の質問に対する答えや適切なチュートリアルは実際には見つかりませんでした。 特定の時間枠内で、特定のハッシュタグを使用してTwitterからツイートを取得し、QGISまたはArcMapのいずれかのマップにそれらのツイートの場所をプロットすることに興味があります。 ツイートにはジオロケーションが関連付けられている可能性があることは知っていますが、そもそもこの情報を抽出するにはどうすればよいですか？

14 python  point  twitter 

私はPyQGISクックブックのジオメトリ処理に関するページにアクセスして再訪しました：http ://documentation.qgis.org/2.0/en/docs/pyqgis_developer_cookbook/geometry.html しかし、 Pythonコンソールから描画するポリゴン。誰でも助けることができますか？

14 python  geometry  qgis-2 

リモートセンシングのスペシャリストとしての長いキャリアにつながることを願って研究を始めています。現在、一部のアプリケーションではArcGISを使用し、他のアプリケーションではENVIを学習しています。私はプログラミング言語を学ぶことが不可欠であることに気付き、IDLとPythonの選択に直面しています。コミュニティが土地表面プロセスのリモートセンシングに推奨するプログラミング言語を聞いてみたいです。

14 python  remote-sensing  software-recommendations  idl 

ジオタグ付きの写真があるので、簡単に参照できるように、また座標からシェープファイルを作成するために、写真の場所にウォーターマークを前面に配置したいと思います。 誰でもこれを行うことができるオープンソースソフトウェアを知っていますか？ このリンクは、Iamが達成しようとしているものの例です

14 python  open-source-gis  geotag  exif