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

私はほとんどすべてのpythonジオプロセシングにESRIのArcpyサイトパッケージを使用しているようです。ESRIの功績として、これらは非常に優れたツールスイートであり、多大な成果を上げることができます。ただし、ESRI Arcpyドメイン外のジオプロセシングスクリプトも作成したいと思います。たとえば、ラスターをポリゴンにクリップする場合、ESRIの次のスクリプトから始めます。 # Import system modules import arcpy from arcpy import env from arcpy.sa import * # Set environment settings env.workspace = "C:/sapyexamples/data" # Set local variables inRaster = "elevation" inMaskData = "mask.shp" # Check out the ArcGIS Spatial Analyst extension license arcpy.CheckOutExtension("Spatial") # Execute ExtractByMask outExtractByMask = ExtractByMask(inRaster, inMaskData) …

69 python  geoprocessing 

これは私にはそのような単純な質問のように思えます（おそらくそうです）が、答えを与える例を見つけることができないようです。PostGISを使用して、ポリゴンの外側にあるポイントを選択したいだけです。最終的に、これはST_Intersectsの逆であり、私が見る限りです。 例：タックスロットレイヤーと住所ポイントレイヤーがあります。ST_Intersectsを使用する必要があると思いますが、逆選択を行うようにするにはどうすればよいですか？次のコードの前にNOTステートメントを追加することも考えられましたが、うまくいきませんでした。 CREATE table t_intersect AS SELECT hp.gid, hp.st_address, hp.city, hp.st_num, hp.the_geom FROM public.parcel as par, public.housepoints as hp WHERE ST_Intersects(hp.the_geom,par.the_geom);

41 postgis  geoprocessing  sql  intersection 

次のワークフローでポイントシェープファイルを取得するように設計されたPythonコードがあります。 ポイントをマージ 互いに1 m以内のポイントが1つのポイントになるようにポイントを統合します フィーチャレイヤーを作成します。z<10のポイントが選択されています バッファポイント ポリゴンからラスター1mの解像度 再分類、ここで1-9 = 1; NoData = 0 各シェープファイルには、約5x7 kmをカバーする約250,000〜350,000ポイントがあります。入力として使用されるポイントデータは、ツリーの位置を表します。各ポイント（つまりツリー）には、クラウン半径を表す「z」値が関連付けられており、バッファプロセスで使用されます。私の目的は、最終的なバイナリ出力を別のプロセスで使用して、天蓋カバーを記述するラスターを作成することです。 4つのシェープファイルでテストを実行したところ、700MBのラスターが生成され、35分かかりました（i5プロセッサと8GB RAM）。このプロセスを3500個のシェープファイルで実行する必要があるので、プロセスを合理化するためのアドバイスをいただければ幸いです（添付コードを参照）。一般的に、ジオプロセシングビッグデータを処理する最良の方法は何ですか？より具体的には、効率を高めるのに役立つ可能性のあるコードまたはワークフローの調整はありますか？ 編集： ジオプロセシングタスクの時間（全体の％）： マージ= 7.6％ 積分= 7.1％ Lyrへの機能= 0 バッファー= 8.8％ ポリゴンからラスター= 74.8％ 再分類= 1.6％ # Import arcpy module import arcpy # Check out any necessary licenses arcpy.CheckOutExtension("spatial") # Script arguments temp4 = arcpy.GetParameterAsText(0) …

34 arcpy  geoprocessing  big-data 

arcpyに基づいていないPythonスクリプトを使用して、属性によるシェープファイルのクエリ、選択からの新しいレイヤーの作成、ポリゴンの面積の計算、ポリゴンのポイントへの変換などを行います。 他のPythonモジュールまたはライブラリを使用するコード例はありますか？私はarcpyを使用してこれを簡単に行うことができますが、他のオプションを調べたいと思いました。

33 python  geoprocessing 

シェープファイルを作成し、すべてのフィーチャの頂点が一定量シフトしているとしましょう。任意のシフトですべてのフィーチャをシフトする最も簡単な方法は何ですか（したがって、頂点の（x、y）位置）？この修正を適用するファイルがたくさんあるので、Bash / OGRの回答が望ましいでしょう:) 最後に、いい機能があるので、私はこのためにSpatialiteを使用することになりましたShiftCoords。しかし、スレッドは非常に有益でした！皆さんありがとう！

33 python  geoprocessing  open-source-gis  fiona 

この投稿はコミュニティWikiです。この投稿を改善するには、既存の回答を編集してください。現在、新しい回答は受け付けていません。 そこに提示されている地理データの一般化/単純化を自動化する既存のソリューション（ソフトウェア、ライブラリ、メソッドなど）を探しています。 地理データを簡素化するためのいくつかのツールを知って使用していますか？ 私はすでにこれらのものを見ました： maptimizeしてKMLマーカーレイヤーの表現を簡素化します。 ESRIによるモデルビルダー レーザースキャンによる半径の明瞭さ 形状のオンラインフィルタリングアルゴリズムを計算するmapshaper 軸システムによる軸 一般化し、複数の表現のウェブサイト上のICA委員会は、このトピックに関する膨大な情報源です。多くの論文では、このトピックに関連する高度な方法とユースケースを紹介しています。 非常に有名なダグラス＆peuckerフィルタリングアルゴリズムは、多くの場合、幾何学的な簡略化のために使用され、ほぼすべてのGISソフトに存在しています。ただし、その効果は、ジオメトリの形状を実際に単純化することなく、ジオメトリを非常に骨の多い/角度のあるものにすることです。たとえば、道路や川を表す線にはあまり適しておらず、滑らかで角張っていない必要があります。

25 geoprocessing  simplify  generalization 

