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