longavailable
commented
1 year ago Firstly, it is a very interesting feedback. Many thanks!
How did it happen?
Voronoi diagram for polygons is based on shapely.ops.voronoi_diagram.
As it mentioned in the mannual:
All vertices of the geometry will be used as the input points to the diagram.
Let us have a look on the geometries from your inputs.
Simple geometries with TOO less vertices!
Then it should be added more vertices along their edges.
I added a few to functions to reach it. See v0.1.3.
When you type:
vd = voronoiDiagram4plg(builtup, boundary)
It return the same output as you posted.
And when:
vd = voronoiDiagram4plg(input, boundary, densify=True)
It will proceed a automatical densifying preprocessing. And it will return the following.
It looks much better!
You can also assign an expected spacing to spacing parameter, for example:
vd = voronoiDiagram4plg(input, boundary, densify=True, spacing=5000)
At last, when you post a specific issue, it 's better to upload the related attachments. That will be time-saving for a reviewer.
I have generated the voronoi regions (in red) for the geometries in blue. As it can be seen, the resulting voronoi regions interset with the original geometries. Any idea?