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)
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 14 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 projectToElements(), and projectToNodes().

evaluatePixel()

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

static

Definition at line 24 of file RasterDataToMesh.cpp.

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

projectToElements()

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

Definition at line 59 of file RasterDataToMesh.cpp.

{
    if (!checkMesh(mesh))
    {
        return false;
    }
 
    auto& elems = mesh.getElements();
    auto& props = mesh.getProperties();
    std::string const name =
        BaseLib::getUniqueName(props.getPropertyVectorNames(), array_name);
    auto vec = props.createNewPropertyVector<double>(
        name, MeshLib::MeshItemType::Cell, 1);
    double const no_data = raster.getHeader().no_data;
    std::transform(elems.cbegin(), elems.cend(), std::back_inserter(*vec),
                   [&](auto const elem)
                   {
                       auto node = getCenterOfGravity(*elem);
                       return evaluatePixel(raster.getValueAtPoint(node),
                                            no_data, default_replacement);
                   });
    return true;
}

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

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

projectToNodes()

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

Definition at line 34 of file RasterDataToMesh.cpp.

{
    if (!checkMesh(mesh))
    {
        return false;
    }
 
    auto& nodes = mesh.getNodes();
    auto& props = mesh.getProperties();
    std::string const name =
        BaseLib::getUniqueName(props.getPropertyVectorNames(), array_name);
    auto vec = props.createNewPropertyVector<double>(
        name, MeshLib::MeshItemType::Node, 1);
    double const no_data = raster.getHeader().no_data;
    std::transform(nodes.cbegin(), nodes.cend(), std::back_inserter(*vec),
                   [&](auto const node)
                   {
                       return evaluatePixel(raster.getValueAtPoint(*node),
                                            no_data, default_replacement);
                   });
    return true;
}

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

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