getHDF: Create or update a local, user-defined subset of the global...

In getMODIS: Utility package for creating or updating a local user-defined subsets of the global MODIS grid data online archive. Pre-processing functionalities for MODIS grid data using MODIS Reporjection Tool (MRT)

Description

Create or update a local user-defined subset of the global MODIS grid data archive (starting from version 0.2.4: ftp://e4ftl01.cr.usgs.gov/,until version 0.2-3: ftp://e4ftl01u.ecs.nasa.gov/). Based on user-specific parameters the function checks in the local archive for available data and downloads missing data from the online MODIS data pool. When run in a schedule job, the function manage the continuous update of local MODIS data archive.

Usage

getHDF(LocalArcPath, HdfName, product, startdate, enddate, tileH, tileV, extent, collection, dlmethod="auto", quiet=FALSE, wait=1, checkSize=FALSE, log=TRUE)

Arguments

LocalArcPath	Character. Local path to look for and to download MODIS files. If missing the default archive '~/MODIS_ARC/' (LINUX), 'C:/Documents and Settings/username/My Documents/MODIS_ARC/' (Windows XP) and 'C:/Users/username/Documents/MODIS_ARC' (Vista/Windows 7) is used.
HdfName	Character (vector or list). Full HDF file name(s) to download. If just few single (AND known!) files should be downloaded use HdfName. If HdfName is provided the other file-related parameters are ignored: startdate, enddate, collection, tileV,.... If many files have to be checked and downloaded it is simpler, safer and faster to use following parameters.
product	Character. (Ignored if HdfName is provided.) MODIS grid product to be downloaded. For 'MODIS TERRA Vegetation Indices at 250 m ground resolution' use: MOD13Q1. For vegetation indices from TERRA and AQUA use MxD13Q1. See https://lpdaac.usgs.gov/lpdaac/products/modis_products_table for more Products.
startdate	Character. (Ignored if HdfName is provided.) Subset start date of MODIS time series (format: "YYYY.MM.DD").
enddate	Character. (Ignored if HdfName is provided.) Subset end date of MODIS time series (format: "YYYY.MM.DD"). Can be a future date, for continuous update in a schedule job.
tileH	Character vector. (Ignored if HdfName is provided.) Horizontal tile number(s) (between 0 and 35). I.e.: tileH = 1:5. Refer to http://modis-land.gsfc.nasa.gov/MODLAND_grid.htm to find your area.
tileV	Character vector. (Ignored if HdfName is provided.) Vertical tile number(s) (between 0 and 17). I.e.: tileV = 9:17.
extent	see 'extent' in ?getTILE.
collection	Character or integer. (Ignored if HdfName is provided.) Refers to the MODIS product collection-version. Often '005', but see MODIS pages for more information.
dlmethod	Default auto. Argument passed to ?download.file
checkSize	Logical. Default FALSE, MODIS-XML files are not downloaded with their MODIS-HDF file and a "checkSum" is not done. If TRUE MODIS-HDF and corresponding MODIS-XML files are checked and downloaded and the filesize is checked. (!keep FALSE for now!)
wait	Numeric. Default 1. Inserts a break (in seconds) after every internal download.file() or getURL() call. It reduces the chance of FTP connection error that often occurs after many requests. Requires the 'audio' package. (Hopefully a temporary solution!)
quiet	Logical. Refers to all internal download.file() actions. Default FALSE, if TRUE the getHDF becomes really boring.
log	Logical. Default TRUE, a log file is generated in LocalArcPath

Value

NONE, a CSV file with reports, and the HDF archive

Author(s)

Matteo Mattiuzzi

See Also

getXML

Examples

## Not run: 
# one spezific file (no regular erpression allowed here for now!)
getHDF(HdfName="MOD13Q1.A2009209.h19v04.005.2009227222141.hdf", wait=0)

# DO NOT RUN THE NEXT EXAMPLE (unless you need this data)!!!
# Get all MODIS TERRA 13Q1 beginning from 01. January 2009 up to today.
# If run in a scheduled batch job it will update the data until end of 2015.
getHDF(product="MOD13Q1",startdate="2009.01.01",enddate="2015.12.31",tileH=19,tileV=4:5) 

# or the same with an (list) Extent
Austria <- list(lat_min=46.12,lat_max=49.3,lon_max=17.47,lon_min=9.2)
getHDF(product="MYD13Q1",startdate="2009.01.01",enddate="2015.12.31",extent=Austria)

# require(mapdata)
getHDF(product="MYD13Q1",startdate="2009.01.01",enddate="2015.12.31",extent="austria")

## End(Not run)

getMODIS documentation built on May 2, 2019, 4:51 p.m.