Coloring the projection of the structures on a PCA chart based on occupancy and time frame

GROMACS version:2019
GROMACS modification: No

Hi all, I was wondering if you know of any way to plot the dynamic evolution of the conformers generated from a MD trajectory on a PCA chart based on colors. I have already obtained the 2D projections according to the gmx covar and gmx anaeig commands. I want to know if there is a way, by using the files generated on Gromacs, to generate:

  1. a density plot of the conformers (from non occupancy to high occupancy)
  2. color de conformers according to the time frame of each conformer (from first to last)