Esta é a minha primeira vez fazendo estatística espacial e estou aprendendo aos poucos. Já vi algumas outras pessoas aqui com problemas semelhantes, mas nenhuma das respostas pareceu funcionar para mim.
Estou tentando fazer um MCP, mas quando vejo os resultados, o tamanho da área é bem menor do que deveria. No entanto, quando ploto os polígonos sobre um mapa de satélite, parece correto.
Meus dados estão configurados assim:
set.seed(1)
animal_sp <- data.frame(ID = rep(c("S1", "S2"), each = 10),
Long = runif(20, 34.2002, 34.55328),
Lat = runif(20, -19.11998, -18.82139))
E o código que usei é o seguinte:
library(sp)
sp::coordinates(animal_sp) <- c("Long", "Lat") # define the variables
proj4string(animal_sp) <- CRS("+init=epsg:4326") # Set the projection
# Create the MCP using the adehabitatHR package:
library(adehabitatHR)
animal_mcp <- mcp(animal_sp,
percent = 100,
unout = "m2") # Make the MCP
animal_mcp # look at the results
Os resultados que obtenho são:
Object of class "SpatialPolygonsDataFrame" (package sp):
Number of SpatialPolygons: 2
Variables measured:
id area
S1 S1 0.02919310
S2 S2 0.03671279
Eu deveria obter uma área de ~260 km².
Qualquer ajuda será muito apreciada!