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

オリジナルのArcGIS pythonインストールを「壊す」ことなく、arcpyを使用できる個別のPythonインストール（AnacondaやPython XYディストリビューションなど）を作成する方法を尋ねるQ＆Aがたくさんあります。 私はこれと同じ問題を抱えていましたが、「普遍的な」解決策を見つけたと思います。これまでのところ、作業中の約8台のコンピューターでこれを実行しましたが、すべて正常に動作します。 ArcGISを使用できる別のPythonディストリビューションをインストールする方法は？

28 arcpy  installation  python-2.7  anaconda 

Shapelyを新しいマシンで実行しようとしています。ソースからのインストール、Webページ上のインストーラーの使用、およびの使用の両方を試みましたpip install shapely。私がライブラリから作成しようとするすべての参照from shapely import *は、エラーをスローしますWindowsError: [Error 126] The specified module could not be found。Shapelyを活用しようとして、同様の経験があり、それを正しく機能させる方法を見つけましたか？ Pythonエラーがスローされました： WindowsError: [Error 126] The specified module could not be found 環境： Windows 7 64ビット Python 2.7

26 python  python-2.7  shapely 

私の質問は、ポリゴンシェープファイルの垂直線の延長です。最初にその質問を参照してください。 表示されるのは、ユーザー定義の間隔で、バウンディングボックスに対して垂直線を生成する方法です。OGR、Fiona、Shapelyなどを使用して次のクリッピングのステップを実行できることは理解していますが、それらの使用法は理解していません。 ポリゴンシェープファイルの1行を読み取るにはどうすればよいですか？Shapelyを使用するすべてのアプリケーションは、LineString、Point、またはPolygonを生成する方法を示しますが、既存のシェープファイルを読み取ることはありません 少なくともスケルトン構造をサポートしてくれますので、その上に構築できます。

23 python  polygon  line  python-2.7  clip 

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 

ラスターをポリゴンに変換するオープンソースのPythonソリューションを探しています（ArcPyは使用しません）。 ラスターをポリゴンに変換するGDAL関数を知っていました。これはマニュアルです：http : //pcjericks.github.io/py-gdalogr-cookbook/raster_layers.html#polygonize-a-raster-band それにもかかわらず、出力はファイルとして保存されず、一時的にメモリ内にある、形の良いポリゴンまたは任意のオブジェクトになる可能性があります。この問題を処理するパッケージまたはコードはありますか？ ラスターがnumpy配列で処理されている場合、アプローチは以下のとおりです。

13 python-2.7  shapely  polygonize 

このhttp://tinyurl.com/odfbanuのような本当に単純なシェープファイルを少しラスタライズする必要があります。これは、米国のシェープファイルクーティン郡です。GDALRasterizeLayerはすべてのポリゴンをラスターに書き込みませんか？しかし、Giopandasやfionaを使用して、tiffの記述部分にrastraioを使用する方法があるかどうか疑問に思っていました。 したがって、私の目標は、ラスタライズして、共通の値（例ではLSAD）を共有するすべてのポリゴンに値を割り当てることです。 だから私はスレッドでshongololoに触発されたコードの最初を書きました： Pythonで属性に基づいてポリゴンをディゾルブします（シェイプリー、フィオナ）？。 from geopandas import GeoDataFrame name_in = 'cb_2013_us_county_20m.shp' #Open the file with geopandas counties = GeoDataFrame.from_file(name_in) #Add a column to the Geodataframe containing the new value for i in range (len(counties)): LSAD = counties.at[i,'LSAD'] if LSAD == 00 : counties['LSAD_NUM'] == 'A' elif LSAD == 03 : …

10 python-2.7  fiona  geopandas 

QGIS Pythonコンソールは初めてです。コンソールからラスターレイヤーを読み込み、qmlファイルからcrsと特定のスタイルを設定しようとしています。 layerPath = os.path.join(raster_dir,raster_name) fileInfo = QtCore.QFileInfo(layerPath) baseName = fileInfo.baseName() layer = QgsRasterLayer(layerPath, baseName) 投影情報を設定してQGISにロードしたい QgsMapLayerRegistry.instance().addMapLayer(layer) 私は窓に灰色のラスターを持っています。ロードするqmlファイルを設定するにはどうすればよいですか？ ご協力いただきありがとうございます

