HE2CF.h

// This file is part of the hdf4 data handler for the OPeNDAP data server.
// Copyright (c) The HDF Group
//
// Author: Hyo-Kyung Lee <hyoklee@hdfgroup.org>
//
// This is free software; you can redistribute it and/or modify it under the
// terms of the GNU Lesser General Public License as published by the Free
// Software Foundation; either version 2.1 of the License, or (at your
// option) any later version.
//
// This software is distributed in the hope that it will be useful, but
// WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
// or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
// License for more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with this software; if not, write to the Free Software Foundation,
// Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
//
// You can contact OPeNDAP, Inc. at PO Box 112, Saunderstown, RI. 02874-0112.
#ifndef _HE2CF_H_
#define _HE2CF_H_
 
#include <mfhdf.h>
#include <hdf.h>
//#include <HdfEosDef.h>
#include <iostream>
#include <sstream>
#include <vector>
#include <map>
#include <algorithm>
#include <libdap/debug.h>
#include <libdap/DAS.h>
#include "HDFCFUtil.h"
 

class HE2CF
{
 
    private:
 
        libdap::DAS* das = nullptr;
 
        // SDStart ID
        int32 sd_id = -1;
 
        // Hopen ID 
        int32 file_id = -1;
 
        // Number of global attributes
        int32 num_global_attributes = -1;
 
        // ECS metadata 
        std::string metadata;
 
        // group name
        std::string gname;
 
        // Store all metadata names.
        // Not an ideal approach, need to re-visit later. KY 2012-6-11
        std::vector<std::string>  eosmetadata_namelist;
    
        // Data field SDS name to SDS id 
        std::map < std::string, int32 > vg_dsd_map;
 
        // Data field vdata name to vdata id
        std::map < std::string, int32 > vg_dvd_map;
 
        // Geolocation field SDS name to SDS id 
        std::map < std::string, int32 > vg_gsd_map;
 
        // Geolocation field vdata name to vdata id
        std::map < std::string, int32 > vg_gvd_map;
 
        // Add metadata_name to the metadata name list.
        void set_eosmetadata_namelist(const std::string &metadata_name)
        {
            eosmetadata_namelist.push_back(metadata_name);
        }
 
        // Is this EOS metadata
        bool is_eosmetadata(const std::string& metadata_name) {
            return (std::find(eosmetadata_namelist.begin(),eosmetadata_namelist.end(),metadata_name) !=eosmetadata_namelist.end());
        }
         
        // Get HDF-EOS2 "Data Fields" and "Geolocation Fields" group object tag and object reference numbers.
        bool get_vgroup_field_refids(const std::string&  _gname, int32* _ref_df, int32* _ref_gf);
 
        // Open SD 
        bool open_sd(const std::string& filename,const int sd_id);
 
        // Open vgroup
        bool open_vgroup(const std::string& filename,const int fileid);
 
        // Combine ECS metadata coremetadata.0, coremetadata.1 etc. into one string.
        bool set_metadata(const std::string& metadataname,std::vector<std::string>&non_num_names, std::vector<std::string>&non_num_data);
 
        // This routine will generate three ECS metadata lists. Note in theory list sl1 and sl2 should be sorted.
        // Since the ECS metadata is always written(sorted) in increasing numeric order, we don't perform this now.
        // Should watch if there are any outliers. 
        void arrange_list(std::list<std::string> & sl1, std::list<std::string>&sl2,std::vector<std::string>&v1,const std::string & name,int& flag) const;
 
        // Obtain SD attribute value
        void obtain_SD_attr_value(const std::string &,std::string&) const;
 
        // Create SDS name to SDS ID and Vdata name to vdata ID maps.
        bool set_vgroup_map(int32 refid,bool isgeo);
 
        // Write the long_name attribute.
        bool write_attr_long_name(const std::string& long_name, 
            const std::string& varname,
            int fieldtype);
        bool write_attr_long_name(const std::string& group_name, 
            const std::string& long_name, 
            const std::string& varname,
            int fieldtype);
 
        // Write the SD attribute.
        bool write_attr_sd(int32 sds_id, const std::string& newfname,int fieldtype);

        short check_scale_offset(int32 sds_id, bool is_scale);
 
 
        // Write the Vdata attribute.
        bool write_attr_vdata(int32 vd_id, const std::string& newfname,int fieldtype);
        void throw_error(const std::string &_error);
    
    public:
        HE2CF();
        virtual ~HE2CF();

        bool   close();

        string get_metadata(const std::string& metadataname,bool&suffix_is_num,std::vector<std::string>&non_num_names, std::vector<std::string>&non_num_data);
        bool set_non_ecsmetadata_attrs();

        bool   open(const std::string& filename,const int sd_id, const int file_id);

        void   set_DAS(libdap::DAS* das);

        bool   write_attribute(const std::string& gname, 
            const std::string& fname, 
            const std::string& newfname,
            int n_groups,
            int fieldtype);

        bool write_attribute_FillValue(const std::string& varname, int type, float val);

        bool write_attribute_coordinates(const std::string& varname, const std::string& coord);

        bool write_attribute_units(const std::string& varname, const std::string& units);
 
 
};
 
 
#endif