data-raw/wq_algorithms_prep.R

In RAJohansen/waterquality: Satellite Derived Water Quality Detection Algorithms

library(tidyverse)
library(waterquality)

wq_algorithms = tribble(
  ~name, ~funs,~worldview2, ~sentinel2, ~landsat8, ~modis, ~meris, ~OLCI, ~type,
  "Al10SABI", Al10SABI, c(7, 5, 2, 3), c(9, 4, 2, 3), c(5, 4, 2, 3), NA, c(13, 7, 2, 4), c(17, 8, 3, 5), "chlorophyll",
  "Am092Bsub", Am092Bsub, NA, c(5, 4), NA, NA, c(8, 7), c(10, 8), "chlorophyll",
  "Am09KBBI", Am09KBBI, NA, c(9,8), NA, NA, c(8,7),c(10, 8), "phycocyanin",
  "Be16FLHblue_WV2", Be16FLHblue_WV2, c(3, 5, 2), NA, NA, NA, NA, NA,"chlorophyll",  
  "Be16FLHblue_S2", Be16FLHblue_S2, NA, c(3,4,2), NA, NA, NA, NA,"chlorophyll",  
  "Be16FLHblue_LS8", Be16FLHblue_LS8, NA, NA, c(3, 4, 2), NA, NA, NA,"chlorophyll",  
  "Be16FLHblue_MERIS", Be16FLHblue_MERIS, NA, NA, NA, NA, c(4, 7, 3), NA,"chlorophyll",  
  "Be16FLHblue_OLCI", Be16FLHblue_OLCI, NA, NA, NA, NA, NA, c(5, 8, 4),"chlorophyll",  
  "Be16FLHviolet_WV2", Be16FLHviolet_WV2, c(3, 5, 1), NA, NA, NA, NA, NA,"chlorophyll",  
  "Be16FLHviolet_S2", Be16FLHviolet_S2, NA, c(3,4,1), NA, NA, NA, NA,"chlorophyll",  
  "Be16FLHviolet_LS8", Be16FLHviolet_LS8, NA, NA, c(3, 4, 1), NA, NA, NA,"chlorophyll",  
  "Be16FLHviolet_MERIS", Be16FLHviolet_MERIS, NA, NA, NA, NA, c(4, 7, 2), NA,"chlorophyll",  
  "Be16FLHviolet_OLCI", Be16FLHviolet_OLCI, NA, NA, NA, NA, NA, c(5, 8, 3),"chlorophyll",  
  "Be16NDPhyI", Be16NDPhyI, c(6, 4), NA, NA, NA, c(9, 6), c(11, 7), "phycocyanin",
  "De933BDA", De933BDA, c(4,5,4), NA, NA, NA, c(6,7,6), c(7,8,7), "chlorophyll",
  "Gi033BDA", Gi033BDA, NA, c(4, 5, 6), NA, NA, c(8, 9, 10), c(10, 11, 12),"chlorophyll", 
  "Go04MCI", Go04MCI, NA, NA, NA, NA, c(9, 8, 10), c(11, 10, 12), "phycocyanin",
  "HU103BDA", HU103BDA, c(4, 4, 6), NA, NA, NA, c(6, 6, 9), c(7, 7, 11), "phycocyanin",
  "Kn07KIVU", Kn07KIVU, c(2, 5, 3), NA, c(1, 4, 3), NA, c(2, 7, 4), c(3, 8, 5), "chlorophyll", 
  "MI092BDA", MI092BDA, c(6,4), NA, NA, NA, c(9, 6), c(11, 7), "phycocyanin", 
  "MM092BDA", MM092BDA, c(6,4), NA, NA, NA, c(9,6), c(11,7), "phycocyanin",
  "MM12NDCI", MM12NDCI, c(6, 5), c(5, 4), NA, NA, c(9, 8), c(11, 10),"chlorophyll",
  "MM143BDAopt", MM143BDAopt, c(4,5, 6), NA, NA, NA, c(6, 7, 9), c(7, 8, 11), "phycocyanin",
  "SI052BDA", SI052BDA, c(6, 4), NA, NA, NA, c(9, 6), c(11, 7),"phycocyanin",
  "SM122BDA", SM122BDA, c(6, 4), NA, NA, NA, c(9, 6), c(11, 7), "phycocyanin",
  "SY002BDA", SY002BDA, c(5, 4), NA, NA, NA, c(7, 6), c(8, 7), "phycocyanin",
  "Be16NDTIblue", Be16NDTIblue, c(5, 2), c(4, 2), c(4, 2), NA, c(7, 3), c(8, 4), "chlorophyll",
  "Be16NDTIviolet", Be16NDTIviolet, c(5, 1), c(4, 1), c(4, 1), NA, c(7, 2), c(8, 3),"chlorophyll",  
  "Be16FLHBlueRedNIR_WV2", Be16FLHBlueRedNIR_WV2, c(5, 7, 2), NA, NA, NA, NA, NA,"phycocyanin",
  "Be16FLHBlueRedNIR_S2", Be16FLHBlueRedNIR_S2, NA, c(4, 9, 2), NA, NA, NA, NA,"phycocyanin",
  "Be16FLHBlueRedNIR_LS8", Be16FLHBlueRedNIR_LS8, NA, NA, c(4, 5, 2), NA, NA, NA,"phycocyanin",
  "Be16FLHBlueRedNIR_MERIS", Be16FLHBlueRedNIR_MERIS, NA, NA, NA, NA, c(7, 13, 3), NA,"phycocyanin",
  "Be16FLHBlueRedNIR_OLCI", Be16FLHBlueRedNIR_OLCI, NA, NA, NA, NA, NA, c(8, 17, 4),"phycocyanin",
  "Be16FLHGreenRedNIR_WV2", Be16FLHGreenRedNIR_WV2, c(5, 7, 3), NA, NA, NA, NA, NA,"phycocyanin",
  "Be16FLHGreenRedNIR_S2", Be16FLHGreenRedNIR_WV2, NA, c(4, 9, 3), NA, NA, NA, NA,"phycocyanin",
  "Be16FLHGreenRedNIR_LS8", Be16FLHGreenRedNIR_WV2, NA, NA, c(4, 5, 3), NA, NA, NA,"phycocyanin",
  "Be16FLHGreenRedNIR_MERIS", Be16FLHGreenRedNIR_WV2, NA, NA, NA, NA, c(7, 13, 5), NA,"phycocyanin",
  "Be16FLHGreenRedNIR_OLCI", Be16FLHGreenRedNIR_WV2, NA, NA, NA, NA, NA, c(8, 17, 6),"phycocyanin",
  "Be16FLHVioletRedNIR_WV2", Be16FLHVioletRedNIR_WV2, c(5, 7, 1), NA, NA, NA, NA, NA,"phycocyanin",
  "Be16FLHVioletRedNIR_S2", Be16FLHVioletRedNIR_S2, NA, c(4, 9, 1), NA, NA, NA, NA,"phycocyanin",
  "Be16FLHVioletRedNIR_LS8", Be16FLHVioletRedNIR_LS8, NA, NA, c(4, 5, 1), NA, NA, NA,"phycocyanin",
  "Be16FLHVioletRedNIR_MERIS", Be16FLHVioletRedNIR_MERIS, NA, NA, NA, NA, c(7, 13, 2), NA,"phycocyanin",
  "Be16FLHVioletRedNIR_OLCI", Be16FLHVioletRedNIR_OLCI, NA, NA, NA, NA, NA, c(8, 17, 3),"phycocyanin",
  "Wy08CI", Wy08CI, NA, NA, NA, NA, c(8, 7, 9), c(10, 8, 11), "phycocyanin",
  "Da052BDA", Da052BDA, c(6, 5), c(5, 4), NA, NA, c(9, 7), c(11, 8), "phycocyanin",
  "Be162B643sub629", Be162B643sub629, c(5, 4), NA, NA, NA, c(7,6), c(8,7), "phycocyanin",
  "Be162B700sub601", Be162B700sub601, c(6, 4), NA, NA, NA, c(9, 6), c(11, 7), "phycocyanin", 
  "Be162BsubPhy", Be162BsubPhy, c(6, 4), NA, NA, NA, c(9, 6), c(11, 7), "phycocyanin",
  "Be16NDPhyI644over615", Be16NDPhyI644over615, c(5, 4), NA, NA, NA, c(7, 6), c(8, 7), "phycocyanin",
  "Be16NDPhyI644over629", Be16NDPhyI644over629, c(5, 4), NA, NA, NA, c(7, 6), c(8, 7), "phycocyanin",
  "Be16Phy2BDA644over629", Be16Phy2BDA644over629, c(5, 4), NA, NA, NA, c(7, 6), c(8, 7), "phycocyanin",
  "MM12NDCIalt", MM12NDCIalt, c(6, 5), NA, NA, NA, c(9, 7), c(11, 8), "phycocyanin",
  "TurbBe16GreenPlusRedBothOverViolet", TurbBe16GreenPlusRedBothOverViolet, c(3, 5, 1), c(3, 4, 1), c(3, 4, 1), NA, c(5, 7, 2), c(6, 8, 3), "turbidity",
  "TurbBe16RedOverViolet", TurbBe16RedOverViolet, c(5, 1), c(4, 1), c(4, 1), NA, c(7, 2), c(8, 2), "turbidity",
  "TurbBow06RedOverGreen", TurbBow06RedOverGreen, c(5, 3), c(4, 3), c(4, 3), NA, c(7, 5), c(8, 6),"turbidity",
  "TurbChip09NIROverGreen", TurbChip09NIROverGreen, c(7, 3), c(9, 3), c(5, 3), NA, c(13, 5), c(17, 6),"turbidity",
  "TurbDox02NIRoverRed", TurbDox02NIRoverRed, c(7, 5), c(9, 4), c(5, 4), c(2, 1), c(13, 7), c(17, 8),"turbidity",
  "TurbFrohn09GreenPlusRedBothOverBlue", TurbFrohn09GreenPlusRedBothOverBlue, c(3, 5, 1), c(3, 4, 2), c(3, 4, 2), NA, c(5, 7, 2), c(6, 8, 3),"turbidity",
  "TurbHarr92NIR", TurbHarr92NIR, c(7), c(9), c(5), c(2), c(13), c(17), "turbidity",
  "TurbLath91RedOverBlue", TurbLath91RedOverBlue, c(5, 2), c(4, 2), c(4, 2), NA, c(7, 2), c(8, 3), "turbidity",
  "TurbMoore80Red", TurbMoore80Red, c(5), c(4), c(4), c(1), c(7), c(8), "turbidity"
)

wq_algorithms = wq_algorithms %>% 
  gather(satellite, bands, worldview2:OLCI) %>% 
  filter(!is.na(bands))

# save data
usethis::use_data(wq_algorithms, overwrite = TRUE)