2016-07-24 1 views
0

J'exécute une analyse MCP (domaine vital) en R sur mes données (appelé Ade) qui est un ensemble de points GPS pour les emplacements d'animaux.Le message d'erreur xy devrait être de classe SpatialPoints

Je courais le script suivant:

Ade <- read.csv("Ade.csv") 
#Get the folder SP that is used to convert the data 
library(sp) 
#Convert 
Ade <- matrix(runif(1564), ncol=2) 
head (Ade) 
Ade <- data.frame(Ade) 
#Convert Ade to spatial points data frame by first ... 
#Tell R that the coordinates are in the column X and column Y 
coords <- SpatialPoints(Ade[, c("X", "Y")]) 
AdeSpat <- SpatialPointsDataFrame(coords, Ade) 
#Specify the coordinate reference system of the data 
proj4string(AdeSpat) <- CRS("+proj=XY +ellps=WGS84") 
#View data 
getClass(AdeSpat) 
class(AdeSpat) 
head (AdeSpat) 
#Two steps in loading AdehabitatHR into R 
install.packages("adehabitatHR") 
library(adehabitatHR) 
#Run MCP 
mcp <- mcp(AdeSpat$coordinates, percent=95, unin = c ("km"), unout = c("km2")) 

mais quand je lance ce que je reçois le message d'erreur suivant:

Error in mcp(AdeSpat$coordinates, percent = 95, unin = c("km"), unout = 
c("km2")) : xy should be of class SpatialPoints 

Quelqu'un peut-il me expliquer ce que cela signifie, et comment je peux résoudre le problème, j'ai cherché sur internet mais je ne trouve aucune réponse pour l'erreur "xy devrait être de classe SpatialPoints".

Quand je vérifie la trame de données, je reçois la sortie suivante: vous

> class(AdeSpat) 
[1] "SpatialPointsDataFrame" 
attr(,"package") 
[1] "sp" 

> head (AdeSpat) 
    coordinates ID ZONE  X  Y 
1 (370984, 9752290) 1 36 M 370984 9752293 
2 (370980, 9752290) 2 36 M 370980 9752288 
3 (370980, 9752290) 3 36 M 370980 9752288 
4 (370995, 9752290) 4 36 M 370995 9752287 
5 (370995, 9752290) 5 36 M 370995 9752287 
6 (371032, 9752310) 6 36 M 371032 9752312 
Coordinate Reference System (CRS) arguments: +proj=XY +ellps=WGS84 

Merci J

Répondre

0

Vous devez mettre directement l'objet qui contient vos données spatiales points.

En vous cas, votre formule sera:

mcp < - mcp (coords, pour cent = 95, UNIN = c ("km"), unout = c ("km2"))

J'espère que ça t'aidera!

Murielle