data/shading.m

   1 % Load in the data
   2 lit = load("-ascii", "shading_lit.adf");
   3 unlit = load("-ascii", "shading_unlit.adf");
   4 % apply acos to column 1 (shading)
   5 lit(:,1) = acos(lit(:,1)')';
   6 unlit(:,1) = acos(unlit(:,1)')';
   7
   8 sizel = size(lit);
   9 sizeu = size(unlit);
  10
  11 lit_b = lit(:,2);
  12 lit_g = lit(:,3);
  13 lit_ = lit(:,4);
  14 lit_mir1 = lit(:,6);
  15 lit_ltir = lit(:,7);
  16 lit_htir = lit(:,8);
  17 lit_mir2 = lit(:,9);
  18
  19 A=(0:0.005:1.55)';
  20
  21 C=140*cos(A);
  22 lit_choice = lit_ltir;
  23
  24 %C=115*cos(A);
  25 %lit_choice = lit_htir;
  26
  27 C=100*cos(A);
  28 lit_choice = (lit_mir1 + lit_mir2);
  29
  30 scatter([lit_choice(1:20:sizel(1),:);C], [lit(1:20:sizel(1),1);A]);
  31
  32 unlit_b = unlit(:,2);
  33 unlit_g = unlit(:,3);
  34 unlit_r = unlit(:,4);
  35 unlit_nir = unlit(:,5);
  36 unlit_mir1 = unlit(:,6);
  37 unlit_ltir = unlit(:,7);
  38 unlit_htir = unlit(:,8);
  39 unlit_mir2 = unlit(:,9);
  40
  41 C=140*cos(A);
  42 unlit_choice = unlit_ltir;
  43
  44 scatter([unlit_choice(1:5:sizeu(1),:);C], [unlit(1:5:sizeu(1),1);A]);