Make trait maps

Here we create the trait maps and export them as GeoTiffs.

Note: In this section we will use R because its provides more convinient raster processing pipelines.

library(raster)

# Load sPlot Data
sPlot <- read.csv("sPlotOpen/cwm_loc.csv")
# Load iNat Data
iNat <- read.csv("iNat_TRY_log.csv")

xy_1 <- cbind(sPlot$Longitude, sPlot$Latitude)
xy_2 <- cbind(iNat$decimalLongitude, iNat$decimalLatitude)

# raster for a 2 degree resolution map

r <- raster(ncols = 180, nrows = 90)

loop.vector <- 5:22 # loop over trait columns in sPlotOpen dataframe

for (i in loop.vector) { # Loop over loop.vector
	vals_1 <- exp(sPlot[, i])
	name1 <- colnames(sPlot[i])
	r1 <- rasterize(xy_1, r, vals_1, fun = mean)
	r1[is.infinite(r1)] <- NA
	crs(r1) <- "+proj=longlat"

	vals_2 <- exp(iNat[name1])
	r2 <- rasterize(xy_2, r, vals_2, fun = mean)
	r2[is.infinite(r2)] <- NA
	crs(r2) <- "+proj=longlat"

	# export as GeoTiff -->  separate file for each trait

	filename1 <- paste("sPlot_", name1, "_", deg, "deg.tif", sep = "")
	writeRaster(r1, filename1, overwrite = TRUE)
	filename2 <- paste("iNat_", name1, "_", deg, "deg.tif", sep = "")
	writeRaster(r2, filename2, overwrite = TRUE)
}