Get MODIS tile id(s)

Description

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

Usage

1
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
# 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)