This paper proposes a method of gamut estimation using data segmentation and 2D surface splines. The device data is first segmented into hue intervals, and then each hue interval is analyzed iteratively to construct a 2D gamut boundary descriptor for that hue interval. The accuracy and smoothness of the gamut boundary estimate can be controlled by the number of hue intervals selected and the sampling within each hue interval.
Mark Shaw, Mark Shaw,
"Gamut estimation using 2D surface splines", Proc. SPIE 6058, Color Imaging XI: Processing, Hardcopy, and Applications, 605807 (16 January 2006); doi: 10.1117/12.660086; https://doi.org/10.1117/12.660086