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

/***************************************************************************
 *   Copyright (C) 2005 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 IIRFILEDATA_H
#define IIRFILEDATA_H

#include "hdffiledata.h"
#include "pixel.h"
#include "satellitefiledata.h"

static const int GET_SHOT_TIME_INDEX_TRACE=0;
static const int BRUTE_FORCE_SEARCH=1;

class IIRFileData : public SatelliteFileData, public HDFFileData {

public:
    enum ProductType{
            UNDEFINED,
            L1,
            L2_TRACK,
            L2_SWATH,
    };
    static const int IIR_swath_size;
    static const int nb_record_max;

private:
    int level;
    ProductType product_type;
    bool day_mode;

// // // //     float *lat_data;
// // // //     float *lon_data;
// // // //     double *time_data;

    vector<PixelType> v_pixel;

    static const string latitude_sds_name;
    static const string longitude_sds_name;
    string time_sds_name;
    static const float colocation_tolerance;

    int lat_lon_index_max[2];

    void init();
    void check_filename( const string& short_filename ) const;
    void parse_filename( const string& short_filename );
    void set_time_coverage();
    void set_lat_lon_index_max();
protected:
    virtual void load_v_pixel() ;

public:
    IIRFileData(const string& _name="", const string& mode="r");
    ~IIRFileData();
    const bool is_day() const;
    void get_latitude_data(float* data);
    void get_longitude_data(float* data);
    void get_time_data(double* data);
    const int get_track_size() const ;
    const int get_swath_size() const ;

    const bool get_index(const float &lat, const float& lon, const double &time, int * near_point_idx, const float coloc_tolerance=IIRFileData::colocation_tolerance);
     const int get_nearest_point_index( const float &lat, const float &lon,const double &time);
    const bool contains_data_at(const float &lat,const float &lon, const double &time) const;
    void load_geolocation_data();
    void free_geolocation_data();
    const bool is_geolocation_data_loaded() const;

    const bool get_index(const float &lat, const float& lon, int &nearest_pix_idx, const float colocation_tolerance=IIRFileData::colocation_tolerance);

    const float get_nearest_point_distance(const float &lat,const float &lon,const float coloc_tolerance=IIRFileData::colocation_tolerance);
    const bool contain_location(const float &lat,const float &lon, const double &tolerance=IIRFileData::colocation_tolerance);
    const bool contain_data(const float &lat, const float &lon, const double &time, const double &colocation_tolerance=IIRFileData::colocation_tolerance) ;
    virtual void get_pixel_coord ( const vector < int > & ipix, float &lat, float &lon, double &time ) {
        UnimplementedMethod e(__FILE__,__LINE__,__PRETTY_FUNCTION__);
        throw e;
    };

     virtual void get_vindex(vector < vector < int > > &v_index, const float &lat, const float& lon,
                                            const float colocation_tolerance ) {
        UnimplementedMethod e(__FILE__,__LINE__,__PRETTY_FUNCTION__);
        throw e;
     };
    virtual void close_data_file() {
        free_hdf_file();
    };
    virtual void open_data_file() {
        load_hdf_file();
    };

};

#endif