Más

Salida de isoplet inesperada de GME

Salida de isoplet inesperada de GME


Convertí un archivo csv en SpatialPointsDataFrame proyectado en R, lo ingresé en bkde2D, convertí el resultado en ráster, guardé el ráster en un TIFF e ingresé ese ráster en la herramienta GME Isopleth.

Aquí está mi código R:

data = read.csv ("_ PolarBearData.csv", na.strings = "NA", header = T) oldproj = "+ proj = longlat + ellps = WGS84 + datum = WGS84" myproj = "+ proj = stere + lat_0 = 90 + lon_0 = -55 + k = 1 + x_0 = 0 + y_0 = 0 + ellps = WGS84 + datum = WGS84 + unidades = m + no_defs + towgs84 = 0,0,0 "coords = cbind (datos $ longitud, datos $ latitude) spdata = SpatialPointsDataFrame (coords, data, proj4string = CRS (oldproj)) # asegúrese de que se proyecte spdata = spTransform (spdata, CRS (myproj)) # parámetros de KDE xy = cbind (coordenadas (spdata) [, 1], coordenadas (spdata [, 2]) bandwidth = 34390 grid_x = 130 grid_y = 288 # Ejecute KDE y guárdelo en ráster TIFF datakernel = bkde2D (xy, bandwidth, gridsize = c (grid_x, grid_y)) ras = raster (list (x = datakernel) $ x1, y = datakernel $ x2, z = datakernel $ fhat)) writeraster (ras, paste ("mydata", '.tif', sep = ")," GTiff ", sobreescritura = T)

A continuación, ejecuté este código GME para producir contornos de volumen del 95%:

isopleth (in = "I:  Data  mydata.tif", out = "I:  Data  Contours.gdb! mylines", cuantiles = c (.95), poly = "I:  Data  Contours.gdb! mypoly ");

Aquí hay una imagen de mi salida: Bad GME Isopleth

La trama de blanco a negro (fondo) representa un KDE generado en R con "bkde2D" del paquete KernSmooth. Los puntos amarillos fueron la entrada a KDE. La línea roja muestra la isopleta del 95% generada en GME. GME está estirando los contornos hacia el sur cuando la entrada es un TIFF sin alterar de R, pero de lo contrario, si el ráster es de Arc, por ejemplo, GME produce contornos correctos. ¿Por qué pasó esto?

ACTUALIZACIÓN: También debo mencionar que agregué el ráster de R directamente en ArcMap, y tenía un CRS desconocido a pesar de que lo especifiqué en R (al menos, para la entrada). Intenté usar Definir proyección para solucionar este problema, pero GME aún produce el mismo resultado. El ráster en sí tiene coordenadas que tienen sentido (coincide con los grupos de puntos en la imagen), pero GME no produce buenas isopletas a partir de él.


No tengo idea de dónde puede estar fallando GME. El software es opaco, sin código fuente para examinar. Como tal, no hay forma de saber qué método se ha implementado o cuál puede ser la fuente del problema.

A continuación, se ofrece una solución para derivar el volumen porcentual en R. Puede ser útil comparar los resultados para asegurarse de que no haya problemas con sus datos.


Ver el vídeo: भगलक चतन पर डरवन क परभव,Impact of Darvin on Geographical Thought For NET,PGT,GEOGRAPHY