"lidR" package - Error in projection of las files as catalog

[JessieElliott]

I am currently trying to read in a catalog of las files (statewide 10TB) and extract around 500 plots from a shapefile of central points. The shapefile reads in correctly, and I have checked it in ArcGIS to be sure. However, when I attempt to display the catalog of las tiles, they appear disjointed and fail to line up with any of the points from the shapefile (which should fall well within the spatial extent). It appears that the northing may be correct, but the easting could be off. I attempted to change the projection of the las files to match the shapefile, but this did not remedy the issue.

Any help at all would be wonderful! Thank you!

Catalog: class : LAScatalog extent : 70000 , 3195000 , 940000 , 2380000 (xmin, xmax, ymin, ymax) coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 area : 997518.4 kunits² points : 379.02 billion points density : 0.4 points/units² num. files : 41353

Shapefile: Object of class SpatialPointsDataFrame Coordinates: min max coords.x1 -87.97372 -84.84951 coords.x2 37.80275 41.72992 Is projected: FALSE proj4string : [+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0] Number of points: 530 Data attributes: PLOT LAT LON
1005 : 1 Min. :37.80 Min. :-87.97
1024 : 1 1st Qu.:38.25 1st Qu.:-86.71
1031 : 1 Median :38.70 Median :-86.34
1035 : 1 Mean :38.82 Mean :-86.39
1042 : 1 3rd Qu.:39.22 3rd Qu.:-86.21
1046 : 1 Max. :41.73 Max. :-84.85
(Other):524

[technocrat]

Hi, and welcome!

Please see the FAQ: What's a reproducible example (`reprex`) and how do I do one? Using a reprex, complete with representative data will attract quicker and more answers. This isn't yet a coding question so it doesn't apply.

Still, it's easy to assume when asking a question the possibility that those reading it might not actually understand the world of ArcGIS and what las files are and what R tools are being used to fetch the catalog object.

So, I can only offer a guess that these are LIDAR files and that what's needed might be

library(rLiDAR)


# Importing LAS file:
LASfile <- system.file("extdata", "LASexample1.las", package="rLiDAR")
# Reading LAS file 
LAS<-readLAS(LASfile,short=TRUE)
head(LAS)
#>             X       Y    Z Intensity ReturnNumber
#> [1,] 525924.4 3377909 0.05       144            1
#> [2,] 525924.0 3377910 0.00       120            1
#> [3,] 525923.6 3377910 0.07       121            1
#> [4,] 525923.2 3377911 0.01       147            1
#> [5,] 525922.8 3377911 0.03       136            1
#> [6,] 525922.4 3377912 0.00       137            1

^{Created on 2020-04-08 by the reprex package (v0.3.0)}