Buenos días, soy nuevo en esta página y no se si estoy realizando en el sitio adecuado la consulta. La expongo de todos modos.
Estoy realizando un código para mi trabajo de fin de master, en una parte del código busco obtener la normalización de los DSM de algunos ayuntamientos de Galicia. Los archivos LiDAR empleados proceden del CNIG, son los datos de la 2ª cobertura (2015-actualidad), los datos de los DTM proceden también del CNIG y son los datos con un paso de maya de 2x2. La información LiDAR por ayuntamientos es muy pesada y he tenido que llevar a cabo la lectura de los archivos .las por partes.
Después de la lectura de los archivos calculo el DSM con la siguiente función:
DSM_estrada_1 <- lidR::grid_canopy(LiDARfiles_estrada_1, res = 2, algorithm = lidR::p2r(subcircle = 1, na.fill = tin()))
Este primer paso funciona perfectamente en cualquiera de los datos seleccionados, el problema surge al normalizar.
Al realizar la función para alguno de los ayuntamientos seleccionados me sale el siguiente error:
LiDARNorm <- lidR::normalize_height(LiDARfiles_estrada_1, InpDTM, na.rm = TRUE)
Error: On-disk rasters not supported in 'raster_as_dataframe()'
In addition: Warning message:
In raster_cell_from_xy(raster, x, y) :
NAs introduced by coercion to integer range
EL tema es que en algunos ayuntamientos si me funciona y en otros me sale ese error. No consigo encontrar por ningún lado una posible solución, no se si alguien podría ayudarme?
Un saludo y muchas gracias de antemano