How to use pymatgen to generate the potential energy surface of NEB calculation?

Dear all,
I’m working on the VASP+VTST calculation of transition state using CI-NEB, I want to plot the energy surface contour for MEP, but can not find a useful method. Can pymatgen implement this function? Any help would be great appreciated.
Thank you