R Lesen und Bearbeiten von LAS-Dateien

Question

Ich habe Probleme mit den & rLiDAR-Pakete zum Lesen in & LAS-Dateien bearbeiten. Ich führe die Beispiele in ihren PDF-Dokumenten durch, bekomme aber immer Fehler. Unten ist für das Rlasspaket.

> 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. 

Unten ist der rLiDAR Fehler ...

> 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 

ich meine LAS-Datei in jedem Paket Ordner in meinem system.file() ähnlich wie die Beispiele auf ihre PDF-Dateien gespeichert haben. Jede Hilfe, um diese Probleme zu lösen oder mich an eine andere Bibliothek zu verweisen, würde sehr geschätzt werden!

Answer 1

Sie rufen den Dateipfad nicht korrekt auf. Wenn Sie die Beispieldatei zugreifen möchten, können Sie system.file Funktion verwenden, aber wie so:

> 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 Siehe zu sehen, wo sich die Datei befindet.

Um Ihre eigene Datei zu importieren, geben Sie einfach den richtigen Pfad an, z.

lasfile <- C:/Users/Paul.Victor/Documents/myproject/myfile.laz