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

ブールラスターがあります。 ラスターの灰色の領域で、特定のサイズのポリゴンを連続した範囲内に収めたいと思います。 基本的に、不規則なポリゴンがあります。可能な限り、不規則なポリゴンの範囲内で既知のポリゴンを「適合」させたいと思います。 多角形の方向は重要ではなく、正方形でもかまいません。グラフィカルにフィットするようにしたいのですが、ポリゴンに数字を付けるだけで（フィットする＃）、同様に機能します。 ArcGIS Desktop 10を使用しています。

25 arcgis-desktop  raster  arcgis-10.0  polygon 

プロットしたいSpatialPolygonsまたはSpatialPolygonsDataFramesがあります。Rの3つのプロットシステム（ベースグラフィックス、ラティスグラフィックス、およびggplot2）を使用して色を変更するにはどうすればよいですか？ サンプルデータ： library(sp) Srs1 = Polygons(list(Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))), "s1") Srs2 = Polygons(list(Polygon(cbind(c(5,4,2,5),c(2,3,2,2)))), "s2") SpDF <- SpatialPolygonsDataFrame( SpatialPolygons(list(Srs1,Srs2)), data.frame( z=1:2, row.names=c("s1","s2") ) ) spplot(SpDF, zcol="z")

25 polygon  r 

ArcMap 10に読み込まれたフィーチャクラスには約12個のポリゴンがあり、すべて地理WGS 1984にあります。 そのフィーチャクラスの各ポリゴンの各頂点に関連付けられた座標を簡単に取得するにはどうすればよいですか？ 理想的には、スプレッドシート形式でうまく表にしたいです。

25 arcgis-desktop  arcgis-10.0  arcmap  polygon  vertices 

無効なポリゴンジオメトリから発生する「TopologyException：Input geom 1 is invalid」自己交差エラーは広く議論されています。しかし、R機能だけに依存する便利な解決策をWeb上で見つけられませんでした。 たとえば、私はここでmap("state", ...) Josh O'Brienの素敵な答えを次の出力から「SpatialPolygons」オブジェクトを作成することができました。 library(maps) library(maptools) map_states = map("state", fill = TRUE, plot = FALSE) IDs = sapply(strsplit(map_states$names, ":"), "[[", 1) spydf_states = map2SpatialPolygons(map_states, IDs = IDs, proj4string = CRS("+init=epsg:4326")) plot(spydf_states) この広く適用されているデータセットの問題は、以下のポイントで自己交差が発生することです。 rgeos::gIsValid(spydf_states) [1] FALSE Warning message: In RGEOSUnaryPredFunc(spgeom, byid, "rgeos_isvalid") : Self-intersection at or …

24 r  polygon  rgeos  self-intersection 

特定のポリゴン境界内に一連のラインフィーチャがあります。各ラインについて、レイヤー内の他のどのラインよりも可能なすべてのポイントが指定されたラインに近いポリゴンを生成したいと思います。過去にDelaunayの三角形分割を使用してポイント入力フィーチャに対してこれを実行しましたが、ラインフィーチャを使用して同様のプロセスを実行した場合、それを見つけることができませんでした。 ETA：Geogeekのソリューションは私に思いつきましたが、入力ラインの頂点の数が少ない直線セクションでは、結果のポリゴンが、本来あるべきではないラインに近づきすぎます（重なっている場合もあります）。ここで、赤い線は私の入力であり、頂点とそこから生成されたティーセンポリゴンを見ることができます。 おそらく、迅速で（非常に）汚い解決策は、各ラインを（ラインの頂点のみではなく）等間隔の多数のポイントセットに変換し、それらからティーセンポリゴンを生成し、元のラインIDに基づいてそれらをディゾルブすることです。

24 qgis  arcgis-desktop  arcgis-10.1  polygon  proximity 

