One more big benefit after adding Lloyd’s algorithm to even the spacing of inhibitory centers was that they also get pushed away from the primary veins. At one point, I thought I would have to do this manually. One possibility was during the dart throwing algorithm, to reject points that were too close to the primary veins but that problem has solved itself.
TopologyException #
Traceback (most recent call last):
File "/home/nmclaren/Code/InsectWings/./main_orthoptera.py", line 203, in <module>
randomize_parameters()
File "/home/nmclaren/Code/InsectWings/./main_orthoptera.py", line 97, in randomize_parameters
vein_renderer = VeinRenderer(parameters)
^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/nmclaren/Code/InsectWings/vein_renderer.py", line 27, in __init__
self._interveinal_regions = self._get_interveinal_regions(self._root_segments, parameters)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/nmclaren/Code/InsectWings/vein_renderer.py", line 135, in _get_interveinal_regions
result.append(InterveinalRegionRenderer(self._get_polygon(prev_segment, segment), parameters))
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/nmclaren/Code/InsectWings/interveinal_region_renderer.py", line 13, in __init__
self._voronoi_polygons = self._get_voronoi_polygons(self._inhibitory_centers, self._polygon)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/nmclaren/Code/InsectWings/interveinal_region_renderer.py", line 80, in _get_voronoi_polygons
return intersection(normalize(voronoi_polygons(inhibitory_centers, extend_to=extent)).geoms, extent)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/nmclaren/Code/InsectWings/venv/lib/python3.12/site-packages/shapely/decorators.py", line 77, in wrapped
return func(*args, **kwargs)
^^^^^^^^^^^^^^^^^^^^^
File "/home/nmclaren/Code/InsectWings/venv/lib/python3.12/site-packages/shapely/set_operations.py", line 131, in intersection
return lib.intersection(a, b, **kwargs)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
shapely.errors.GEOSException: TopologyException: side location conflict at 351.24287691041724 346.75142461791654. This can occur if the input geometry is invalid.
What I really need to do is separate the generation of Voronoi lines from generation of the primary veins. Then I can see what the primary veins look like before the error tosses me out of the program.
What’s the situation with a debugger for Python?
Hey, here we go: https://marketplace.visualstudio.com/items?itemName=ms-python.debugpy
The purpose of that extension is more about compatibility issues. It sounds like the main Python extension can do what I’m looking for.
The main intent of this extension is to offer:
1. **Independence and Compatibility:** The Python Debugger extension aims to separate the debugging functionality from the main Python extension to prevent compatibility issues.
Actually. I’ve already got both extensions installed.
I just stopped the code from generating the cross veins if it detects a collision between primary veins. That solved the problem.
Now I can look at the cross veins for different wing shapes:
I’m still getting the occasional error as I generate different wings using parameter randomization.
Traceback (most recent call last):
File "/home/nmclaren/Code/InsectWings/./main_orthoptera.py", line 203, in <module>
randomize_parameters()
File "/home/nmclaren/Code/InsectWings/./main_orthoptera.py", line 97, in randomize_parameters
vein_renderer = VeinRenderer(parameters)
^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/nmclaren/Code/InsectWings/vein_renderer.py", line 29, in __init__
self._interveinal_regions = self._get_interveinal_regions(self._root_segments, parameters)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/nmclaren/Code/InsectWings/vein_renderer.py", line 139, in _get_interveinal_regions
result.append(InterveinalRegionRenderer(self._get_polygon(prev_segment, segment), parameters))
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/nmclaren/Code/InsectWings/interveinal_region_renderer.py", line 13, in __init__
self._inhibitory_centers = self._lloyds_algorithm(inhibitory_centers, 50)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/nmclaren/Code/InsectWings/interveinal_region_renderer.py", line 69, in _lloyds_algorithm
voronoi_polygons = self._get_voronoi_polygons(inhibitory_centers, self._polygon)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/nmclaren/Code/InsectWings/interveinal_region_renderer.py", line 64, in _get_voronoi_polygons
return intersection(normalize(voronoi_polygons(inhibitory_centers, extend_to=extent)).geoms, extent)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/home/nmclaren/Code/InsectWings/venv/lib/python3.12/site-packages/shapely/decorators.py", line 77, in wrapped
return func(*args, **kwargs)
^^^^^^^^^^^^^^^^^^^^^
File "/home/nmclaren/Code/InsectWings/venv/lib/python3.12/site-packages/shapely/set_operations.py", line 131, in intersection
return lib.intersection(a, b, **kwargs)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
shapely.errors.GEOSException: TopologyException: side location conflict at 299.64486182880285 356.84611036182139. This can occur if the input geometry is invalid.
Says the geometry might be invalid. I can just skip wings with invalid polygons..?
That didn’t work at all. I still saw the same error. So instead I’m catching the GEOSException but I haven’t been able to verify because I’ve been unlucky lucky with my parameter randomization RNG.
Next Steps #
I want to randomize other primary vein parameters. Right now all I’m randomizing are the curves.
Also, I want to filter out any results where the primary veins extend outside a certain rectangle. Just ignore anything that goes outside my viewport.
And finally, I want to filter out results where the interveinal space goes from very narrow to very very wide.
Those two should take care of a lot of the displeasing wings I’m seeing. Reduce the number of rejections I have to make. (Again, it would be pretty cool to map the parameter space and see if there’s a pattern to this that I can exploit…)
One more thing. This process can be improved if I generate primary veins and cross veins separately. That way I don’t have to wait the extra several hundred milliseconds while the cross veins are generated before hitting “next.”
Separate Primary and Cross Vein Generation #
Okay, I had enough time to do the last one at least. Much more convenient.
Plus I can try different variations of cross veins for a given set of primary veins.
