search

caiohamamura / rGEDIsimulator

6 stars 1 forks source link

About GEDI Simulator #1

Open lisiqing981025 opened 1 year ago

lisiqing981025 commented 1 year ago

hi, When I want to use the GEDI waveform simulator, it tells me that the spot number is not 0 or anything else.

library(rGEDI) outdir <- "G:/data/simulate GEDI/ALS/test" lasfile_amazon <- file.path(outdir, "Amazon.las") lasfile_savanna <- file.path(outdir, "Savanna.las") library(lidR) lidR 4.0.4 using 4 threads. Help on . Bug report on <github.com/r-lidar/lidR>. library(lidR) library(plot3D) las_amazon<-readLAS(lasfile_amazon) las_savanna<-readLAS(lasfile_savanna) xcenter_amazon = mean(bbox(las_amazon)[1,]) ycenter_amazon = mean(bbox(las_amazon)[2,]) xcenter_savanna = mean(bbox(las_savanna)[1,]) ycenter_savanna = mean(bbox(las_savanna)[2,]) wf_amazon<-gediWFSimulator(input=lasfile_amazon,output=file.path(getwd(),"gediWF_amazon_simulation.h5"),coords = c(xcenter_amazon, ycenter_amazon)) Wave 0 of 1 Waveforms written to C:/Users/QDY.DESKTOP-9847DTA/Documents/gediWF_amazon_simulation.h5 wf_savanna<-gediWFSimulator(input=lasfile_savanna,output=file.path(getwd(),"gediWF_savanna_simulation.h5"),coords = c(xcenter_savanna, ycenter_savanna)) Wave 0 of 1 Waveforms written to C:/Users/QDY.DESKTOP-9847DTA/Documents/gediWF_savanna_simulation.h5

png("gediWf.png", width = 8, height = 6, units = 'in', res = 300) par(mfrow=c(2,2), mar=c(4,4,0,0), oma=c(0,0,1,1),cex.axis = 1.2) scatter3D(las_amazon@data$X,las_amazon@data$Y,las_amazon@data$Z,pch = 16,colkey = FALSE, main="",

  cex = 0.5,bty = "u",col.panel ="gray90",phi = 30,alpha=1,theta=45,
  col.grid = "gray50", xlab="UTM Easting (m)", ylab="UTM Northing (m)", zlab="Elevation (m)")

shot_number = 0 simulated_waveform_amazon = getLevel1BWF(wf_amazon, shot_number) Error in getLevel1BWF(wf_amazon, shot_number) : Shot number 0 was not found within the dataset!. Please try another shot number shot_number = 0 simulated_waveform_amazon = getLevel1BWF(wf_amazon, shot_number) Error in getLevel1BWF(wf_amazon, shot_number) : Shot number 0 was not found within the dataset!. Please try another shot number shot_number = 1 simulated_waveform_amazon = getLevel1BWF(wf_amazon, shot_number) Error in getLevel1BWF(wf_amazon, shot_number) : Shot number 1 was not found within the dataset!. Please try another shot number

How can I solve it? thanks!