Abstract: The method includes identifying a set of points defining a gamut for a color image or image forming device, selecting a parameter &agr;, and computing the gamut boundary of the color image from the set of points and &agr; using Delaunay triangulation. The system includes a processing unit, a memory, and software operative on the processing unit. The software is capable of identifying a gamut boundary for the image, identifying a gamut boundary for the image forming unit, and calculating the percentage of colors in the image capable of being rendered by the image forming unit from the gamut boundary for the image and the image forming unit. The software is also capable of identifying an alpha shape associated with a color gamut and calculating a color gamut volume. The alpha shape includes a plurality of tetrahedra and each of the plurality of tetrahedra has a volume. The plurality of volumes are summed to generate the color gamut volume.