psy*_*sky 5 r geospatial lidar las lidr
激光雷达数据只是 3D 坐标，通常采用las文件格式。\xd0\xa1内容示例
library(rgdal)\nlibrary(raster)\nlibrary(tmap)\nlibrary(tmaptools)\nlibrary(lidR)\nlibrary(RStoolbox)\nlas_cat <- readLAScatalog("C:/1/078-638.las")\nsummary(las_cat)\nopt_chunk_size(las_cat) <- 500\nplot(las_cat, chunk_pattern = TRUE)\nlas_cat\n#> class : LAScatalog (v1.2 format 1)\n#> extent : 637999, 638240.5, 6077999, 6079999 (xmin, xmax, ymin, ymax)\n#> coord. ref. : NA \n#> area : 483081.1 units\xc2\xb2\n#> points : 3.68 million points\n#> density : 7.6 points/units\xc2\xb2\n#> density : 5.6 pulses/units\xc2\xb2\n#> num. files : 1\n有没有办法获得data.frameR 的典型点坐标？\n作为示例，我们可以使用来自此http://data.wvgis.wvu.edu/elevation/的数据。\n此外，有没有办法获得角度反射，data.frame激光雷达文件中的反射数量？
请查找以下一种可能的解决方案以获取data.table,data.frame包含所有信息的解决方案。你可以用来as.data.frame()得到一个 puredata.frame但 adata.table 是adata.frame
雷普莱克斯
\n注意：我使用了库.las中内置的数据集lidR，因为它更方便。
lidRlibrary(lidR)\n\nLASfile <- system.file("extdata", "example.laz", package="rlas")\nlas <- readLAS(LASfile)\n\nsummary(las)\n#> class : LAS (v1.0 format 1)\n#> memory : 21.2 Kb \n#> extent : 339002.9, 339015.1, 5248000, 5248001 (xmin, xmax, ymin, ymax)\n#> coord. ref. : NAD83 / UTM zone 17N \n#> area : 16 m\xc2\xb2\n#> points : 30 points\n#> density : 1.88 points/m\xc2\xb2\n#> density : 1.62 pulses/m\xc2\xb2\n#> File signature: LASF \n#> File source ID: 0 \n#> Global encoding:\n#> - GPS Time Type: GPS Week Time \n#> - Synthetic Return Numbers: no \n#> - Well Know Text: CRS is GeoTIFF \n#> - Aggregate Model: false \n#> Project ID - GUID: 00000000-0000-0000-0000-000000000000 \n#> Version: 1.0\n#> System identifier: LAStools (c) by rapidlasso GmbH \n#> Generating software: las2las (version 201011) \n#> File creation d/y: 343/2011\n#> header size: 227 \n#> Offset to point data: 323 \n#> Num. var. length record: 1 \n#> Point data format: 1 \n#> Point data record length: 28 \n#> Num. of point records: 30 \n#> Num. of points by return: 26 4 0 0 0 \n#> Scale factor X Y Z: 0.001 0.001 0.001 \n#> Offset X Y Z: 6e+05 6500000 0 \n#> min X Y Z: 339002.9 5248000 973.145 \n#> max X Y Z: 339015.1 5248001 978.345 \n#> Variable Length Records (VLR):\n#> Variable Length Record 1 of 1 \n#> Description: by LAStools of rapidlasso GmbH \n#> Tags:\n#> Key 1024 value 1 \n#> Key 3072 value 26917 \n#> Key 3076 value 9001 \n#> Key 4099 value 9001 \n#> Extended Variable Length Records (EVLR): void\npayload(las)\n#> Z gpstime Intensity ReturnNumber NumberOfReturns ScanDirectionFlag\n#> 1 975.589 269347.3 82 1 1 1\n#> 2 974.778 269347.3 54 1 1 1\n#> 3 974.471 269347.3 27 2 2 1\n#> 4 974.025 269347.3 55 2 2 1\n#> 5 974.298 269347.3 117 1 1 0\n#> 6 974.985 269347.3 81 1 1 0\n#> 7 975.182 269347.3 84 1 1 1\n#> 8 974.434 269347.3 104 1 1 1\n#> 9 974.159 269347.3 91 1 1 1\n#> 10 973.145 269347.3 99 1 1 1\n#> 11 976.739 269347.5 87 1 1 1\n#> 12 976.823 269347.5 83 1 1 1\n#> 13 977.227 269347.5 87 1 1 1\n#> 14 975.873 269347.5 87 1 1 1\n#> 15 975.782 269347.5 107 1 1 1\n#> 16 975.353 269347.5 76 1 1 1\n#> 17 974.704 269347.5 113 1 1 1\n#> 18 977.170 269347.5 64 1 1 0\n#> 19 977.757 269347.5 89 1 1 0\n#> 20 978.212 269347.5 98 1 1 0\n#> 21 978.309 269347.5 50 1 2 0\n#> 22 974.816 269347.5 31 2 2 0\n#> 23 978.345 269347.5 51 1 2 1\n#> 24 974.824 269347.5 32 2 2 1\n#> 25 978.014 269347.5 85 1 1 1\n#> 26 977.781 269347.5 94 1 1 1\n#> 27 976.455 269347.5 78 1 1 1\n#> 28 976.313 269347.7 71 1 1 1\n#> 29 975.735 269347.7 75 1 1 1\n#> 30 975.674 269347.7 106 1 1 1\n#> EdgeOfFlightline Classification Synthetic_flag Keypoint_flag Withheld_flag\n#> 1 1 1 FALSE FALSE FALSE\n#> 2 0 1 FALSE FALSE FALSE\n#> 3 0 1 FALSE FALSE FALSE\n#> 4 0 1 FALSE FALSE FALSE\n#> 5 0 1 FALSE FALSE FALSE\n#> 6 0 1 FALSE FALSE FALSE\n#> 7 1 1 FALSE FALSE FALSE\n#> 8 0 1 FALSE FALSE FALSE\n#> 9 0 1 FALSE FALSE FALSE\n#> 10 0 1 FALSE FALSE FALSE\n#> 11 1 1 FALSE FALSE FALSE\n#> 12 0 1 FALSE FALSE FALSE\n#> 13 0 1 FALSE FALSE FALSE\n#> 14 0 1 FALSE FALSE FALSE\n#> 15 0 1 FALSE FALSE FALSE\n#> 16 0 1 FALSE FALSE FALSE\n#> 17 0 1 FALSE FALSE FALSE\n#> 18 0 1 FALSE FALSE FALSE\n#> 19 0 1 FALSE FALSE FALSE\n#> 20 0 1 FALSE FALSE FALSE\n#> 21 0 1 FALSE FALSE FALSE\n#> 22 0 2 FALSE FALSE FALSE\n#> 23 1 1 FALSE FALSE FALSE\n#> 24 1 2 FALSE FALSE FALSE\n#> 25 0 1 FALSE FALSE FALSE\n#> 26 0 1 FALSE FALSE FALSE\n#> 27 0 1 FALSE FALSE FALSE\n#> 28 1 1 FALSE FALSE FALSE\n#> 29 0 2 FALSE FALSE FALSE\n#> 30 0 1 FALSE FALSE FALSE\n#> ScanAngleRank UserData PointSourceID X Y\n#> 1 -21 32 17 339002.9 5248001\n#> 2 -21 32 17 339003.0 5248000\n#> 3 -21 32 17 339002.9 5248000\n#> 4 -21 32 17 339002.9 5248000\n#> 5 -21 32 17 339003.6 5248000\n#> 6 -21 32 17 339003.5 5248000\n#> 7 -21 32 17 339003.6 5248000\n#> 8 -21 32 17 339003.7 5248000\n#> 9 -21 32 17 339003.6 5248000\n#> 10 -21 32 17 339003.7 5248000\n#> 11 -22 32 17 339009.6 5248001\n#> 12 -22 32 17 339009.5 5248001\n#> 13 -22 32 17 339009.2 5248000\n#> 14 -22 32 17 339009.4 5248001\n#> 15 -22 32 17 339009.3 5248000\n#> 16 -22 32 17 339009.3 5248000\n#> 17 -22 32 17 339009.3 5248000\n#> 18 -22 32 17 339009.5 5248000\n#> 19 -22 32 17 339009.5 5248000\n#> 20 -22 32 17 339009.5 5248000\n#> 21 -22 32 17 339009.6 5248000\n#> 22 -22 32 17 339010.7 5248001\n#> 23 -22 32 17 339009.6 5248000\n#> 24 -22 32 17 339010.6 5248001\n#> 25 -22 32 17 339009.5 5248000\n#> 26 -22 32 17 339009.4 5248000\n#> 27 -22 32 17 339009.7 5248000\n#> 28 -22 32 17 339015.1 5248000\n#> 29 -22 32 17 339015.1 5248000\n#> 30 -22 32 17 339015.0 5248000\n由reprex 包于 2022 年 3 月 18 日创建(v2.0.1)
\n