Richard Tran has written an excellent Jupyter notebook on how you can use pymatgen to automatically generate surface slabs and analyzing calculated surface energies to construct the Wulff shape. Check it out here.
matgenb – Generating surface slabs and Wulff shapes
