5

パッケージでSpatialPolygonsDataFrame使用するシェープファイルを読み込んで作成したがreadOGRあります。パッケージを使用して、その地域で収集された調査データからの補間に使用されるrgdalサンプリンググリッドを生成するために使用しようとしています。ただし、SpatialPointsDataFrame は調査よりもはるかに広い領域を網羅しているため、内挿は調査が実施された場所から遠く離れた値を予測しています。調査データとシェープファイルは両方とも、同じ proj4string を使用して投影されました。spsamplesp

測量ステーションによって設定された座標を使用して SpatialPolygonsDataFrame をサブセット化したいのですが、関連する値がオブジェクトのどこに格納されているかわかりません。

シェープファイルがオンラインでホストされていないため、関連データを提供できないことを残念に思います。ただし、オランダ向けのこの投稿に対する Paul Hiemstra の応答からいくつかのコードを借ります。

library(ggplot2)
library(sp)
library(automap)
library(rgdal)

#get the spatial data for the Netherlands
con <- url("http://gadm.org/data/rda/NLD_adm0.RData")
print(load(con))
close(con)

class(gadm)
bbox(gadm)
>         min      max
>r1  3.360782  7.29271
>r2 50.755165 53.55458

調査がこの分野で行われたとしましょう：

bbox(surveys)
    >         min      max
    >r1     4.000    7.000
    >r2    51.000   53.000

SpatialPolygonsDataFrame のその領域を切り取るにはどうすればよいですか?

編集：この質問は私の答えのようです。十分に検索できなかったことをお詫びします (ただし、コメントは、rgeos をどこに向けるべきかについての感覚を与えてくれました)。ただし、gIntersection Rがクラッシュします...

4

1 に答える 1

-1

ポリゴンのサイズに応じて、次のようなことができます

range = cbind(c(4,7), c(51,53))
centroids <- coordinates(spdf)

spdf.subset <- spdf[centroids[,1] > range[1,1] &
                    centroids[,1] < range[2,1] &
                    centroids[,2] > range[1,2] &
                    centroids[,2] < range[2,2],]
于 2012-06-19T22:09:01.037 に答える