В настоящее время я пытаюсь прочитать в каталоге файлов las (по всему штату 10 ТБ) и извлечь около 500 графиков из шейп-файла центральных точек. Шейп-файл правильно читается, и я проверил его в ArcGIS, чтобы быть уверенным. Однако, когда я пытаюсь отобразить каталог плиток las, они кажутся несвязанными и не могут совпасть ни с одной из точек из шейп-файла (который должен находиться в пределах пространственного экстента). Похоже, что север может быть правильным, но восток может быть отключен. Я попытался изменить проекцию файлов las, чтобы она соответствовала шейп-файлу, но это не помогло решить проблему.
Любая помощь вообще была бы замечательной! Спасибо!
Каталог:
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
Шейп-файл:
Объект класса 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
Вот код, который я использовал.
'''shp <- readOGR("D:/All LAS/Shapefile/plot4.shp")
'''ctg <- catalog(folder = "D:/All LAS")
'''plot(ctg)
'''EABplots <- lasclip(ctg, shp, radius = 7.3)