LayeredMeshGenerator.h

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#pragma once
 
#include <limits>
#include <memory>
#include <string>
#include <vector>
 
namespace GeoLib
{
class Raster;
}
 
namespace MeshLib
{
class Mesh;
class Node;
class Element;
}  // namespace MeshLib

class LayeredMeshGenerator
{
public:
    virtual bool createLayers(MeshLib::Mesh const& mesh,
                              std::vector<GeoLib::Raster const*> const& rasters,
                              double minimum_thickness,
                              double noDataReplacementValue = 0.0) final;

    virtual bool createRasterLayers(
        MeshLib::Mesh const& mesh,
        std::vector<GeoLib::Raster const*> const& rasters,
        double minimum_thickness,
        double noDataReplacementValue) = 0;

    std::unique_ptr<MeshLib::Mesh> getMesh(std::string const& mesh_name) const;
 
protected:
    LayeredMeshGenerator() = default;
    virtual ~LayeredMeshGenerator() = default;

    virtual void addLayerToMesh(MeshLib::Mesh const& mesh_layer,
                                unsigned layer_id,
                                GeoLib::Raster const& raster) = 0;

    MeshLib::Node* getNewLayerNode(MeshLib::Node const& dem_node,
                                   MeshLib::Node const& last_layer_node,
                                   GeoLib::Raster const& raster,
                                   std::size_t new_node_id) const;

    static double calcEpsilon(GeoLib::Raster const& low,
                              GeoLib::Raster const& high);
 
    double _elevation_epsilon{0.0001};
    double _minimum_thickness{std::numeric_limits<float>::epsilon()};
    std::vector<int> _materials;
    std::vector<MeshLib::Node*> _nodes;
    std::vector<MeshLib::Element*> _elements;
};