Extraction scripts for Xue and Zhang datasets

.gitignore

@@ -39,3 +39,4 @@ vignettes/*.pdf
 .Renviron
 
 *.Rproj
+data/

README.md

@@ -17,3 +17,5 @@ Climate data were provided by Dr. Abby Frazier in 2022. Descriptions of the data
 * Only 1 Cf
 * Rain and T2
 
+- Right now Xue doesn't add anything from Zhang. Try getting daily files from NCAR
+

Xue_extraction.R

@@ -0,0 +1,14 @@
+library(tidyr)
+library(dplyr)
+library(ncdf4)
+library(stars)
+library(raster)
+library(here)
+
+DataDir <- "D:/HI_Data/Xue/"
+
+Xue <- read_ncdf(paste0(DataDir,"hist_mean_rain.nc"))
+
+nc<-nc_open(paste0(DataDir,"hist_mean_rain.nc"))
+nc$dim
+

Zhang_extraction.R

@@ -0,0 +1,25 @@
+library(tidyr)
+library(dplyr)
+library(ncdf4)
+library(stars)
+library(raster)
+library(here)
+
+DataDir <- "D:/HI_Data/Zhang/"
+
+test.tif <- raster("D:/HI_Data/Zhang/Processed/HI_Downscaling_Share/DynDS_RainfallTIF_AllFolders/DynDS_7MonthlyTIF_byIsland_250m/Maui/present/ma_rf_pres_1990_01_250m.tif")
+plot(test.tif)
+crs(test.tif)
+
+boundary <- st_read('./data/HALE/HALE_boundary.shp')
+boundary <- st_transform(boundary, st_crs(test.tif))
+
+plot(test.tif)
+plot(boundary[1])
+
+t<-mask(test.tif,boundary)
+plot(t)
+
+Zhang <- read_stars("D:/HI_Data/Zhang/Processed/HI_Downscaling_Share/DynDS_RainfallTIF_AllFolders/DynDS_7MonthlyTIF_byIsland_250m/Maui/present/ma_rf_pres_1990_01_250m.tif")
+Zhang
+plot(Zhang)