VGT

Created Tuesday 30 April 2019

Status Map (quality flags)

"SMlegend" <- function()
{
desc = rep("o",7)
res = matrix("0",ncol=8,nrow=256)
for (i in 1:256){
  a = as.integer(intToBits(i-1)[1:8])
  #print(a)
  if(a[1]==0 & a[2]==0) desc[1] = "clear"
  if(a[1]==0 & a[2]==1) desc[1] = "shadow"
  if(a[1]==1 & a[2]==0) desc[1] = "uncert"
  if(a[1]==1 & a[2]==1) desc[1] = "cloud"
  if(a[3]==0) desc[2] = "no ice"
  if(a[3]==1) desc[2] = "ice"
  if(a[4]==0) desc[3] = "water"
  if(a[4]==1) desc[3] = "land"
  if(a[5]==0) desc[4] = "MIR bad"
  if(a[5]==1) desc[4] = "MIR good"
  if(a[6]==0) desc[5] = "B3 bad"
  if(a[6]==1) desc[5] = "B3 good"
  if(a[7]==0) desc[6] = "B2 bad"
  if(a[7]==1) desc[6] = "B2 good"
  if(a[8]==0) desc[7] = "B0 bad"
  if(a[8]==1) desc[7] = "B0 good"
  b = paste(a, collapse="")
  res[i,] = c(b,desc)
  #print(desc)
}
dimnames(res)[[1]] = 0:255
res
}