arcStats: Get Summary of Local MODIS Data

In MatMatt/MODIS: Acquisition and Processing of MODIS Products

Get Summary of Local MODIS Data

Description

In the same manner as getHdf(), this function quantifies the availability of local MODIS .hdf data and gives you a visual and/or tabular overview.

Usage

arcStats(
  product,
  collection = NULL,
  extent = "global",
  begin = "2000.01.01",
  end = format(Sys.Date(), "%Y.%m.%d"),
  asMap = TRUE,
  outName = NULL,
  ...
)

Arguments

product	character, see getProduct(). MODIS grid product to be checked.
collection	character or integer, see getCollection(). MODIS product version.
extent	Extent information, defaults to "global". See getTile().
begin	character. Begin date of MODIS time series, see transDate().
end	character. End date, defaults to Sys.Date().
asMap	Controls output type. Possible options are TRUE (.png), FALSE (.csv) or "both".
outName	character. Name of output file, defaults to "product.collection.YYYYMMDDHHMMSS.png" (or .csv) of the function call or, if applicable, "product.collection.extent.YYYYMMDDHHMMSS.png" (or .csv).
...	Arguments passed to MODISoptions(), most importantly 'outProj' and 'outDirPath'.

Value

An invisible NULL (probably this will change to a matrix-like object similar to the .csv output). If asMap = TRUE, a "table.csv" and a "image.png" file in 'outDirPath'.

Author(s)

Matteo Mattiuzzi

Examples

## Not run: 
# The following examples expect that you have some data stored locally!
########################################################### 
# generates 2 png's and 2 csv's, one for TERRA one for AQUA
arcStats(product="M.D13Q1")

# -"- with the specified countries
arcStats(product="M.D13Q1",extent=c("austria","germany","italy"))

# generates 1 png and 1 csv for AQUA
arcStats(product="MYD13Q1",begin="2005001",outName="MyDataStart2005")

# generates 1 png for AQUA for the selected area and plots it in 'Sinusoidal'
arcStats(product="MYD13Q1",begin="2005001",asMap=TRUE, outName="InteractiveSelection2005",
         extent=getTile(), outProj="asIn")

# -"- and plots it in 'Geographic' Coordinates.
arcStats(product="MYD13Q1",begin="2005001",asMap=TRUE, outName="InteractiveSelection2005",
         extent=getTile(), outProj="GEOGRAPHIC")

## End(Not run)

MatMatt/MODIS documentation built on Feb. 1, 2023, 12:39 a.m.