Looking for adjacency municipalities.

municipality_adjacency.py

import geopandas as gpd

df = gpd.read_file('municipalities.shp')

df["adjacents"] = None

for index, municipality in df.iterrows():
    # get 'not disjoint' countries
    # CVE_ENT is the ID of the Municipality
    adjacents = df[~df.geometry.disjoint(municipality.geometry)].CVE_ENT.tolist()
    # remove own name from the list
    adjacents = [name for name in neighbors if municipality.CVE_ENT != name]
    # add names of adjacents as adjacents value
    df.at[index, "adjacents"] = ", ".join(adjacents)