For the prototype / dev version, we treat all obstacles as cutlines. In a later stage, we may either introduce a new layer in the schematisation (clone_cell_cutline) or add a boolean attribute to the obstacle layer (is_clone_cell_cutline).
Implement the following in threedigrid-builder:
- [ ] cut the cells into clone cells w/ the cutlines
- "snap to pixel grid" the cutlines
- [ ] rasterize the cells ids into a raster ("area mask")
- minimum 1 pixels per clone cell domain
- [ ] rework the cells and lines make the flowlines (fortran) based on the area mask
- every clone cell domain needs to contain at least 1 data pixel
- reorder the 2D nodes and lines
- create flowline (type) between cloned cells (needed info: exchange level, cutline length). the kcu type of this flowline will be a new kcu type.
- [ ] write geometry of cell domains to gridadmin [extend the data format, not change it]
- [ ] write geometry of clone cell domain centerpoint (point-on-surface) to the gridadmin
- [ ] write geometry of clone cell / clone cell flowlines to gridadmin (linestring between the clone domain points (point-on-surface))
- [ ] rasterize the node id's and write this as a separate "node id mask" GeoTIFF (to be used in threedi-tables later)
For the prototype / dev version, we treat all obstacles as cutlines. In a later stage, we may either introduce a new layer in the schematisation (clone_cell_cutline) or add a boolean attribute to the obstacle layer (is_clone_cell_cutline).
Implement the following in threedigrid-builder: