/home/pascal/depot/filedata/src/dardarfiledata.h

/***************************************************************************
 *   Copyright (C) 2011 by Nicolas PASCAL   *
 *   nicolas.pascal@icare.univ-lille1.fr   *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 *   This program 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 General Public License for more details.                          *
 *                                                                         *
 *   You should have received a copy of the GNU General Public License     *
 *   along with this program; if not, write to the                         *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/
#ifndef DARDARFILEDATA_H
#define DARDARFILEDATA_H

#include "cloudsatfiledata.h"

class DARDARFileData : public SatelliteFileData, public HDFFileData {
    static const float colocation_tolerance;

    /***** constants *****/
    string latitude_name;
    string longitude_name;
    string time_name;
    enum Product_ID {
        UNDEFINED = 0,
        DARDAR_CLOUD,
        DARDAR_MASK,
    };

    /***** Product characteristics *****/
    Product_ID product_id;
    string product;
    string version;
    int lat_lon_index_max;
    short orbit_nb;

    void parse_filename(const string & short_filename);
    bool check_filename(const string & short_filename);
    void init();
    void set_lat_lon_index_max();
    void set_time_coverage();
public:
    static const int nb_dardar_bin=436;
    static const int16 dardar_height [];
    DARDARFileData(const string &name, const string &mode="r");
    ~DARDARFileData();
    void load_geolocation_data();
    void free_geolocation_data();
    void load_v_pixel() ;
    void close_data_file() {
        free_hdf_file();
    };
    void open_data_file() {
        load_hdf_file();
    };
    const int get_nb_geo_points() { return lat_lon_index_max + 1; }
    string get_product() const { return product; }
    string get_version() const { return version; }

    const bool contain_location(const float &lat,const float &lon, const double &tolerance=0.01);
    const bool contain_data(const float &lat, const float &lon, const double &time, const double &colocation_tolerance=0.01) ;
    const float get_nearest_point_distance(const float &lat,const float &lon,const float coloc_tolerance=DARDARFileData::colocation_tolerance);
     virtual const bool get_index(const float &lat, const float& lon, int &nearest_pix_idx, const float colocation_tolerance=0.01);
     void get_vindex(vector < vector < int > > &v_index, const float &lat, const float& lon,
                                            const float colocation_tolerance=0.01 ) {
        UnimplementedMethod e(__FILE__,__LINE__,__PRETTY_FUNCTION__);
        throw e;
     };
    void get_pixel_coord ( const vector < int > & ipix, float &lat, float &lon, double &time ) {
        get_pixel_coord (ipix[0], lat, lon, time);
    };
    void get_pixel_coord ( const int ipix, float &lat, float &lon, double &time ) {
        if (!is_geolocation_data_loaded())
            load_geolocation_data();
        lat  = lat_data  [ ipix ];
        lon  = lon_data  [ ipix ];
        time = time_data [ ipix ];
    };

    static const int16* get_dardar_height() {
        return DARDARFileData::dardar_height;
    }

};

#endif