[发明专利]基于平移评价优选的光谱图像数据融合方法

摘要

本发明公开了一种基于平移评价优选的光谱图像数据融合方法，属于超光谱图像数据融合技术领域。本发明首先对数据矩阵的像素向量进行分类，得到类矩阵；采用类矩阵主成分变换的降维投影矩阵，对数据矩阵进行投影，获得降维数据矩阵；再分配给对色空间黑白、红绿和黄蓝通道，转换至sRGB颜色空间，获得R、G、B通道颜色分量数据；通过平移映射颜色分量数据，综合评价系列数字图像，获取每个通道的最优平移量；最后利用三个最优平移量分别平移映射各通道颜色分量数据，得到融合结果图像。本发明采用最优平移量融合的方法，获得了高亮度的伪彩色融合图像，提高了图像的对比度和色彩区分度，便于人眼从复杂背景中快速识别目标。

主权项

一种基于平移评价优选的光谱图像数据融合方法，其特征在于：包括以下步骤：步骤1：载入数据矩阵：读取光谱数据立方体，并整合成二维数据矩阵X(i,j)，其中i＝1、2、…、P,j＝1、2、…、L，P为图像的像素数，其值等于图像高度H与图像宽度W的乘积，L为光谱波段数；步骤2：数据矩阵像素分类获取类矩阵：对数据矩阵X(i,j)的像素向量X(i,:)＝[X(i,1) X(i,2) ... X(i,L)]进行分类，从每类像素向量中提取一个像素向量构成类矩阵Y(k,j)，其中i＝1、2、…、P,j＝1、2、…、L，k＝1、2、…、Pc，Pc为像素向量分类后的类数；步骤3：类矩阵主成分变换获取降维数据矩阵：利用类矩阵Y(k,j)光谱维主成分变换得到的降维投影矩阵Tn(j,n)，投影变换数据矩阵X(i,j)的光谱向量X(:,j)＝[X(1,j) X(2,j) ... X(P,j)]T，得到降维数据矩阵Xn(i,n)，其中k＝1、2、…、Pc；j＝1、2、…、L，n＝1、2、3，i＝1、2、…、P；步骤4：降维数据矩阵经过对色空间转换到sRGB颜色空间，得到R通道颜色分量数据Rs(:)、G通道颜色分量数据Gs(:)、B通道颜色分量数据Bs(:)；步骤5：采用R通道颜色分量数据Rs(:)和R通道平移步长Δmr，通过平移评价优选方法，得到R通道的最优平移量MR；采用G通道颜色分量数据Gs(:)和G通道平移步长Δmg，通过平移评价优选方法，得到G通道的最优平移量MG；采用B通道颜色分量数据Bs(:)和B通道平移步长Δmb，通过平移评价优选方法，得到B通道的最优平移量MB；所述平移评价优选方法包括以下步骤：步骤5.1：将输入本方法的颜色分量数据存入数组V(j)，j＝1、2、…、P；求出颜色分量数据数组的最大值Vmax和最小值Vmin，计算平移总步数Nm＝[(Vmax‑Vmin)/Δm]+1，其中Δm为输入本方法的平移步长；初始化以下数组为0：平移量M(i)，标准差Ed(i)，熵Ee(i)，平均梯度Eg(i)，标准差分值Sd(i)，熵分值Se(i)，平均梯度分值Sg(i)，综合评价值Sa(i)，其中i＝1、2、…、Nm；步骤5.2：循环标记i等于1，开始平移循环；步骤5.3：判断循环标记i是否小于等于平移总步数Nm，如果为假，跳出平移循环，执行步骤5.6；如果为真，执行步骤5.4；步骤5.4：采用“分段平移、数字映射、单项评价”方法，得到平移量M(i)、标准差Ed(i)、熵Ee(i)、平均梯度Eg(i)；所述“分段平移、数字映射、单项评价”方法包括以下步骤：步骤5.4.1：计算平移量M(i)＝Vmin+(i‑1)×Δm，得到平移后的颜色分量数据Vm(j)＝V(j)‑M(i)，且j＝1、2、…、P；步骤5.4.2：数字映射平移后的颜色分量数据Vm(:)的每个像素，得到平移后的数字图像Vb(:)，即Vb(j)＝Map(Vm(j))，且j＝1、2、…、P，其中，函数y＝Map(x)为数字映射函数；步骤5.4.3：计算平移后的数字图像Vb(:)的单项评价值：标准差Ed(i)、熵Ee(i)、平均梯度Eg(i)：计算标准差，其计算公式为：其中V‾b=Σj=1PVb(j)/P;]]>计算熵，其计算公式为：q(j)表示灰度值为j的像素数nj与图像的像素数P之比；计算平均梯度，首先将平移后的数字图像Vb(:)转换为二维数组形式V2b(m,n)＝Vb(m×W+n)，其中m＝1、2、…、H，n＝1、2、…、W；然后计算平均梯度，公式为：Eg(i)=1(H-1)(W-1)Σm=1H-1Σn=1W-1(V2b(m,n+1)-V2b(m,n))2+(V2b(m+1,n)-V2b(m,n))2;]]>步骤5.5：循环标记i加1，返回执行5.3步骤；步骤5.6：采用单项打分方法，分别对标准差Ed(:)、熵Ee(:)和平均梯度Eg(:)进行单项打分，得到标准差分值Sd(:)、熵分值Se(:)和平均梯度分值Sg(:)；所述单项打分方法包括以下步骤：步骤5.6.1：初始化Nm行3列的二维数组Buf(i,j)和Out(i,j)，其中i＝1、2、…、Nm，j＝1、2、3,二维数组Buf(i,j)和Out(i,j)的第一列Buf(:,1)、Out(:,1)存储单项评价值的序号，第二列Buf(:,2)、Out(:,2)存储单项评价值，第三列Buf(:,3)、Out(:,3)存储单项评价值的分值；步骤5.6.2：对单项评价值即Buf(:,2)进行排序，并且单项评价值的序号Buf(:,1)随着Buf(:,2)顺序的变化而变化；步骤5.6.3：取Buf(1,3)等于1，循环标记i等于2，开始循环；步骤5.6.4：判断i是否小于等于Nm，如果为假执行步骤5.6.6；如果为真，再判断Buf(i,2)是否等于Buf(i‑1,2)，如果为真，取Buf(i,3)等于Buf(i‑1,3)，如果为假Buf(i,3)等于i；步骤5.6.5：循环标记i加1，并返回执行步骤5.6.4；步骤5.6.6：循环标记i等于1；步骤5.6.7：判断i是否小于等于Nm，如果为真，Out(Buf(i,1),1)等于Buf(i,1)，Out(Buf(i,1),2)等于Buf(i,2)，Out(Buf(i,1),3)等于Buf(i,3)；如果为假，执行5.6.9步骤；步骤5.6.8：循环标记i加1，并返回执行步骤5.6.7；步骤5.6.9：输出Out数组的第3列数据Out(:,3)作为对应单项评价值的分值；步骤5.7：计算综合评价值Sa(i)＝Sd(i)+Se(i)+Sg(i),其中i＝1、2、…、Nm；步骤5.8：将最大综合评价值Smax初始化为Sa(1)，最优平移量Mo初始化为M(1)，循环标记i等于1，开始循环；步骤5.9：判断i是否小于等于Nm，如果为假，跳出循环，执行步骤5.11；如果为真，再判断Sa(i)是否大于等于Smax，如果为真，最大综合评价值Smax等于Sa(i)，最优平移量Mo等于M(i)，如果为假，执行步骤5.10；步骤5.10：循环标记i加1，并返回执行5.9步骤；步骤5.11：输出最优平移量Mo；步骤6：融合结果图像：根据R通道的最优平移量MR，平移R通道颜色分量数据Rs(:)；根据G通道的最优平移量MG，平移G通道颜色分量数据Gs(:)；根据B通道的最优平移量MB，平移B通道颜色分量数据Bs(:)，并分别数字映射，得到最终的融合图像Im(i,1:3)，其计算公式如下：Im(i,1)＝Map(Rs(i)‑MR)；Im(i,2)＝Map(Gs(i)‑MG)；Im(i,3)＝Map(Bs(i)‑MB)；其中i＝1、2、…、P，y＝Map(x)为数字映射函数。