Merge pull request #583 from gdsfactory/remove-gdstk-dependency

removing gdstk dependencies anywhere other than the tidy3d module

Author: tvt173

gplugins/common/utils/get_component_with_local_layers.py

@@ -1,9 +1,9 @@
 from functools import lru_cache
 
 import gdsfactory as gf
-import gdstk
 import numpy as np
 import shapely.geometry as sg
+import shapely.ops as so
 from natsort import natsorted
 
 
@@ -225,10 +225,12 @@ def get_icon_poly(name):
     polys = min(polys_priority.items(), key=lambda x: x[0])[1]
     if not polys:
         return default
-    polys = gdstk.boolean(polys, [], operation="or")
+    polys = [sg.Polygon(p) for p in polys]
+    polys = so.unary_union(polys)
+    polys = _extract_vertices(polys)
     if not polys:
         return default
-    poly = rdp(polys[0].points)
+    poly = rdp(polys[0])
     if (poly.shape[0] < 3) or (poly.shape[0] > 100):
         return default
     poly = (poly - c.bbox_np()[0:1]) / (c.bbox_np()[1:2] - c.bbox_np()[0:1])
@@ -242,6 +244,14 @@ def get_icon_poly(name):
     return poly
 
 
+def _extract_vertices(geometry):
+    if geometry.geom_type == 'Polygon':
+        return [list(geometry.exterior.coords)]
+    elif geometry.geom_type == 'MultiPolygon':
+        return [list(p.exterior.coords) for p in geometry.geoms]
+    else:
+        raise TypeError(f"Unhandled geometry type: {geometry.geom_type}")
+
 def rdp(poly, eps=0.1):
     poly = np.asarray(poly)
     if not poly.shape:

gplugins/gfviz/netlist.py

@@ -164,7 +164,7 @@ def write_sparameters_meep_mpi(
     with open(parameters_file, "wb") as outp:
         pickle.dump(settings, outp, pickle.HIGHEST_PROTOCOL)
 
-    # Save component to disk through gds for gdstk compatibility
+    # Save component to disk through gds
     component_file = tempfile.with_suffix(".gds")
     component.write_gds(component_file, with_metadata=True)
 

gplugins/gmeep/write_sparameters_meep_mpi.py

@@ -15,7 +15,6 @@ classifiers = [
 dependencies = [
   "gdsfactory>=9.3.5",
   "pint",
-  "gdstk",
   "tqdm",
   "numpy",
   "xarray==2025.1.2"
@@ -71,7 +70,7 @@ sax = [
   "ray"
 ]
 schematic = ["bokeh", "ipywidgets", "natsort"]
-tidy3d = ["tidy3d>=2.8.2,<2.9", "meshio", "numpy", "meshwell~=1.0.7"]
+tidy3d = ["tidy3d>=2.8.2,<2.9", "meshio", "numpy", "meshwell~=1.0.7", "gdstk"]
 vlsir = ["vlsir", "vlsirtools"]
 
 [tool.codespell]