Package Agriwater': January 30, 2019

Package ‘agriwater’
January 30, 2019

Title Evapotranspiration and Energy Fluxes Spatial Analysis
Version 1.0.0
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). Teixeira (2010) <doi:10.3390/rs0251287>.
Teixeira et al. (2015) <doi:10.3390/rs71114597>.
Silva, Manzione, and Albuquerque Filho (2018) <doi:10.3390/horticulturae4040044>.
Maintainer Cesar de Oliveira Ferreira Silva <cesaroliveira.f.silva@gmail.com>
Depends R (>= 3.2.0)
License MIT + file LICENSE
Encoding UTF-8
LazyData true
BugReports https://github.com/cesarofs/agriwater/issues
Imports raster, sp, rgdal
RoxygenNote 6.1.0
Suggests knitr, rmarkdown
VignetteBuilder knitr
NeedsCompilation no
Author Cesar de Oliveira Ferreira Silva [cre]
(<https://orcid.org/0000-0002-5152-6497>),
Antonio Heriberto de Castro Teixeira [aut],
Rodrigo Lilla Manzione [ctb] (<https://orcid.org/0000-0002-0754-2641>)
Repository CRAN
Date/Publication 2019-01-30 22:34:00 UTC
1
2 albedo_l8
R topics documented:
albedo_l8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
albedo_modis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
albedo_s2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
evapo_l8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
evapo_l8t . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
evapo_l8t_grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
evapo_l8_grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
evapo_modis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
evapo_modis_grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
evapo_s2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
evapo_s2_grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
kc_l8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
kc_l8t . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
kc_l8t_grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
kc_l8_grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
kc_modis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
kc_modis_grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
kc_s2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
kc_s2_grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
radiation_l8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
radiation_l8t . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
radiation_l8t_grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
radiation_l8_grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
radiation_modis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
radiation_modis_grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
radiation_s2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
radiation_s2_grid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
reflectance_l8 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
Index 22
albedo_l8 Surface Albedo using Landsat-8 images.
Description
Surface Albedo using Landsat-8 images.
Usage
albedo_l8(doy)
Arguments
doy is the Day of Year (DOY)
albedo_modis 3
Value
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24").
albedo_modis Surface Albedo using MODIS images.
Description
Surface Albedo using MODIS images.
Usage
albedo_modis()
Value
albedo_s2 Surface Albedo using Sentinel-2 images.
Description
Surface Albedo using Sentinel-2 images.
Usage
albedo_s2()
Value
4 evapo_l8t
evapo_l8 Actual evapotranspiration (ETa) using Landsat-8 images with single

agrometeorological data.
Description
Actual evapotranspiration (ETa) using Landsat-8 images with single agrometeorological data.
Usage
evapo_l8(doy, RG, Ta, ET0, a, b)
Arguments

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
Value
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Tem-
perature ("LST"), net radiation ("Rn_MJ"), Crop Coefficient ("kc") and Actual Evapotranspiration
(evapo).
evapo_l8t Actual evapotranspiration (ETa) using Landsat-8 (including thermal

bands) images with single agrometeorological data.
Description
Actual evapotranspiration (ETa) using Landsat-8 (including thermal bands) images with single
Usage
evapo_l8t(doy, RG, Ta, ET0, a, b)

evapo_l8t_grid 5
Arguments

Value
(evapo).
evapo_l8t_grid Actual evapotranspiration (ETa) using Landsat-8 (including thermal

bands) images with a grid of agrometeorological data.
Description
Actual evapotranspiration (ETa) using Landsat-8 (including thermal bands) images with a grid of
Usage
evapo_l8t_grid(doy, a, b)
Arguments

Value
(evapo).
6 evapo_modis
evapo_l8_grid Actual evapotranspiration (ETa) using Landsat-8 images with a grid

of agrometeorological data.
Description
Actual evapotranspiration (ETa) using Landsat-8 images with a grid of agrometeorological data.
Usage
evapo_l8_grid(doy, a, b)
Arguments
Value
(evapo).
evapo_modis Actual evapotranspiration (ETa) using MODIS with single agromete-

orological data.
Description
Actual evapotranspiration (ETa) using MODIS with single agrometeorological data.
Usage
evapo_modis(doy, RG, Ta, ET0, a, b)
Arguments
evapo_modis 7
Value
(evapo).
Examples
library(agriwater)
# dependencies of package 'agriwater'

library(sp)
library(raster)
library(rgdal)
# 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
wd <- tempdir()
xy <- matrix(rnorm(100, mean = 0.07, sd = 0.01), 10, 10)
rast <- raster(xy)
extent(rast) <- c(-40.5,-40.45,-9.5,-9.45)
projection(rast) <- CRS("+proj=longlat +datum=WGS84")
writeRaster(rast, file.path(wd, "B2"),format = "GTiff", overwrite=TRUE)
xy <- matrix(rnorm(100, mean = 0.05, sd = 0.015),10, 10)
rast <- raster(xy)
extent(rast) <- c(-40.5,-40.45,-9.5,-9.45)
mask <- as(extent(rast), 'SpatialPolygons')

projection(mask) <- CRS("+proj=longlat +datum=WGS84")
shapefile(mask, file.path(wd,"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()
DOY = 134
a = 1.8
b = -0.008
RG = 17.6
Ta = 27.9
ET0 = 3.8
evapo_modis(DOY, RG, Ta, ET0, a, b)
#Exiting temporary folder and returning to previous workspace

8 evapo_s2
setwd(initial)
evapo_modis_grid Actual evapotranspiration (ETa) using MODIS with a grid of agrome-

teorological data.
Description
Actual evapotranspiration (ETa) using MODIS with a grid of agrometeorological data.
Usage
evapo_modis_grid(doy, a, b)
Arguments
Value
(evapo).
evapo_s2 Actual evapotranspiration (ETa) using Sentinel-2 images with single

Description
Actual evapotranspiration (ETa) using Sentinel-2 images with single agrometeorological data.
Usage
evapo_s2(doy, RG, Ta, ET0, a, b)
Arguments
evapo_s2_grid 9
Value
(evapo).
evapo_s2_grid Actual evapotranspiration (ETa) using Sentinel-2 images with a grid

of agrometeorological data.
Description
Actual evapotranspiration (ETa) using Sentinel-2 images with a grid of agrometeorological data.
Usage
evapo_s2_grid(doy, a, b)
Arguments
Value
(evapo).
kc_l8 Crop coefficient (ETa / ET0) using Landsat-8 images with single
Description
Crop coefficient (ETa / ET0) using Landsat-8 images with single agrometeorological data.
Usage
kc_l8(doy, RG, Ta, a, b)
10 kc_l8t
Arguments

Value
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Temper-
ature ("LST"), Crop Coefficient ("kc") and net radiation ("Rn_MJ").
kc_l8t Crop coefficient (ETa / ET0) using Landsat-8 images (including ther-
mal bands) with single agrometeorological data.
Description
Crop coefficient (ETa / ET0) using Landsat-8 images (including thermal bands) with single agrom-
eteorological data.
Usage
kc_l8t(doy, RG, Ta, a, b)
Arguments

Value
kc_l8t_grid 11
kc_l8t_grid Crop coefficient (ETa / ET0) using Landsat-8 images (including ther-
mal bands) with a grid of agrometeorological data.
Description
Crop coefficient (ETa / ET0) using Landsat-8 images (including thermal bands) with a grid of
Usage
kc_l8t_grid(doy, a, b)
Arguments
Value
kc_l8_grid Crop coefficient (ETa / ET0) using Landsat-8 images with a grid of
Description
Crop coefficient (ETa / ET0) using Landsat-8 images with a grid of agrometeorological data.
Usage
kc_l8_grid(doy, a, b)
Arguments
Value
12 kc_modis
kc_modis Crop coefficient (ETa / ET0) using MODIS with single agrometeoro-
logical data.
Description
Crop coefficient (ETa / ET0) using MODIS with single agrometeorological data.
Usage
kc_modis(doy, RG, Ta, a, b)
Arguments
Value
Examples
library(agriwater)

library(sp)
library(raster)
library(rgdal)

wd <- tempdir()
initial = getwd()
setwd(wd)
# creating raster which simulate MODIS reflectances - for using


rast <- raster(xy)
extent(rast) <- c(-40.5,-40.45,-9.5,-9.45)
kc_modis_grid 13

rast <- raster(xy)
extent(rast) <- c(-40.5,-40.45,-9.5,-9.45)


# kc_l8(), kc_l8t(), kc_s2(), kc_l8_grid(), kc_l8t_grid(),
# kc_s2(), kc_s2_grid() and kc_modis_grid()
DOY = 134
a = 1.8
b = -0.008
RG = 17.6
Ta = 27.9
kc_modis(DOY, RG, Ta, a, b)

setwd(initial)
kc_modis_grid Crop coefficient (ETa / ET0) using MODIS with a grid of agrometeo-
rological data.
Description
Crop coefficient (ETa / ET0) using MODIS with a grid of agrometeorological data.
Usage
kc_modis_grid(doy, a, b)
Arguments
Value
14 kc_s2_grid
kc_s2 Crop coefficient (ETa / ET0) using Sentinel-2 images with single
Description
Crop coefficient (ETa / ET0) using Sentinel-2 images with single agrometeorological data.
Usage
kc_s2(doy, RG, Ta, a, b)
Arguments
Value
kc_s2_grid Crop coefficient (ETa / ET0) using Sentinel-2 images with a grid of
Description
Crop coefficient (ETa / ET0) using Sentinel-2 images with a grid of agrometeorological data.
Usage
kc_s2_grid(doy, a, b)
Arguments
Value
radiation_l8 15
radiation_l8 Energy balance using Landsat-8 images with single agrometeorologi-

cal data.
Description
Energy balance using Landsat-8 images with single agrometeorological data.
Usage
radiation_l8(doy, RG, Ta, ET0, a, b)
Arguments

Value
It returns in raster format (.tif) the Surface Albedo at 24h scale ("Alb_24"), NDVI, Surface Tempera-
ture ("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").
radiation_l8t Energy balance using Landsat-8 images (including thermal bands)

with single agrometeorological data.
Description
Energy balance using Landsat-8 images (including thermal bands) with single agrometeorological
data.
Usage
radiation_l8t(doy, RG, Ta, ET0, a, b)

16 radiation_l8t_grid
Arguments

Value
radiation_l8t_grid Energy balance using Landsat-8 images (including thermal bands)

with a grid of agrometeorological data.
Description
Energy balance using Landsat-8 images (including thermal bands) with a grid of agrometeorological
data.
Usage
radiation_l8t_grid(doy, a, b)
Arguments

Value
radiation_l8_grid 17
radiation_l8_grid Energy balance using Landsat-8 images with a grid of agrometeoro-

logical data.
Description
Energy balance using Landsat-8 images with a grid of agrometeorological data.
Usage
radiation_l8_grid(doy, a, b)
Arguments
Value
radiation_modis Energy balance using Landsat-8 images with single agrometeorologi-

cal data.
Description
Energy balance using Landsat-8 images with single agrometeorological data.
Usage
radiation_modis(doy, RG, Ta, ET0, a, b)
Arguments
18 radiation_modis
Value
Examples
library(agriwater)

library(sp)
library(raster)
library(rgdal)

wd <- tempdir()
initial = getwd()
setwd(wd)
# creating raster which simulate Sentinel-2 reflectances - for using

wd <- tempdir()
setwd <- tempdir()
rast <- raster(xy)
extent(rast) <- c(-40.5,-40.45,-9.5,-9.45)
rast <- raster(xy)
extent(rast) <- c(-40.5,-40.45,-9.5,-9.45)
# creating mask of study area


# radiation_l8(), radiation_l8t(), radiation_s2(),
# radiation_l8_grid(), radiation_l8t_grid(),
# radiation_s2_grid(), radiation_s2() and radiation_modis_grid()
DOY = 134
a = 1.8
b = -0.008
RG = 17.6
Ta = 27.9
ET0 = 3.8
radiation_modis(DOY, RG, Ta, ET0, a, b)
radiation_modis_grid 19

setwd(initial)
radiation_modis_grid Energy balance using Landsat-8 images with a grid of agrometeoro-

logical data.
Description
Energy balance using Landsat-8 images with a grid of agrometeorological data.
Usage
radiation_modis_grid(doy, a, b)
Arguments

Value
radiation_s2 Energy balance using Sentinel-2 images with single agrometeorologi-

cal data.
Description
Energy balance using Sentinel-2 images with single agrometeorological data.
Usage
radiation_s2(doy, RG, Ta, ET0, a, b)

20 radiation_s2_grid
Arguments

Value
radiation_s2_grid Energy balance using Sentinel-2 images with a grid of agrometeoro-

logical data.
Description
Energy balance using Sentinel-2 images with a grid of agrometeorological data.
Usage
radiation_s2_grid(doy, a, b)
Arguments

Value
reflectance_l8 21
reflectance_l8 Reflectancies from Landsat-8 images.
Description
Reflectancies from Landsat-8 images.
Usage
reflectance_l8(doy)
Arguments
Value
Index
albedo_l8, 2
albedo_modis, 3
albedo_s2, 3
evapo_l8, 4
evapo_l8_grid, 6
evapo_l8t, 4
evapo_l8t_grid, 5
evapo_modis, 6
evapo_modis_grid, 8
evapo_s2, 8
evapo_s2_grid, 9
kc_l8, 9
kc_l8_grid, 11
kc_l8t, 10
kc_l8t_grid, 11
kc_modis, 12
kc_modis_grid, 13
kc_s2, 14
kc_s2_grid, 14
radiation_l8, 15
radiation_l8_grid, 17
radiation_l8t, 15
radiation_l8t_grid, 16
radiation_modis, 17
radiation_modis_grid, 19
radiation_s2, 19
radiation_s2_grid, 20
reflectance_l8, 21
22

Package Agriwater': January 30, 2019

Uploaded by

Copyright:

Available Formats

Package Agriwater': January 30, 2019

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Package Agriwater': January 30, 2019

Uploaded by

Copyright:

Available Formats

Package ‘agriwater’

January 30, 2019

albedo_l8 Surface Albedo using Landsat-8 images.

albedo_modis Surface Albedo using MODIS images.

Surface Albedo using MODIS images.

albedo_s2 Surface Albedo using Sentinel-2 images.

Surface Albedo using Sentinel-2 images.

evapo_l8 Actual evapotranspiration (ETa) using Landsat-8 images with single

evapo_l8(doy, RG, Ta, ET0, a, b)

doy is the Day of Year (DOY)

evapo_l8t Actual evapotranspiration (ETa) using Landsat-8 (including thermal

evapo_l8t(doy, RG, Ta, ET0, a, b)

doy is the Day of Year (DOY)

evapo_l8t_grid Actual evapotranspiration (ETa) using Landsat-8 (including thermal

doy is the Day of Year (DOY)

evapo_l8_grid Actual evapotranspiration (ETa) using Landsat-8 images with a grid

evapo_modis Actual evapotranspiration (ETa) using MODIS with single agromete-

# dependencies of package 'agriwater'

# Using a temporary folder to run example

# creating raster which simulate Sentinel-2 reflectances - for using

mask <- as(extent(rast), 'SpatialPolygons')

# using "agriwater" - it's the same procedure as the used for

#Exiting temporary folder and returning to previous workspace

evapo_modis_grid Actual evapotranspiration (ETa) using MODIS with a grid of agrome-

evapo_s2 Actual evapotranspiration (ETa) using Sentinel-2 images with single

evapo_s2_grid Actual evapotranspiration (ETa) using Sentinel-2 images with a grid

doy is the Day of Year (DOY)

kc_l8t(doy, RG, Ta, a, b)

doy is the Day of Year (DOY)

# dependencies of package 'agriwater'

# Using a temporary folder to run example

# creating raster which simulate MODIS reflectances - for using

xy <- matrix(rnorm(100, mean = 0.07, sd = 0.01), 10, 10)

writeRaster(rast, file.path(wd, "B2"),format = "GTiff", overwrite=TRUE)

mask <- as(extent(rast), 'SpatialPolygons')

# using "agriwater" - it's the same procedure as the used for

#Exiting temporary folder and returning to previous workspace

radiation_l8 Energy balance using Landsat-8 images with single agrometeorologi-

Energy balance using Landsat-8 images with single agrometeorological data.

radiation_l8(doy, RG, Ta, ET0, a, b)

doy is the Day of Year (DOY)

radiation_l8t Energy balance using Landsat-8 images (including thermal bands)

radiation_l8t(doy, RG, Ta, ET0, a, b)

doy is the Day of Year (DOY)

radiation_l8t_grid Energy balance using Landsat-8 images (including thermal bands)

doy is the Day of Year (DOY)

radiation_l8_grid Energy balance using Landsat-8 images with a grid of agrometeoro-

radiation_modis Energy balance using Landsat-8 images with single agrometeorologi-

# dependencies of package 'agriwater'

# Using a temporary folder to run example

# creating raster which simulate Sentinel-2 reflectances - for using

# creating mask of study area

# using "agriwater" - it's the same procedure as the used for

#Exiting temporary folder and returning to previous workspace

radiation_modis_grid Energy balance using Landsat-8 images with a grid of agrometeoro-

Energy balance using Landsat-8 images with a grid of agrometeorological data.

doy is the Day of Year (DOY)

radiation_s2 Energy balance using Sentinel-2 images with single agrometeorologi-

Energy balance using Sentinel-2 images with single agrometeorological data.