root/climate/research/LST_missingdata.R @ b74bf08c
| 1 | f693d9e4 | Adam M. Wilson | ### 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() |