Удаление третье измерение в SpatialPointsDataFrame

Question

У меня есть SpatialPoints объект с 3-х измерениях:

x <- c(1,1,1,2,2,2,3,3,3) 
y <- c(1,2,3,1,2,3,1,2,3) 
z <- c(1,3,5,2,1,2,1,2,3) 
xyz <- cbind(x,y,z) 
ss <- SpatialPoints(xyz) 
dimensions(ss) 

И в raster объекта:

rr <- raster(matrix(round(runif(49,0,10)),7,7), xmn=0, xmx=4, ymn=0, ymx=4, crs=NA, template=NULL) 

Я хочу, чтобы извлечь значения растра с использованием SpatialPoints объекта:

extract(rr,ss) 
#Error in .xyValues(x, coordinates(y), ..., df = df) : 
# xy should have 2 columns only. 
#Found these dimensions: 9, 3 

Вы можете визуализировать данные, если вы нт:

plot(rr) 
plot(ss, add=T) 

Таким образом, проблема в том, что функция extract пакета растрового требует 2 объекта измерения SpatialPoints. Моя (по моим реальным данным) в 3-х мерных. Есть ли способ отказаться от 3-го измерения моей формы фигуры? Я пробовал:

coordinates(ss) <- coordinates(ss)[,-3] 
#Error in `coordinates<-`(`*tmp*`, value = c(1, 1, 1, 2, 2, 2, 3, 3, 3, : 
# setting coordinates cannot be done on Spatial objects, where they have #already been set 

Я не хочу, чтобы перестроить форму с нуля.

Answer 1

Просто перезаписать coords слот объекта S4:

ss@coords <- ss@coords[, 1:2] 

---

Я не знаю, как создается ваш объект SpatialPoints, но если вы используете rgdal::readOGR есть pointDropZ аргумент (по умолчанию FALSE)

Answer 2

@rcs ответ лучше, но это работает, а также:

ss <- SpatialPoints(coordinates(ss)[,-3]) 

, и если у вас есть SpatialPointsDataFrame:

ss <- SpatialPointsDataFrame(coordinates(ss)[,-3], ss@data,proj4string=CRS(proj4string(ss)))