1
|
#Correct last row -------------------------------------------------------------------------------------------
|
2
|
library (raster)
|
3
|
|
4
|
path<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/ClearDay_Original_IMG_Extracts/ByDate"
|
5
|
|
6
|
#Get h08v04 tiles:
|
7
|
h08v04_infiles <- list.files(path=path, pattern="*.h08v04.005.ClearDay.img$")
|
8
|
for (i in 1: length(h08v04_infiles)){
|
9
|
h08v04_infiles[i]<- paste (path,"/",h08v04_infiles[i],sep="")
|
10
|
}
|
11
|
h08v04_infiles
|
12
|
|
13
|
#raster(h09v04_infiles[311])[1200,1]
|
14
|
#raster(h09v04_infiles[311])[1200,2]
|
15
|
|
16
|
rastfunction<- function (h08v04_infiles){
|
17
|
r<-raster(h08v04_infiles)
|
18
|
val<-raster(h08v04_infiles)[1199,]
|
19
|
r[1200,]<-val
|
20
|
name1<-sub("_Extracts/ByDate","_Extracts/ByDate/BadPixelsCorrected_Tiles",h08v04_infiles)
|
21
|
name<-sub(".ClearDay.img",".ClearDay_corrected.img",name1)
|
22
|
writeRaster(r,filename=name,format="HFA")
|
23
|
}
|
24
|
|
25
|
for (i in 1:length (h08v04_infiles)){
|
26
|
rastfunction(h08v04_infiles[i])
|
27
|
}
|