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

ファイルから印刷コンポーザーをロードし、アトラスを生成して画像にエクスポートするプラグインを構築しようとしています。これまでのところ、テンプレートの読み込みと画像へのエクスポートに成功しています。 凡例のレイヤー（tocにもあります）をエクスポートしたマップに追加できなかったため、マップが空白になり、フィールド式が機能しなくなりました。 # Get layers in the legend and append, must be a cleaner way to do this? layers = self.iface.legendInterface().layers() layerStringList = [] for layer in layers: layerID = layer.id() layerStringList.append(layerID) # Add layer to map render myMapRenderer = QgsMapRenderer() myMapRenderer.setLayerSet(layerStringList) myMapRenderer.setProjectionsEnabled(False) # Load template myComposition = QgsComposition(myMapRenderer) myFile …

10 python  qgis-plugins  pyqgis  qgis-2  print-composer 

ベクターレイヤーを追加しながら色を変更することはできますか？水の領域を持つベクターレイヤーがあり、青い色のスタイルで読み込みたいのですが。レイヤーは次のように追加されます。 QgsMapLayerRegistry.instance().addMapLayer(self.vlayer) ロード後またはロード中に色を変更できますか？

10 qgis  python  pyqgis  color 

USGSが提供するNational Elevation Datasetから高精度の標高（1/3-1/9アーク秒の解像度）データを見つけました。これは、int IMG、Arcgrid、およびGridFloat形式で提供されます。ArcGridが有料のソフトウェアパッケージに対応していることは知っていますが、無料で入手できるツールを使用することに専念しています。 地上レベルと相関させようとしているGPSデータがあります。 IMG、ArcGrid、またはGridFloat形式のデータファイルを使用してgpsデータを高度に変換できるpythonライブラリはありますか？

10 python  elevation  esri-grid-format 

次の行を使用して、ユーザーにステータスを通知します。 iface.mainWindow().statusBar().showMessage("Status:" + str(i)) プラグインが私のデータセットで実行されるのに約2分かかりますが、Windowsはプラグインを「応答なし」として検出し、ステータスの更新の表示を停止します。新しいユーザーにとって、プログラムがクラッシュしたように見えるので、これはあまり良くありません。 プラグインのステータスに関してユーザーが暗闇に放置されないようにするための回避策はありますか？

10 qgis  python 

時折、残念ながら残念ながら、私のAGSサービスは理由もなく停止します。ユーザーが私に電話して、地図にデータがない理由を尋ねるまで、私は通常これに気づきません。 サービスが停止したときにAGSから応答を受け取る方法（私はPythonを考えていました）はありますか？他のAGSユーザーはこの問題にどのように対処しますか？

10 python  arcgis-server  administration 

Pythonでcsvモジュールを使用して、pythonスクリプトと同じフォルダーにあるcsvファイルを開き、シェープファイルモジュールpyshpを使用してシェープファイルを作成する方法を理解しようとしています。 csvファイルは次のようになりますが、数千行のレコードを持つことができます。 id_nr;date;target;start_lat;start_lon 1;2012-05-21;navpoint 25x;55.123654;13.456954 1;2012-05-23;navpoint 11f;55.143654;12.456954

10 python  shapefile  csv  pyshp 

ワークスペースをループしてシェイプファイルの名前を変更し、フィーチャデータセットにコピーする便利なスクリプトツールがあります。ただし、ワークスペースのどこかに破損したシェープファイルがある場合、スクリプトは失敗し、処理が停止します。 このようなエラーをどのように処理しますか？エラーファイルを印刷し、forループの次のシェイプファイルの処理を続行して完了する方法はありますか？ import arcpy from arcpy import env # Allow overwriting of output env.overwriteOutput = True # Parameters env.workspace = arcpy.GetParameterAsText(0) state = arcpy.GetParameterAsText(1) gdb = arcpy.GetParameterAsText(2) # Get a list of shapefiles in folder fcs = arcpy.ListFeatureClasses() # Find the total count of shapefiles in list fcCount = len(fcs) …

10 python  arcpy 

