1
|
#Correct 1st column-------------------------------------------------------------------------------------------------
|
2
|
library (raster)
|
3
|
|
4
|
path<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/ClearDay_Original_IMG_Extracts/ByDate"
|
5
|
|
6
|
#Get h09v04 tiles:
|
7
|
h09v04_infiles <- list.files(path=path, pattern="*.h09v04.005.ClearDay.img$")
|
8
|
for (i in 1: length(h09v04_infiles)){
|
9
|
h09v04_infiles[i]<- paste (path,"/",h09v04_infiles[i],sep="")
|
10
|
}
|
11
|
h09v04_infiles
|
12
|
|
13
|
#raster(h09v04_infiles[311])[1200,1]
|
14
|
#raster(h09v04_infiles[311])[1200,2]
|
15
|
|
16
|
rastfunction<- function (h09v04_infiles){
|
17
|
r<-raster(h09v04_infiles)
|
18
|
val<-raster(h09v04_infiles)[,2]
|
19
|
r[,1]<-val
|
20
|
name1<-sub("_Extracts/ByDate","_Extracts/ByDate/BadPixelsCorrected_Tiles",h09v04_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 (h09v04_infiles)){
|
26
|
rastfunction(h09v04_infiles[i])
|
27
|
}
|
28
|
|
29
|
#Correct last row, overwrite corrected files-------------------------------------------------------------------------------------------
|
30
|
path1<- "/data/project/organisms/MODIS_LST_Oregon/ClearDayGDAL/ClearDay_Original_IMG_Extracts/ByDate/BadPixelsCorrected_Tiles"
|
31
|
|
32
|
#Get h09v04 tiles:
|
33
|
newh09v04_infiles <- list.files(path=path1, pattern="*.h09v04.005.ClearDay_corrected.img$")
|
34
|
for (i in 1: length(newh09v04_infiles)){
|
35
|
newh09v04_infiles[i]<- paste (path1,"/",newh09v04_infiles[i],sep="")
|
36
|
}
|
37
|
newh09v04_infiles
|
38
|
|
39
|
#raster(h09v04_infiles[311])[1200,1]
|
40
|
#raster(h09v04_infiles[311])[1200,2]
|
41
|
|
42
|
newrastfunction<- function (newh09v04_infiles){
|
43
|
r<-raster(newh09v04_infiles)
|
44
|
val<-raster(newh09v04_infiles)[1199,]
|
45
|
r[1200,]<-val
|
46
|
writeRaster(r,filename=newh09v04_infiles,format="HFA",overwrite=TRUE)
|
47
|
}
|
48
|
|
49
|
for (i in 1:length (newh09v04_infiles)){
|
50
|
newrastfunction(newh09v04_infiles[i])
|
51
|
}
|