pytuflow.CATCHJson.data_point

pytuflow.CATCHJson.data_point#

CATCHJson.data_point(locations, data_types, time, averaging_method=None)#

Extracts the data value for the given point locations and data types at the specified time.

The locations can be a single point in the form of a tuple (x, y) or in the Well Known Text (WKT) format. It can also be a list of points, or a dictionary of points where the key will be used in the column name in the resulting DataFrame.

The locations argument can also be a single GIS file path e.g. Shapefile or GPKG (but any format supported by GDAL is also supported). GPKG’s should follow the TUFLOW convention if specifying the layer name within the database database.gpkg >> layer. If the GIS layer has a field called name, label, or ID then this will be used as the column name in the resulting DataFrame.

The returned value will be a single float if a single location and data type is provided, or a tuple if the data type is a vector result type. If multiple locations and/or data types are provided, a DataFrame will be returned with the data types as columns and the point locations as the index.

Parameters:

  • locations (Point | list[Point] | dict[str, Point] | PathLike) – The location to extract the data for.

  • data_types (str | list[str]) – The data types to extract the data for.

  • time (TimeLike) – The time to extract the data for.

  • averaging_method (str, optional) –

    The depth-averaging method to use. Only applicable for 3D results.

    The averaging methods are:

    • singlelevel

    • multilevel

    • depth

    • height

    • elevation

    • sigma

    The averaging method parameters can be adjusted by building them into the method string in a URI style format. The format is as follows:

    <method>?dir=<dir>&<value1>&<value2>

    Where

    • <method> is the averaging method name

    • <dir> is the direction, top or bottom (i.e. from top or from bottom) - only used by certain averaging methods

    • <value1>, <value2>… are the values to be used in the averaging method (the number required to be passed depends on the averaging method)

    e.g. 'singlelevel?dir=top&1' uses the single level averaging method and takes the first vertical layer from the top. Or 'sigma&0.1&0.9' uses the sigma averaging method and averages values located between the 10th and 90th water column depth.

Returns:

The data value(s) for the given location(s) and data type(s).

Return type:

float | tuple[float, float] | pd.DataFrame

Examples

Get the water level data for a given point defined as (x, y):

>>> res = ... # Assume mesh is a loaded CATCHJson result
>>> res.data_point((293250, 6178030), 'water level', 1.5)
42.723076

Get velocity vector data for a given point defined as (x, y):

>>> res.data_point((293250, 6178030), 'vector velocity', 1.5)
(0.282843, 0.154213)

Get the maximum water level and depth for multiple points defined in a shapefile. Time is passed as -1 since it is a static dataset (it could be any time value since it won’t affect the result):

>>> res.data_point('/path/to/points.shp', ['max water level', 'max depth'], -1)
      max water level  max depth
pnt1        40.501997   2.785571
pnt2        43.221862   3.450053
Contents