Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

about GEDI simulate #58

Closed
lisiqing981025 opened this issue Nov 2, 2023 · 1 comment
Closed

about GEDI simulate #58

lisiqing981025 opened this issue Nov 2, 2023 · 1 comment

Comments

@lisiqing981025
Copy link

lisiqing981025 commented Nov 2, 2023

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 <gis.stackexchange.com>. 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!

@caiohamamura
Copy link
Collaborator

Hi @lisiqing981025 , we are splitting up the rGEDIsimulator, we will not be maintaining it in this repository, you can post the issue in the new repository https://github.com/caiohamamura/rGEDIsimulator

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants