extra/PlanetScope_RSR_SatID_0f_10.R

In belasi01/HyperocR: Hyperspectral OCR data processing

library(plotly)

##### convert RSR to RData

PlanetScope_0F_10=read.table("extra/PlanetScope_RSR_SatID_0f_10.csv", skip=1, sep=",")
names(PlanetScope_0F_10) <- c("waves", "B1", "B2", "B3", "B4")
PlanetScope_0F_10$waves <- PlanetScope_0F_10$waves*1000

p <- plot_ly(x=~PlanetScope_0F_10$waves,y=~PlanetScope_0F_10$B1, type = 'scatter', mode = 'lines', fill = 'tozeroy', color = I("blue"), name ="B1") %>%
  add_trace(x=~PlanetScope_0F_10$waves,y=~PlanetScope_0F_10$B2, name = 'B2', color = I("green"), fill = 'tozeroy') %>%
  add_trace(x=~PlanetScope_0F_10$waves,y=~PlanetScope_0F_10$B3, name = 'B3', color = I("red"), fill = 'tozeroy') %>%
  add_trace(x=~PlanetScope_0F_10$waves,y=~PlanetScope_0F_10$B4, name = 'B4', color = I("black"), fill = 'tozeroy') %>%
  layout(xaxis = list(title = 'Wavelenght', range=c(400,900)),
         yaxis = list(title = 'RSR'))
print(p)

save(PlanetScope_0F_10, file = "data/PlanetScope_0F_10_RSRs.RData")