NASA Ocean Color

ocssw V2022
 /*---------------------------------------------------------------------*/
 /* get_smoke.c - Eric Vermote's smoke index                            */
 /*                                                                     */
 /* Inputs:                                                             */
 /*     l2rec - level-2 structure containing one complete scan after    */
 /*             atmospheric correction.                                 */
 /* Outputs:                                                            */
 /*     smoke - smoke index, 1 value per pixel.                         */
 /*                                                                     */
 /* Written by: Bryan Franz, SAIC-GSC, February 2000                    */
 /*                                                                     */
 /*---------------------------------------------------------------------*/
  
 #include <stdlib.h>
 #include <math.h>
 #include "l12_proto.h"
  
 static float fill = 0.0;
  
 void get_smoke(l2str *l2rec, float smoke[]) {
     int32_t ip, ipb;
  
     float rho1;
     float rho2;
     float rho4;
     float rho5;
     float rho6;
  
     l1str *l1rec = l2rec->l1rec;
  
     for (ip = 0; ip < l1rec->npix; ip++) {
  
         ipb = l1rec->l1file->nbands*ip;
  
         rho1 = l1rec->rhos[ipb + 0];
         rho2 = l1rec->rhos[ipb + 1];
         rho4 = l1rec->rhos[ipb + 3];
         rho5 = l1rec->rhos[ipb + 4];
         rho6 = l1rec->rhos[ipb + 5];
  
         if (rho1 < 0.0 || rho6 > 0.4)
             smoke[ip] = fill;
  
         else if (rho2 < rho4 || rho5 < rho6)
             smoke[ip] = 0.0;
  
         else
             smoke[ip] = rho1 * (rho2 - rho4) * (rho5 - rho6) * 10000.0;
  
     }
  
 }
  
get_smoke
void get_smoke(l2str *l2rec, float smoke[])
Definition: get_smoke.c:20
l1rec
read l1rec
Definition: HOWTO_Add_a_sensor.txt:72
l12_proto.h