gsvd-process/process-gsvd-fluxnet.R

In jmzobitz/BRDF: A collection of files to compute BRDF using GSVD

### 1. Call a site name
### NEED TO LOAD UP TIDYVERSE AND GEIGEN FIRST?
library(tidyverse)
library(devtools)
library(BRDF)
# Pull up the data
site_list <- fluxnet %>% split(.$site)


# # ### Compute the least squares solutions
#  ls_process <- function(site_list) {
#
#    site_process <- site_list %>%
#      gather(key=band,value=value,band1:band7) %>%
#      split(.$band) %>%
#      map(.f=ls_solution) %>%
#      bind_rows(.id="band")
#
#    return(site_process)
#
#  }
#
#  # Pull up the data
#  ls_results <- fluxnet %>%
#    split(.$site) %>%
#    map(.f=ls_process) %>%
#    bind_rows(.id="site")
#
#
# save(ls_results,file='../Saved-Files-12.08/fluxnet-ls-results.Rda')
load('../Saved-Files-12.08/fluxnet-ls-results.Rda')

#######



### 2.1 Compute the kernel matrices

# Compute the kernel matrices, as a list by site
data_list <- map(site_list,~kernel_matrix(.x$site))

# Compute the GSVD, using the package geigen, save it for later
#gsvd_list <-map(data_list,gsvd_compute)
#save(gsvd_list,file='../Saved-Files-12.08/fluxnet-gsvd-nosnow-flag.Rda')

#load('../Saved-Files-12.08/fluxnet-gsvd.Rda')
load('../Saved-Files-12.08/fluxnet-gsvd-nosnow-flag.Rda')



### 4.1 Analyze the results to optimize across all bands (solution and norm)

# # Compute the smoothing lambda's, saving them for later
# lambda_list <- map2(gsvd_list,data_list,~gsvd_optimize(.x,.y$rho))
#
#
# ### Select out a band
#
# solution_list <- pmap(list(x=gsvd_list,y=lambda_list,z=data_list),.f=function(x,y,z){ gsvd_solution(x,y,z) })
#
# save(solution_list,file='../Saved-Files-12.08/solution-list-gsvd-nosnow-flag.Rda')

load('../Saved-Files-12.08/solution-list-gsvd-nosnow-flag.Rda')
# Process white sky albedo
albedo_list <- map(solution_list,albedo_compute)
use_data(albedo_list,overwrite = TRUE)


# Make our plots
source("gsvd-process/map-plot.R")
source("gsvd-process/data-plot.R")
source("gsvd-process/lCurve-plot.R")
source("gsvd-process/kernel-plot.R")
source("gsvd-process/kernel-taylor-plot.R")
source("gsvd-process/albedo-plot.R")
source("gsvd-process/albedo-taylor-plot.R")
source("gsvd-process/reflectance-taylor-plot.R")
source("gsvd-process/lambda-boxplot.R")

source("gsvd-process/negative-brdf-plot.R")
source('gsvd-process/sample-site.R')