getTILE: Get MODIS tile id(s)

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

Get MODIS tile id(s) for a geographic area.

Usage

getTILE(tileH, tileV, extent)

Arguments

tileH	Integer. Horizontal tile number(s) (between 0 and 35). I.e.: tileH = 19:25
tileV	Integer. Vertical tile number(s) (between 0 and 17). I.e.: tileV = 9:12.
extent	Extent Informations see details.

Details

extent:
If 'list': myExtent <- list(lat_min=mumeric, lat_max=mumeric, lon_max=mumeric, lon_min=mumeric).
If 'character', the country name from mapdata.
Or one of the following objects 'extent','rasterLayer','rasterStack' or 'rasterBrick' from package(raster).

Value

List

Author(s)

Matteo Mattiuzzi

See Also

getHDF, raster package: ?extent,?raster,?stack or ?brick

Examples

# with tileH and tileV
getTILE(tileH=18:19,tileV=4)

# with Extent of class 'list'
Austria <- list(lat_min=46.12,lat_max=49.3,lon_min=9.2,lon_max=17.47)
getTILE(extent=Austria)

# with Extent or raster* object from raster package
require(raster)

# rasterObject
rasterObject <- raster(xmn=9.2,xmx=17.47,ymn=46.12,ymx=49.3)
getTILE(extent=rasterObject)

# raster extent
ext <- extent(rasterObject)
getTILE(extent=ext)

# require(mapdata)
# if you pass a character name to the extent option a standard call to 'map('worldHires',extent,plot=FALSE)' is performed.
# Use with caution! The following example includes much more areas than the main USA. Look at "map('worldHires','usa')" result.  
getTILE(extent="usa")
# so you better always check befor applying an extent of the class 'character' that defines a 'map' object
# if you want subsets or more detailed maps use ?map arguments 
ext <-  map('state', region = c('new york', 'new jersey', 'penn'))
getTILE(extent=ext)

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