1
|
library (raster)
|
2
|
|
3
|
path<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/ClearDay_Original_IMG_Extracts/ByDate/BadPixelsCorrected_Tiles"
|
4
|
|
5
|
#Find monthly means for h09v04 tiles:
|
6
|
h09v04_inFiles <- list.files(path=path, pattern="*.h09v04.005.ClearDay_corrected.img$")
|
7
|
for (i in 1: length(h09v04_inFiles)){
|
8
|
h09v04_inFiles[i]<- paste (path,"/",h09v04_inFiles[i],sep="")
|
9
|
}
|
10
|
h09v04_inFiles
|
11
|
|
12
|
#Separate out monthly groupings
|
13
|
Jan_h09v04<- c(1:length (h09v04_inFiles))*NA
|
14
|
for (i in 1: length(h09v04_inFiles)){
|
15
|
if (length(grep("\\.1\\.1",h09v04_inFiles[i]) !=0)||
|
16
|
length(grep("\\.1\\.2",h09v04_inFiles[i]) !=0)||
|
17
|
length(grep("\\.1\\.3",h09v04_inFiles[i]) !=0)||
|
18
|
length(grep("\\.1\\.4",h09v04_inFiles[i]) !=0)||
|
19
|
length(grep("\\.1\\.5",h09v04_inFiles[i]) !=0)||
|
20
|
length(grep("\\.1\\.6",h09v04_inFiles[i]) !=0)||
|
21
|
length(grep("\\.1\\.7",h09v04_inFiles[i]) !=0)||
|
22
|
length(grep("\\.1\\.8",h09v04_inFiles[i]) !=0)||
|
23
|
length(grep("\\.1\\.9",h09v04_inFiles[i]) !=0)){
|
24
|
Jan_h09v04[i]<- h09v04_inFiles[i]
|
25
|
}
|
26
|
}
|
27
|
Jan_h09v04<-Jan_h09v04[!is.na(Jan_h09v04)]
|
28
|
Jan_h09v04
|
29
|
|
30
|
#Create rasters stack, calc means, scale
|
31
|
Jan_h09v04_stack<- stack(Jan_h09v04) #checked nlayers(), 310: correct
|
32
|
Jan_h09v04_mn<- mean(Jan_h09v04_stack, na.rm=TRUE)
|
33
|
Jan_h09v04_mnScaled<- Jan_h09v04_mn*0.0005
|
34
|
|
35
|
#Write results to new file
|
36
|
Out_Jan<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/Jan_h09v04_AVG_Scaled.img"
|
37
|
writeRaster (Jan_h09v04_mnScaled, filename=Out_Jan,format="HFA",overwrite=TRUE)
|
38
|
|
39
|
Feb_h09v04<- c(1:length (h09v04_inFiles))*NA
|
40
|
for (i in 1: length(h09v04_inFiles)){
|
41
|
if (length(grep("\\.2\\.1",h09v04_inFiles[i]) !=0)||
|
42
|
length(grep("\\.2\\.2",h09v04_inFiles[i]) !=0)||
|
43
|
length(grep("\\.2\\.3",h09v04_inFiles[i]) !=0)||
|
44
|
length(grep("\\.2\\.4",h09v04_inFiles[i]) !=0)||
|
45
|
length(grep("\\.2\\.5",h09v04_inFiles[i]) !=0)||
|
46
|
length(grep("\\.2\\.6",h09v04_inFiles[i]) !=0)||
|
47
|
length(grep("\\.2\\.7",h09v04_inFiles[i]) !=0)||
|
48
|
length(grep("\\.2\\.8",h09v04_inFiles[i]) !=0)||
|
49
|
length(grep("\\.2\\.9",h09v04_inFiles[i]) !=0)){
|
50
|
Feb_h09v04[i]<- h09v04_inFiles[i]
|
51
|
}
|
52
|
}
|
53
|
Feb_h09v04<-Feb_h09v04[!is.na(Feb_h09v04)]
|
54
|
Feb_h09v04
|
55
|
|
56
|
Feb_h09v04_stack<- stack(Feb_h09v04)
|
57
|
nlayers(Feb_h09v04_stack)
|
58
|
Feb_h09v04_mn<- mean(Feb_h09v04_stack, na.rm=TRUE)
|
59
|
Feb_h09v04_mnScaled<- Feb_h09v04_mn*0.0005
|
60
|
|
61
|
Out_Feb<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/Feb_h09v04_AVG_Scaled.img"
|
62
|
writeRaster (Feb_h09v04_mnScaled, filename=Out_Feb,format="HFA",overwrite=TRUE)
|
63
|
|
64
|
Mar_h09v04<- c(1:length (h09v04_inFiles))*NA
|
65
|
for (i in 1: length(h09v04_inFiles)){
|
66
|
if (length(grep("\\.3\\.1",h09v04_inFiles[i]) !=0)||
|
67
|
length(grep("\\.3\\.2",h09v04_inFiles[i]) !=0)||
|
68
|
length(grep("\\.3\\.3",h09v04_inFiles[i]) !=0)||
|
69
|
length(grep("\\.3\\.4",h09v04_inFiles[i]) !=0)||
|
70
|
length(grep("\\.3\\.5",h09v04_inFiles[i]) !=0)||
|
71
|
length(grep("\\.3\\.6",h09v04_inFiles[i]) !=0)||
|
72
|
length(grep("\\.3\\.7",h09v04_inFiles[i]) !=0)||
|
73
|
length(grep("\\.3\\.8",h09v04_inFiles[i]) !=0)||
|
74
|
length(grep("\\.3\\.9",h09v04_inFiles[i]) !=0)){
|
75
|
Mar_h09v04[i]<- h09v04_inFiles[i]
|
76
|
}
|
77
|
}
|
78
|
Mar_h09v04<-Mar_h09v04[!is.na(Mar_h09v04)]
|
79
|
Mar_h09v04
|
80
|
|
81
|
Mar_h09v04_stack<- stack(Mar_h09v04)
|
82
|
nlayers(Mar_h09v04_stack)
|
83
|
Mar_h09v04_mn<- mean(Mar_h09v04_stack, na.rm=TRUE)
|
84
|
Mar_h09v04_mnScaled<- Mar_h09v04_mn*0.0005
|
85
|
|
86
|
Out_Mar<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/Mar_h09v04_AVG_Scaled.img"
|
87
|
writeRaster (Mar_h09v04_mnScaled, filename=Out_Mar,format="HFA",overwrite=TRUE)
|
88
|
|
89
|
Apr_h09v04<- c(1:length (h09v04_inFiles))*NA
|
90
|
for (i in 1: length(h09v04_inFiles)){
|
91
|
if (length(grep("\\.4\\.1",h09v04_inFiles[i]) !=0)||
|
92
|
length(grep("\\.4\\.2",h09v04_inFiles[i]) !=0)||
|
93
|
length(grep("\\.4\\.3",h09v04_inFiles[i]) !=0)||
|
94
|
length(grep("\\.4\\.4",h09v04_inFiles[i]) !=0)||
|
95
|
length(grep("\\.4\\.5",h09v04_inFiles[i]) !=0)||
|
96
|
length(grep("\\.4\\.6",h09v04_inFiles[i]) !=0)||
|
97
|
length(grep("\\.4\\.7",h09v04_inFiles[i]) !=0)||
|
98
|
length(grep("\\.4\\.8",h09v04_inFiles[i]) !=0)||
|
99
|
length(grep("\\.4\\.9",h09v04_inFiles[i]) !=0)){
|
100
|
Apr_h09v04[i]<- h09v04_inFiles[i]
|
101
|
}
|
102
|
}
|
103
|
Apr_h09v04<-Apr_h09v04[!is.na(Apr_h09v04)]
|
104
|
Apr_h09v04
|
105
|
|
106
|
Apr_h09v04_stack<- stack(Apr_h09v04)
|
107
|
nlayers(Apr_h09v04_stack)
|
108
|
Apr_h09v04_mn<- mean(Apr_h09v04_stack, na.rm=TRUE)
|
109
|
Apr_h09v04_mnScaled<- Apr_h09v04_mn*0.0005
|
110
|
|
111
|
Out_Apr<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/Apr_h09v04_AVG_Scaled.img"
|
112
|
writeRaster (Apr_h09v04_mnScaled, filename=Out_Apr,format="HFA",overwrite=TRUE)
|
113
|
|
114
|
May_h09v04<- c(1:length (h09v04_inFiles))*NA
|
115
|
for (i in 1: length(h09v04_inFiles)){
|
116
|
if (length(grep("\\.5\\.1",h09v04_inFiles[i]) !=0)||
|
117
|
length(grep("\\.5\\.2",h09v04_inFiles[i]) !=0)||
|
118
|
length(grep("\\.5\\.3",h09v04_inFiles[i]) !=0)||
|
119
|
length(grep("\\.5\\.4",h09v04_inFiles[i]) !=0)||
|
120
|
length(grep("\\.5\\.5",h09v04_inFiles[i]) !=0)||
|
121
|
length(grep("\\.5\\.6",h09v04_inFiles[i]) !=0)||
|
122
|
length(grep("\\.5\\.7",h09v04_inFiles[i]) !=0)||
|
123
|
length(grep("\\.5\\.8",h09v04_inFiles[i]) !=0)||
|
124
|
length(grep("\\.5\\.9",h09v04_inFiles[i]) !=0)){
|
125
|
May_h09v04[i]<- h09v04_inFiles[i]
|
126
|
}
|
127
|
}
|
128
|
May_h09v04<-May_h09v04[!is.na(May_h09v04)]
|
129
|
May_h09v04
|
130
|
|
131
|
May_h09v04_stack<- stack(May_h09v04)
|
132
|
nlayers(May_h09v04_stack)
|
133
|
May_h09v04_mn<- mean(May_h09v04_stack, na.rm=TRUE)
|
134
|
May_h09v04_mnScaled<- May_h09v04_mn*0.0005
|
135
|
|
136
|
Out_May<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/May_h09v04_AVG_Scaled.img"
|
137
|
writeRaster (May_h09v04_mnScaled, filename=Out_May,format="HFA",overwrite=TRUE)
|
138
|
|
139
|
Jun_h09v04<- c(1:length (h09v04_inFiles))*NA
|
140
|
for (i in 1: length(h09v04_inFiles)){
|
141
|
if (length(grep("\\.6\\.1",h09v04_inFiles[i]) !=0)||
|
142
|
length(grep("\\.6\\.2",h09v04_inFiles[i]) !=0)||
|
143
|
length(grep("\\.6\\.3",h09v04_inFiles[i]) !=0)||
|
144
|
length(grep("\\.6\\.4",h09v04_inFiles[i]) !=0)||
|
145
|
length(grep("\\.6\\.5",h09v04_inFiles[i]) !=0)||
|
146
|
length(grep("\\.6\\.6",h09v04_inFiles[i]) !=0)||
|
147
|
length(grep("\\.6\\.7",h09v04_inFiles[i]) !=0)||
|
148
|
length(grep("\\.6\\.8",h09v04_inFiles[i]) !=0)||
|
149
|
length(grep("\\.6\\.9",h09v04_inFiles[i]) !=0)){
|
150
|
Jun_h09v04[i]<- h09v04_inFiles[i]
|
151
|
}
|
152
|
}
|
153
|
Jun_h09v04<-Jun_h09v04[!is.na(Jun_h09v04)]
|
154
|
Jun_h09v04
|
155
|
|
156
|
Jun_h09v04_stack<- stack(Jun_h09v04)
|
157
|
nlayers(Jun_h09v04_stack)
|
158
|
Jun_h09v04_mn<- mean(Jun_h09v04_stack, na.rm=TRUE)
|
159
|
Jun_h09v04_mnScaled<- Jun_h09v04_mn*0.0005
|
160
|
|
161
|
Out_Jun<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/Jun_h09v04_AVG_Scaled.img"
|
162
|
writeRaster (Jun_h09v04_mnScaled, filename=Out_Jun,format="HFA",overwrite=TRUE)
|
163
|
|
164
|
Jul_h09v04<- c(1:length (h09v04_inFiles))*NA
|
165
|
for (i in 1: length(h09v04_inFiles)){
|
166
|
if (length(grep("\\.7\\.1",h09v04_inFiles[i]) !=0)||
|
167
|
length(grep("\\.7\\.2",h09v04_inFiles[i]) !=0)||
|
168
|
length(grep("\\.7\\.3",h09v04_inFiles[i]) !=0)||
|
169
|
length(grep("\\.7\\.4",h09v04_inFiles[i]) !=0)||
|
170
|
length(grep("\\.7\\.5",h09v04_inFiles[i]) !=0)||
|
171
|
length(grep("\\.7\\.6",h09v04_inFiles[i]) !=0)||
|
172
|
length(grep("\\.7\\.7",h09v04_inFiles[i]) !=0)||
|
173
|
length(grep("\\.7\\.8",h09v04_inFiles[i]) !=0)||
|
174
|
length(grep("\\.7\\.9",h09v04_inFiles[i]) !=0)){
|
175
|
Jul_h09v04[i]<- h09v04_inFiles[i]
|
176
|
}
|
177
|
}
|
178
|
Jul_h09v04<-Jul_h09v04[!is.na(Jul_h09v04)]
|
179
|
Jul_h09v04
|
180
|
|
181
|
Jul_h09v04_stack<- stack(Jul_h09v04)
|
182
|
nlayers(Jul_h09v04_stack)
|
183
|
Jul_h09v04_mn<- mean(Jul_h09v04_stack, na.rm=TRUE)
|
184
|
Jul_h09v04_mnScaled<- Jul_h09v04_mn*0.0005
|
185
|
|
186
|
Out_Jul<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/Jul_h09v04_AVG_Scaled.img"
|
187
|
writeRaster (Jul_h09v04_mnScaled, filename=Out_Jul,format="HFA",overwrite=TRUE)
|
188
|
|
189
|
Aug_h09v04<- c(1:length (h09v04_inFiles))*NA
|
190
|
for (i in 1: length(h09v04_inFiles)){
|
191
|
if (length(grep("\\.8\\.1",h09v04_inFiles[i]) !=0)||
|
192
|
length(grep("\\.8\\.2",h09v04_inFiles[i]) !=0)||
|
193
|
length(grep("\\.8\\.3",h09v04_inFiles[i]) !=0)||
|
194
|
length(grep("\\.8\\.4",h09v04_inFiles[i]) !=0)||
|
195
|
length(grep("\\.8\\.5",h09v04_inFiles[i]) !=0)||
|
196
|
length(grep("\\.8\\.6",h09v04_inFiles[i]) !=0)||
|
197
|
length(grep("\\.8\\.7",h09v04_inFiles[i]) !=0)||
|
198
|
length(grep("\\.8\\.8",h09v04_inFiles[i]) !=0)||
|
199
|
length(grep("\\.8\\.9",h09v04_inFiles[i]) !=0)){
|
200
|
Aug_h09v04[i]<- h09v04_inFiles[i]
|
201
|
}
|
202
|
}
|
203
|
Aug_h09v04<-Aug_h09v04[!is.na(Aug_h09v04)]
|
204
|
Aug_h09v04
|
205
|
|
206
|
Aug_h09v04_stack<- stack(Aug_h09v04)
|
207
|
nlayers(Aug_h09v04_stack)
|
208
|
Aug_h09v04_mn<- mean(Aug_h09v04_stack, na.rm=TRUE)
|
209
|
Aug_h09v04_mnScaled<- Aug_h09v04_mn*0.0005
|
210
|
|
211
|
Out_Aug<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/Aug_h09v04_AVG_Scaled.img"
|
212
|
writeRaster (Aug_h09v04_mnScaled, filename=Out_Aug,format="HFA",overwrite=TRUE)
|
213
|
|
214
|
Sep_h09v04<- c(1:length (h09v04_inFiles))*NA
|
215
|
for (i in 1: length(h09v04_inFiles)){
|
216
|
if (length(grep("\\.9\\.1",h09v04_inFiles[i]) !=0)||
|
217
|
length(grep("\\.9\\.2",h09v04_inFiles[i]) !=0)||
|
218
|
length(grep("\\.9\\.3",h09v04_inFiles[i]) !=0)||
|
219
|
length(grep("\\.9\\.4",h09v04_inFiles[i]) !=0)||
|
220
|
length(grep("\\.9\\.5",h09v04_inFiles[i]) !=0)||
|
221
|
length(grep("\\.9\\.6",h09v04_inFiles[i]) !=0)||
|
222
|
length(grep("\\.9\\.7",h09v04_inFiles[i]) !=0)||
|
223
|
length(grep("\\.9\\.8",h09v04_inFiles[i]) !=0)||
|
224
|
length(grep("\\.9\\.9",h09v04_inFiles[i]) !=0)){
|
225
|
Sep_h09v04[i]<- h09v04_inFiles[i]
|
226
|
}
|
227
|
}
|
228
|
Sep_h09v04<-Sep_h09v04[!is.na(Sep_h09v04)]
|
229
|
Sep_h09v04
|
230
|
|
231
|
Sep_h09v04_stack<- stack(Sep_h09v04)
|
232
|
nlayers(Sep_h09v04_stack)
|
233
|
Sep_h09v04_mn<- mean(Sep_h09v04_stack, na.rm=TRUE)
|
234
|
Sep_h09v04_mnScaled<- Sep_h09v04_mn*0.0005
|
235
|
|
236
|
Out_Sep<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/Sep_h09v04_AVG_Scaled.img"
|
237
|
writeRaster (Sep_h09v04_mnScaled, filename=Out_Sep,format="HFA",overwrite=TRUE)
|
238
|
|
239
|
Oct_h09v04<- c(1:length (h09v04_inFiles))*NA
|
240
|
for (i in 1: length(h09v04_inFiles)){
|
241
|
if (length(grep("\\.10\\.1",h09v04_inFiles[i]) !=0)||
|
242
|
length(grep("\\.10\\.2",h09v04_inFiles[i]) !=0)||
|
243
|
length(grep("\\.10\\.3",h09v04_inFiles[i]) !=0)||
|
244
|
length(grep("\\.10\\.4",h09v04_inFiles[i]) !=0)||
|
245
|
length(grep("\\.10\\.5",h09v04_inFiles[i]) !=0)||
|
246
|
length(grep("\\.10\\.6",h09v04_inFiles[i]) !=0)||
|
247
|
length(grep("\\.10\\.7",h09v04_inFiles[i]) !=0)||
|
248
|
length(grep("\\.10\\.8",h09v04_inFiles[i]) !=0)||
|
249
|
length(grep("\\.10\\.9",h09v04_inFiles[i]) !=0)){
|
250
|
Oct_h09v04[i]<- h09v04_inFiles[i]
|
251
|
}
|
252
|
}
|
253
|
Oct_h09v04<-Oct_h09v04[!is.na(Oct_h09v04)]
|
254
|
Oct_h09v04
|
255
|
|
256
|
Oct_h09v04_stack<- stack(Oct_h09v04)
|
257
|
nlayers(Oct_h09v04_stack)
|
258
|
Oct_h09v04_mn<- mean(Oct_h09v04_stack, na.rm=TRUE)
|
259
|
Oct_h09v04_mnScaled<- Oct_h09v04_mn*0.0005
|
260
|
|
261
|
Out_Oct<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/Oct_h09v04_AVG_Scaled.img"
|
262
|
writeRaster (Oct_h09v04_mnScaled, filename=Out_Oct,format="HFA",overwrite=TRUE)
|
263
|
|
264
|
Nov_h09v04<- c(1:length (h09v04_inFiles))*NA
|
265
|
for (i in 1: length(h09v04_inFiles)){
|
266
|
if (length(grep("\\.11\\.1",h09v04_inFiles[i]) !=0)||
|
267
|
length(grep("\\.11\\.2",h09v04_inFiles[i]) !=0)||
|
268
|
length(grep("\\.11\\.3",h09v04_inFiles[i]) !=0)||
|
269
|
length(grep("\\.11\\.4",h09v04_inFiles[i]) !=0)||
|
270
|
length(grep("\\.11\\.5",h09v04_inFiles[i]) !=0)||
|
271
|
length(grep("\\.11\\.6",h09v04_inFiles[i]) !=0)||
|
272
|
length(grep("\\.11\\.7",h09v04_inFiles[i]) !=0)||
|
273
|
length(grep("\\.11\\.8",h09v04_inFiles[i]) !=0)||
|
274
|
length(grep("\\.11\\.9",h09v04_inFiles[i]) !=0)){
|
275
|
Nov_h09v04[i]<- h09v04_inFiles[i]
|
276
|
}
|
277
|
}
|
278
|
Nov_h09v04<-Nov_h09v04[!is.na(Nov_h09v04)]
|
279
|
Nov_h09v04
|
280
|
|
281
|
Nov_h09v04_stack<- stack(Nov_h09v04)
|
282
|
nlayers(Nov_h09v04_stack)
|
283
|
Nov_h09v04_mn<- mean(Nov_h09v04_stack, na.rm=TRUE)
|
284
|
Nov_h09v04_mnScaled<- Nov_h09v04_mn*0.0005
|
285
|
|
286
|
Out_Nov<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/Nov_h09v04_AVG_Scaled.img"
|
287
|
writeRaster (Nov_h09v04_mnScaled, filename=Out_Nov,format="HFA",overwrite=TRUE)
|
288
|
|
289
|
Dec_h09v04<- c(1:length (h09v04_inFiles))*NA
|
290
|
for (i in 1: length(h09v04_inFiles)){
|
291
|
if (length(grep("\\.12\\.1",h09v04_inFiles[i]) !=0)||
|
292
|
length(grep("\\.12\\.2",h09v04_inFiles[i]) !=0)||
|
293
|
length(grep("\\.12\\.3",h09v04_inFiles[i]) !=0)||
|
294
|
length(grep("\\.12\\.4",h09v04_inFiles[i]) !=0)||
|
295
|
length(grep("\\.12\\.5",h09v04_inFiles[i]) !=0)||
|
296
|
length(grep("\\.12\\.6",h09v04_inFiles[i]) !=0)||
|
297
|
length(grep("\\.12\\.7",h09v04_inFiles[i]) !=0)||
|
298
|
length(grep("\\.12\\.8",h09v04_inFiles[i]) !=0)||
|
299
|
length(grep("\\.12\\.9",h09v04_inFiles[i]) !=0)){
|
300
|
Dec_h09v04[i]<- h09v04_inFiles[i]
|
301
|
}
|
302
|
}
|
303
|
Dec_h09v04<-Dec_h09v04[!is.na(Dec_h09v04)]
|
304
|
Dec_h09v04
|
305
|
|
306
|
Dec_h09v04_stack<- stack(Dec_h09v04)
|
307
|
nlayers(Dec_h09v04_stack)
|
308
|
Dec_h09v04_mn<- mean(Dec_h09v04_stack, na.rm=TRUE)
|
309
|
Dec_h09v04_mnScaled<- Dec_h09v04_mn*0.0005
|
310
|
|
311
|
Out_Dec<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/NewClearDay_MonthlyAvgs_Tiles/Dec_h09v04_AVG_Scaled.img"
|
312
|
writeRaster (Dec_h09v04_mnScaled, filename=Out_Dec,format="HFA",overwrite=TRUE)
|