Build a hexagon grid over a shapely geometry

grid.py

import geopandas as gpd
from shapely.geometry import shape

CB = None # replace by a shapely geometry object.

r=0.2 # radius of the hexagon
h = (r * np.sqrt(3))/2
j = r/2


hex_array=[]

y= ymin
flip = True
while y <ymax:
    x=xmin
    if flip:
        x+=r+r/2
    while x < xmax:
        hexagon = shape(
            {
                "type": "Polygon",
                "coordinates": [
                    [
                        [x +j, y - h],
                        [x + r, y],
                        [x+j,y+h],
                        [x-j, y + h],
                        [x -r, y],
                        [x - h / 2, y-h]
                    ]
                ],
            }
        )
        x+=r*3
        
        
        if CB.intersects(hexagon):
            hex_array.append(hexagon)
    y+=h
    flip=not flip
        
gpd.GeoDataFrame({"geometry":hex_array}).boundary.plot()