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

Python用のgdalライブラリのインストールに問題があります。私はWindows 7を使用していて、Python 2.7.9を実行するWinPython 2.7.9.2を持っています。ラスター操作にgdal + Pythonを使用したいと思います。 私がしまし読ん PIPは、Windows上でパッケージをインストールするための好ましい方法であり、そのは、Python 2.7.9を詰めたということ。PATH環境変数を編集して、pipとpythonを含めました。 ;c:\~\WinPython-32bit-2.7.9.2\python-2.7.9\Scripts\;c:\~WinPython-32bit-2.7.9.2\python-2.7.9\ コマンドプロンプトから直接pythonとpipの両方を実行できるようになりました。 入力するとエラーが表示される c:\>pip install gdal それは与える： C:\~\AppData\Local\Programs\Common\Microsoft\Visual C++ for Python\9.0\VC\Include\xlocale(342) : warning C4530: C++ exception handler used, but unwind semantics are not enabled. Specify /EHsc extensions/gdal_wrap.cpp(2855) : fatal error C1083: Cannot open include file: 'cpl_port.h': No such file or directory …

9 python  gdal  installation 

スタンドアロンアプリケーションでopenlayersプラグインをロードしようとすると、パスが設定されていてプラグインが利用可能であると表示されていても、常にFalseで失敗します。 （Ubuntu 14.04、QGIS 2.4、python-2.7でテスト済み） #!/usr/bin/python # -*- coding: utf-8 -*- import sys, os import qgis import qgis.gui, qgis.utils from qgis.core import * from PyQt4 import QtGui # open an app app = QtGui.QApplication(sys.argv) # supply path to where is your qgis installed QgsApplication.setPrefixPath(u'/usr', True) # load providers QgsApplication.initQgis() # set …

9 qgis  python  qgis-openlayers-plugin  standalone 

Nodata値を管理する算術ルールは何ですか？これらのnull値の定義方法に関係なく？ArcMap Raster Calculatorを使用しています。 1 * Nodata（または1 * null）、1 / Nodata、1 + Nodata、1-Nodata、Nodata-1などの操作について具体的に考えています...

9 arcgis-desktop  python  arcgis-10.1  raster-calculator 

PythonスクリプトでGDALPolygonize（）を使用してラスターをポリゴン化しようとしています。スクリプトは、昨日の午後5時にポリゴン化を開始し、現在も午前9時30分にポリゴン化しています。どれだけ進んでいるかはわかりませんが、Windowsエクスプローラーを更新すると、出力シェープファイルのファイルサイズの変化を確認できるので、まだ進んでいることがわかります。 私のラスターはかなり大きいですが、それでもそれほど長くかかるとは思いません。私のラスターは35,486列、23,682行で、セルサイズは1メートルです。これは、値1がデータを表し、0がNoDataのバイナリラスタです。 ArcGISで、Conversion ToolboxのRaster to Polygonを使用してポリゴン化した場合、56秒かかりました。結果のシェープファイルは200MBですが、GDALPolygonizeによってまだ作成されているシェープファイルはまだ100MBしかありません。そのため、GDALは終夜実行した後、半分ほどの作業が完了したと思います。 仕様：Windows 7 64ビット、8 GB RAM、GDAL 1.10 64ビット、ArcGIS Desktop 10.2、64ビットバックグラウンドジオプロセシングfor ArcGISデスクトップ、Python 2.7.3 64ビット 更新 2日目-GDALPolygonizeはまだ実行中です。それは一晩で2晩続けて、終わらなかった。ArcGISは56秒かかりました。

9 python  gdal  polygonize 

QGISプラグインを開発したところ、[レイヤー]パネルでレイヤーが選択されていない場合、すべてのレイヤーを反復処理して何でもできるが、「アクティブ化」できることがわかりました。 誰かが私にヒントを教えてくれませんか？

9 qgis  python  qgis-plugins  pyqgis 

私はこのようなOpenLayersをウェブサイト持っている1ユーザーにポリゴンを描画する機能を提供し、その後にGeoJSONファイルを作成します。このGeoJSONをPythonスクリプトで処理したいと思います。Pythonスクリプトは、最終的にWebサイトに返され、情報ボックスとして表示される文字列を返します。 このプロジェクトを達成するには、どのようなフレームワーク/アプリ/ステップを検討する必要がありますか？

9 python  openlayers-2  geojson 

既存のラスターの空間投影に基づいてシェープファイルレイヤーを作成しようとしています。（最終的に、これはgdal.Polygonizeで使用され、ラスターデータから形状を取得します。）ただし、CreateLayerを使用しようとすると、エラーが発生します。 raster_uri-シェープファイルにしたいラスターファイルの場所。私はこれが存在することを知っています。 out_uri-場所新しいシェープファイルを印刷する場所。 layer_name-文字列のレイヤー名入力。 filed_name-文字列フィールド名入力。 これは私が使っていたコードです。 raster = gdal.Open(raster_uri) driver = ogr.GetDriverByName("ESRI Shapefile") ds = driver.CreateDataSource(out_uri) spat_ref = osr.SpatialReference() proj = raster.GetProjectionRef() spat_ref.ImportFromWkt(proj) layer = ds.CreateLayer(layer_name, spat_ref, ogr.wkbPolygon) 次のエラーメッセージが表示されます。 TypeError: in method 'DataSource_CreateLayer', argument 2 of type 'char const *' spat_refのタイプの出力が表示される<class 'osgeo.osr.SpatialReference'>ので、どのchar const *が不満を言っているのかわかりません。 興味深いことに、同じファイルを使用してpythonシェルからこれをすべて実行しようとすると、レイヤーは正常に作成されます。何が問題になるのかについてのアイデアはありますか？

