0

特定の緯度/経度座標が特定の米国の州に含まれるかどうかを識別しようとしています。パッケージを使用して、maps州境の緯度/経度座標を取得しました。座標が各州に含まれる緯度と経度の値の範囲内にあるかどうかを尋ねる方法があると思いますが、エレガントな解決策を見つけていません。

基本的に、これは、地域の名前を知らなくてもポイントが地域内にあるかどうかを識別する方法です (mapsパッケージは、地域の名前付けに関してかなり特異であるため、たとえば、ソ連など)。ただし、地域名を抽出する機能も備えています。後で(シェーディングなどのために)。

#Example coordinates
points=data.frame(lat=42,long=-100)
#Load US map
library(maps)
states=map_data("state")
#Summarize value of lat and long values for each region
library(plyr)
ddply(states,c("region"),function(x) return(c( range(x$long), range(x$lat) )) )

最後に、ポイントがその地域で発生するかどうかdata.frame statesを示す最後の列を用意したいと思います。true/falseしたがって、上記の例では、すべてのエントリregion=="nebraska"true.

4

0 に答える 0