1
|
library (raster)
|
2
|
|
3
|
MA_path<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles"
|
4
|
img_files<-list.files(path=MA_path, pattern="*.img$")
|
5
|
|
6
|
#Create empty vectors to hold correct tiles to be mosaiced
|
7
|
JanTiles_MonthlyAvgs<- c(1:length(img_files))*NA
|
8
|
FebTiles_MonthlyAvgs<- c(1:length(img_files))*NA
|
9
|
MarTiles_MonthlyAvgs<- c(1:length(img_files))*NA
|
10
|
AprTiles_MonthlyAvgs<- c(1:length(img_files))*NA
|
11
|
MayTiles_MonthlyAvgs<- c(1:length(img_files))*NA
|
12
|
JunTiles_MonthlyAvgs<- c(1:length(img_files))*NA
|
13
|
JulTiles_MonthlyAvgs<- c(1:length(img_files))*NA
|
14
|
AugTiles_MonthlyAvgs<- c(1:length(img_files))*NA
|
15
|
SepTiles_MonthlyAvgs<- c(1:length(img_files))*NA
|
16
|
OctTiles_MonthlyAvgs<- c(1:length(img_files))*NA
|
17
|
NovTiles_MonthlyAvgs<- c(1:length(img_files))*NA
|
18
|
DecTiles_MonthlyAvgs<- c(1:length(img_files))*NA
|
19
|
|
20
|
#Separate pairs of tiles (h08v04 and h09v04) into vectors
|
21
|
for (i in 1:length(img_files)){
|
22
|
if (length(grep("^Jan.h.",img_files[i]) !=0)){
|
23
|
img_files[i]-> JanTiles_MonthlyAvgs[i]
|
24
|
}else if (length(grep("^Feb.h.",img_files[i]) !=0)){
|
25
|
img_files[i]-> FebTiles_MonthlyAvgs[i]
|
26
|
}else if (length(grep("^Mar.h.",img_files[i]) !=0)){
|
27
|
img_files[i]-> MarTiles_MonthlyAvgs[i]
|
28
|
}else if (length(grep("^Apr.h.",img_files[i]) !=0)){
|
29
|
img_files[i]-> AprTiles_MonthlyAvgs[i]
|
30
|
}else if (length(grep("^May.h.",img_files[i]) !=0)){
|
31
|
img_files[i]-> MayTiles_MonthlyAvgs[i]
|
32
|
}else if (length(grep("^Jun.h.",img_files[i]) !=0)){
|
33
|
img_files[i]-> JunTiles_MonthlyAvgs[i]
|
34
|
}else if (length(grep("^Jul.h.",img_files[i]) !=0)){
|
35
|
img_files[i]-> JulTiles_MonthlyAvgs[i]
|
36
|
}else if (length(grep("^Aug.h.",img_files[i]) !=0)){
|
37
|
img_files[i]-> AugTiles_MonthlyAvgs[i]
|
38
|
}else if (length(grep("^Sep.h.",img_files[i]) !=0)){
|
39
|
img_files[i]-> SepTiles_MonthlyAvgs[i]
|
40
|
}else if (length(grep("^Oct.h.",img_files[i]) !=0)){
|
41
|
img_files[i]-> OctTiles_MonthlyAvgs[i]
|
42
|
}else if (length(grep("^Nov.h.",img_files[i]) !=0)){
|
43
|
img_files[i]-> NovTiles_MonthlyAvgs[i]
|
44
|
}else if (length(grep("^Dec.h.",img_files[i]) !=0)){
|
45
|
img_files[i]-> DecTiles_MonthlyAvgs[i]
|
46
|
}
|
47
|
}
|
48
|
|
49
|
#Get rid of NA's and add full path name to each file
|
50
|
JanTiles_MonthlyAvgs<-JanTiles_MonthlyAvgs[!is.na(JanTiles_MonthlyAvgs)]
|
51
|
JanTiles_MonthlyAvgs<- paste("/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/",JanTiles_MonthlyAvgs,sep="")
|
52
|
FebTiles_MonthlyAvgs<-FebTiles_MonthlyAvgs[!is.na(FebTiles_MonthlyAvgs)]
|
53
|
FebTiles_MonthlyAvgs<- paste("/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/",FebTiles_MonthlyAvgs,sep="")
|
54
|
MarTiles_MonthlyAvgs<-MarTiles_MonthlyAvgs[!is.na(MarTiles_MonthlyAvgs)]
|
55
|
MarTiles_MonthlyAvgs<- paste("/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/",MarTiles_MonthlyAvgs,sep="")
|
56
|
AprTiles_MonthlyAvgs<-AprTiles_MonthlyAvgs[!is.na(AprTiles_MonthlyAvgs)]
|
57
|
AprTiles_MonthlyAvgs<- paste("/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/",AprTiles_MonthlyAvgs,sep="")
|
58
|
MayTiles_MonthlyAvgs<-MayTiles_MonthlyAvgs[!is.na(MayTiles_MonthlyAvgs)]
|
59
|
MayTiles_MonthlyAvgs<- paste("/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/",MayTiles_MonthlyAvgs,sep="")
|
60
|
JunTiles_MonthlyAvgs<-JunTiles_MonthlyAvgs[!is.na(JunTiles_MonthlyAvgs)]
|
61
|
JunTiles_MonthlyAvgs<- paste("/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/",JunTiles_MonthlyAvgs,sep="")
|
62
|
JulTiles_MonthlyAvgs<-JulTiles_MonthlyAvgs[!is.na(JulTiles_MonthlyAvgs)]
|
63
|
JulTiles_MonthlyAvgs<- paste("/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/",JulTiles_MonthlyAvgs,sep="")
|
64
|
AugTiles_MonthlyAvgs<-AugTiles_MonthlyAvgs[!is.na(AugTiles_MonthlyAvgs)]
|
65
|
AugTiles_MonthlyAvgs<- paste("/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/",AugTiles_MonthlyAvgs,sep="")
|
66
|
SepTiles_MonthlyAvgs<-SepTiles_MonthlyAvgs[!is.na(SepTiles_MonthlyAvgs)]
|
67
|
SepTiles_MonthlyAvgs<- paste("/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/",SepTiles_MonthlyAvgs,sep="")
|
68
|
OctTiles_MonthlyAvgs<-OctTiles_MonthlyAvgs[!is.na(OctTiles_MonthlyAvgs)]
|
69
|
OctTiles_MonthlyAvgs<- paste("/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/",OctTiles_MonthlyAvgs,sep="")
|
70
|
NovTiles_MonthlyAvgs<-NovTiles_MonthlyAvgs[!is.na(NovTiles_MonthlyAvgs)]
|
71
|
NovTiles_MonthlyAvgs<- paste("/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/",NovTiles_MonthlyAvgs,sep="")
|
72
|
DecTiles_MonthlyAvgs<-DecTiles_MonthlyAvgs[!is.na(DecTiles_MonthlyAvgs)]
|
73
|
DecTiles_MonthlyAvgs<- paste("/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/",DecTiles_MonthlyAvgs,sep="")
|
74
|
|
75
|
#Create mosaics
|
76
|
Jan_name<-"/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Mosaics/JanMosaic_MonthlyAvg_Scaled.img"
|
77
|
mosaic(raster(JanTiles_MonthlyAvgs[1]),raster(JanTiles_MonthlyAvgs[2]), fun=mean,na.rm=TRUE, filename=Jan_name, datatype="HFA")
|
78
|
Feb_name<-"/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Mosaics/FebMosaic_MonthlyAvg_Scaled.img"
|
79
|
mosaic(raster(FebTiles_MonthlyAvgs[1]),raster(FebTiles_MonthlyAvgs[2]), fun=mean,na.rm=TRUE, filename=Feb_name, datatype="HFA")
|
80
|
Mar_name<-"/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Mosaics/MarMosaic_MonthlyAvg_Scaled.img"
|
81
|
mosaic(raster(MarTiles_MonthlyAvgs[1]),raster(MarTiles_MonthlyAvgs[2]), fun=mean,na.rm=TRUE, filename=Mar_name, datatype="HFA")
|
82
|
Apr_name<-"/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Mosaics/AprMosaic_MonthlyAvg_Scaled.img"
|
83
|
mosaic(raster(AprTiles_MonthlyAvgs[1]),raster(AprTiles_MonthlyAvgs[2]), fun=mean,na.rm=TRUE, filename=Apr_name, datatype="HFA")
|
84
|
May_name<-"/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Mosaics/MayMosaic_MonthlyAvg_Scaled.img"
|
85
|
mosaic(raster(MayTiles_MonthlyAvgs[1]),raster(MayTiles_MonthlyAvgs[2]), fun=mean,na.rm=TRUE, filename=May_name, datatype="HFA")
|
86
|
Jun_name<-"/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Mosaics/JunMosaic_MonthlyAvg_Scaled.img"
|
87
|
mosaic(raster(JunTiles_MonthlyAvgs[1]),raster(JunTiles_MonthlyAvgs[2]), fun=mean,na.rm=TRUE, filename=Jun_name, datatype="HFA")
|
88
|
Jul_name<-"/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Mosaics/JulMosaic_MonthlyAvg_Scaled.img"
|
89
|
mosaic(raster(JulTiles_MonthlyAvgs[1]),raster(JulTiles_MonthlyAvgs[2]), fun=mean,na.rm=TRUE, filename=Jul_name, datatype="HFA")
|
90
|
Aug_name<-"/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Mosaics/AugMosaic_MonthlyAvg_Scaled.img"
|
91
|
mosaic(raster(AugTiles_MonthlyAvgs[1]),raster(AugTiles_MonthlyAvgs[2]), fun=mean,na.rm=TRUE, filename=Aug_name, datatype="HFA")
|
92
|
Sep_name<-"/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Mosaics/SepMosaic_MonthlyAvg_Scaled.img"
|
93
|
mosaic(raster(SepTiles_MonthlyAvgs[1]),raster(SepTiles_MonthlyAvgs[2]), fun=mean,na.rm=TRUE, filename=Sep_name, datatype="HFA")
|
94
|
Oct_name<-"/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Mosaics/OctMosaic_MonthlyAvg_Scaled.img"
|
95
|
mosaic(raster(OctTiles_MonthlyAvgs[1]),raster(OctTiles_MonthlyAvgs[2]), fun=mean,na.rm=TRUE, filename=Oct_name, datatype="HFA")
|
96
|
Nov_name<-"/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Mosaics/NovMosaic_MonthlyAvg_Scaled.img"
|
97
|
mosaic(raster(NovTiles_MonthlyAvgs[1]),raster(NovTiles_MonthlyAvgs[2]), fun=mean,na.rm=TRUE, filename=Nov_name, datatype="HFA")
|
98
|
Dec_name<-"/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Mosaics/DecMosaic_MonthlyAvg_Scaled.img"
|
99
|
mosaic(raster(DecTiles_MonthlyAvgs[1]),raster(DecTiles_MonthlyAvgs[2]), fun=mean,na.rm=TRUE, filename=Dec_name, datatype="HFA")
|