2D-hoogtebepaling

def mesh2d_altitude_from_raster(
    network,
    rasterpath,
    where: RasterStatPosition = "face",
    stat="mean",
    fill_option: FillOption = "fill_value",
    fill_value=None,
):
    """
    Method to determine level of nodes

    This function works faster for large amounts of cells, since it does not
    draw polygons but checks for nearest neighbours (voronoi) based
    on interpolation.

    Note that the raster is not clipped. Any values outside the bounds are
    also taken into account.
    
    Arguments:
      network: network object containing the mesh
      rasterpath: path to the rasterfile containing elevations
      where: 'face' or 'node'; default = face
      stat: rasterstat statistic to be used (default = mean)
      fill_option: "fill-value", 'nearest', or 'interpolate'. Fill missing cells with a constant value, by nearest-neighbour interpolation, or by linear interpolation.
      fill_value: value to fill missing cells with (also needed for interpolation in case that fails)
    """