The following Sage code is working perfectly. It generates the polyhedron from a vertex list of interest and computes the Ehrhart polynomial.
p = Polyhedron(vertices = vertex_list)
p = p.ehrhart_polynomial(engine = 'latte')
How can I now compute the number of integer lattice points inisde the convex hull of the polyhedron? Moreover, can I enumerate them?