| Title: | Evapotranspiration and Energy Fluxes Spatial Analysis |
|---|---|
| Description: | Spatial modeling of energy balance and actual evapotranspiration using satellite images and meteorological data. Options of satellite are: Landsat-8 (with and without thermal bands), Sentinel-2 and MODIS. Respectively spatial resolutions are 30, 100, 10 and 250 meters. User can use data from a single meteorological station or a grid of meteorological stations (using any spatial interpolation method). Silva, Teixeira, and Manzione (2019) <doi:10.1016/j.envsoft.2019.104497>. |
| Authors: | Cesar de Oliveira Ferreira Silva [aut, cre], Antonio Heriberto de Castro Teixeira [ctb], Rodrigo Lilla Manzione [aut] |
| Maintainer: | Cesar de Oliveira Ferreira Silva <[email protected]> |
| License: | MIT + file LICENSE |
| Version: | 1.0.2 |
| Built: | 2024-10-09 05:00:03 UTC |
| Source: | https://github.com/cesarofs/agriwater |
Surface Albedo using Landsat-8 images.
albedo_l8(doy)albedo_l8(doy)
doy |
is the Day of Year (DOY) |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24").
Surface Albedo using MODIS images.
albedo_modis()albedo_modis()
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24").
library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate Sentinel-2 reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.01),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B1.tif"),filetype = "GTiff", overwrite=TRUE) # creating mask of study area mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" albedo_modis() #Exiting temporary folder and returning to previous workspace setwd(initial)library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate Sentinel-2 reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.01),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B1.tif"),filetype = "GTiff", overwrite=TRUE) # creating mask of study area mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" albedo_modis() #Exiting temporary folder and returning to previous workspace setwd(initial)
Surface Albedo using Sentinel-2 images.
albedo_s2()albedo_s2()
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24").
library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate Sentinel-2 reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.07, sd = 0.01), 2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B3.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.03, sd = 0.018),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B4.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B8.tif"),filetype = "GTiff", overwrite=TRUE) mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" albedo_s2() #Exiting temporary folder and returning to previous workspace setwd(initial)library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate Sentinel-2 reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.07, sd = 0.01), 2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B3.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.03, sd = 0.018),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B4.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B8.tif"),filetype = "GTiff", overwrite=TRUE) mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" albedo_s2() #Exiting temporary folder and returning to previous workspace setwd(initial)
Actual evapotranspiration (ETa) using Landsat-8 images with single agrometeorological data.
evapo_l8(doy, RG, Ta, ET0, a, b)evapo_l8(doy, RG, Ta, ET0, a, b)
doy |
is the Day of Year (DOY) |
RG |
is the global solar radiation |
Ta |
is the average air temperature |
ET0 |
is the reference evapotranspiration |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), net radiation ("Rn_MJ"), Crop Coefficient ("kc") and Actual Evapotranspiration (evapo).
Actual evapotranspiration (ETa) using Landsat-8 images with a grid of agrometeorological data.
evapo_l8_grid(doy, a, b)evapo_l8_grid(doy, a, b)
doy |
is the Day of Year (DOY) |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), net radiation ("Rn_MJ"), Crop Coefficient ("kc") and Actual Evapotranspiration (evapo).
Actual evapotranspiration (ETa) using Landsat-8 (including thermal bands) images with single agrometeorological data.
evapo_l8t(doy, RG, Ta, ET0, a, b)evapo_l8t(doy, RG, Ta, ET0, a, b)
doy |
is the Day of Year (DOY) |
RG |
is the global solar radiation |
Ta |
is the average air temperature |
ET0 |
is the reference evapotranspiration |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), net radiation ("Rn_MJ"), Crop Coefficient ("kc") and Actual Evapotranspiration (evapo).
Actual evapotranspiration (ETa) using Landsat-8 (including thermal bands) images with a grid of agrometeorological data.
evapo_l8t_grid(doy, a, b)evapo_l8t_grid(doy, a, b)
doy |
is the Day of Year (DOY) |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), net radiation ("Rn_MJ"), Crop Coefficient ("kc") and Actual Evapotranspiration (evapo).
Actual evapotranspiration (ETa) using MODIS with single agrometeorological data.
evapo_modis(doy, RG, Ta, ET0, a, b)evapo_modis(doy, RG, Ta, ET0, a, b)
doy |
is the Day of Year (DOY) |
RG |
is the global solar radiation |
Ta |
is the average air temperature |
ET0 |
is the reference evapotranspiration |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), net radiation ("Rn_MJ"), Crop Coefficient ("kc") and Actual Evapotranspiration (evapo).
library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate Sentinel-2 reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B1.tif"),filetype = "GTiff", overwrite=TRUE) mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" - it's the same procedure as the used for # evapo_l8(), evapo_l8t(), evapo_modis_grid(), evapo_l8_grid(), # evapo_l8t_grid(), evapo_s2() and evapo_s2_grid() evapo_modis(doy = 134, RG = 17.6, Ta = 27.9, ET0 = 3.8, a = 1.8, b = -0.008) #Exiting temporary folder and returning to previous workspace setwd(initial)library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate Sentinel-2 reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B1.tif"),filetype = "GTiff", overwrite=TRUE) mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" - it's the same procedure as the used for # evapo_l8(), evapo_l8t(), evapo_modis_grid(), evapo_l8_grid(), # evapo_l8t_grid(), evapo_s2() and evapo_s2_grid() evapo_modis(doy = 134, RG = 17.6, Ta = 27.9, ET0 = 3.8, a = 1.8, b = -0.008) #Exiting temporary folder and returning to previous workspace setwd(initial)
Actual evapotranspiration (ETa) using MODIS with a grid of agrometeorological data.
evapo_modis_grid(doy, a, b)evapo_modis_grid(doy, a, b)
doy |
is the Day of Year (DOY) |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), net radiation ("Rn_MJ"), Crop Coefficient ("kc") and Actual Evapotranspiration (evapo).
Actual evapotranspiration (ETa) using Sentinel-2 images with single agrometeorological data.
evapo_s2(doy, RG, Ta, ET0, a, b)evapo_s2(doy, RG, Ta, ET0, a, b)
doy |
is the Day of Year (DOY) |
RG |
is the global solar radiation |
Ta |
is the average air temperature |
ET0 |
is the reference evapotranspiration |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), net radiation ("Rn_MJ"), Crop Coefficient ("kc") and Actual Evapotranspiration (evapo).
library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate Sentinel-2 reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.07, sd = 0.01), 2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B3.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.03, sd = 0.018),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B4.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B8.tif"),filetype = "GTiff", overwrite=TRUE) mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" evapo_s2(doy = 134, RG = 17.6, Ta = 27.9, ET0 = 3.8, a = 1.8, b = -0.008) #Exiting temporary folder and returning to previous workspace setwd(initial)library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate Sentinel-2 reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.07, sd = 0.01), 2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B3.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.03, sd = 0.018),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B4.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B8.tif"),filetype = "GTiff", overwrite=TRUE) mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" evapo_s2(doy = 134, RG = 17.6, Ta = 27.9, ET0 = 3.8, a = 1.8, b = -0.008) #Exiting temporary folder and returning to previous workspace setwd(initial)
Actual evapotranspiration (ETa) using Sentinel-2 images with a grid of agrometeorological data.
evapo_s2_grid(doy, a, b)evapo_s2_grid(doy, a, b)
doy |
is the Day of Year (DOY) |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), net radiation ("Rn_MJ"), Crop Coefficient ("kc") and Actual Evapotranspiration (evapo).
Crop coefficient (ETa / ET0) using Landsat-8 images with single agrometeorological data.
kc_l8(doy, RG, Ta, a, b)kc_l8(doy, RG, Ta, a, b)
doy |
is the Day of Year (DOY) |
RG |
is the global solar radiation |
Ta |
is the average air temperature |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc") and net radiation ("Rn_MJ").
Crop coefficient (ETa / ET0) using Landsat-8 images with a grid of agrometeorological data.
kc_l8_grid(doy, a, b)kc_l8_grid(doy, a, b)
doy |
is the Day of Year (DOY) |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc") and net radiation ("Rn_MJ").
Crop coefficient (ETa / ET0) using Landsat-8 images (including thermal bands) with single agrometeorological data.
kc_l8t(doy, RG, Ta, a, b)kc_l8t(doy, RG, Ta, a, b)
doy |
is the Day of Year (DOY) |
RG |
is the global solar radiation |
Ta |
is the average air temperature |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc") and net radiation ("Rn_MJ").
Crop coefficient (ETa / ET0) using Landsat-8 images (including thermal bands) with a grid of agrometeorological data.
kc_l8t_grid(doy, a, b)kc_l8t_grid(doy, a, b)
doy |
is the Day of Year (DOY) |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc") and net radiation ("Rn_MJ").
Crop coefficient (ETa / ET0) using MODIS with single agrometeorological data.
kc_modis(doy, RG, Ta, a, b)kc_modis(doy, RG, Ta, a, b)
doy |
is the Day of Year (DOY) |
RG |
is the global solar radiation |
Ta |
is the average air temperature |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc") and net radiation ("Rn_MJ").
library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate MODIS reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.07, sd = 0.01), 2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"), filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B1.tif"), filetype = "GTiff", overwrite=TRUE) mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" kc_modis(doy = 134, RG = 17.6, Ta = 27.9, a = 1.8, b = -0.008) #Exiting temporary folder and returning to previous workspace setwd(initial)library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate MODIS reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.07, sd = 0.01), 2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"), filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B1.tif"), filetype = "GTiff", overwrite=TRUE) mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" kc_modis(doy = 134, RG = 17.6, Ta = 27.9, a = 1.8, b = -0.008) #Exiting temporary folder and returning to previous workspace setwd(initial)
Crop coefficient (ETa / ET0) using MODIS with a grid of agrometeorological data.
kc_modis_grid(doy, a, b)kc_modis_grid(doy, a, b)
doy |
is the Day of Year (DOY) |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc") and net radiation ("Rn_MJ").
Crop coefficient (ETa / ET0) using Sentinel-2 images with single agrometeorological data.
kc_s2(doy, RG, Ta, a, b)kc_s2(doy, RG, Ta, a, b)
doy |
is the Day of Year (DOY) |
RG |
is the global solar radiation |
Ta |
is the average air temperature |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc") and net radiation ("Rn_MJ").
library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate Sentinel-2 reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.07, sd = 0.01), 2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B3.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.03, sd = 0.018),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B4.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B8.tif"),filetype = "GTiff", overwrite=TRUE) mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" kc_s2(doy = 134, RG = 17.6, Ta = 27.9, a = 1.8, b = -0.008) #Exiting temporary folder and returning to previous workspace setwd(initial)library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate Sentinel-2 reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.07, sd = 0.01), 2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B3.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.03, sd = 0.018),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B4.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B8.tif"),filetype = "GTiff", overwrite=TRUE) mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" kc_s2(doy = 134, RG = 17.6, Ta = 27.9, a = 1.8, b = -0.008) #Exiting temporary folder and returning to previous workspace setwd(initial)
Crop coefficient (ETa / ET0) using Sentinel-2 images with a grid of agrometeorological data.
kc_s2_grid(doy, a, b)kc_s2_grid(doy, a, b)
doy |
is the Day of Year (DOY) |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc") and net radiation ("Rn_MJ").
Energy balance using Landsat-8 images with single agrometeorological data.
radiation_l8(doy, RG, Ta, ET0, a, b)radiation_l8(doy, RG, Ta, ET0, a, b)
doy |
is the Day of Year (DOY) |
RG |
is the global solar radiation |
Ta |
is the average air temperature |
ET0 |
is the reference evapotranspiration |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc"), Actual Evapotranspiration (evapo), latent heat flux "LE_MJ"), net radiation ("Rn_MJ"), ground heat flux ("G_MJ") and the sensible heat flux ("H_MJ").
Energy balance using Landsat-8 images with a grid of agrometeorological data.
radiation_l8_grid(doy, a, b)radiation_l8_grid(doy, a, b)
doy |
is the Day of Year (DOY) |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc"), Actual Evapotranspiration (evapo), latent heat flux "LE_MJ"), net radiation ("Rn_MJ"), ground heat flux ("G_MJ") and the sensible heat flux ("H_MJ").
Energy balance using Landsat-8 images (including thermal bands) with single agrometeorological data.
radiation_l8t(doy, RG, Ta, ET0, a, b)radiation_l8t(doy, RG, Ta, ET0, a, b)
doy |
is the Day of Year (DOY) |
RG |
is the global solar radiation |
Ta |
is the average air temperature |
ET0 |
is the reference evapotranspiration |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc"), Actual Evapotranspiration (evapo), latent heat flux "LE_MJ"), net radiation ("Rn_MJ"), ground heat flux ("G_MJ") and the sensible heat flux ("H_MJ").
Energy balance using Landsat-8 images (including thermal bands) with a grid of agrometeorological data.
radiation_l8t_grid(doy, a, b)radiation_l8t_grid(doy, a, b)
doy |
is the Day of Year (DOY) |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc"), Actual Evapotranspiration (evapo), latent heat flux "LE_MJ"), net radiation ("Rn_MJ"), ground heat flux ("G_MJ") and the sensible heat flux ("H_MJ").
Energy balance using Landsat-8 images with single agrometeorological data.
radiation_modis(doy, RG, Ta, ET0, a, b)radiation_modis(doy, RG, Ta, ET0, a, b)
doy |
is the Day of Year (DOY) |
RG |
is the global solar radiation |
Ta |
is the average air temperature |
ET0 |
is the reference evapotranspiration |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc"), Actual Evapotranspiration (evapo), latent heat flux "LE_MJ"), net radiation ("Rn_MJ"), ground heat flux ("G_MJ") and the sensible heat flux ("H_MJ").
library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate Sentinel-2 reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B1.tif"),filetype = "GTiff", overwrite=TRUE) # creating mask of study area mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" - it's the same procedure as the used for # radiation_l8(), radiation_l8t(), radiation_s2(), # radiation_l8_grid(), radiation_l8t_grid(), # radiation_s2_grid(), radiation_s2() and radiation_modis_grid() radiation_modis(doy = 134, RG = 17.6, Ta = 27.9, ET0 = 3.8, a = 1.8, b = -0.008) #Exiting temporary folder and returning to previous workspace setwd(initial)library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate Sentinel-2 reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B1.tif"),filetype = "GTiff", overwrite=TRUE) # creating mask of study area mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" - it's the same procedure as the used for # radiation_l8(), radiation_l8t(), radiation_s2(), # radiation_l8_grid(), radiation_l8t_grid(), # radiation_s2_grid(), radiation_s2() and radiation_modis_grid() radiation_modis(doy = 134, RG = 17.6, Ta = 27.9, ET0 = 3.8, a = 1.8, b = -0.008) #Exiting temporary folder and returning to previous workspace setwd(initial)
Energy balance using Landsat-8 images with a grid of agrometeorological data.
radiation_modis_grid(doy, a, b)radiation_modis_grid(doy, a, b)
doy |
is the Day of Year (DOY) |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc"), Actual Evapotranspiration (evapo), latent heat flux "LE_MJ"), net radiation ("Rn_MJ"), ground heat flux ("G_MJ") and the sensible heat flux ("H_MJ").
Energy balance using Sentinel-2 images with single agrometeorological data.
radiation_s2(doy, RG, Ta, ET0, a, b)radiation_s2(doy, RG, Ta, ET0, a, b)
doy |
is the Day of Year (DOY) |
RG |
is the global solar radiation |
Ta |
is the average air temperature |
ET0 |
is the reference evapotranspiration |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc"), Actual Evapotranspiration (evapo), latent heat flux "LE_MJ"), net radiation ("Rn_MJ"), ground heat flux ("G_MJ") and the sensible heat flux ("H_MJ").
library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate Sentinel-2 reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.07, sd = 0.01), 2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B3.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.03, sd = 0.018),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B4.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B8.tif"),filetype = "GTiff", overwrite=TRUE) mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" radiation_s2(doy = 134, RG = 17.6, Ta = 27.9, ET0 = 3.8, a = 1.8, b = -0.008) #Exiting temporary folder and returning to previous workspace setwd(initial)library(agriwater) # dependencies of package 'agriwater' library(terra) # Using a temporary folder to run example wd <- tempdir() initial = getwd() setwd(wd) # creating raster which simulate Sentinel-2 reflectances - for using # real data, please download: # https://drive.google.com/open?id=14E1wHNLxG7_Dh4I-GqNYakj8YJDgKLzk xy <- matrix(rnorm(4, mean = 0.07, sd = 0.01), 2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B2.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B3.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.03, sd = 0.018),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B4.tif"),filetype = "GTiff", overwrite=TRUE) xy <- matrix(rnorm(4, mean = 0.05, sd = 0.015),2, 2) rast <- rast(xy, crs="+proj=longlat +datum=WGS84") ext(rast) <- c(-40.5,-40.45,-9.5,-9.45) writeRaster(rast, file.path(wd, "B8.tif"),filetype = "GTiff", overwrite=TRUE) mask <- as.polygons(rast) writeVector(mask, file.path(getwd(),"mask.shp"), overwrite=TRUE) # using "agriwater" radiation_s2(doy = 134, RG = 17.6, Ta = 27.9, ET0 = 3.8, a = 1.8, b = -0.008) #Exiting temporary folder and returning to previous workspace setwd(initial)
Energy balance using Sentinel-2 images with a grid of agrometeorological data.
radiation_s2_grid(doy, a, b)radiation_s2_grid(doy, a, b)
doy |
is the Day of Year (DOY) |
a |
is one of the regression coefficients of SAFER algorithm |
b |
is one of the regression coefficients of SAFER algorithm |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temperature ("LST"), Crop Coefficient ("kc"), Actual Evapotranspiration (evapo), latent heat flux "LE_MJ"), net radiation ("Rn_MJ"), ground heat flux ("G_MJ") and the sensible heat flux ("H_MJ").
Reflectancies from Landsat-8 images.
reflectance_l8(doy)reflectance_l8(doy)
doy |
is the Day of Year (DOY) |
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24").