search

DHI / mikeio

Read, write and manipulate dfs0, dfs1, dfs2, dfs3, dfsu and mesh files.
https://dhi.github.io/mikeio
BSD 3-Clause "New" or "Revised" License
136 stars 53 forks source link

Fix GeometryFM.contains (optional use of shapely) #700

Open ecomodeller opened 1 month ago

jsmariegaard commented 4 weeks ago

I guess it is not possible to use _point_in_polygon() ?

ecomodeller commented 4 weeks ago

I guess it is not possible to use _point_in_polygon() ?

I suppose we can, but that is only part of the solution, we still need to construct a proper polygon with holes for the entire domain, unless we loop over all elements, which we could do in a clever way similar to how find_index works (when #685 has been fixed) .