这篇专利文献(CN102385745A)提出了一种单源伪双能探测器安检数据的伪彩配色方法。
需要首先明确的是:这篇专利中的"颜色表"并不是一张预先画好的静态图片(如传统的256色查找表 LUT) ,而是一套基于 HSV 色彩空间、非线性映射以及概率隶属度加权的"动态颜色生成算法"。系统通过实时计算每个像素的物理属性,动态"渲染"出最终的 RGB 颜色。
以下是该文献中"制作/生成颜色表"的核心逻辑和具体步骤解析:
一、 核心思想:从"硬边界"到"软过渡"
传统的双能X射线配色通常是"一刀切"的:算出来是有机物就涂橙色,是无机物就涂蓝色。这会导致物质边界处出现生硬的色块,且无法表现混合物。
该专利的"颜色表"制作核心在于:先分别计算出有机物(橙)、混合物(绿)、无机物(蓝)三种基础色的HSV值,然后根据该像素属于这三类物质的"隶属度(概率)",将三种颜色的RGB值进行加权融合。
二、 具体配色规则与映射步骤
系统对每个像素点生成颜色的过程分为以下四个关键步骤:
1. 亮度(V值)的计算:防止过曝
首先对原始的高低能信号进行归一化对数变换,得到亮度值 VVV(即 normLogV):
- 公式 :V=log(x+1)/log(3951)V = \log(x+1) / \log(3951)V=log(x+1)/log(3951)
- 作用:X射线信号动态范围极大,对数变换可以压缩高亮区域,防止图像局部过曝,同时保留低密度区域的细节。
- 极值处理 :如果原始信号超过满度值(>3950),直接记为黑色 ;低于本底值(<20),记为白色。
2. 色相(H值)的精细化:同色系分"深、中、浅"
为了让图像更有层次感,专利没有给有机物/无机物只分配一个固定的H值,而是根据亮度(V值)将每种颜色分为深、中、浅三个层级:
- 橙色(代表有机物,如爆炸物、毒品、塑料) :
- 深橙:H = 21 (条件:V < 0.73)
- 中橙:H = 25 (条件:0.73 ≤ V ≤ 0.85)
- 浅橙:H = 30 (条件:V > 0.85)
- 绿色(代表混合物) :
- 深绿:H = 115 (条件:V < 0.73)
- 中绿:H = 126 (条件:0.73 ≤ V ≤ 0.85)
- 浅绿:H = 140 (条件:V > 0.85)
- 蓝色(代表无机物,如金属、玻璃) :
- 深蓝:H = 213 (条件:V < 0.25)
- 中蓝:H = 217 (条件:0.25 ≤ V ≤ 0.5)
- 浅蓝:H = 221 (条件:V > 0.5)
3. 饱和度(S值)的非线性拉伸:凸显危险品
这是该专利的一个巧妙设计。低密度的违禁品(如炸药、毒品)在X射线下信号很弱,如果按线性映射,它们的颜色会灰蒙蒙的,看不清。因此,专利设计了非线性映射函数来强行提高特定区间的饱和度(色彩鲜艳度):
- 对于橙色和蓝色 :采用分段非线性函数(
NaLEOB),在接近0和接近1的区域进行拉伸,使得低密度和高密度物质的色彩对比更强烈。 - 对于绿色 :采用 S=1−1−x2S = 1 - \sqrt{1-x^2}S=1−1−x2 的映射函数,将饱和度拉伸集中在中间区域(0.5附近)。
4. 核心:基于"隶属度"的 RGB 色彩融合
经过前三步,系统已经算出了纯橙色、纯绿色、纯蓝色 在当前像素下的 HSV 值,并通过公式将它们转换为了 RGB 值 (Rorg,Gorg,BorgR_{org}, G_{org}, B_{org}Rorg,Gorg,Borg 等)。
接下来,利用 EM 算法(混合高斯模型)计算出的隶属度向量 (Porg,Pgreen,Pblue)(P_{org}, P_{green}, P_{blue})(Porg,Pgreen,Pblue)(即该像素属于有机物、混合物、无机物的概率,三者之和为1),进行最终的 RGB 加权融合:
- 最终红色分量 :Rfinal=Porg×Rorg+Pgreen×Rgreen+Pblue×RblueR_{final} = P_{org} \times R_{org} + P_{green} \times R_{green} + P_{blue} \times R_{blue}Rfinal=Porg×Rorg+Pgreen×Rgreen+Pblue×Rblue
- 最终绿色分量 :Gfinal=Porg×Gorg+Pgreen×Ggreen+Pblue×GblueG_{final} = P_{org} \times G_{org} + P_{green} \times G_{green} + P_{blue} \times G_{blue}Gfinal=Porg×Gorg+Pgreen×Ggreen+Pblue×Gblue
- 最终蓝色分量 :Bfinal=Porg×Borg+Pgreen×Bgreen+Pblue×BblueB_{final} = P_{org} \times B_{org} + P_{green} \times B_{green} + P_{blue} \times B_{blue}Bfinal=Porg×Borg+Pgreen×Bgreen+Pblue×Bblue
(注:专利原文OCR识别中的公式存在乱码,如 Grinaal=PorgGovg...,但其数学本质就是上述的线性加权融合。)
三、 这种"动态颜色表"的优势
通过上述算法"制作"出来的伪彩图像,具有传统静态颜色查找表无法比拟的优势:
- 消除硬边界(平滑过渡) :如果一个像素处于有机物和无机物的交界处,它的 PorgP_{org}Porg 可能是 0.5,PblueP_{blue}Pblue 也是 0.5。最终的颜色就是橙色和蓝色的完美混合色(可能会呈现出过渡的灰褐色或暗黄色),这在视觉上非常自然,帮助安检员识别被遮挡或形状不规则的物品边缘。
- 凸显低密度危险品:通过饱和度的非线性拉伸,即使是非常薄的塑料炸药(信号极弱),也能在图像上呈现出非常鲜艳的橙色,而不是融入背景的灰暗色中。
- 层次分明 :同一种物质(比如一块厚塑料和一块薄塑料),因为厚度(亮度V值)不同,会呈现出"深橙"到"浅橙"的渐变,安检员可以通过颜色的深浅直观判断物品的厚度和密度。

总结
这篇文献中的"颜色表"本质上是一个 "物理信号 -> 隶属度概率 -> HSV非线性映射 -> RGB概率加权" 的数学模型。它摒弃了死板的查表法,利用概率融合和非线性视觉增强,极大地提升了单源伪双能X光机对复杂行李中危险品的辨识度。