Classify_poi : Only spatial polygons are supported

[MaximeAuffret]

Hello,
I use lidR to produce DEM from lidar data.
I recently updated all my packages and now the function classify_poi, that I use to classify the buildings, doesn't work anymore. here's what I get :

Usually it worked without any problem.
I tried to change the chunk size, the output DEM resolution, the sf/sfc or sp format of the building file.
Do you know where it come from ?

Usually it worked without any problem.
I tried to change the chunk size, the output DEM resolution, the sf/sfc or sp format of the building file.
Do you know where it come from ?

# The library that I use
library(lidR)
library(knitr)
library(shiny)
library(manipulateWidget)
library(raster)
library(ggplot2)
library(sf)

# The script
BAT <- st_read(dir_bat, quiet = T)    # Building file (.shp)
ZI_LAS = readLascatalog(x)              # Las file from IGN (Test urban area of Manosque (04) in France) the file is declassified for tests

mycfs <- csf(
    sloop_smooth      <- FALSE,
    class_threshold   <- 0.05,
    cloth_resolution  <- 0.5,
    rigidness         <- 1,
    iterations        <- 500,
    time_step         <- 0.65
  ) 

opt_laz_compression(ZI_LAS) = T
opt_output_files(ZI_LAS) <- paste0(dir_sol,"/Sol_{XCENTER}_{YCENTER}_{ID}")

ZI_SOL = classify_ground(ZI_LAS,algorithm = mycfs)
  
opt_chunk_buffer(ZI_SOL) <- 20
opt_chunk_size(ZI_SOL)   <- 250
plot(ZI_SOL, chunk = TRUE)

dir.create(paste0(dir_sol,"/BAT"))
dir_bat = paste0(dir_sol,"/BAT")

opt_output_files(ZI_SOL) <- paste0(dir_bat,"/Bat_{XCENTER}_{YCENTER}_{ID}")

st_crs(BAT) <- st_crs(ZI_SOL)
ZI_BAT <- classify_poi(ZI_SOL, class = as.integer(6), roi = BAT )

UPDATE : I use my scipt on two computer, it work in one of them so I think the problem comes from a package but i don't find wich one.

[Jean-Romain]

Please provide a minimal reproducible example in the form of text copy pasted (no image).

[Jean-Romain]

I guess a minimal reproducible example would be

library(lidR)
library(sf)

x = "..."

BAT <- st_read(dir_bat, quiet = T)    # Building file (.shp)
ZI_SOL = readLAScatalog(x)            # Las file from IGN (Test urban area of Manosque (04) in France) the file is declassified for tests

opt_output_files(ZI_SOL) <- paste0(tmpdir(),"/Bat_{XCENTER}_{YCENTER}_{ID}")

st_crs(BAT) <- st_crs(ZI_SOL)
ZI_BAT <- classify_poi(ZI_SOL, class = as.integer(6), roi = BAT )

Now I need some data to reproduce.

[MaximeAuffret]

You can download the data here : https://wxs.ign.fr/mct8y415ma63eudk50qaf2bi/telechargement/prepackage/JEU-TEST_LIDARHD_PACK_17-05-2022$LIDARHD_1-0_HYB_JEU-TEST-DALLE-0923_6308-2021/file/LIDARHD_1-0_HYB_JEU-TEST-DALLE-0923_6308-2021.7z

The building file is attched here.
BAT.zip

I removed the classification by using the extent (raster package) of the data and classify_poi (before the update).

[Jean-Romain]

The following worked fine. Again please provide a minimal reproducible code. Do not let met guess what code I should run to reproduce.

library(lidR)
library(sf)

x = "lidR issues/issue 630/LIDARHD_1-0_HYB_JEU-TEST-DALLE-0923_6308-2021/"
dir_bat = "lidR issues/issue 630/BAT"

BAT <- st_read(dir_bat, quiet = T)    # Building file (.shp)
ZI_SOL = readLAScatalog(x)            # Las file from IGN (Test urban area of Manosque (04) in France) the file is declassified for tests

opt_chunk_buffer(ZI_SOL) <- 20
opt_chunk_size(ZI_SOL)   <- 250
opt_laz_compression(ZI_LAS) = T
opt_output_files(ZI_SOL) <- paste0(tempdir(),"/Sol_{XCENTER}_{YCENTER}_{ID}")

plot(ZI_SOL, chunk = TRUE)
plot(BAT, add = T)

st_crs(BAT) <- st_crs(ZI_SOL)
ZI_BAT <- classify_poi(ZI_SOL, class = as.integer(6), roi = BAT )

[Jean-Romain]

Almost mimaly reproducible but I can't copy paste it to reproduce and I need to make modifications to make it working. The following seems minimal but not reproducible because it works for me. Can you confirm if it works for you?

I'm using 4.0.3 but I don't remember I changed anything in classify_poi() since 4.0.0. There is nothing in the changelog.

library(lidR)
library(sf)

LAS = readLAScatalog("issue 630/DONNEES_NON_CLASSIFIES/")
bat = st_read("issue 630/BAT/BAT.shp")

opt_output_files(LAS) = paste0(tempdir(),"/*")

st_crs(bat) <- st_crs(LAS)
classify_poi(LAS, class = as.integer(6), roi = bat)

[MaximeAuffret]

It is not working.
I am re-installing Rtools 4 and updating all the packages again to see.

[Jean-Romain]

You can only install 4.0.2 from CRAN

[MaximeAuffret]

It seems to be working now. I don't know what happened but it was certainly the packages.
Maybe the version of Rtools ?

Now I've got :
sf= 1.0.8 (1.07 before)
lidR = 4.0.1 (4.0.0 before)

Thank you for your help. And I apologize for the inconvenience.

[Jean-Romain]

Ok, I don't know either. I'm glad it works