1 |
7526fb1c
|
Jim Regetz
|
#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 |
|
|
}
|