getHdf: Create or update a local, user-defined subset of the global...
In aubreyd/modis4rwrfhydro: MODIS acquisition and processing

Description Usage Arguments Value Author(s) References Examples

Create or update a local user-defined subset of the global MODIS grid data archive (starting from version MODIS 0.2.4: ftp://e4ftl01.cr.usgs.gov/, until version MODIS 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 the local MODIS data archive.

1 2	getHdf(product, begin=NULL, end=NULL, tileH=NULL, tileV=NULL, extent=NULL, collection=NULL, HdfName, quiet=FALSE, wait=0.5, checkIntegrity=FALSE, ...)

`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: `begin,end,collection,tileV,...`. If many files have to be checked and downloaded it is simpler, safer and faster to use the 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 `M.D13Q1`. Run `getProduct()` for an overview of supported products.
`begin`	Character. (Ignored if `HdfName` is provided.) Subset begin date of MODIS time series. See: `?transDate` for formatting.
`end`	Character. (Ignored if `HdfName` is provided.) Subset end date of MODIS time series. See: `?transDate` for formatting. Can even be a future date, for continuous update in schedule jobs.
`tileH`	Character vector. (Ignored if `HdfName` is provided.) Horizontal tile number(s) (for MODIS 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) (for MODIS 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 or run `?getCollection` for more information.
`checkIntegrity`	Logical. Default `FALSE`, no file integrity check is performed. If `TRUE` the HDF files size is checked and the function tries to re-download the file. Currently it uses GDAL (gdalinfo), that supports HDF4, run `MODISoptions()` to check availability).
`wait`	Numeric. 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. (Hopefully a temporary solution!)
`quiet`	Logical. Refers to all internal `download.file()` actions. Default `FALSE`, if `TRUE` the `getHdf()` becomes really boring.
`...`	Arguments found in MODISoptions, sections 'STORAGE' and 'DOWNLOAD'

An invisible vector of downloaded data+path, if log=TRUE, a CSV file with reports

Matteo Mattiuzzi

MODIS data is obtained through the online Data Pool at the NASA Land Processes Distributed Active Archive Center (LP DAAC), USGS/Earth Resources Observation and Science (EROS) Center, Sioux Falls, South Dakota https://lpdaac.usgs.gov/get_data.

SRTM data is obtained through CGIAR-CSI and mirror servers. For the use this data please read http://srtm.csi.cgiar.org/SELECTION/SRT_disclaimer.htm. Jarvis A., H.I. Reuter, A. Nelson, E. Guevara, 2008, Hole-filled seamless SRTM data V4, International Centre for Tropical Agriculture (CIAT), available from http://srtm.csi.cgiar.org

## Not run: 
# one specific file (no regular erpression allowed here for now!)
# ca 2 mB!!!
a <- getHdf(HdfName="MYD11A1.A2009001.h18v04.005.2009013145459.hdf", wait=0)
a

# DO NOT RUN THE NEXT EXAMPLE (unless you need this data)!!!
# Get all MODIS TERRA and AQUA 11A1 beginning from 18. October 2011 up to today.
# (Can be ran in a sceduled job, for daily archive update)
# getHdf(product="M.D11A1",begin="2011.10.18",tileH=18:19,tileV=4) 

# same Tiles with a LIST extent.
# approximatley 21 mB!
Austria <- list(xmax=17.47,xmin=9.2,ymin=46.12,ymax=49.3)
b <- getHdf(product="M.D11A1",begin="2009001",end="2009-01-02",extent=Austria)
b

# require(mapdata)
getHdf(product="M.D11A1",begin="2009.01.01",end="2009002",extent="austria")
# without specification of the extent... interactive selection see: "getTile()"
c <- getHdf(product="M.D11A1",begin="2009.01.01",end="2009002")
c

# SRTM data
# SRTM server limits downloads! (It seams to be blocked, but (normally) it continues after a while,
# you have to be patient!).
# The files are zips, this function only performs the download!
d <- getHdf(product="SRTM",extent="austria")  
d

## End(Not run)