Wie man XY und ID zum Lesen in adehabitatRS eingibt

Question

Ich habe versucht, adehabitat laufen zu lassen, aber nachdem meine Spezies/Präsenz-Koordinaten transformiert wurden, fehlen XY-Komponenten.

CSV-Datei mit species/Präsenz- Koordinaten

> registros <- read.csv("~/Desktop/datos habitat.csv") 
> head(registros,3) 
     Especie   X  Y 
1  C.thous -71.42153 5.98031 
2 H.hydrochaeris -71.42153 5.98031 
3  L.pardalis -71.42153 5.98031 

Projektionskoordinatenform lat-long

> coordinates(registros) <-c("X","Y") 
> proj4string(registros) <- CRS("+init=epsg:4326") 
> registros <- spTransform(registros, CRS("+init=epsg:21898")) 

XY-Koordinaten fehlen gegangen, nur Spezies Variable, wenn links.

Also wenn ich versuche, eine Liste zu erstellen, um adehabitat nur SP laufen zu lassen. Variable ist in den Loks vorhanden.

$map 
class  : SpatialPixelsDataFrame 
dimensions : 2614, 3084, 8061576, 2 (nrow, ncol, npixels, nlayers) 
resolution : 9.259e-05, 9.259e-05 (x, y) 
extent  : -71.59446, -71.30891, 5.876562, 6.118592 (xmin, xmax, ymin, ymax) 
coord. ref. : +init=epsg:21898 +proj=tmerc +lat_0=4.599047222222222 +lon_0=-71.08091666666667 +k=1 +x_0=1000000 +y_0=1000000 +ellps=intl +towgs84=307,304,-318,0,0,0,0 +units=m +no_defs 
names  :   v_eco.agua, v_eco.bosque 
min values : -2.60905458751601e-05,   0 
max values : 0.0236267629555045,   376 

$relocs 
class  : SpatialPointsDataFrame 
features : 270 
extent  : 949058.7, 966743.8, 1147074, 1165429 (xmin, xmax, ymin, ymax) 
coord. ref. : +init=epsg:21898 +proj=tmerc +lat_0=4.599047222222222 +lon_0=-71.08091666666667 +k=1 +x_0=1000000 +y_0=1000000 +ellps=intl +towgs84=307,304,-318,0,0,0,0 +units=m +no_defs 
variables : 1 
names  :  Especie 
min values :   C.paca 
max values : T.tetradactyla 

Irgendwelche Vorschläge, wie kann ich dieses Problem lösen? Danke!

Answer 1

Sie sind nicht weg, es ist nur, weil die head Methode kommt frmo raster und druckt sie nicht. Wenn Ihre R-Sitzung sp lädt, aber raster lädt, zeigt der Befehl head (von sp) Ihre Koordinaten an.