如何判断XYZ三元组是否为有效物理颜色？CIE 1931色度图生成咨询

嘿，针对你提到的两个关于CIE XYZ颜色空间和1931色度图的问题，我来分享一些实用的解决方案：

真实颜色的核心定义是：它可以对应物理世界中存在的某种光谱分布，换句话说，这个XYZ必须能通过CIE标准观察者的颜色匹配函数，由非负的光谱功率分布（SPD）积分得到。具体有两种常用判断方法：

方法1：通过色度坐标xy判断（最简便） #

首先把XYZ转换为CIE 1931色度坐标：

x = X / (X + Y + Z)
y = Y / (X + Y + Z)

（注意：如果X+Y+Z=0，这代表纯黑，属于有效颜色）

接下来只需要检查(x,y)是否落在CIE 1931色度图的马蹄形色域边界内部（包括边界）。判断方法有两种：

• 预定义边界点集：获取CIE官方提供的380~780nm单色光的xy坐标列表，加上紫边界（连接780nm和380nm的直线），构成一个闭合多边形，然后用点-in-polygon算法判断目标点是否在多边形内。
• 解析拟合公式：利用已有的色域边界拟合方程（比如CIE推荐的多项式拟合），直接代入x计算对应的y上限/下限，判断目标y是否在范围内。

方法2：线性规划验证（更严谨） #

如果你需要最严谨的验证，可以把问题转化为线性规划：是否存在非负的光谱功率分布S(λ)（λ∈380~780nm），满足：

X = ∫₃₈₀⁷⁸⁰ S(λ) · x̄(λ) dλ
Y = ∫₃₈₀⁷⁸⁰ S(λ) · ȳ(λ) dλ
Z = ∫₃₈₀⁷⁸⁰ S(λ) · z̄(λ) dλ

其中x̄, ȳ, z̄是CIE 1931标准观察者的颜色匹配函数。不过这种方法计算量较大，一般只用于学术场景，日常用第一种方法足够。

你提到的随机高斯叠加光谱是可行的，但还有更高效、准确的方法：

方法1：直接绘制色域边界（最推荐） #

CIE 1931色度图的边界本身就是有明确定义的：

• 马蹄形的弧形部分：380~780nm所有单色光的xy坐标轨迹
• 直线部分（紫边界）：连接780nm（x=0.7347, y=0.2653）和380nm（x=0.1741, y=0.0050）的直线

具体步骤：

• 下载或计算380~780nm每个波长对应的xy坐标（可以通过标准观察者函数计算单色光的XYZ，再转xy）
• 按波长顺序连接这些弧形点，再连接首尾点形成紫边界，构成闭合的色域多边形
• 对多边形进行描边和内部填充，就得到完整的CIE 1931色度图

方法2：系统生成内部真实颜色点 #

如果需要填充色域内部的采样点，可以用以下方法：

• 凸组合采样：随机选取色域边界上的两个点，按任意非负比例混合（比如(x1*a + x2*(1-a), y1*a + y2*(1-a))，a∈[0,1]），得到的点一定是真实颜色的坐标
• 多单色光混合：随机选取3~5个不同波长的单色光，分配非负权重（权重和为1），计算对应的XYZ后转xy，这些点也会落在色域内部
• 你的随机高斯方法：叠加多个高斯函数模拟自然光谱（注意高斯的中心波长在380~780nm，振幅非负），再转XYZ和xy，这种方法能生成接近自然光照的真实颜色点

• 不要直接遍历所有XYZ组合：XYZ的取值范围理论上是无限的，而且大部分组合对应非真实颜色，效率极低
• 绘制时注意：CIE 1931色度图的x轴范围一般是0_{0.8，y轴是0}0.9，这样能完整展示马蹄形色域

内容的提问来源于stack exchange，提问作者Chris_F