28 const std::vector<GeoLib::Raster const*>& rasters,
29 double minimum_thickness,
30 double noDataReplacementValue)
46 std::unique_ptr<MeshLib::Mesh> top(
50 top = std::make_unique<MeshLib::Mesh>(mesh);
54 *top, *rasters.back(), noDataReplacementValue))
59 std::unique_ptr<MeshLib::Mesh> bottom(
new MeshLib::Mesh(*top));
70 ERR(
"The materials vector is not empty.");
76 const std::size_t nRasters(rasters.size());
77 for (std::size_t i = 1; i < nRasters; ++i)
94 const std::vector<MeshLib::Node*>& nodes(dem_mesh.
getNodes());
95 const std::size_t node_id_offset(
_nodes.size());
96 const std::size_t last_layer_node_offset(node_id_offset - nNodes);
98 for (std::size_t i = 0; i < nNodes; ++i)
101 *
_nodes[last_layer_node_offset + i],
106 const std::vector<MeshLib::Element*>& layer_elements(
112 std::array<MeshLib::Node*, 3> tri_nodes = {
121 std::array<MeshLib::Node*, 4> quad_nodes = {
135 const unsigned nLayerBoundaries(nLayers - 1);
137 const std::vector<MeshLib::Element*>& layer_elements(layer.
getElements());
140 const std::size_t nElemNodes(elem->getNumberOfBaseNodes());
141 for (
unsigned i = 0; i < nElemNodes; ++i)
143 if (elem->getNeighbor(i) ==
nullptr)
145 for (
unsigned j = 0; j < nLayerBoundaries; ++j)
147 const std::size_t offset(j * nNodes);
159 Eigen::Map<Eigen::Vector3d const>(n0->
getCoords());
161 Eigen::Map<Eigen::Vector3d const>(n1->
getCoords());
163 Eigen::Map<Eigen::Vector3d const>(n2->
getCoords());
165 Eigen::Map<Eigen::Vector3d const>(n3->
getCoords());
166 double const eps = std::numeric_limits<double>::epsilon();
168 if ((v2 - v1).
norm() > eps)
170 const std::array tri_nodes = {n0, n2, n1};
174 if ((v3 - v0).
norm() > eps)
176 const std::array tri_nodes = {n0, n3, n2};
187 std::size_t nElementsPerLayer)
189 for (std::size_t i = nLayers - 1; i > 0; --i)
191 const std::size_t lower_offset((i - 1) * nElementsPerLayer);
192 const std::size_t upper_offset(i * nElementsPerLayer);
193 for (std::size_t j = 0; j < nElementsPerLayer; ++j)
200 for (std::size_t k = 0; k < nElemNodes; ++k)
209 if (count == nElemNodes)
Definition of Duplicate functions.
Definition of the LayeredVolume class.
void ERR(char const *fmt, Args const &... args)
Definition of the MeshLayerMapper class.
Definition of the class Properties that implements a container of properties.
Definition of the Quad class.
Definition of the GeoLib::Raster class.
Definition of the Tri class.
Class Raster is used for managing raster data.
double _minimum_thickness
std::vector< int > _materials
double _elevation_epsilon
static double calcEpsilon(GeoLib::Raster const &low, GeoLib::Raster const &high)
Calculates a data-dependent epsilon value.
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
std::vector< MeshLib::Node * > _nodes
std::vector< MeshLib::Element * > _elements
void removeCongruentElements(std::size_t nLayers, std::size_t nElementsPerLayer)
Removes duplicate 2D elements (possible due to outcroppings)
std::vector< MeshLib::Node > _attribute_points
void addLayerToMesh(const MeshLib::Mesh &dem_mesh, unsigned layer_id, GeoLib::Raster const &raster) override
Adds another layer to the subsurface mesh.
bool createRasterLayers(const MeshLib::Mesh &mesh, const std::vector< GeoLib::Raster const * > &rasters, double minimum_thickness, double noDataReplacementValue=0.0) override
void addLayerBoundaries(const MeshLib::Mesh &layer, std::size_t nLayers)
Creates boundary surfaces between the mapped layers to make the volumes watertight.
const T * getCoords() const
std::size_t searchByElementType(MeshElemType eleType)
Marks all elements of the given element type.
const std::vector< std::size_t > & getSearchedElementIDs() const
return marked elements
virtual unsigned getNumberOfBaseNodes() const =0
virtual void setNode(unsigned idx, Node *node)=0
static bool layerMapping(MeshLib::Mesh const &mesh, const GeoLib::Raster &raster, double nodata_replacement=0.0, bool const ignore_nodata=false)
std::vector< Node * > const & getNodes() const
Get the nodes-vector for the mesh.
unsigned getDimension() const
Returns the dimension of the mesh (determined by the maximum dimension over all elements).
std::vector< Element * > const & getElements() const
Get the element-vector for the mesh.
std::size_t getNumberOfNodes() const
Get the number of nodes.
double norm(MatrixOrVector const &x, MathLib::VecNormType type)
MeshLib::Mesh * removeElements(const MeshLib::Mesh &mesh, const std::vector< std::size_t > &removed_element_ids, const std::string &new_mesh_name)
std::vector< Node * > copyNodeVector(const std::vector< Node * > &nodes)
Creates a deep copy of a Node vector.
std::size_t getNodeIndex(Element const &element, unsigned const idx)
std::vector< Element * > copyElementVector(std::vector< Element * > const &elements, std::vector< Node * > const &new_nodes, std::vector< std::size_t > const *const node_id_map)
MeshLib::Node getCenterOfGravity(Element const &element)
Calculates the center of gravity for the mesh element.