|
1
|
### script to return bounding box for a list of tiles
|
|
2
|
|
|
3
|
## import commandline arguments
|
|
4
|
library(getopt)
|
|
5
|
## get options
|
|
6
|
opta <- getopt(matrix(c(
|
|
7
|
'tile', 't', 1, 'character'
|
|
8
|
), ncol=4, byrow=TRUE))
|
|
9
|
if(is.null(opta$tile)) stop("Please provide a list of tiles. For example, -t \"h11v08,h12v09\"")
|
|
10
|
|
|
11
|
tiles=sub(" ","",tolower(do.call(c,strsplit(opta$tile,split=","))))
|
|
12
|
|
|
13
|
### get tile boundaries
|
|
14
|
tb=read.table("http://landweb.nascom.nasa.gov/developers/sn_tiles/sn_bound_10deg.txt",skip=6,nrows=648,header=T)
|
|
15
|
tb$tile=paste("h",sprintf("%02d",tb$ih),"v",sprintf("%02d",tb$iv),sep="")
|
|
16
|
|
|
17
|
|
|
18
|
## get minmax for all tiles
|
|
19
|
tb2=tb[tb$tile%in%tiles,]
|
|
20
|
|
|
21
|
print(tb2)
|
|
22
|
|
|
23
|
## print summary
|
|
24
|
print(paste("Getting bounding box for tiles: ",paste(tb2$tile,collapse=", ")))
|
|
25
|
print(paste("lon[",min(tb2$lon_min),",",max(tb2$lon_max),"] lat[",min(tb2$lat_min),",",max(tb2$lat_max),"]"))
|
|
26
|
|
|
27
|
|
