Visualize FEM results in julia

Hello everyone. I wrote a spectral element solver in julia, but cant visualize my results for unstructured domains.
I have x,y coordinates from a classical finite element mesh with corresponding values. I would like to plot these with a heatmap, but all plot packages i have seen so far only plot heatmaps in rectangular shaped matrices. Is there any way, i can plot the heatmaps for all elements in a unstructered grid?

Check the viz function in Meshes.jl. The GDSJL book is a good resource to learn more advanced features.