Description Usage Arguments Details Value Author(s) See Also Examples
Compute basic statistics (min, max, mean and standard deviation) of NDVI maps between Ystart and Yend.
1 2 3 4 |
ndvidirectory |
path to the ndvi images folder (GeoTIFF files). |
region |
name of the region, same as in the name of the ndvi images. |
Ystart |
starting year of the analysis (in four-digit format). |
Yend |
ending year of the analysis (in four-digit format). |
xlim |
a two dimensions vector containing the longitude limits of the area of interest. |
ylim |
a two dimensions vector containing the latitude limits of the area of interest. |
type |
type of nomenclature followed by the NDVI images : "VITO_CLIP", "VITO_VGT" or "GIMMS" (see details section for more information). |
outname |
prefix of the name of the file where the maps containing the min, max, mean and sd are saved. |
outext |
extension of the created maps. |
org |
a two dimensions vector indicating the layout of the maps. |
shapefile |
name of the shape/kml file with objects (points or polygons) to be displayed on the maps. |
shapedir |
path to the shape/kml file. |
shapeext |
extension of the shape/kml file : "shp" for a Esri shapefile or "kml" for a kml file from Google earth . |
pal |
name of the palette of colors used to make the map from the package RColorBrewer. |
maplocalstat
load the maps between Ystart
and Yend
. Then it computes basics statistics (max, min, mean, sd) discarding NA values.
More about the variable type :
type = GIMMS : nomenclature from Global Land Cover Facility
[RE][YY][MMM][DDD].[SAT]-VIg.tif
RE : identifies the continent (AF for Africa,AZ for Australia and New Zealand, EA for Eurasia, NA for North America and SA for South America and Central America)
YY : the two last digits of the year
MMM : three-letter abbreviation of the month
DDD : identifies the composite period, 15a denotes the days 1-15 of the month, and 15b denotes the days from 16 to the end of the month.
SAT : identifies the satellite, NOAA-7, 9, 11, 14, 16, or 17 from which the data originated.
type = VITO_VGT : FreeVGT VITO's nomenclature
NDV_[YYYYMMDD]_[REGION]_Extract.tif
YYYY : identifies the year (four digits)
MM : identifies the months (two digits)
DD : identifies the date of the period : "01", "11" or "21" (10-day composite).
REGION : idendifies the region selected while extracting maps
type = VITO_CLIP : nomenclature after clipping the images in a particular region.
[REGION][YYYY]M[MM]P[D].tif
REGION : the region as entered while clipping data.
YYYY : identifies the year (four digits)
MM : identifies the months (two digits)
D : identifies the period (0, 1 or 2).
The function returns 4 maps containing the basic statistics :
mean |
The mean of the NDVI values between |
max |
The maximum of the NDVI values between |
min |
The minimum of the NDVI values between |
sd |
The standard deviation of the NDVI values between |
Romain Frelat and Bruno Gerard
mapmaxyear
,PeriodAnomaly
,MaxAnomaly
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 | #local variables
ndvidirectory=paste(system.file("extdata/VITO_Mzimba",
package="ndvits"), "/", sep="")
region="Mzimba"
Ystart=2004
Yend=2009
#executing maplocalstat
res=maplocalstat(ndvidirectory, region, Ystart, Yend,
outname="Mzimba", outext="show")
#adding points to the map
shape="SLP_Mzimba"
shapedir=paste(system.file("extdata/shape", package="ndvits"),
"/", sep="")
res=maplocalstat(ndvidirectory, region, Ystart, Yend,
outname="Mzimba", outext="show", shapefile=shape,
shapedir=shapedir)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.