modif path

interpolation.R

+library(tidyverse)
+require(fasterize)
+require(raster)
+library(terra)
+library(sfExtras)
+library(sf)
+library(gstat)
+
+## pour Edith
+
+load("~/recherche/beyond/Camille/data camille/final_data.RData")
+
+## pour Camille
+
+setwd("C:/Users/BioSP/Documents/Th?se/data/qgis_output/resultats final")
+load("final_data.RData")
+
+final_data$Altitude <- as.numeric(final_data$Altitude)
+
+################################################################################
+
+xy = st_centroid_coords(final_data)
+mydata = final_data %>%  st_drop_geometry()
+mydata=cbind(mydata,xy)
+
+coordinates(mydata)<- c("longitude","latitude")
+crs(mydata) = crs(final_data)
+
+sb<-final_data
+sb_ras <- raster(sb, res = 500,  crs = st_crs(final_data))
+sb_ras_ID<-fasterize(sb,sb_ras,field = "ID")
+
+################################################################################
+
+interpola <- function(var){
+  myvar <- deparse(substitute(var))
+  sel=!is.na(final_data[[myvar]])
+  mg <- gstat(id = myvar,
+              formula =var~1,
+              data=mydata[sel,],
+              nmax=12,
+              set=list(idp = 2))
+  
+  sb_ras_v<-fasterize(sb,sb_ras,field = myvar)
+  
+  z <- interpolate(sb_ras_v, mg, xyNames=c("longitude","latitude"))
+  idp <- mask(z, sb_ras_ID)
+  return(idp)
+}
+
+interpola(var = drainc_q)