SpatialPixelsDataFrame オブジェクトに一連の集水域ポリゴンを (シェープファイルから) オーバーレイするときに奇妙な問題が発生します。70 の流域がありますがNA
、60 の流域で s を取得しています。
データは Dropbox で入手できます。次のスクリプトは、ファイルのダウンロード、解凍などを行い、再現可能な例を作成します。
## load packages required
require(rgdal)
require(sp)
## proj4 string
BNG <- CRS("+init=epsg:27700")
## load in the grid
download.file("https://www.dropbox.com/s/7dvcf6ddo6r3043/depogrid.rds?dl=1",
"./depogrid.rds", method = "wget")
myGrid <- readRDS("depogrid.rds")
## convert to spatial pixels
depo.sp <- SpatialPixelsDataFrame(data.matrix(myGrid[, 1:2]),
data = myGrid[, 3, drop = FALSE],
proj4string = BNG)
## load the shapefile
download.file("https://www.dropbox.com/s/bwzmeumnrja3c0o/dart.zip?dl=1",
"./dart.zip", method = "wget")
unzip("./dart.zip")
shp <- readOGR(".", "DARTall")
## project shp to OSGB36
shp <- spTransform(shp, BNG)
## overlay shapefile on depo grid
depo.dart <- over(shp, depo.sp)
問題はNA
次のとおりです。
> head(depo.dart)
Deposition
0 NA
1 NA
2 NA
3 NA
4 NA
5 0.2604089
> sum(is.na(depo.dart[,1]))
[1] 60
> nrow(depo.dart)
[1] 70
グリッドが英国地域全体をカバーし、集水域が小さいため、 sはないはずですNA
。おそらく、何か間違ったことをしているでしょう! もしそうなら、私は何を間違っていますか?どうすれば修正できますか? の 70 の集水域ごとにDeposition
(データ)の平均値を求めたいと思います。SpatialPixelsDataFrame
shp
[グリッドはダミー データであることに注意してください。実際の堆積グリッドを配布することはできませんが、データとオブジェクトは基本的に (正確に) と同じdepo.sp
です。
追加情報
> over(shp, depo.sp, returnList = TRUE)
[[1]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[2]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[3]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[4]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[5]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[6]]
Deposition
3143 0.2604089
[[7]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[8]]
Deposition
3289 0.1139663
[[9]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[10]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[11]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[12]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[13]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[14]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[15]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[16]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[17]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[18]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[19]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[20]]
Deposition
3585 0.8610678
[[21]]
Deposition
3288 0.2242067
[[22]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[23]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[24]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[25]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[26]]
Deposition
3141 0.3296081
[[27]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[28]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[29]]
Deposition
3438 0.3609788
[[30]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[31]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[32]]
Deposition
3287 0.4812801
[[33]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[34]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[35]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[36]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[37]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[38]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[39]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[40]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[41]]
Deposition
3436 0.3223757
[[42]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[43]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[44]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[45]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[46]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[47]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[48]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[49]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[50]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[51]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[52]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[53]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[54]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[55]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[56]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[57]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[58]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[59]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[60]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[61]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[62]]
Deposition
2580 0.2820617
[[63]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[64]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[65]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[66]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[67]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[68]]
Deposition
3141 0.3296081
3287 0.4812801
3288 0.2242067
3438 0.3609788
[[69]]
[1] Deposition
<0 rows> (or 0-length row.names)
[[70]]
[1] Deposition
<0 rows> (or 0-length row.names)
> str(over(shp, depo.sp, returnList = TRUE))
List of 70
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 1 obs. of 1 variable:
..$ Deposition: num 0.26
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 1 obs. of 1 variable:
..$ Deposition: num 0.114
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 1 obs. of 1 variable:
..$ Deposition: num 0.861
$ :'data.frame': 1 obs. of 1 variable:
..$ Deposition: num 0.224
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 1 obs. of 1 variable:
..$ Deposition: num 0.33
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 1 obs. of 1 variable:
..$ Deposition: num 0.361
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 1 obs. of 1 variable:
..$ Deposition: num 0.481
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 1 obs. of 1 variable:
..$ Deposition: num 0.322
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 1 obs. of 1 variable:
..$ Deposition: num 0.282
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 4 obs. of 1 variable:
..$ Deposition: num [1:4] 0.33 0.481 0.224 0.361
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)
$ :'data.frame': 0 obs. of 1 variable:
..$ Deposition: num(0)