NASA Ocean Color

ocssw V2022

 /*
  * To change this license header, choose License Headers in Project Properties.
  * To change this template file, choose Tools | Templates
  * and open the template in the editor.
  */
  
 #include <cfloat>
 #include <cmath>
 #include <cstdio>
 #include <cstdlib>
 #include <cstring>
 #include <pugixml.hpp>
 #include <vector>
 #include <boost/foreach.hpp>
  
 #include <dfutils.h>
 #include <genutils.h>
 #include <sensorDefs.h>
 #include <timeutils.h>
  
 #include <boost/algorithm/string.hpp>
  
 #include "filetype.h"
  
  
 extern "C" file_format chk_safe_xml(char *filename) {
     file_format ret = {FT_INVALID, -1, -1};
     pugi::xml_document rootNode;
     pugi::xml_node dataNode;
  
     const char *productName;
     const char *platformName;
  
     if (!rootNode.load_file(filename)) {
         return ret;
     }
  
     productName = rootNode.first_element_by_path("xfdu:XFDU/informationPackageMap/xfdu:contentUnit").attribute("textInfo").value();
     if (strstr(productName, "OLCI") != NULL && strstr(productName, "SENTINEL-3") != NULL && strstr(productName, "Level 1") != NULL) {
         dataNode = rootNode.first_element_by_path("xfdu:XFDU/metadataSection").find_child_by_attribute("metadataObject", "ID", "platform");
         platformName = dataNode.first_element_by_path("metadataWrap/xmlData/sentinel-safe:platform/sentinel-safe:number").child_value();
         if(strcmp("A", platformName) == 0) {
             if (want_verbose != 0){
                 printf("Found OLCI Sentinel-3A product: %s\n", productName);
             }
             ret.type = FT_OLCI;
             ret.sensor_id = OLCIS3A;
             ret.subsensor_id = OLCI_S3A;
         } else if(strcmp("B", platformName) == 0) {
             if (want_verbose != 0){
                 printf("Found OLCI Sentinel-3B product: %s\n", productName);
             }
             ret.type = FT_OLCI;
             ret.sensor_id = OLCIS3B;
             ret.subsensor_id = OLCI_S3B;
         } else {
             printf("Unknown platform: %s  for product: %s\n", platformName, productName);
             exit(EXIT_FAILURE);
         }
     } else if (strstr(productName, "MSI") != NULL && strstr(productName, "SENTINEL-2") != NULL && strstr(productName, "Level-1C") != NULL) {
         dataNode = rootNode.first_element_by_path("xfdu:XFDU/metadataSection").find_child_by_attribute("metadataObject", "ID", "platform");
         platformName = dataNode.first_element_by_path("metadataWrap/xmlData/safe:platform/safe:number").child_value();
         if(!strcmp(platformName, "2A")) {
             if (want_verbose != 0){
                 printf("Found MSI Sentinel-2A product: %s\n", productName);
             }
             ret.type = FT_MSIL1C;
             ret.sensor_id = MSIS2A;
             ret.subsensor_id = MSI_S2A;
         } else if(!strcmp(platformName, "2B")) {
             if (want_verbose != 0){
                 printf("Found MSI Sentinel-2B product: %s\n", productName);
             }
             ret.type = FT_MSIL1C;
             ret.sensor_id = MSIS2B;
             ret.subsensor_id = MSI_S2B;
         } else {
             printf("Unknown platform: %s  for product: %s\n", platformName, productName);
             exit(EXIT_FAILURE);
         }
     } else if (strstr(productName, "MERIS") != NULL && strstr(productName, "ENVISAT") != NULL && strstr(productName, "Level 1") != NULL) {
         if (want_verbose != 0){
             printf("Found MERIS L1B SAFE product: %s\n", productName);
         }
         ret.type = FT_MERISL1BSAFE;
         ret.sensor_id = MERIS;
         ret.subsensor_id = -1;
     }
     return ret;
 }
  