Hola Andrés, El argumento x dice que debe ser una conductanceMatrix no un raster.
Manuel *Manuel Spínola, Ph.D.* Instituto Internacional en Conservación y Manejo de Vida Silvestre Universidad Nacional Apartado 1350-3000 Heredia COSTA RICA [email protected] <[email protected]> [email protected] Teléfono: (506) 8706 - 4662 Sitio web institucional: ICOMVIS <http://www.icomvis.una.ac.cr/index.php/manuel> Sitio web personal: Sitio personal <https://mspinola-sitioweb.netlify.app> Blog sobre Ciencia de Datos: Blog de Ciencia de Datos <https://mspinola-ciencia-de-datos.netlify.app> El El lun, 14 oct 2024 a la(s) 09:20, Andrés Ocaña Carretón via R-help-es < [email protected]> escribió: > ¡Buenas tardes! > > Antes de nada, dar las gracias por permitirme participar en esta lista. > > Soy un usuario recién llegado a R de la mano del análisis geoespacial. Soy > un arqueólogo que se ha topado con un package para R creado por un > arqueólogo para la ejecución de procesos para a creación de rutas de menor > coste cuyo nombre es leastcostpath y que se encuentra en el repositorio de > GitHub. Contiene diferentes scripts relacionados con esta cuestión. > > En concreto, estoy intentando realizar un cálculo de tipo FETE (from > everywhere to everywhere approach). Para ello he utilizado una capa de tipo > ráster, en este caso de costes, y otra de tipo vectorial con puntos. Ambas > las he creado en QGIS. La de costes es un mapa de pendientes en grados > reclasificado. He usado una capa de costes por ser la que se utiliza en > algunos ejemplos de uso que he visto, si bien en ellos se utiliza otra > generada a partir de un mdt con la aplicación de algún algoritmo, por el > ejemplo el de Tobler. Esta se puede crear con otro script incluido en el > paquete, pero que no he sido capaz de ejecutar porque me da error. > > Tras cargar todos los packages necesarios en RStudio, ejecuto el > siguiente código: > > r <- raster("pendiente_grados_reclas.tif") > locs <- st_read("puntos.geojson") > > lcps <- create_FETE_lcps(x = r, locations = locs, cost_distance = FALSE, > ncores = 1) > > Y me da el siguiente mensaje de error: > "Error in .local(x, ...) invalid layer names" > > El código usado es igual al que viene en la ayuda del script y los nombres > de las capas son iguales a las creadas, por lo que no comprendo el error. > He estado buscando por internet pero no encuentro nada parecido. > > ¿Alquien me podría ayudar a solucionarlo ? ¿O indicarme dónde puede estar > el problema? > > Incluyo el texto de la ayuda en cuanto al uso del script por si resulta de > interés. > > Usage > create_FETE_lcps(x, locations, cost_distance = FALSE, ncores = 1) > > Arguments > x conductanceMatrix > > locations sf 'POINT' or 'MULTIPOINT', SpatVector, data.frame or matrix > containing the locations coordinates > > cost_distance logical if TRUE computes total accumulated cost from origin > to destination. FALSE (default) > > ncores numeric Number of cores used when calculating least-cost paths > from-everywhere-to-everywhere. 1 (default) > > Un saludo > > Andrés Ocaña Carretón > > > > > > [[alternative HTML version deleted]] > > _______________________________________________ > R-help-es mailing list > [email protected] > https://stat.ethz.ch/mailman/listinfo/r-help-es > [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
