I'm trying to reproject a LAZ file from EPSG 2226 to EPSG 6418.
The reprojection works fine (I followed this: Can I re-project an LAS file in LidR) but I cannot store the information onto the LAS header and subsequently, as well as I save my LAS file (and get out of R environment), the CRS is lost.
I used the code below:
library(lidR) las_tmp <- readLAS("./USGS_LPC_CA_Sonoma_2013_SOCO_0054_108_LAS_2017.laz") las_tmp_reproj <- spTransform(las_tmp, sp::CRS("+init=epsg:6418")) #> Warning messages: #> 1: In showSRID(uprojargs, format = "PROJ", multiline = "NO") : #> Discarded datum NAD83_National_Spatial_Reference_System_2011 in CRS definition #> 2: In sp::proj4string(x) : CRS object has comment, which is lost in output #> 3: EPSG code not found: header not updated. Try to use the function epsg() manually to ensure CRS will be written in file. #> 4: In showSRID(uprojargs, format = "PROJ", multiline = "NO") : #> Discarded datum Unknown based on GRS80 ellipsoid in CRS definition
As said by
warnings() number 3, the epsg code is not found and header not updated. I followed the tips and tried the following command:
epsg(las_tmp_reproj) <- as.integer(6418) #> Warning message: #> In showSRID(SRS_string, format = "PROJ", multiline = "NO") : #> Discarded datum NAD83_National_Spatial_Reference_System_2011 in CRS definition
Unfortunately, after this step, the access to the epsg information
epsg(las_tmp_reproj) throws me the following error:
epsg(ctg_tmp_reproj) #> Error in if (tags[[i]]$key == 3072) return(i) : #> l'argument est de longueur nulle
Is it related with the datum itself? The error message is not really explicit and I'm pretty stuck