LANに50台のコンピューターがあるとします。各コンピューターには、米国の特定の州にあるすべてのパーセルポリゴン用のジオデータベースがあります。 z $ / acre 未満の値を持つ別のパーセルのyフィート以内にあるx $ / acre 以上の値を持つすべてのパーセルを検出するジオプロセシングタスクを作成します。 データが50台のコンピューターに分散されていることを知らず、気にせずに、このクエリを作成して実行したいと思います。境界条件に留意してください。また、ある状態の高価な区画が別の状態の安価な区画に近いケースを返すクエリも必要です。 このような分散ジオプロセシングをサポートするアーキテクチャはありますか？ アーキテクチャは、抽象的に記述することも、AzureまたはAmazon Web Servicesに固有の実装として記述することもできます。または、できれば、豊富なArcGISデスクトップライセンスを使用して、コンピューターが夜間にアイドル状態になる典型的なオフィスとして。

24 geoprocessing  cloud-gis  amazon-web-services 

私が使用しているFUSIONとレーザー測量データを処理するためのコマンドラインFUSIONライダーツールキット（LTK）を。Googleの広範な検索（ "Lidar Python"）では、libLASおよびpyLASがPython LiDARライブラリとして生成されましたが、これらはLASデータへの読み取りおよび書き込みアクセスのみを提供するようです。 特に、点群からの天蓋表面モデルに加えて、強度と密度の画像を作成することに興味があります。FUSION LTKで可能な同じ種類のタスクを実行できる、Pythonで一般に受け入れられているツールのセットはありますか？

22 python  geoprocessing  lidar  las  fusion 

GISの中心概念は、データセットに関する質問に答えることです。データベースの観点から。空間拡張を伴うSQLは、このような質問をする方法です。機械可読テキストベースのフォームで質問を表現できる他の方法は何ですか？さまざまなアプローチの利点は何ですか？

22 geoprocessing  analysis 

他の多くの分野には、GISツールで潜在的に解決できる視覚的/「空間的」問題があります。私はGISツールを奇妙な非地理的タスクに使用したことがあり、これがどのくらい起こるのか、そしてどこでこれが行われたかの良い例があるかどうか疑問に思っていました。 私が興味を持っている種類の例： 地理情報技術を使用して、海洋ストロマトライトの堆積物粒子内の内石ボーリングプロセスの空間構造を定量化する [PDF] 方法論を要約するには： ストロマトライト（岩の一種）の挿し木の画像を撮影するために、顕微鏡が使用されました。これらの画像はBMPとしてエクスポートされ、Erdas Imagineが教師付き分類に使用されました。分類はArcView 3.xにエクスポートされ、クリーンアップされ、適切な分類がポリゴンに変換されました。その後、いくつかの空間選択と密度マッピングを行いました。 論文の図： ご覧のとおり、非地理データを使用しているが、まだ空間要素がある問題を探しています。GISの側面は、地図作成、視覚化、分析など、何でもかまいません。 科学論文/研究は特に興味深いですが、必須ではありません。

18 geoprocessing  cartography  analysis  visualisation  gis-principle 

現在、LANDSATデータを使用して、リモートセンシングコースに登録しています。 私のコースの一部は、放射輝度値に基づいてシーンを標準化することです。これは退屈な作業であるため、誰かがすでにそれを自動化する方法を見つけていると思います。 LANDSATシーンを標準化するツール、ESRIなどがありますか？

17 geoprocessing  grass  remote-sensing 

GISで時間属性を処理する方法について何か方法やアイデアはありますか？たとえば、オブジェクトが特定の時間間隔の場所にあり、それらが別の時間間隔の別の場所にあることを指定する方法です。たとえば、政治的境界が変更された場合などです。

16 geoprocessing  data  time  modelling  spatio-temporal-data 

これは、他の質問に対する明確な前兆であったはずです（質問しませんでした ）。PostGISでスパイダーダイアグラム（ハブライン）を作成する方法は？ レイヤーAのポイント（店舗）とレイヤーBのポイント（顧客）の関係がわからない場合は、一般的に「顧客1は最寄りの店舗からサービスを受けています」と言いたいと思います。私はこの事実が真実ではないかもしれないことを理解していますが、それはまともな代理になる可能性があります。 PostGISを使用して、レイヤーA（店舗）の最も近いポイントのIDをレイヤーB（顧客）の各ポイントに割り当てる最も効率的な方法は何ですか。私が探している出力は以下のようなものです。 Customer | Store 1 | A 2 | A 3 | B 4 | C

15 postgis  geoprocessing  business  nearest-neighbor 

PostGISにインポートする必要があるesriパーソナルジオデータベースがあります。 最終的には、WebサーバーにMDBをアップロードし、レイヤーをPostGISにインポートできるようにします。 Windowsでは、FWツールを使用して、以下が正常に機能します。 ogr2ogr -f "PostgreSQL" PG:"host=10.10.10.10 user=someuser dbname=poi password=somepassword port=5432" S:\GISData\Test.mdb -a_srs EPSG:26986 しかし、私はこれを自動化したい（そしてプログラム的にそれをしたい）ので、LinuxではLOTがもっと簡単になるだろう[と思う]。

15 postgis  open-source-gis  geoprocessing  esri-geodatabase 

FMEでは、複数のリーダーを単一のトランスフォーマーに接続する方法はありますか？私はFME 2015を使用しており、多数のフィーチャクラスを1つのフィーチャにクリップしたいのですが、これを行う簡単な方法が見つからないようです。 編集：入力データセットから複数のフィーチャタイプを選択し、それらをすべてフィーチャタイプの個々の出力矢印をクリックしてトランスフォーマの入力矢印にリンクすることなく単一のトランスフォーマにリンクしようとします。 編集2：これはスクリーンショットです、うまくいけばこれが役立つでしょう：

14 geoprocessing  fme  clip