*2D-hoogtebepaling

De hoogte van het 2D-grid kan met twee functies bepaald:

altitude_constant: Ken een constante hoogte aan de rekencellen. Dit kan zowel op de faces (cellen) als nodes (hoekpunten) gedaan worden, al zal dat in praktijk niet uitmaken bij een constante hoogte.
altitude_from_raster: Bepaal de hoogte uit een hoogteraster. Op basis van het grid, het gekozen hoogteraster een aantal instellingen wordt de hoogte van het rekenrooster bepaald. De docstring van de functie verduidelijkt dit (Engels):

def altitude_from_raster(self, rasterpath, where='face', stat='mean', missing='default'):
    """
    Method to determine level within cell or on nodes. The values are determined by
    applying a statistic to the pixels within the cell bounds or around the node.

    In case of the option 'node' Voronoi polygons are drawn around the nodes. These cells
    are cut of at the edges of the the grid. This option might take a bit longer, since
    all the polygons need to be drawn and clipped at the edges.

    In case of msising values, which can occur:
    - due to no-data parts in the grid that is sampled
    - when the cell sizes within which the altitude is determined is smaller than a raster pixel.
    The missing data can be filled.

    Parameters
    ----------
    rasterpath : str
        Path to raster
    where : str
        Locations where the altitude is determined. Can be on the faces, so within the
        cell boundaries, or node, on the cell edged. Default: 'face'
    stat : str
        Statistic to determined from values within polygon bounds. A string is required that
        describes a function that is known by numpy, such as 'mean' or 'max', without any
        further arguments (quantile is not possible, since we'd need to specify which quantile).
        Default: 'mean'
    missing : str
        How to fill the missing values.
        - default: No filling, the missing values will have a NaN-value in the grid
        - nearest: Fill the missing data with the nearest cell value that has a value
        - interpolation: Interpolate the missing data with cell values that are present
            Default: 'default'
    """

Beide functies zijn een ‘method’ van de class Mesh2D.

Met de functie set_missing_z_value kan een waarde worden gespecificeerd voor de cellen waar geen waarde bekend is. Deze waarde wordt in de mdu bij BedLevUni ingevuld.

Merk op dat de keuze voor de hoogtebepaling op de nodes of faces niet los gezien kan worden van de rekeninstellingen in het mdu-bestand. Zie (op het moment van schrijven) paragraaf 8.4.2.3 - “Conveyance in 2D” en paragraaf 8.8.1 - “Definitions” in de handleiding.