とが NE lat lng
ありSW lat lng
ます。私の目標は単純です。RStudioを使用しており、上記の2つの緯度と経度から境界ボックスを作成し、ポリゴンを生成して、正方形ポリゴンにランダムポイントを追加します。data.frameにインポートできるcsvファイルにNELat / LngとSW-Lat / Lngがあります。
始めるには少しガイダンスが必要です。R言語は非常に新しいもので、これまで慣れていたものとはまったく異なります。
これまでのところ私はこれを持っています:
coords = cbind(78.46801, 19.53407)
coordsmax = cbind(78.83157, 19.74557 )
sp = SpatialPoints(coordsmax)
sp2 = SpatialPoints(coords)
r1 = rbind(coords, coordsmax[1, ]) # join
P1 = Polygon(r1)
Ps1 = Polygons(list(P1), ID = "a")
plot(Ps1)
他の場所でいくつかの例をみましたが、4つの角のポリゴンをプロットできません。
r1には単一の重複した座標しかありませんが、それは 'r1 = rbind(coords、coordsmax) `であることを意味していますか??
—
mdsumner
はい、修正しました。それでも動作しません
—
-user134611
問題を1つずつ慎重に処理することは、解決策を得るための良い方法です。2つまたは1つの座標のみを指定した場合、Polygon()がエラーにならないのはバグだと考えています。spの問題でこれを指摘しました。あなたq、library(raster); SPs1 <-as(extent(r1)、 "SpatialPolygons"); plot(SPs1)
—
mdsumner
stackexchangeの回答機能を使用する必要があることを残念に思いますが、私の評判はコメントするほど高くありません。@ dof1985(x_min、y_min)、(x_max、y_min)、(x_max、y_max)、(x_max、y_min)、(x_min、y_min)を使用してポリゴンを作成します。私は同じことを試みましたが、失敗しました。あなたの例を詳しく見て、(x_min、y_min)、(x_max、y_min)、...の代わりに使用することを調べて、(x_min、y_min)、(x_min、y_max)、...の可能性があります私のエラーの理由は?あなたの例を修正してください:)
—
ExploreR