root/climate/research/LST_missingdata.R @ f41365c2
| 1 |
### Short script to make a plot showing missing LST data in tiles from Alberto
|
|---|---|
| 2 |
library(rasterVis) |
| 3 |
|
| 4 |
setwd("~/Downloads/nasa/") |
| 5 |
|
| 6 |
|
| 7 |
f=data.frame(full=T,path=list.files(pattern="tif$"),stringsAsFactors=F) |
| 8 |
f$month=as.numeric(do.call(rbind,strsplit(f$path,"_|[.]"))[,7]) |
| 9 |
f=f[order(f$month),] |
| 10 |
f$mn=month.name[f$month] |
| 11 |
|
| 12 |
d=stack(f$path) |
| 13 |
names(d)=f$mn |
| 14 |
|
| 15 |
|
| 16 |
colramp=colorRampPalette(c("blue","orange","red")) |
| 17 |
|
| 18 |
png("Climatologies.png",width=1500,height=600,pointsize=22) |
| 19 |
levelplot(d,col.regions=c("grey",colramp(99)),at=c(-0.5,0.5,seq(1,70,len=99)),main="Land Surface Temperature - Monthly Climatologies",sub="Tile H08v05 (California and Northern Mexico) \n Grey indicates missing data") |
| 20 |
dev.off() |