1 |
7526fb1c
|
Jim Regetz
|
# Mosaics two tiles of daily average clear day coverage for the clear day
|
2 |
|
|
# MODIS layer extracted from the MODIS LST files
|
3 |
|
|
#
|
4 |
|
|
# Reads a directory of .img files and mosaics the h08 and h09 tiles together
|
5 |
|
|
#
|
6 |
|
|
# Developed by John Donoghue
|
7 |
|
|
# Created: 20 October 2010
|
8 |
|
|
# Last Updated:
|
9 |
|
|
#
|
10 |
|
|
# For NCEAS Working Group Environment and Organisms
|
11 |
|
|
#
|
12 |
|
|
|
13 |
|
|
#setwd("/data/project/organisms/R")
|
14 |
|
|
#source("Mosaic Clear Day Average Grids.r", echo=TRUE, print.eval=TRUE)
|
15 |
|
|
|
16 |
|
|
setwd("/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL")
|
17 |
|
|
library(raster)
|
18 |
|
|
|
19 |
|
|
for(i in 1:366){
|
20 |
|
|
# get the two tiles that comprise each day of the year
|
21 |
|
|
myPattern <- paste("Day_", toString(i), "_.*img$", sep="")
|
22 |
|
|
print(myPattern)
|
23 |
|
|
inFiles <- list.files(pattern=myPattern)
|
24 |
|
|
nFiles <- length(inFiles)
|
25 |
|
|
print (nFiles)
|
26 |
|
|
|
27 |
|
|
# should only be two tiles at this point
|
28 |
|
|
h08Raster <- raster()
|
29 |
|
|
h08Raster <- raster(inFiles[1])
|
30 |
|
|
print(sprintf("Reading file: %s", inFiles[1]))
|
31 |
|
|
h09Raster <- raster()
|
32 |
|
|
h09Raster <- raster(inFiles[2])
|
33 |
|
|
print(sprintf("Reading file: %s", inFiles[2]))
|
34 |
|
|
|
35 |
|
|
newRaster <- merge(h08Raster, h09Raster)
|
36 |
|
|
|
37 |
|
|
# Check if file exists
|
38 |
|
|
outFile <- paste("Day_",toString(i),"_ClearDay_Average",".img",sep="")
|
39 |
|
|
if (file.exists(outFile) == FALSE) {
|
40 |
|
|
print(sprintf("...Writing file: %s",outFile))
|
41 |
|
|
writeRaster(newRaster,filename=outFile,format="HFA",overwrite=TRUE)
|
42 |
|
|
}
|
43 |
|
|
}
|
44 |
|
|
|
45 |
|
|
# finished
|
46 |
|
|
print("Finished")
|