我无法使用rlas& amp; rLiDAR用于阅读和阅读的套餐编辑LAS文件。我在他们的PDF文档上运行了示例,但我一直在收到错误。以下是rlas包。
> lasfile <- system.file("LAS", "42389364730000.las", package = "rlas")
> lasdata <- rlas::readlasdata(lasfile)
ERROR: wrong file signature '~Version Information
VERS. 2.0:
WRAP. YES:
END. PE'
ERROR: cannot open lasreaderlas with file name 'C:\Users\Paul.Victor\Documents\R\R-3.4.3\library\rlas\LAS\42389364730000.las'
Error: LASlib internal error. See message above.
以下是rLiDAR错误......
> lasfile <- system.file("LAS", "42389364730000.LAS", package = "rLiDAR")
>
> lasdata <- readLAS(lasfile, short = TRUE)
Error in readLAS(lasfile, short = TRUE) :
The LASfile input is not a valid LAS file
我已将我的LAS文件保存在system.file()中的每个包的文件夹中,类似于其PDF上的示例。任何帮助解决这些问题或指导我到不同的图书馆将不胜感激!
答案 0 :(得分:2)
问题是您正在使用的库是为Lidar LAS文件设计的,而您要读取的文件是测井LAS文件格式...具有相同文件扩展名的两种不同学科/数据类型。 / p>
答案 1 :(得分:1)
您没有正确调用文件路径。如果要访问示例文件,可以使用system.file
函数,但如下所示:
> library(rlas)
> lasfile <- system.file("extdata", "example.laz", package = "rlas")
> lasdata <- rlas::readlasdata(lasfile)
> str(lasdata)
Classes ‘data.table’ and 'data.frame': 30 obs. of 13 variables:
$ X : num 339003 339003 339003 339003 339004 ...
$ Y : num 5248001 5248000 5248000 5248000 5248000 ...
$ Z : num 976 975 974 974 974 ...
$ gpstime : num 269347 269347 269347 269347 269347 ...
$ Intensity : int 82 54 27 55 117 81 84 104 91 99 ...
$ ReturnNumber : int 1 1 2 2 1 1 1 1 1 1 ...
$ NumberOfReturns : int 1 1 2 2 1 1 1 1 1 1 ...
$ ScanDirectionFlag: int 1 1 1 1 0 0 1 1 1 1 ...
$ EdgeOfFlightline : int 1 0 0 0 0 0 1 0 0 0 ...
$ Classification : int 1 1 1 1 1 1 1 1 1 1 ...
$ ScanAngle : int -21 -21 -21 -21 -21 -21 -21 -21 -21 -21 ...
$ UserData : int 32 32 32 32 32 32 32 32 32 32 ...
$ PointSourceID : int 17 17 17 17 17 17 17 17 17 17 ...
- attr(*, ".internal.selfref")=<externalptr>
请参阅here以查看文件的位置。
要导入您自己的文件,只需指定正确的路径,例如
lasfile <- C:/Users/Paul.Victor/Documents/myproject/myfile.laz