次のコード（pyqgis クックブックの例から要約）を使用してベクターレイヤーを反復処理しているときに、機能が選択されているかどうかを確認する方法はありますか？ provider = vlayer.dataProvider() feat = QgsFeature() allAttrs = provider.attributeIndexes() provider.select(allAttrs) while provider.nextFeature(feat): geom = feat.geometry() attrs = feat.attributeMap() for (k,attr) in attrs.iteritems(): print "%d: %s" % (k, attr.toString()) または、選択した機能のリストをを使用して作成することもできますが、vlayer.selectedFeatures()各機能を直接確認する方法があることを期待しています。

10 qgis  python  pyqgis  select 

各点で不確実性の範囲が異なる方位角データを示す図を作成しようとしています。1991年の論文からのこのオールドスクールフィギュアは、私が目指している「ボウタイプロット」のアイデアを捉えています。 同様の図を作成する方法について何か提案はありますか？私はGISの初心者ですが、大学を通じてArcGISにアクセスできます。私のArcの経験は地質図の作成に限定されているので、あまりエキゾチックなことをする必要はありませんでした。 ArcとQGISのシンボルオプションをいろいろ調べてみましたが、その仕事をすると思われる設定は見ていません。これは、蝶ネクタイの形をしたシンボルを方位角で回転するだけの問題ではないことに注意してください。各「ボウタイ」の角度範囲は異なる必要があります。 私のPythonスキルは「強い中級」、Rのスキルは「中級」と評価しているので、必要に応じてmatplotlib、mpl_toolkits.basemapまたは類似のライブラリと一緒に何かをハッキングするのを嫌いません。しかし、聞いたことのないGISランドからのより簡単な解決策がある場合に備えて、その道を進む前にここでアドバイスを求めたいと思いました。

10 python  algorithm  r  visualisation  spherical-geometry 

クエリが失敗した後にOGRでPostgreSQL出力エラーメッセージをキャッチすることは可能ですか？ 例えば conn.ExecuteSQL('SELECT * FROM non_existing_table;') Postgresではこれはエラーメッセージを返します： "Error: relation 'non_existing_table' does not exist." ogrでは標準を取得します。 "Error: current transaction is aborted, commands ignored until end of transaction block" Postgresのエラーメッセージをキャッチ/解析できれば、おまけです。私はSWIG / CPLについて読んだことを考えると、これは少し長いショットだと思いますが、誰かがアイデアを持っているかどうかを確認するように頼む価値があります。理想的にはこれはPythonで行われますが、C / C ++コメントも使用します。 ありがとう。

10 python  postgis  gdal  postgresql  ogr 

ArcGISを使用せずにラスターから標高データを取得した経験があり、情報をpython listまたはdict？ タプルのリストとしてXYデータを取得します。 xy34 =[perp_obj[j].CalcPnts(float(i.dist), orientation) for j in range (len(perp_obj))] リストをループするか、それを関数またはクラスメソッドに渡して、xyペアの対応する高度を取得します。 私はこのトピックについていくつかの調査を行い、gdal APIは有望に聞こえます。誰もが私に物事、落とし穴、サンプルコードに取り組む方法をアドバイスできますか？ 作業中のマシンでシステムパス変数を編集できないため、GDALはオプションではありません。 誰かが別のアプローチについて知っていますか？

10 python  raster  elevation 

いくつかのテストを行うために、Ceroからマルチスペクトル画像を作成したいと考えています。塩とコショウのノイズが乗った5つの完全に均一なバンド、または中央に異なる値の2乗のような本当にシンプルなもの。明らかに、これは行列のスタック、つまり多次元配列であり、生成は非常に簡単です。私はpythonとgdalを使用してこれを実現したいのですが、gdalはかなり気密性が高く、コツがつかめません。geotiffファイルを作成するのが理想的です。誰かがこれを手伝ってくれませんか？いくつかのポインタまたは非常に穏やかなgdalチュートリアル？皆さん、ありがとうございました。

10 python  raster  gdal  imagery 