北米の国境をいくつかの変数を表すラスターイメージ上にプロットし、Rを使用してプロットの上に等高線を重ねたいと思います。ベースグラフィックスとラティスを使用してこれを行うことに成功しましたが、遅すぎる！私はまだggplot2でこれを行っていませんが、速度の点でより良くなるとは思いません。 gribファイルから作成されたnetcdfファイルにデータがあります。今のところ、私はからRDATAファイルで使用可能だったカナダ、米国、メキシコ、のために国境をダウンロードGADM SpatialPolygonsDataFrameがオブジェクトとしてRに読み込みます。 コードは次のとおりです。 # Load packages library(raster) #library(ncdf) # If you cannot install ncdf4 library(ncdf4) # Read in the file, get the 13th layer # fn <- 'path_to_file' r <- raster(fn, band=13) # Set the projection and extent p4 <- "+proj=lcc +lat_1=50.0 +lat_2=50.0 +units=km +x_0=32.46341 +y_0=32.46341 +lon_0=-107 +lat_0=1.0" projection(r) …

24 polygon  r 

重複するラスターレイヤー（深さ、電流、波などの物理レイヤー）の束のクラスター分析に基づいて、ポリゴンからn個のサブ領域を作成しようとしています。 現在、ポリゴン全体に通常のグリッドを作成し、物理ラスターレイヤーから物理属性を抽出し（例：Gridspotまたは同等のツール）、n個のクラスターに制限されたクラスター分析を実行できます（Rまたはその他の統計パッケージ内）。 次に、各クラスターグループを特定し、それらをGIS（QGISまたはArcMap）にプロットします。ただし、一部のクラスターグループは分散され（空間的外れ値であり、意味のない）、一部はクラスター化される（サブリージョンとしての価値がある）ことを想定しています。 次に、代表的なクランプを手動で描画して、n個のサブ領域を作成できます。 QGISで実行できるArcMap 10.1 Grouping Analysisのようなツールはありますか？10.0のみを実行しています。 n個のサブ領域（バイオ領域）を作成するために、複数のラスターレイヤーのクラスター分析を行うより良い方法の提案はありますか？

23 qgis  raster  polygon  clustering 

2日前に、ESRIシェープファイルのポリゴンの頂点の内部保存順序について質問しました。その質問は回答されました（ポリゴンはシェイプファイルに時計回りまたは反時計回りに保存されていますか？）、また古い投稿で回答されました（ポリゴン作成（時計回りの回転かどうか）） しかし、今では私の質問はより一般的であり、独自の答えがあるかどうかはわかりません。時計回りの順序はESRIシェープファイルのみですか、それとも一般的なGIS形式ですか？また、GISソフトウェアの内部表現についてはどうですか？たとえば、QGISを使用し、ポリゴンを含む* .shpを読んだ場合、外側の境界の内部表現は元のシェープファイルのように時計回りであると仮定しますが、QGISでサポートされているすべてのファイル形式はどうでしょうか？ArcGISの場合はどうですか？また、ポリゴンが反時計回りに保存されているファイル形式が存在する場合、これらのファイルがQGIS、ArcGISなどにロードされると、内部で方向が変更されるため、たとえばPyQGISを使用してデータを読み取ると、ポリゴンは時計回りになります順序付けられました？ 私の目的はQGISのプラグインを作成することですが、データのソースはESRIシェープファイルまたはその他の形式にすることができます。方位角を使用してポリゴンの連続する側面間の角度を確認する必要があるため、順序が時計回りであるかどうかを知る必要があります。1つの解決策は、各ポリゴンの面積を計算することです。正しく覚えていれば、正の場合は順序は時計回りになり、負の場合は順序は反時計回りになります。 面積計算は集中的なタスクではないため、プラグインをそれほど遅くすることはありません。しかし、QGISの特殊なケースでは、元のソースの順序に関係なく、ポリゴンを時計回りに保存するか反時計回りに保存するかを誰もが知っていますか？これまでにESRIシェープファイルを使用しており、layer.getFeatures（）。geometry（）。asPolygon（）の座標は、外側の境界線では時計回りに、穴の場合は反時計回りに、つまり元の* .shpのように保存されます。

23 shapefile  polygon  vertices  specification 

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