9 python  gdal  ogr  osgeo 

私はpythonを使用してQgisのプラグインの開発を始めています。これは非常に新しいので、QGISにメッセージを表示して変数の結果を表示したいと考えています。これどうやってするの？

9 qgis  python  pyqgis 

地上レーザースキャナーを使用して、いくつかの森林区画の非常に高密度の点群を作成しました。次に、1.3メートルを超えるポイントを削除して、粗い木質の破片（枯れた枯れた木）を確認しました。添付されているのは、赤い楕円の内側に粗い木質の破片があるサンプルプロットの陰影付きDEMです。 プロットはまた、1.3メートル未満の木の幹の一部である小さな木、地面、小さな岩で構成されています。画像から、木質の破片はその連続的な形状で識別できます。この画像から木質の破片を抽出するツールを探しています。Arcmap、Envi、またはその他のオープンソースソフトウェアは完璧です。コーディングが役立つ場合は、Pythonの基本的な知識も持っています。

8 python  dem  point-cloud  forest-ecology 

PythonまたはPyQGISを使用して、多くのラスターをいくつかのピクセルサイズ、範囲、および投影システムで整列させたい。 最初に考えるのは、GDALを使用することです。 gdalinfo（参照ラスターから情報を検索するため） gdalwarp（投影システムとピクセルサイズを変換するため） gdal_translate（範囲を変換するため） それらは機能しますが、Pythonでこれらのツールを一緒に使用するのは簡単ではなく、この作業を完了するにはかなりの時間が必要です。また、QGISおよびAlign Rasters Tool を使用して簡単に行うことができます。 PyQGISまたはPythonを使用してこれを行うツールはありますか？（より多くのプログラミングを自動的に行いたい） 更新 私はこの質問からこのコードを見つけます： from osgeo import gdal, gdalconst inputfile = #Path to input file input = gdal.Open(inputfile, gdalconst.GA_ReadOnly) inputProj = input.GetProjection() inputTrans = input.GetGeoTransform() referencefile = #Path to reference file reference = gdal.Open(referencefile, gdalconst.GAReadOnly) referenceProj = reference.GetProjection() referenceTrans = reference.GetGeoTransform() bandreference …

8 qgis  python  raster  pyqgis  gdal 

どのようにして楕円を描くのshapelyですか？

8 python  shapely  ellipse 

複数の離散値を含むラスターを再分類できるGDALのツールを探しています。これまでのところ、GDALとPythonを使用してラスタを再分類しています。ただし、ここでは単一の値しか処理できないようです。 スイートのネイティブツールのようなものはありますか？

8 python  raster  gdal  reclassify 

問題は、下の図のように、一連の障害物を通過する信号（光や音など）の伝播をモデル化することです。信号は底面（地形）を通過できませんが、障害物を通過できます。通過した障害物の数を数えたい。 地形と障害物は2Dの数の多い配列（x、y、z）にあります。これが私がすることです： output = numpy.zeros(terrain.shape) obstacles = terrain + obstacle_heights for i in xrange (obstacles.shape[0]): for j in xrange (obstacles.shape[1]): mask = obstacles[i,j] > terrain[i,j:] output[i,j:][mask] +=1 結果は、[0, 0, 0, 1, 1, 1, 2, 3, 4, 4, 4 ...]行ごとのようなものになります。 この方法は正常に機能します（提供された地形の谷はを使用して埋められますnumpy.maximum.accumulate）。さて、ベクトル化されたソリューションを使用してスピードアップすることは可能でしょうか？

8 python  numpy 

意味のない名前のpdfが何千もあり、いくつかのフォルダに分散しています。名前を変更する必要があります。幸い、それらのほとんどは空間レイヤーのリンクを介して参照されます。属性テーブルにはファイルに関する意味のある情報があり、「hotlinks」という列にはpdfを指すディレクトリパスがあります。 属性テーブルのパスをたどり、レイヤーの属性に基づいてそのパス内のPDFの名前を変更できるQGISのメソッドはありますか？新しい名前を含む新しいフィールドがあります。concat（ "Field1"、 "Field2"、 "Field3"）なので、それに基づいてPDFの名前を変更します。 キャッチ... pdfは散らばっていて、同じフォルダーにないので、スクリプトは1つのフォルダーを指すだけではありません。

8 qgis  python  pdf 

Ubuntu 16.04にlandsat-utilをインストールしようとしています。UbuntuをWindows内の仮想マシンとして実行しています。依存関係の競合を回避するために、クリーンインストール（パッケージをインストールしない）を使用しています。 このページの指示に従いましたが、このコマンドを実行しようとすると毎回エラーが発生します。 $: pip install landsat-util これは出力からの抜粋です： Running setup.py bdist_wheel for pycurl ... error Complete output from command /usr/bin/python -u -c "import setuptools, tokenize;__file__='/tmp/pip-build-oH9kCk/pycurl/setup.py';exec(compile(getattr(tokenize, 'open', open)(__file__).read().replace('\r\n', '\n'), __file__, 'exec'))" bdist_wheel -d /tmp/tmpgxvEEBpip-wheel- --python-tag cp27: Using curl-config (libcurl 7.47.0) running bdist_wheel running build running build_py creating build creating build/lib.linux-x86_64-2.7 creating …

8 python  landsat  landsat-8