GeoPyを使用して住所を緯度、経度にジオコーディングしています。また、各住所の項目別住所コンポーネント（ストリート、市、州、郵便番号）を抽出したいと思います。 GeoPyはアドレスを含む文字列を返しますが、各コンポーネントを分離する信頼できる方法が見つかりません。例えば： {street: '123 Main Street', city: 'Los Angeles', state: 'CA', zip: 90034, country: 'USA'} GoogleジオコーディングAPIはこれらの個々のコンポーネントを返します...これらをGeoPyから取得する方法はありますか？（または別のジオコーディングツール？）

10 python  geocoding  google  address-parsing 

私はpython地理空間プログラミングを初めて使用します。次のスクリプトを実行して、対応するエラーメッセージを取得しました >>> import osgeo >>> import osgeo.ogr >>> shapefile = osgeo.ogr.Open("tl_2009_us_state.shp") >>> numLayers = shapefile.GetLayerCount() Traceback (most recent call last): File "<pyshell#5>", line 1, in <module> numLayers = shapefile.GetLayerCount() AttributeError: 'NoneType' object has no attribute 'GetLayerCount'

10 python  gdal 

特定の業界のすべての雇用主のオフィスの場所を表すポイントフィーチャクラスがあります。フィーチャクラスには、各オフィスで働いている従業員の数を格納するための属性があります。誰かがこのデータを使用するように要求しました。この場合、可能な限り最小の地理単位（国勢調査区）に空間的に結合されています。ただし、プライバシー契約により、データがそのままリリースされることはありません。代わりに、次の2つの条件を満たすように抑制する必要があります。 ポリゴンには少なくとも3つの雇用者（ポイント）が含まれている必要があります。 ポリゴン内の総雇用の80％以下は、単一の雇用者によるものです。 ポイントを国勢調査区に空間的に結合し、それぞれの合計と最大雇用を維持するスクリプトを作成できました。抑制基準を満たさないものにはフラグが立てられます。（ポイントを含まないポリゴンには、抑制するデータがないため、フラグが付けられません。）次に、各ブロックグループをチェックして、フラグが付けられたブロックが含まれているかどうかを確認します。フラグの付いていないブロックのみを含むブロックグループは、ブロックで置き換えられます。次に、結果のフィーチャクラスが抑制基準に対してチェックされ、ブロックグループがデータを適切に抑制しているかどうかがチェックされます。 同じプロセスがTractsでも繰り返され、Tracts（フラグ付きとフラグなし）、ブロックグループとブロック（すべてフラグなし）で構成されるデータセットが残ります。ただし、地理的階層の次の進歩は郡であり、このデータを要求する人には役に立たない。 次に、私の質問は次のとおりです。ポリゴンをできるだけ多くのグループに集約して、すべてがいくつかの最小基準を満たすように一般的に受け入れられている方法はありますか？ 集約に適用したいいくつかのルールを以下に示します。 フラグが付けられたトラクトは、可能な限り、他のフラグが付けられたトラクトとのみ集約する必要があります。 他と隣接していないフラグ付きトラクト（またはまだ基準を満たさない孤立したグループ）の場合、それらはすでに基準を満たしているトラクトと結合できますが、その間に雇用主がいないトラクトも存在する場合があります。含める必要があります。 絶対に不可能な場合を除いて、郡の境界をそのままにしたいと思います（入力フィーチャをそれぞれの郡に分離してから処理することで、これを行うことを期待しています）。 ソリューションは、ArcGISツールまたはオープンソースのPythonライブラリを使用して、Pythonで作成する必要があります。 理想的には、誰かがこの集約を実装する既存の手段を私に指摘することができます。そうでない場合は、自分でアルゴリズムをコーディングしてもかまいませんが、特定の手順/ツールのリストをいただければ幸いです。この問題は、（不連続なポリゴンを使用した）再区画の特殊なケースとして私を襲い、この目的のために、PySALの地域化アルゴリズムの使用を検討しましたが、これらを使用する全従業員の最大雇用者の割合を確認する方法はわかりません。

10 python  redistricting  regionalization  aggregation