fits_read_attributes.cc

// fits_read_attributes.cc
 
// This file is part of fits_handler, a data handler for the OPeNDAP data
// server. 
 
// Copyright (c) 2004,2005 University Corporation for Atmospheric Research
// Author: Patrick West <pwest@ucar.edu> and Jose Garcia <jgarcia@ucar.edu>
//
// This library 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 library 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 library; if not, write to the Free Software
// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
//
// You can contact University Corporation for Atmospheric Research at
// 3080 Center Green Drive, Boulder, CO 80301
 
// (c) COPYRIGHT University Corporation for Atmostpheric Research 2004-2005
// Please read the full copyright statement in the file COPYRIGHT_UCAR.
//
// Authors:
//      pwest       Patrick West <pwest@ucar.edu>
//      jgarcia     Jose Garcia <jgarcia@ucar.edu>
 
// #include <stdlib.h>
// #include <ctype.h>
 
#include <string>
#include <sstream>
 
#include <fitsio.h>
 
#include <libdap/AttrTable.h>
#include <libdap/DAS.h>
 
#include "fits_read_attributes.h"
 
using namespace libdap;
using namespace std;
 
bool fits_handler::fits_read_attributes(DAS &das, const string &filename, string &error)
{
    fitsfile *fptr;
    int status = 0;
    if (fits_open_file(&fptr, filename.c_str(), READONLY, &status)) {
        error = "Can not open fits file ";
        error += filename;
        return false;
    }
 
    int hdutype;
    for (int ii = 1; !(fits_movabs_hdu(fptr, ii, &hdutype, &status)); ii++) {
        AttrTable *at = new AttrTable();
 
        /* get no. of keywords */
        int nkeys, keypos;
        if (fits_get_hdrpos(fptr, &nkeys, &keypos, &status)) {
            fits_close_file(fptr, &status);
            return false;
        }
 
        for (int jj = 1; jj <= nkeys; jj++) {
            char name[FLEN_KEYWORD];
            char value[FLEN_VALUE];
            char comment[FLEN_COMMENT];
 
            if (fits_read_keyn(fptr, jj, name, value, comment, &status)) {
                fits_close_file(fptr, &status);
                return false;
            }
#if 0
            // name, ..., will always be true
            string s_name = "";
            if (name)
                s_name = name;
 
            string s_value = "";
            if (value)
                s_value = value;
 
            string s_comment = "";
            if (comment)
                s_comment = comment;
#endif
            string s_name = name;
            if (s_name.empty()/* == ""*/) {
                ostringstream oss;
                oss << "key_" << jj;
                s_name = oss.str();
#if 0
                char tmp[100];
                ltoa(jj, tmp, 10);
                s_name = (string) "key_" + tmp;
#endif
            }
 
            {
                string com ;
                // No need to add quote here. 
                // libdap4 will add the quote when generating DAS. 
                // KY 2022-08-25
#if 0
                //string com = "\"";
#endif
                com.append(value);
                com.append(" / ");
                com.append(comment);
#if 0
                //com.append("\"");
#endif
                at->append_attr(s_name, "String"/*type*/, com);
            }
 
        }
#if 0
        string str = "HDU_";
        char tmp[100];
        ltoa(ii, tmp, 10);
        str.append(tmp);
#endif
        ostringstream oss;
        oss << "HDU_" << ii;
        das.add_table(oss.str(), at);
    }
 
    if (status == END_OF_FILE) /* status values are defined in fitsioc.h */
        status = 0; /* got the expected EOF error; reset = 0  */
    else {
        /* got an unexpected error                */
        fits_close_file(fptr, &status);
        return false;
    }
 
    if (fits_close_file(fptr, &status))
        return false;
 
    return true;
}