MeshToolsLib::RasterDataToMesh Namespace Reference

Detailed Description

Adding pixel values from a raster onto nodes or cells of a mesh.

Functions
static bool	checkMesh (MeshLib::Mesh const &mesh)
static double	evaluatePixel (double const value, double const no_data, double const replacement)
static MathLib::Point3d const &	evaluationPoint (MeshLib::Node const &node)
	The point at which the raster is sampled for a mesh item.
static MathLib::Point3d	evaluationPoint (MeshLib::Element const &element)
template<typename MeshItems>
static bool	projectToMeshItems (MeshLib::Mesh &mesh, GeoLib::Raster const &raster, double const default_replacement, std::string const &array_name, MeshLib::MeshItemType const item_type, MeshItems const &mesh_items)
bool	projectToNodes (MeshLib::Mesh &mesh, GeoLib::Raster const &raster, double const default_replacement, std::string const &array_name)
bool	projectToElements (MeshLib::Mesh &mesh, GeoLib::Raster const &raster, double const default_replacement, std::string const &array_name)

Function Documentation

checkMesh()

bool MeshToolsLib::RasterDataToMesh::checkMesh ( MeshLib::Mesh const & mesh )

static

Definition at line 17 of file RasterDataToMesh.cpp.

{
    if (mesh.getDimension() > 2)
    {
        ERR("This functionality is currently only available for 2D meshes.");
        return false;
    }
    return true;
}

References ERR(), and MeshLib::Mesh::getDimension().

Referenced by projectToMeshItems().

evaluatePixel()

double MeshToolsLib::RasterDataToMesh::evaluatePixel ( double const value, double const no_data, double const replacement )

static

Definition at line 27 of file RasterDataToMesh.cpp.

{
    if (std::abs(value - no_data) < std::numeric_limits<double>::epsilon())
    {
        return replacement;
    }
    return value;
}

evaluationPoint() [1/2]

MathLib::Point3d MeshToolsLib::RasterDataToMesh::evaluationPoint ( MeshLib::Element const & element )

static

Definition at line 43 of file RasterDataToMesh.cpp.

{
    return getCenterOfGravity(element);
}

evaluationPoint() [2/2]

MathLib::Point3d const & MeshToolsLib::RasterDataToMesh::evaluationPoint ( MeshLib::Node const & node )

static

The point at which the raster is sampled for a mesh item.

Definition at line 38 of file RasterDataToMesh.cpp.

{
    return node;
}

projectToElements()

bool MeshToolsLib::RasterDataToMesh::projectToElements	(	MeshLib::Mesh &	mesh,
		GeoLib::Raster const &	raster,
		double const	default_replacement,
		std::string const &	array_name )

Definition at line 96 of file RasterDataToMesh.cpp.

{
    return projectToMeshItems(mesh, raster, default_replacement, array_name,
                              MeshLib::MeshItemType::Cell, mesh.getElements());
}

References MeshLib::Cell, MeshLib::Mesh::getElements(), and projectToMeshItems().

Referenced by main(), and MeshView::openRasterDataToMeshDialog().

projectToMeshItems()

template<typename MeshItems>

bool MeshToolsLib::RasterDataToMesh::projectToMeshItems ( MeshLib::Mesh & mesh, GeoLib::Raster const & raster, double const default_replacement, std::string const & array_name, MeshLib::MeshItemType const item_type, MeshItems const & mesh_items )

static

Creates a new property vector of the given item_type on mesh and fills it with the raster values sampled at the evaluation point of each of the mesh_items.

Definition at line 52 of file RasterDataToMesh.cpp.

{
    if (!checkMesh(mesh))
    {
        return false;
    }
 
    auto& props = mesh.getProperties();
    std::string const name =
        BaseLib::getUniqueName(props.getPropertyVectorNames(), array_name);
    auto* const values = props.createNewPropertyVector<double>(
        name, item_type, mesh_items.size(), 1);
    if (values == nullptr)
    {
        // MeshLib::Properties keys its property vectors by name alone, and the
        // name above is unique among the existing ones, so this is a broken
        // invariant rather than malformed input.
        OGS_FATAL("Could not create the property vector '{:s}'.", name);
    }
 
    double const no_data = raster.getHeader().no_data;
    ranges::transform(mesh_items, values->begin(),
                      [&](auto const item)
                      {
                          return evaluatePixel(
                              raster.getValueAtPoint(evaluationPoint(*item)),
                              no_data, default_replacement);
                      });
    return true;
}

References checkMesh(), GeoLib::Raster::getHeader(), MeshLib::Mesh::getProperties(), BaseLib::getUniqueName(), GeoLib::RasterHeader::no_data, and OGS_FATAL.

Referenced by projectToElements(), and projectToNodes().

projectToNodes()

bool MeshToolsLib::RasterDataToMesh::projectToNodes	(	MeshLib::Mesh &	mesh,
		GeoLib::Raster const &	raster,
		double const	default_replacement,
		std::string const &	array_name )

Definition at line 88 of file RasterDataToMesh.cpp.

{
    return projectToMeshItems(mesh, raster, default_replacement, array_name,
                              MeshLib::MeshItemType::Node, mesh.getNodes());
}

References MeshLib::Mesh::getNodes(), MeshLib::Node, and projectToMeshItems().

Referenced by main(), and MeshView::openRasterDataToMeshDialog().