1 |
7526fb1c
|
Jim Regetz
|
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")
|