Reading marker of cross sections

[hendryliu]

Is there any functionality to read the markers of cross sections, e.g. left levee bank and lowest point?

[gedaskir]

Here is an example script to do that:

import mikeio1d
from mikeio1d.res1d import Res1D

filename = "MyResults.res1d"
res1d = Res1D(filename)

reaches = list(res1d.data.Reaches)

for reach in reaches:
    grid_points = list(reach.GridPoints)

    for grid_point in grid_points:
        try:
            cross_section = grid_point.CrossSection

            points = list(cross_section.Points)
            left = points[cross_section.LeftLeveeBank]
            lowest = points[cross_section.LowestPoint]
            right = points[cross_section.RightLeveeBank]

            print("Gridpoint:", reach.Name, grid_point.Chainage)
            print("  Cross section ID:", cross_section.ID)
            print("    Left levee bank coordinates:", left.X, left.Z)
            print("    Lowest point coordinates:", lowest.X, lowest.Z)
            print("    Right levee bank coordinates:", right.X, right.Z)

        except:
            pass

[hendryliu]

The script doesn't work now with this error:

---> 14 cross_section = grid_point.CrossSection
AttributeError: 'IRes1DGridPoint' object has no attribute 'CrossSection'

Is there any solution?
Thanks.

[JesperGr]

Not all grid points have cross sections, only every second grid point has cross sections. I guess that is what the try: Except: pass tries to handle. There is a few other alternatives:

• Check if it is an H-gridpoint, alternatively if it has the CrossSection property.
• Loop over odd grid points, in range(0,count,2).

Also, some reaches do not have cross sections at alle, so do check the CrossSection property for null values.

[ryan-kipawa]

MIKEIO 1D depends on the Python.NET package that had a major release since this script was written (v3.0 and above). I think this is relevant here, since the error is noting an Interface ('IRes1DGridPoint') instead of an actual object. In addition to what @JesperGr mentioned, you could test this as a fix:

It's a bit messy but should work temporarily. We'll improve this in future versions of MIKEIO 1D.

from mikeio1d.dotnet import pythonnet_implementation as impl

reaches = list(impl(res1d.data.Reaches))

for reach in reaches:
    reach = impl(reach)
    grid_points = list(impl(reach.GridPoints))

    for grid_point in grid_points:
        grid_point = impl(grid_point)
        try:
            print(type(grid_point))
            cross_section = impl(grid_point.CrossSection)

            points = list(impl(cross_section.Points))
            left = points[impl(cross_section.LeftLeveeBank)]
            lowest = points[impl(cross_section.LowestPoint)]
            right = points[impl(cross_section.RightLeveeBank)]

            print("Gridpoint:", impl(reach.Name), impl(grid_point.Chainage))
            print("  Cross section ID:", cross_section.ID)
            print("    Left levee bank coordinates:", left.X, left.Z)
            print("    Lowest point coordinates:", lowest.X, lowest.Z)
            print("    Right levee bank coordinates:", right.X, right.Z)

        except Exception as e:
            pass