Início / coding / Perguntas / 79557649
Accepted
haddonchris
Asked: 2025-04-06 07:39:49 +0800 CST

Por que meu MCP está me dando um tamanho de área muito pequeno?

  • 772

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!

1 respostas

  1. Best Answer

    Como @r2evans mencionou, você precisa projetar seu objeto sp para o CRS da Zona UTM correspondente antes de executar mcp(). Atualmente, mcp()está tratando suas coordenadas como metros literais, daí os valores de área pequenos.

    Veja como corrigir seu problema. Observe que você tem "sorte" de todos os seus dados estarem dentro da mesma Zona UTM. Se seus dados abrangerem várias Zonas UTM, você precisará pesquisar online como lidar com esse problema/escolher o CRS correto em relação à extensão dos seus dados.

    library(sp)
library(adehabitatHR)

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))

# Function to return UTM Zone/s
utm_epsg <- \(lon, lat) {

  utm_zone <- floor((lon + 180) / 6) + 1
  if(lat >= 0) {
    epsg <- 32600 + utm_zone
  } else {
    epsg <- 32700 + utm_zone
  }
  
}

# Get UTM Zone/s for animal_sp
unique(mapply(utm_epsg, animal_sp$Long, animal_sp$Lat))
# [1] 32736

# Define coordinates
coordinates(animal_sp) <- c("Long", "Lat")

# Set intial projection
proj4string(animal_sp) <- CRS("EPSG:4326")

# Project to relevant UTM Zone CRS
animal_sp_utm <- spTransform(animal_sp, "EPSG:32736")

# Create the MCP
animal_mcp <- mcp(animal_sp_utm, 
                  percent = 100,
                  unout = "m2") # Make the MCP

animal_mcp
# Object of class "SpatialPolygonsDataFrame" (package sp):
# 
# Number of SpatialPolygons:  2
# 
# Variables measured:
#    id      area
# S1 S1 601430654
# S2 S2 461876993
    • 2

relate perguntas