R读取和编辑LAS文件

时间:2017-12-13 16:11:07

标签: r

我无法使用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上的示例。任何帮助解决这些问题或指导我到不同的图书馆将不胜感激!

2 个答案:

答案 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