Realización de un contorno de volumen del 95% en ArcMap sin herramientas Hawth o GME?

Question

Tengo algunos puntos que representan ubicaciones de animales a lo largo de unos meses. Estoy emulando el trabajo de investigadores anteriores en este campo. Aún soy novato en el uso de SIG, por lo que antes pensé que quería pasar las ubicaciones de mis animales por un Kernel Density Estimate. Sin embargo, releyendo la literatura anterior y comparando los mapas, creo que necesito un tipo particular de KDE. La investigación anterior se publicó en 2011, por lo que el autor utilizó las herramientas de Hawth para realizar una KDE fija con contornos de volumen (referenciada aquí y más abajo). Por desgracia, Hawth no funciona con Arc10 y su análogo, GME, no funciona con ArcMap 10.6, que está en todos los ordenadores de la universidad.

¿Existe alguna otra herramienta que pueda descargar para realizar este contorno de volumen porcentual?

He aquí una foto de la KDE del informe citado

Gestión de los bosques de pinos y enebros en la Base Aérea de Kirtland: Pinyon Jay summer and winter home ranges and habitat use Informe final de 2009 :

Para el análisis kernel utilizamos el Too de Hawth Estimador. El análisis kernel es un método estadístico no paramétrico para estimar densidades de probabilidad a partir de un conjunto de puntos. Cuando se utiliza para para analizar datos sobre el área de campeo, los métodos de densidad kernel describen la probabilidad de encontrar un animal en el área de campeo. probabilidad de encontrar un animal en cualquier lugar. El kernel fijo calcula una densidad produce líneas de contorno que representan el límite del área que contiene un porcentaje especificado del volumen de una densidad de probabilidad de probabilidad. Un contorno de volumen del 95%, por ejemplo, suele contener 95% de los puntos utilizados para generar la estimación de la densidad del kernel. Kernel se establecieron como sigue: sc normal bivariante, factor de suavizado de parámetro único (h)-1000, raster y contornos de volumen porcentual: 50, 90 y 95.

Answer 1

Aunque esta respuesta es probablemente un poco tarde, hay una manera muy simple de superar este problema utilizando R que no requiere conocimientos y puede ser útil para cualquier persona con un problema similar . El código siguiente calculará las curvas de nivel necesarias y luego las escribirá en un shapefile que se puede proyectar en ArcGIS.

library(sp)
library(adehabitatHR)

# read in a .csv which contains lat and long for your data points

mydata<-read.csv("mydata.csv")

# create spatial points object

datapoints<-SpatialPoints(cbind(mydata$lon,mydata$lat))

kud <- kernelUD(datapoints, h=0.2) #kernel density utilization distribution (h=smoothing factor)

homerange75 <- getverticeshr(kud, percent=75) #convert into vector object, use 75,50,25 for quartic kernel

# set coordinates and write to shapfile which can be projected in ArcGIS

WGScoor<-  homerange75
proj4string(WGScoor)<- CRS("+proj=longlat +datum=WGS84")
raster::shapefile(WGScoor, "homerange75_new.shp")

Answer 2

ArcGIS 10.6 es la última versión, y probablemente la última, de ArcGIS Desktop antes de que ESRI vuelque todos sus esfuerzos en ArcGIS Pro. Como tal, las herramientas como herramientas de Hawth y GME probablemente no se actualizará para ArcGIS 10.6.

Así que la solución a mi problema, salvo convertirme en competente en modelización espacial en R esta semana, era encontrar a alguien con ArcMap 10.3 y utilizar GME para crear isopletas a partir de mis anteriores estimaciones de densidad del núcleo. ¡A veces hay que bajar antes de poder subir!