Tengo un gran conjunto de datos con 36.000 puntos que representan usos comerciales del suelo, cada uno con un campo que contiene los metros cuadrados. He ejecutado un análisis de densidad de núcleo en este conjunto de datos, produciendo una trama que muestra la densidad de metros cuadrados comerciales en toda el área metropolitana. Tengo que dividir esta trama en regiones correspondientes a los máximos locales, a los que llamo "centro". Ya he determinado la ubicación de los centros, y ahora tengo que hacer una de estas dos cosas:
-
utilizar una herramienta de agrupación de puntos, como "partitioning around medoids", para agrupar los puntos en clusters alrededor de los centros que he identificado. El problema de este método es que es computacionalmente intenso, y más aún si intento utilizar una matriz de disimilitud para ponderar los puntos por tamaño.
-
dividir de algún modo el raster de densidad del núcleo (que se asemeja aproximadamente a un raster de terreno) en "colinas" individuales alrededor de cada centro. Pero no se me ocurre ninguna herramienta para hacerlo.
Este problema me atormenta desde hace tiempo, y esperaba poder realizar el método de agrupación en R, pero lleva mucho tiempo y me estoy quedando sin tiempo. ¿Alguien conoce un método sencillo para dividir los rásters de densidad en vecindarios de intensidad o para agrupar rápidamente grandes conjuntos de datos?