23 python  polygon  line  python-2.7  clip 

IDが重複する場合にSPDFをrbindingするための適切なRイディオムは何ですか？ここで（多くの場合そうであるように）IDは基本的に無意味であるため、rbindにそれらを無視させることはできません。 library(sp) library(UScensus2000) library(UScensus2000tract) data(state) # for state names states <- gsub( " ", "_", tolower(state.name) ) datanames <- paste(states,"tract", sep=".") data( list=datanames ) lst <- lapply(datanames,get) nation <- do.call( rbind, lst ) Error in validObject(res) : invalid class “SpatialPolygons” object: non-unique Polygons ID slot values # This non-exported function …

22 polygon  r 

したがって、Rでシェープファイルを読み取る通常の方法は、次のようにmaptoolsパッケージを使用することです。 sfdata <- readShapeSpatial("/path/to/my/shapefile.shp", proj4string=CRS("+proj=longlat")) ただし、shapefile.shpがなく、代わりに一連のポリゴン座標があるユースケースがあります 16.484375 59.736328125,17.4951171875 55.1220703125,24.74609375 55.0341796875,22.5927734375 61.142578125,16.484375 59.736328125 およびその対応する投影 coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 このデータからsfdata（「ポリゴンオブジェクト」）を直接「インスタンス化」するにはどうすればよいですか？（これらのデータを使用してシェープファイルを作成し、新しく作成されたシェープファイルから読み取るというラウンドアバウト方法を使用することなく）

22 shapefile  polygon  r 

PostGISの特定のポリゴン内にあることが保証されているポイントを見つけるにはどうすればよいですか？ 私はST_Centroid機能を知っています。ただし、重心は常にポリゴン内にあるとは限りません。以下を参照してください。 さらに、多角形の境界上にある点の使用を避けたいのですが、境界内にある（ドーナツ型の多角形の穴内にはない）点が必要です。

22 postgis  polygon  point  geometry  point-in-polygon 

対応するLocal Enterprise Partnershipsテーブルデータを、英国の地方自治体を含むポリゴンシェープファイルに結合しました。 スクリーンショットに示されているように、LEP_NAMEによってLocal Authorityポリゴンを溶解またはマージするにはどうすればよいですか？ 理想的には、結果のマージされたポリゴンでLocal Authority属性を保持したいと思います。スクリプトに関する知識がないため、ツールベースのソリューションを探しています。どんな助けやポインタも大歓迎です。 追加の注意：ディゾルブツールを使用しようとすると、上記のようにLEP_NAMEフィールドは使用できません。LEP属性データはExcelで作成され、CSV形式として保存され、Local Authorityポリゴンシェープファイルに結合されました。他のスクリーンショットに示されているように、IDENTIFI0フィールドが重複しています-これが問題の原因である可能性がありますか？

22 qgis  polygon  fields-attributes  merge  dissolve 

架空のマップを作成しているので、たくさんのポイント、ライン、そして確かにポリゴンを作成する必要があります。その後、データをgeojsonとしてエクスポートします。しかしその前に、私は離れてすべての要素に一意のIDを与えなければなりません。 最大のポリゴンが最小のIDを取得するなど、特別な並べ替えは必要ありません。私は必要なすべてのポリゴン とIDを、終わりにせずに やって、それを手動で私が今しなければならないように。 誰かがそれを行う方法を知っていれば素晴らしいでしょう。

22 qgis  polygon  unique-id 

この投稿では、凸多角形内で最大面積長方形を見つける方法に関するアルゴリズム / アイデアを探しています。 次の図で、数字はフィットされた長方形の面積です。示されているように、望ましい長方形は各次元で異なり、任意の角度にすることができます。 編集： 言及された問題に対処する方法が明確にわからないため、ここで質問しました。それでも、maximum-area-rectangle は、ポリゴンのエッジ上に1つのエッジが（必ずしも同じ長さのエッジである必要はありませんが）配置されているものの1つであると推測されます。

21 polygon  geometry  algorithm