How exactly does Emmeans calculate Estimated Means and how to replicate that in Python?

I’m looking for/want to implement a Python-only version of pymer4, especially Linear Mixed-Effects Models.