9 qgis  python-2.7 

多くのShapely LineStringと交差するShapely Polygon / MultiPolygonを決定するために使用しているコードがあります。この質問への回答を通じて、コードは次のようになっています。 import fiona from shapely.geometry import LineString, Polygon, MultiPolygon, shape # Open each layer poly_layer = fiona.open('polygon_layer.shp') line_layer = fiona.open('line_layer.shp') # Convert to lists of shapely geometries the_lines = [shape(line['geometry']) for line in line_layer] the_polygons = [(poly['properties']['GEOID'], shape(poly['geometry'])) for poly in poly_layer] # Check for Polygons/MultiPolygons …

9 intersection  python-2.7  shapely  fiona  rtree 

USGS Earth Explorer Webサイト（http://earthexplorer.usgs.gov/）にPython でログインし、Landsat Archiveコレクションからユーザー名、パスワード、パス/行（または緯度/経度）、日付を入力するだけで画像をダウンロードしたい、および雲量のしきい値。私はすでにこれを試しました：https : //github.com/olivierhagolle/LANDSAT-Download/wiki うまくいきませんでした。それから私はこれを調べました：https : //github.com/developmentseed/landsat-util しかし、私は仕事の制限のためにツールを追加できなかったのでそれを使用できません。 そう。私はPythonの経験が限られているため、それを使用してWebサイトにアクセスしたことがありません。Windows 7でPython 2.7.8を使用していて、Chromeを介してWebサイトにアクセスしようとしています（バージョン51.0.2704.106 mが役立つ場合）。 Webサイトへのアクセス、ログイン（すでにログイン資格情報を持っている）、画像の検索、フォルダーへのダウンロードのための簡単なスクリプトのアイデアはありますか？Pythonでできることのように思えますが、私にとってオンラインソリューションはすべて高度すぎるようです。必要なのは、さまざまな座標や日付で繰り返し使用できるスタンドアロンスクリプトだけです。

9 python  python-2.7  landsat  web  earth-explorer 

pyhton、gdal、および/またはnumpyを使用して、10クラスのラスターから8クラスのラスターにラスターファイルを再分類したいと思います。クラスは整数として表されます。この投稿の手順に従って、GDALとPython、numpy.equal doc、およびgdal_calc doc を使用してラスタを再分類しました。しかし、役に立たない。 再分類されるラスターファイルには、0〜11の範囲の整数値が含まれ、値100および255も含まれます。次に、再分類（from value：to value）を示します。 nodata：4、0：4、1：1、2：2、3：3、4：3、5：4、6：5、7：5、8：6、9：7、10：8、100： nodata、255：nodata、 私ができることは、tkinter.FileDialogを使用して再分類するラスターファイルを選択し、retrans = gdal.Open（raster、GA_ReadOnly）でgeotransformやピクセルサイズなどのラスター情報を取得することです。 上記を解決するにはどうすればよいですか？ 場合によっては、再分類されるラスターがかなり大きくなる可能性がある（500MBから5GB）ことに注意してください。

9 gdal  python-2.7  numpy  reclassify 

ArcGISで1つの問題に直面しています。私はナビゲーションデータベースに取り組んでいます。私たちのデータベースでは、単一車線の道路は1本の線で表され、複数車線の道路（中央に仕切りのある道路）は2本の平行線（図では赤い色の線）で表されます。 私はいくつかのポイントがマルチレーンストリートの内側と外側にあるポイントシェープファイルを持っています。 マルチレーンストリート内にあるポイントを見つけるArcPyスクリプトを作成したいと思います。つまり、これらの平行線の間（写真でマーク）。 これを達成する方法がわかりません。誰かが私を助けてくれますか？ 私はそれにいくつかの練習をしました、そして私はラインの片側にバッファを作成することはマルチレーンポリゴン（写真で示される）の中に作成することができることを発見しました。 しかし今問題は、ポリゴンが実際にラインと交差している（つまり、マルチレーン境界に重なっている）ことです。不要なポイントをキャッチします。このポリゴンをストリートラインに合わせる方法はありますか？ 注：ストリートラインも移動するため、ここでは統合は機能しません。ストリートラインに沿ってポリゴンを揃えるだけです。

