functions/loadtileGMSFromLas.R

loadTileGMSFromLas <- function(path, coordX, coordY){
  
  coordX<-str_pad(as.integer(floor(coordX/1000)*1000), 6, pad = "0")
  coordY<-str_pad(as.integer(floor(coordY/1000)*1000), 6, pad = "0")
  
  #uuid<-UUIDgenerate()
  multifileFlag<-checkIfMultiTile(path, coordX, coordY)
  multifiles <- NULL
  if(multifileFlag){
    multifiles<-list.files(path = path, pattern = paste0("ahn_", coordX,"_", coordY,"_*.*las"), full.names = T, recursive = T)
  }
  #dir.create(paste0(temp_dir,uuid,"/"))
  centralFile<-list.files(path = path, paste0("ahn_", coordX,"_", coordY,".las"), full.names = T, recursive = T)
  files<-c(centralFile,multifiles)
  if(length(files)!=0){
  #lapply(files,file.copy,to=paste0(temp_dir,uuid,"/"))
  #currentFiles<-list.files(path = paste0(temp_dir, uuid,"/"), full.names = TRUE)
  #print(paste(currentFiles))
  #lapply(paste(lasZipLocation, currentFiles), system)
  ##system("/usr/people/pagani/opt/testFile/LAStools/bin/laszip .laz")
  #system(paste0("rm ", temp_dir, uuid,"/*.laz"))
  files_las<-files
  out.matrix<-lapply(files_las, readLAS)
  outDF<-do.call(rbind, out.matrix)
  outDF<-data.frame(outDF)
  #out<-data.frame(out.matrix)
  #system(paste0("rm ", temp_dir, uuid,"/*.las"))
  
  rm(out.matrix,multifiles,centralFile,files,files_las)
  gc()
  return(outDF)
  }
  else {
    return(NULL)
  }

}
KNMI-DataLab/SkyViewFactor documentation built on May 30, 2019, 8 a.m.