WeatherData.h

Go to the documentation of this file.

#ifndef GlobeEngine_WeatherData_h
#define GlobeEngine_WeatherData_h

#include <vmmlib/vmmlib.hpp>
#include <vector>
#include "FlowDataBlueprint.h"

#include "ncFile.h"

namespace geFlow {

    struct ClimateGridInfo {
        float* longCoords;
        float* latCoords;
        float* heightCoords;
    };

    struct ClimateDataInfo {
        vmml::Vector4ui sizes;
        std::vector<std::string> names;
    };

    struct ClimateDataLayer2D {
        std::vector<vmml::Vector2f> directions;
    };

    struct ClimateDataLayer1D {
        vmml::Vector2f minMax;
        std::vector<float> data;
    };

    class WeatherData {
            
    public:
        WeatherData();
        ~WeatherData();
        void clear();
        void loadFromNetCDF(std::shared_ptr<geFlow::WeatherDataLoadingInfo> _loadinginfo);
    
        const std::vector<vmml::Vector2f>& getWindDirections(int _timestep = 0) const;
        const vmml::Vector2f& getWindDirectionAtIdx(int _longIdx, int _latIdx, int _timestep = 0) const;

        unsigned int getTotalDataSize() const;
        unsigned int getTimeSteps() const;
        vmml::Vector2s getLongLatCellExtent() const;
        vmml::Vector3s getLongHeightLatCellExtent() const;
        const vmml::Vector4ui& getCellExtent() const;
        const vmml::Vector2d& getGridLongLatMidPoint() const;
        const vmml::Vector4d& getGridLongLatBounds() const;
        const vmml::Vector2d& getGridLongLatBoundSize() const;
        const ClimateGridInfo& getGridInfo() const;
        ClimateDataLayer1D& getCloudInfo(int _timestep = 0);
        const double getGridHeight() const;
    private:
        void setDataArraySize();
        void loadWindInformation(const netCDF::NcFile& _dataFile, 
            std::shared_ptr<geFlow::WeatherDataLoadingInfo> _info);
        void loadCloudInformation(const netCDF::NcFile& _dataFile,
            std::shared_ptr<geFlow::WeatherDataLoadingInfo> _info);
        void loadVolumeInformation(const netCDF::NcFile& _dataFile,
            std::shared_ptr<geFlow::WeatherDataLoadingInfo> _info);
        // usually assumes ordered data as in 
        // the variable: dataInfo
        ClimateLoadingRange getLoadingRange(
            std::vector<std::string> _orderingStrings,
            vmml::Vector4ui _ranges);
        int getIdxForName(std::string _name) const;
        // returns 0=totaldataperlayer and 1=totaldatatoload
        static vmml::Vector2ui computeDataSizeForVariable(vmml::Vector4ui _dataSize);

        ClimateGridInfo grid;
        std::vector<ClimateDataLayer2D> wind;
        std::vector<ClimateDataLayer1D> layer1D;
        // usually assumes ordered data:
        // longitude = x-axis
        // height = y-axis
        // latitude = z-axis
        // timesteps = w-axis
        ClimateDataInfo dataInfo;
        //vmml::Vector4ui dataSize;
        // long min, long max, lat min, lat max 
        vmml::Vector2d gridMidpoint;
        vmml::Vector4d gridbounds;
        vmml::Vector2d gridboundsLenght;
        // top and bottom bound
        vmml::Vector2d volumeHeightbounds;
        double volumeHeight;

        unsigned int totDataSize;
    };
}
#endif