Download MODIS images

-----Original Message-----
From: r-sig-geo-bounces at stat.math.ethz.ch [mailto:r-sig-geo-
bounces at stat.math.ethz.ch] On Behalf Of Rosca Bogdan
Sent: Saturday, July 03, 2010 7:04 PM
To: r-sig-geo at stat.math.ethz.ch
Subject: [R-sig-Geo] Download MODIS images

Hi!
I am trying to work with T. Hengl's script downloaded from:
http://spatial-analyst.net/book/Rscripts in order to get modis images for my
country. I have only modify the script in order to sample data according to
my country border (Romania, which is not so far away from Croatia) and
everything run smooth until these lines:

# download the blocks in a loop (year 2006):
for (i in 268:313){
getlist <- strsplit(getURL(paste(MOD11A2, dates$dirname[[i]], "/", sep=""),
.opts=curlOptions(ftplistonly=
TRUE)), "\r\n")[[1]]
BLOCK1 <- getlist[grep(getlist,
pattern="MOD11A2.********.h18v04.*************.hdf")[1]]
BLOCK2 <- getlist[grep(getlist,
pattern="MOD11A2.********.h19v04.*************.hdf")[1]]

# write up the file names back to the dates.txt:
for(j in 2:3){
  dates[i,j] <- get(paste("BLOCK", j-1, sep=""))
}

# Download all blocks from the list to a local drive:
# while(!is.na(dates[i,2])&!is.na(dates[i,3])&!is.na(dates[i,4])&!is.na
(dates[i,5])&!is.na(dates[i,6])&!is.na(dates[i,7])&!is.na(dates[i,8])&!is.na
(dates[i,9])&!is.na(dates[i,10])){
download.file(paste(MOD11A2a, dates$dirname[[i]], "/", BLOCK1,sep=""),
destfile=paste(getwd(), "/", BLOCK1, sep=""), mode='wb', method='wget',
quiet=T, cacheOK=FALSE)
download.file(paste(MOD11A2, dates$dirname[[i]], "/", BLOCK2,sep=""),
destfile=paste(getwd(), "/", BLOCK2,sep=""), mode='wb', method='wget',
quiet=T, cacheOK=FALSE)

# remove "." from the file name:
dirname1 <- sub(sub(pattern="\\.", replacement="_", dates$dirname[[i]]),
pattern="\\.", replacement="_", dates$dirname[[i]])
# mosaic the blocks:
mosaicname <- file(paste(MRT, "TmpMosaic.prm", sep=""), open="wt")
write(paste(workd, BLOCK1, sep=""), mosaicname)
write(paste(workd, BLOCK2, sep=""), mosaicname, append=T)
close(mosaicname)
# generate temporary mosaic:
shell(cmd=paste(MRT, 'mrtmosaic -i ', MRT, 'TmpMosaic.prm -s "0 0 0 0 1 0 0
0 0 0 0 0" -o ', workd, 'TmpMosaic.hdf', sep=""))
close(mosaicname)

# resample to UTM 33N:
filename <- file(paste(MRT, "mrt", dirname1, ".prm", sep=""), open="wt")
write(paste('INPUT_FILENAME = ', workd, 'TmpMosaic.hdf', sep=""), filename)
write('  ', filename, append=TRUE)
# write('SPECTRAL_SUBSET = ( 0 1 0 0 0 0 0 0 0 0 0 )', filename,
append=TRUE)
write('SPECTRAL_SUBSET = ( 1 )', filename, append=TRUE)
write('  ', filename, append=TRUE)
write('SPATIAL_SUBSET_TYPE = OUTPUT_PROJ_COORDS', filename, append=TRUE)
write('  ', filename, append=TRUE)
write(paste('SPATIAL_SUBSET_UL_CORNER = (', XUL, YUL, ')'), filename,
append=TRUE)
write(paste('SPATIAL_SUBSET_LR_CORNER = (', XLR, YLR, ')'), filename,
append=TRUE)
write('  ', filename, append=TRUE)
write(paste('OUTPUT_FILENAME = ', workd, 'LST', dirname1, '.tif', sep=""),
filename, append=TRUE)
write('  ', filename, append=TRUE)
write('RESAMPLING_TYPE = NEAREST_NEIGHBOR', filename, append=TRUE)
write('  ', filename, append=TRUE)
write('OUTPUT_PROJECTION_TYPE = UTM', filename, append=TRUE)
write('  ', filename, append=TRUE)
write('OUTPUT_PROJECTION_PARAMETERS = ( ', filename, append=TRUE)
write(paste(lon.c, lat.c, '0.0'), filename, append=TRUE)
write(' 0 0.0 0.0', filename, append=TRUE)
write(' 0.0 0.0 0.0', filename, append=TRUE)
write(' 0.0 0.0 0.0', filename, append=TRUE)
write(' 0.0 0.0 0.0 )', filename, append=TRUE)
write('  ', filename, append=TRUE)
write('UTM_ZONE = 35', filename, append=TRUE)
write('DATUM = WGS84', filename, append=TRUE)
write('  ', filename, append=TRUE)
write('OUTPUT_PIXEL_SIZE = 1000', filename, append=TRUE)
write('  ', filename, append=TRUE)
close(filename)

# Mosaic the images to get the whole area:
shell(cmd=paste(MRT, 'resample -p ', MRT, 'mrt', dirname1, '.prm', sep=""))
# delete all hdf files!
# unlink(paste(getwd(), '/', BLOCK1, sep=""))
# unlink(paste(getwd(), '/', BLOCK2, sep=""))
}

I keep getting an error like this:
Error in grep(getlist, pattern =
"MOD11A2.********.h18v04.*************.hdf") :
 invalid regular expression 'MOD11A2.********.h18v04.*************.hdf'
In addition: Warning message:
In grep(getlist, pattern = "MOD11A2.********.h18v04.*************.hdf") :
 regcomp error:  'Invalid use of repetition operators'

BLOCK1 <- getlist[grep(getlist, pattern="MOD11A2.*.h18v04.*.hdf")[1]]

Please help.
Thank's,
Rosca Bogdan

	[[alternative HTML version deleted]]