8 arcpy  arcgis-10.2  python-2.7 

2014年4月11日更新 機能の削除ツールでスクリプトがハングアップしていたようですので、以下の回答で提案されているように、テーブルの切り捨てに切り替えました。追加ツールから未使用の変数も削除しました。 2014年4月10日更新 このスクリプトを同僚のコンピューターで実行し（彼のコンピューターにはより多くのメモリがあり、ArcGIS 10.0 / Python26が含まれています）、すばやく実行されました。やったー！テクニカルサポートがArcGIS 10.0 CDを見つけたら、インストールしてテストし、コンピューターの速度が向上するかどうかを確認します。明確にするために、同じスクリプトを実行しており、ネットワークドライブとデータベース接続は同じようにマップされており、印刷ステートメントは同じです。それが起こったら、私はここにアップデートを投稿します。 更新を終了 Oracleデータベースで更新を実行する一部のPythonスクリプトの速度を上げる必要があります。スクリプトを開始するためのスケジュールされたタスクとバッチファイルを介して、これらのPythonスクリプトを1年以上実行しました。先週、XPからWindows 7マシンに移行し、ArcGIS 10.0-> 10.1に移行しました。それ以来、スクリプトはひどく遅くなっています。小さなフィーチャクラス（〜20のフィーチャを含む）を使用してこのスクリプトを実行すると、30秒で実行されます。中程度のフィーチャクラス（約80,000レコード）を使用すると、15分で実行されます。すぐに転送できるようにする必要があるフィーチャクラスには、約1,000,000レコードが含まれています。スクリプトは、ファイルが存在するかどうかを確認するための印刷ステートメント（以下のコードのifステートメント）までしか機能しません。このプロセスは、私のXP / ArcGIS 10.0コンピューターで完了するのに35分しかかかりません。 以下は、私がテストしてきた簡略化されたコードです。速度を上げるために私ができることについて誰かが提案をしていますか？ありがとう、パティ import arcpy, os, sys from arcpy import env arcpy.env.overwriteOutput = True from datetime import datetime import smtplib import string import urllib #Define variables inWorkspace = "O:/LANDING_PAD/BOE/example.gdb" lpFeatures = inWorkspace + os.sep + …

8 arcgis-10.1  python-2.7  windows-7 

Windows 7を実行しているマシンにGRASS 6.4.3 svnをインストールしています。また、マシンにWindows用のPython 2.7.3をインストールしました。次に、GRASS GUI自体のPythonシェル内ではなく、PythonからGRASSライブラリを呼び出す必要があります。env.batこのために環境変数を設定しました。しかし、接続は失敗します。これはenv.bat私が持っているものです： set GISBASE=C:\Program Files\GRASS GIS 6.4.3svn set HOME=%USERPROFILE% set LD_LIBRARY_PATH = %GISBASE%\lib set PATH= %GISBASE%\etc;%PATH% set PATH=%GISBASE%\Python27;%PATH% set PATH=%GISBASE%\lib;%GISBASE%\bin;%GISBASE%\extralib;%GISBASE%\msys\bin;%PATH% set PATH=C:\Python27;%PATH% set PYTHONLIB= C:\Python27;%PYTHONLIB% set PYTHONPATH= %GISBASE%\etc\python;%PYTHONPATH% rem Environmental variables for GRASS stand-alone installer set GRASS_WISH=%GISBASE%\extrabin\wish.exe set GRASS_PYTHON=python set GRASS_PROJSHARE=%GISBASE%\proj set GRASS_HTML_BROWSER=explorer set GRASS_SH=%GISBASE%\msys\bin\sh.exe set …

8 grass  python-2.7