Revision f693d9e4
Added by Adam M. Wilson almost 11 years ago
| climate/research/LST_missingdata.R | ||
|---|---|---|
| 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() |
|
Also available in: Unified diff
Adding Missing data plot for Alberto's data