基于接触式图像传感器的多谱段图像重建方法及系统的利记博彩app
【专利摘要】本发明提供一种基于接触式图像传感器的多谱段图像重建方法及系统，主要包括前去噪阶段、重建阶段和后去噪阶段，所述前去噪阶段和后去噪阶段，分别对每一个谱段的图像进行切块处理，非局部自相似块搜索，形成图像矩阵，根据特征矩阵与保留的特征值得到去噪后的图像矩阵，恢复得到完整的图像；所述重建阶段包括对各个谱段图像纵向插值放大，根据传感器发光顺序对n个谱段的图像进行排序，对各个谱段图像切分为小块，对每一个小块求其直方图分布，得到直方图分布集合，进行直方图拉伸，将同谱段的处理后的图像小块按原顺序拼在一起，得到所有谱段的处理后的完整图像。本发明可以充分利用各谱段信息，重建过程仅需进行线性映射，效率高，适用性好。
【专利说明】
基于接触式图像传感器的多谱段图像重建方法及系统
技术领域
[0001] 本发明涉及线阵相机图像处理领域，具体涉及基于接触式图像传感器的多谱段图 像处理及重建方法及系统。
【背景技术】
[0002] 随着信息技术的发展，机器视觉技术在检测领域的应用越来越广、可靠性越来越 高、检测速度越来越快，已经成为众多领域不可或缺的技术手段。其中在表面检测领域中， 接触式图像传感器得到了广泛的应用。然而，由于表面检测对待测的图像分辨率要求很高， 分辨率的高低直接影响了表面检测技术的效果，但分辨率的提高也同时意味着成本的上 升，而这显然与采用接触式图像传感器控制成本的初衷相悖。为了解决这一问题，本专利探 寻了一种结合接触式图像传感器成像方式进行多谱段图像重建从而获得更高分辨率图像 的方法，力求在同等的分辨率情况下获得更高分辨率图像。现有的对图像进行重建的专利 (包括发明授权专利和发明公布专利)如下：
[0003] 申请号201010238674.7, "一种高分辨率视频图像重建方法"，此方法的思想是对 视频中的图像进行插值，并利用运动矢量集对被插值点进行运动补偿，从而获得重建后的 高分辨率视频图像。此方法需要多帧图像提供互补信息，从而实现单帧图像的重建，然而在 表面检测的应用中无法获取单一被测表面的多帧视频图像，每一被测表面仅有一张图像， 因此该方法不适用于接触式图像传感器采集的多谱段图像的重建。
[0004] 申请号201310522407.6, "图像重建方法及装置"，此方法的思想是采用压缩感知 的思想，利用观测矩阵与原图像之间的关系来重建原图像，获取更高分辨率的图像。此方法 需使用特定的装置采集图像，不适用于接触式图像传感器采集的多谱段图像的重建。
[0005] 申请号201410212391.0, "低分辨率红外图像重建高分辨率红外图像的方法"，此 方法的思想是利用同场景高分辨率的可见光图像的高频信息重建对应的红外图像，使得重 建后的红外图像拥有高频细节及清晰的边缘。此方法需要低分辨率的红外图像拥有对应的 高分辨率可见光图像，然而接触式图像传感器所获取的所有谱段图像均为低分辨率图像， 无法提供高频细节，因此该方法不适用于接触式图像传感器采集的多谱段图像的重建。
[0006] 因此，探寻一种结合接触式图像传感器成像方式进行多谱段重建从而获得更高分 辨率图像的软件方法，必将推动表面检测领域的发展。
【发明内容】
[0007] 本发明提供一种基于接触式图像传感器多谱段图像重建技术方案，能够结合各谱 段的信息，获得纵向分辨率更高的多谱段图像，应用于接触式图像传感器所采集的多谱段 图像。
[0008] 本发明技术方案提供一种基于接触式图像传感器的多谱段图像重建方法，依次包 括前去噪阶段、重建阶段和后去噪阶段，
[0009] 所述前去噪阶段和后去噪阶段，分别包括如下步骤，
[0010]步骤1.1，首先获取η个谱段图像，前去噪阶段是利用接触式图像传感器获取η个谱 段图像，后去噪阶段是利用重建阶段所得所有谱段的处理后的完整图像；
[0011]步骤1.2,对每一个谱段的图像分别进行处理，将单谱段图像切分成mXm像素的小 块，同时相邻小块间有s像素的重叠区域，m和s为预设的取值；
[0012] 步骤1.3,非局部自相似块搜索，包括对每一个小块寻找同谱段中以该小块为中心 的W X W像素区域中的k个最相似的小块，k为预设的取值；
[0013] 步骤1.4,将k个最相似小块的列向量形式^…处串接起来形成图像矩阵^己为 XkXM，对图像矩阵XkXM进行奇异值分解得到特征值λ、特征向量μ以及特征矩阵U;对特征值进 行降序排列，保留前r个较大的特征值使得特征根的能量积值参数大于预设阈值;将特征矩 阵与保留的特征值相乘得到去噪后的图像矩阵VkXM;
[0014] 步骤1.5,将图像矩阵X^xm中的小块恢复至对应位置，并拼成完整的图像，对于重 叠部分的图像取各小块重叠部位像素灰度值的平均值；
[0015] 所述重建阶段包括如下步骤，
[0016] 步骤2.1，设图像纵向是接触式图像传感器与被摄物之间发生相对运动的方向，首 先输入接触式图像传感器多谱段图像，然后对各个谱段图像纵向插值放大，根据传感器发 光顺序对η个谱段的图像进行排序；
[0017] 步骤2.2,将各个谱段图像切分为小块图像，对每一个谱段的图像进行分别处理， 对单谱段图像切分成a X a像素的小块，a为η的整数倍，
[0018] 步骤2.3,将η个谱段图像的同一位置的小块形成一个集合，记为C，则
[0019] C={ci，C2，."，cn}
[0020] 所述Cl为第i个谱段的图像小块；
[0021 ]对每一个小块求其直方图分布，得到直方图分布集合，记为H，则
[0022] H={hi,h2,---,hn}
[0023] 所述lu为第i个谱段图像小块的直方图分布向量；
[0024]步骤2.4,对η个谱段同一位置的图像块依次处理，处理第i个谱段的图像块(^时， 将其它η-I个谱段的图像块进行直方图拉伸，使其直方图范围与第i个谱段的图像块的直方 图范围一致;直方图拉伸使用如下公式，
[0025]
[0026]所述Cj为其余η-I个谱段中的第j个谱段的小块，j为Cj的拉伸结果，maxi和maxj分 另ll为ci和cj中的最大灰度级，miru和minj分别为ci和cj中的最小灰度级;最终得到关于第i个 谱段的处理后图像块集合，记为C i，则
[0027] C7 i = {c71,C7 2, . . . ,Ci, . . .
[0028] 步骤2.5,对处理后图像块集合1进行处理，取第1个图像小块中的第（1+nXd) 行，将取出的各行按顺序拼接在一起，得到处理后的第i个谱段的图像小块，1 = 1，2,…，n，d 为使（1+n X d)不大于a的任意整数；
[0029]步骤2.6,将同谱段的处理后的图像小块按原顺序拼在一起，得到所有谱段的处理 后的完整图像。
[0030]而且，n = 3。
[0031] 而且，步骤2.1中，对各谱段图像每两行之间进行线性插值处理，其纵向分辨率增 大为原图的η倍，并保证第i个谱段的图像的原始数据位于放大图像的第(i+nXd)行，d=l， 2,…，N，N为原图行像素数。
[0032] 本发明相应提供一种基于接触式图像传感器的多谱段图像重建系统，包括前去噪 模块、重建模块和后去噪模块，
[0033] 所述前去噪模块和后去噪模块，分别包括如下单元，
[0034]第一单元，用于首先获取η个谱段图像，前去噪模块的第一单元是利用接触式图像 传感器获取η个谱段图像，后去噪模块的第一单元是利用重建模块所得所有谱段的处理后 的完整图像；
[0035] 第二单元，用于对每一个谱段的图像分别进行处理，将单谱段图像切分成mXm像 素的小块，同时相邻小块间有s像素的重叠区域，m和s为预设的取值；
[0036] 第三单元，用于非局部自相似块搜索，包括对每一个小块寻找同谱段中以该小块 为中心的W X W像素区域中的k个最相似的小块，k为预设的取值；
[0037] 第四单元，用于将k个最相似小块的列向量形式Χ1、Χ2···Χ1?串接起来形成图像矩阵， 记为XkXM，对图像矩阵XkXM进行奇异值分解得到特征值λ、特征向量μ以及特征矩阵U;对特征 值进行降序排列，保留前r个较大的特征值使得特征根的能量积值参数大于预设阈值;将特 征矩阵与保留的特征值相乘得到去噪后的图像矩阵X' kxM;
[0038] 第五单元，用于将图像矩阵X'kxM中的小块恢复至对应位置，并拼成完整的图像，对 于重叠部分的图像取各小块重叠部位像素灰度值的平均值；
[0039]所述重建模块包括如下子单元，
[0040] 第一子单元，用于设图像纵向是接触式图像传感器与被摄物之间发生相对运动的 方向，首先输入接触式图像传感器多谱段图像，然后对各个谱段图像纵向插值放大，根据传 感器发光顺序对η个谱段的图像进行排序；
[0041] 第二子单元，用于将各个谱段图像切分为小块图像，对每一个谱段的图像进行分 别处理，对单谱段图像切分成a X a像素的小块，a为η的整数倍，
[0042]第三子单元，用于将η个谱段图像的同一位置的小块形成一个集合，记为C，则
[0043] C={ci，C2，...，cn}
[0044] 所述Cl为第i个谱段的图像小块；
[0045] 对每一个小块求其直方图分布，得到直方图分布集合，记为H，则
[0046] H={hi，h2，H.，hn}
[0047] 所述hi为第i个谱段图像小块的直方图分布向量；
[0048] 第四子单元，用于对η个谱段同一位置的图像块依次处理，处理第i个谱段的图像 块(^时，将其它η-I个谱段的图像块进行直方图拉伸，使其直方图范围与第i个谱段的图像 块的直方图范围一致;直方图拉伸使用如下公式，
[0049]
[0050]所述Cj为其余η-I个谱段中的第j个谱段的小块，j为Cj的拉伸结果，maxi和maxj分 另ll为ci和cj中的最大灰度级，miru和minj分别为ci和cj中的最小灰度级;最终得到关于第i个 谱段的处理后图像块集合，记为C i，则
[0051 ] C' i = {V 1，C' 2, · · ·，Ci，· · ·，C' η}
[0052]第五子单元，用于对处理后图像块集合(/1进行处理，取第1个图像小块中的第（1+ nXd)行，将取出的各行按顺序拼接在一起，得到处理后的第i个谱段的图像小块，1 = 1， 2,…，n，d为使（Ι+nXd)不大于a的任意整数；
[0053]第六子单元，用于将同谱段的处理后的图像小块按原顺序拼在一起，得到所有谱 段的处理后的完整图像。
[0054] 而且，n = 3。
[0055] 而且，第一子单元中，对各谱段图像每两行之间进行线性插值处理，其纵向分辨率 增大为原图的η倍，并保证第i个谱段的图像的原始数据位于放大图像的第(i+nXd)行，d = 1，2,…，N，N为原图行像素数。
[0056] 本发明为一种关于接触式图像传感器多谱段图像的重建技术方案，对接触式图像 传感器多个谱段的图像进行重建，最终获得所处理谱段的更高分辨率图像，计算量小，便于 硬件快速实现，图像重建效果好，实时性好。
[0057]本发明的有益效果是：
[0058] 1、可以充分利用接触式图像传感器多谱段图像的各谱段信息，并提高所有谱段图 像的纵向分辨率；
[0059] 2、重建过程仅需进行线性映射，效率高；
[0060] 3、适用性好。
【附图说明】
[0061] 图1为本发明实施例的基于接触式图像传感器多谱段图像重建方法流程图。
[0062] 图2为本发明实施例的重建过程中当n = 3，j = 2的情况下第1谱段和第3谱段的示 意图。
【具体实施方式】
[0063] 为了能更清晰地理解本发明，以下结合附图和具体实施例对本发明进行详细说 明。
[0064]本发明实施例的整个重建过程包括前去噪阶段、重建阶段、后去噪阶段三个步骤： 其中，前去噪阶段和后去噪阶段，两去噪阶段过程相同，去噪过程包括:图像分块、非局部自 相似块搜索、奇异值分解去噪、图像块拼合;重建过程包括:图像分块、对各谱段图像块进行 直方图拉伸、图像块拼合。
[0065]如图1，本发明实施例首先获取接触式图像传感器采集的包含红绿蓝3个谱段的图 像，利用SVD方法对其进行去噪处理，得到无噪的3谱段图像;然后对各段图像切块，对于同 一位置的3谱段图像小块，利用直方图范围拉伸分别重建各小块，将所有小块按原顺序拼合 在一起后得到重建后的3谱段图像;最后利用SVD方法对3个谱段的图像分别进行去噪处理， 消除重建中引入的噪声，得到无噪的重建后3谱段图像。
[0066]实施例具体的实施方式如下：
[0067] 一、前去噪阶段：
[0068]前去噪阶段包括如下步骤：
[0069] 步骤1.1:首先利用接触式图像传感器获取η个谱段图像:实施例中n = 3,即利用接 触式图像传感器获取3个谱段图像；
[0070] 具体实施时，η的取值根据情况而定，例如还可能有红外、紫外等谱段；
[0071] 步骤1.2:对每一个谱段的图像分别进行处理，将单谱段图像切分成mXm像素的小 块，同时相邻小块间有s像素的重叠区域;具体实施时，本领域技术人员可自行设定m和s的 取值；
[0072] 步骤1.3:非局部自相似块搜索:对每一个小块寻找其同谱段中以该小块为中心的 W X W像素区域中的k个最相似的小块，相似度通过欧几里得距离进行度量，其计算方法如 下：
[0073]
[0074] ^^，1个^力」而安/艾里仰1认/艾即^丨、1^像小块的列向量形式七为1小块的第1个 像素点，yi为y小块的第i个像素点，Μ为一个小块的总像素数，M=m X m;具体实施时，本领域 技术人员可自行设定k的取值；
[0075] 步骤1.4:将k个最相似小块的列向量形式^、^…处串接起来形成图像矩阵^己为 XkXM，则
[0076] XkXM= [XlX2."Xk]
[0077] 对图像矩阵XkxMX进行SVD分解(奇异值分解)得到特征值λ、特征向量μ以及特征矩 阵U;对特征值进行降序排列，保留前r个较大的特征值使得特征根的能量积值参数大于预 设阈值，具体实施时，本领域技术人员可自行设定取值，建议取值范围在80 %到95 %之间， 例如90%;将特征矩阵与保留的r个特征值相乘得到去噪后的图像矩阵X^xm;
[0078] 步骤1.5:将图像矩阵X^xm中的小块恢复至对应位置，并拼成完整的图像，对于重 叠部分的图像取各小块重叠部位像素灰度值的平均值。
[0079]二、重建阶段
[0080]重建阶段包括如下步骤：
[0081]步骤2.1:首先设图像纵向是接触式图像传感器与被摄物之间发生相对运动的方 向，输入接触式图像传感器多谱段图像，参见图1中S101，如果将相对运动方向设为横向，处 理方式等同，后续进行横向插值和列合并即可；
[0082]然后对各个谱段图像纵向插值放大，参见图1中S102:
[0083] 根据传感器发光顺序对η个谱段的图像进行排序，对η个谱段的图像进行分别处 理，对单谱段图像每两行之间进行线性插值处理，使其纵向分辨率增大为原图的η倍，对第i 个谱段的图像，原图的第1行位于处理后图像的第j行，若j>i，则删除前(j_i)行，并在最后1 行之后复制（j_i)次最后1行，若j〈i，则在第1行之前复制（i _ j)次第一行并删除最后（i_j) 行，使得第i个谱段的图像的原始数据位于处理后图像的第（i+nXd)行，d = l，2,…，D，D为 原图行像素数，从而保证重建图像时每一行图像都能对应一个谱段的原始数据;见图2当η =3, j = 2的情况下第1谱段和第3谱段的示意图。
[0084] 实施例根据传感器发光顺序对3个谱段的图像进行排序，本实施例中传感器发光 顺序为红、蓝、绿，对3个谱段的图像进行分别处理，对单谱段图像每两行之间进行线性插值 处理，使其纵向分辨率增大为原图的3倍，并保证第i个谱段的图像的原始数据位于放大图 像的第（i+nXd)行，d = l，2,…，N，N为原图行像素数；
[0085]参见图2,设n = 3，j = 2的情况下第1谱段和第3谱段的示意图，N行经插值处理后变 为3N行，以第1谱段为例，第一幅为原图，第二幅是插值后的图像，第三幅是做复制和删除操 作以后的图像，原图的第1行位于处理后图像的第2行，因为j>i，则删除第1行，并在最后1行 之后复制（j-i)次最后1行，使得第1个谱段的图像的原始数据位于处理后图像的第（l+3d) 行，d=l，2,…，N，N为原图行像素数，如第1谱段处理后处于，从而保证重建图像时每一行图 像都能对应一个谱段的原始数据。
[0086]步骤2.2:将各个谱段图像切分为小块图像，参见图1中S103:对每一个谱段的图像 进行分别处理，对单谱段图像切分成a Xa像素的小块，a为η的整数倍：X [0087]实施例中a = 6,即对每一个谱段的图像进行分别处理，对单谱段图像切分成6X6 像素的小块。
[0088] 步骤2.3:对每个小块图像计算其直方图分布，参见图1中S104:将η个谱段图像的 同一位置的小块形成一个集合，记为C，则
[0089] C={ci,C2,---,Cn}
[0090] 所述Cl为第i个谱段的图像小块;对每一个小块求其直方图分布，得到直方图分布 集合，记为Η，则
[0091] H={hi，h2，."，hn}
[0092] 所述lu为第i个谱段图像小块的直方图分布向量。
[0093] 实施例中，将3个谱段图像的同一位置的小块形成一个集合，记为C，为便于标记起 见，下标直接采用R、G、B，则
[0094] C={cr，cb，cg}
[0095] 其中CR为红光谱段的图像小块，CB为蓝光谱段的图像小块，CC为绿光谱段的图像小 块;对每一个小块求其直方图分布，得到直方图分布集合，记为Η，则
[0096] H={hR，hB，hc}
[0097] 其中hR为红光谱段图像小块的直方图分布向量，hB为蓝光谱段图像小块的直方图 分布向量，he为绿光谱段图像小块的直方图分布向量。
[0098]步骤2.4:对η个谱段的图像块同一位置的依次处理，处理第i个谱段的图像块(^ 时，将其它η-I个谱段的图像块进行直方图拉伸，使其直方图范围与第i个谱段的图像块的 直方图范围一致;直方图拉伸使用如下公式：
[0099]
[0?00]所述Cj为其余n-ι个谱段中的第j个谱段的小块，j为Cj的拉伸结果，maxi和maxj分 另ll为ci和cj中的最大灰度级，miru和minj分别为ci和cj中的最小灰度级;最终得到关于第i个 谱段的处理后图像块集合，记为C i，则
[0101 ] C7 i = {c7 1,07 2, ··· ,Ci, ··· ,C7 η}
[0102]实施例对3个谱段的图像块依次处理，以红光谱段图像块CR的重建为例，将蓝光谱 段的图像块cb和绿光谱段的图像块cc分别进行直方图拉伸，使其直方图范围与红光谱段图 像块CR的直方图范围一致;直方图拉伸使用如下公式：
[0103]
[0104]其中Cj为其余n-1个谱段中的第j个谱段的小块，maxi和maxj分别为Ci和Cj中的最大 灰度级，miru和minj分别为ci和cj中的最小灰度级；
[0105] 其中，j = l，2,3,分别代表红光谱段、蓝光谱段、绿光谱段；以j = l为例，最终得到 关于红光谱段的处理后图像块集合，记为C'r，则
[0106] C/r={cr,c/g,c/b}
[0107] 其中，所述CR为红光谱段的小块，b为绿光谱段的小块、蓝光谱段的小块的拉 伸结果。
[0108] 步骤2.5:图像块拼合:对处理后图像块集合(:^进行处理，取第1个图像小块中的 第（1+nXd)行（1 = 1，2,…，n，d为使（1+nXd)不大于a的任意整数），将取出的各行按顺序拼 接在一起，即得处理后的第i个谱段的图像小块：
[0109] 实施例中，以红光谱段图像块CR的重建为例，对处理后图像块集合C'R进行处理，取 cr中的第1、4行，6中的第2、5行，A中的第3、6行，将6行按原顺序拼接在一起，即得到处理 后的红光谱段的图像小块。
[0110] 步骤2.6:将同谱段的处理后的图像小块按原顺序拼在一起，得到所有谱段的处理 后的完整图像。
[0111] 具体实施时，步骤2.4~2.6可针对各频段依次进行，参见图1中S105~S109。包括 以下步骤：
[0112] 选择一个谱段进行重建，处理该谱段的某图像块(^时，将其它η-I个谱段的图像块 进行直方图拉伸，即根据像素灰度做线性映射，使其直方图范围与第i个谱段的待处理图像 块的直方图范围一致；
[0113] 对同一位置的所有谱段小块图像，取其反应真实数据的一行，按谱段的发光顺序 拼接在一起，得到重建后的待处理谱段图像小块；
[0114] 将重建后的图像小块按原顺序拼接在一起，得到重建后的特定谱段图像；
[0115] 判断是否所有谱段全部重建完毕，是则输出重建后的接触式传感器多谱段高分辨 率图像，否则返回选下一个谱段进行重建。
[0116] 三、后去噪阶段
[0117] 再次执行步骤1.1~步骤1.5,即可完成后去噪阶段。唯一区别在于，步骤1.1是根 据步骤2所得所有谱段的处理后的完整图像。
[0118] 具体实施时，本发明所提供方法可基于软件技术实现自动运行流程，也可采用模 块化方式实现相应系统。
[0119] 本发明相应提供一种基于接触式图像传感器的多谱段图像重建系统，包括前去噪 模块、重建模块和后去噪模块，
[0120] 所述前去噪模块和后去噪模块，分别包括如下单元，
[0121] 第一单元，用于首先获取η个谱段图像，前去噪模块的第一单元是利用接触式图像 传感器获取η个谱段图像，后去噪模块的第一单元是利用重建模块所得所有谱段的处理后 的完整图像；
[0122] 第二单元，用于对每一个谱段的图像分别进行处理，将单谱段图像切分成mXm像 素的小块，同时相邻小块间有S像素的重叠区域，m和S为预设的取值；
[0123] 第三单元，用于非局部自相似块搜索，包括对每一个小块寻找同谱段中以该小块 为中心的W X W像素区域中的k个最相似的小块，k为预设的取值；
[0124] 第四单元，用于将k个最相似小块的列向量形式Χ1、Χ2···Χ1?串接起来形成图像矩阵， 记为XkXM，对图像矩阵XkXM进行奇异值分解得到特征值λ、特征向量μ以及特征矩阵U;对特征 值进行降序排列，保留前r个较大的特征值使得特征根的能量积值参数大于预设阈值;将特 征矩阵与保留的特征值相乘得到去噪后的图像矩阵X' kxM;
[0125] 第五单元，用于将图像矩阵X'kxM中的小块恢复至对应位置，并拼成完整的图像，对 于重叠部分的图像取各小块重叠部位像素灰度值的平均值；
[0126] 所述重建模块包括如下子单元，
[0127] 第一子单元，用于设图像纵向是接触式图像传感器与被摄物之间发生相对运动的 方向，首先输入接触式图像传感器多谱段图像，然后对各个谱段图像纵向插值放大，根据传 感器发光顺序对η个谱段的图像进行排序；
[0128] 第二子单元，用于将各个谱段图像切分为小块图像，对每一个谱段的图像进行分 别处理，对单谱段图像切分成a X a像素的小块，a为η的整数倍，
[0129] 第三子单元，用于将η个谱段图像的同一位置的小块形成一个集合，记为C，则
[0130] C={ci，C2，."，cn}
[0131] 所述Cl为第i个谱段的图像小块；
[0132] 对每一个小块求其直方图分布，得到直方图分布集合，记为H，则
[0133] H={hi,h2,---,hn}
[0134] 所述lu为第i个谱段图像小块的直方图分布向量；
[0135] 第四子单元，用于对η个谱段同一位置的图像块依次处理，处理第i个谱段的图像 块(^时，将其它η-I个谱段的图像块进行直方图拉伸，使其直方图范围与第i个谱段的图像 块的直方图范围一致;直方图拉伸使用如下公式，
[0136]
[0137]所述Cj为其余η-I个谱段中的第j个谱段的小块，j为Cj的拉伸结果，maxi和maxj分 另ll为ci和cj中的最大灰度级，miru和minj分别为ci和cj中的最小灰度级;最终得到关于第i个 谱段的处理后图像块集合，记为C i，则
[0138] C7 i = {c71,〇7 2, ··· ,Ci, ··· ,c7 η}
[0139] 第五子单元，用于对处理后图像块集合1进行处理，取第1个图像小块中的第（1+ nXd)行，将取出的各行按顺序拼接在一起，得到处理后的第i个谱段的图像小块，1 = 1， 2,…，n，d为使（Ι+nXd)不大于a的任意整数；
[0140] 第六子单元，用于将同谱段的处理后的图像小块按原顺序拼在一起，得到所有谱 段的处理后的完整图像。
[0141] 各模块具体实现可参见相应步骤，本发明不予赘述。
[0142] 上述实施例并不代表本发明的全部应用。凡在本发明的精神和原则之内，所做的 各种变化、修改，均在本发明的保护范围之内。
【主权项】
1. 一种基于接触式图像传感器的多谱段图像重建方法，其特征在于：依次包括前去噪 阶段、重建阶段和后去噪阶段， 所述前去噪阶段和后去噪阶段，分别包括如下步骤， 步骤1.1，首先获取η个谱段图像，前去噪阶段是利用接触式图像传感器获取η个谱段图 像，后去噪阶段是利用重建阶段所得所有谱段的处理后的完整图像； 步骤1.2,对每一个谱段的图像分别进行处理，将单谱段图像切分成mXm像素的小块， 同时相邻小块间有S像素的重叠区域，m和S为预设的取值； 步骤1.3,非局部自相似块捜索，包括对每一个小块寻找同谱段中W该小块为中屯、的W XW像素区域中的k个最相似的小块，k为预设的取值； 步骤1.4,将k个最相似小块的列向量形式xi、X2'''xk串接起来形成图像矩阵，记为XkxM， 对图像矩阵XkXM进行奇异值分解得到特征值λ、特征向量μΚ及特征矩阵U;对特征值进行降 序排列，保留前r个较大的特征值使得特征根的能量积值参数大于预设阔值;将特征矩阵与 保留的特征值相乘得到去噪后的图像矩阵X^kXM; 步骤1.5,将图像矩阵X^kxM中的小块恢复至对应位置，并拼成完整的图像，对于重叠部 分的图像取各小块重叠部位像素灰度值的平均值； 所述重建阶段包括如下步骤， 步骤2.1，设图像纵向是接触式图像传感器与被摄物之间发生相对运动的方向，首先输 入接触式图像传感器多谱段图像，然后对各个谱段图像纵向插值放大，根据传感器发光顺 序对η个谱段的图像进行排序； 步骤2.2,将各个谱段图像切分为小块图像，对每一个谱段的图像进行分别处理，对单 谱段图像切分成a X a像素的小块，a为η的整数倍， 步骤2.3，将η个谱段图像的同一位置的小块形成一个集合，记为C，则 C= kl,C2, ... ,Cn} 所述Cl为第i个谱段的图像小块； 对每一个小块求其直方图分布，得到直方图分布集合，记为Η，则 Η=化l，h2，...，hn} 所述hi为第i个谱段图像小块的直方图分布向量； 步骤2.4,对η个谱段同一位置的图像块依次处理，处理第i个谱段的图像块Cl时，将其它 n-1个谱段的图像块进行直方图拉伸，使其直方图范围与第i个谱段的图像块的直方图范围 一致;直方图拉伸使用如下公式，所述cj为其余n-1个谱段中的第j个谱段的小块，c'j为cj的拉伸结果，maxi和maxj分别为 Ci和Cj中的最大灰度级，mini和minj分别为Ci和Cj中的最小灰度级;最终得到关于第i个谱段 的处理后图像块集合，记为C/ 1，则 C' i= k' 1, C' 2 , ,Ci , , C' η} 步骤2.5,对处理后图像块集合(/1进行处理，取第1个图像小块中的第（1+nXd)行，将取 出的各行按顺序拼接在一起，得到处理后的第i个谱段的图像小块，1 = 1，2,…，n，d为使(1+ nXd)不大于a的任意整数； 步骤2.6,将同谱段的处理后的图像小块按原顺序拼在一起，得到所有谱段的处理后的 完整图像。2. 根据权利要求1所述基于接触式图像传感器的多谱段图像重建方法，其特征在于:n =3。3. 根据权利要求1或2所述基于接触式图像传感器的多谱段图像重建方法，其特征在 于:步骤2.1中，对各谱段图像每两行之间进行线性插值处理，其纵向分辨率增大为原图的η 倍，并保证第i个谱段的图像的原始数据位于放大图像的第（i+nXd)行，d = l，2,…，Ν，Ν为 原图行像素数。4. 一种基于接触式图像传感器的多谱段图像重建系统，其特征在于:包括前去噪模块、 重建模块和后去噪模块， 所述前去噪模块和后去噪模块，分别包括如下单元， 第一单元，用于首先获取η个谱段图像，前去噪模块的第一单元是利用接触式图像传感 器获取η个谱段图像，后去噪模块的第一单元是利用重建模块所得所有谱段的处理后的完 整图像； 第二单元，用于对每一个谱段的图像分别进行处理，将单谱段图像切分成mXm像素的 小块，同时相邻小块间有S像素的重叠区域，m和S为预设的取值； 第Ξ单元，用于非局部自相似块捜索，包括对每一个小块寻找同谱段中W该小块为中 屯、的WXW像素区域中的k个最相似的小块，k为预设的取值； 第四单元，用于将k个最相似小块的列向量形式串接起来形成图像矩阵，记为 XkXM，对图像矩阵XkXM进行奇异值分解得到特征值λ、特征向量μΚ及特征矩阵U;对特征值进 行降序排列，保留前r个较大的特征值使得特征根的能量积值参数大于预设阔值;将特征矩 阵与保留的特征值相乘得到去噪后的图像矩阵X^kXM; 第五单元，用于将图像矩阵χ/kxM中的小块恢复至对应位置，并拼成完整的图像，对于重 叠部分的图像取各小块重叠部位像素灰度值的平均值； 所述重建模块包括如下子单元， 第一子单元，用于设图像纵向是接触式图像传感器与被摄物之间发生相对运动的方 向，首先输入接触式图像传感器多谱段图像，然后对各个谱段图像纵向插值放大，根据传感 器发光顺序对η个谱段的图像进行排序； 第二子单元，用于将各个谱段图像切分为小块图像，对每一个谱段的图像进行分别处 理，对单谱段图像切分成a X a像素的小块，a为η的整数倍， 第Ξ子单元，用于将η个谱段图像的同一位置的小块形成一个集合，记为C，则 C= {>1,02, ... ,Cn} 所述Cl为第i个谱段的图像小块； 对每一个小块求其直方图分布，得到直方图分布集合，记为Η，则 Η=化I,h2,... ,hn} 所述hi为第i个谱段图像小块的直方图分布向量； 第四子单元，用于对η个谱段同一位置的图像块依次处理，处理第i个谱段的图像块Cl 时，将其它n-1个谱段的图像块进行直方图拉伸，使其直方图范围与第i个谱段的图像块的 直方图范围一致;直方图拉伸使用如下公式，所述Cj为其余n-1个谱段中的第j个谱段的小块，c'j为Cj的拉伸结果，maxi和maxj分别为 Ci和Cj中的最大灰度级，mini和minj分别为Ci和Cj中的最小灰度级;最终得到关于第i个谱段 的处理后图像块集合，记为C/ 1，则 C' i= k' 1, C' 2 , ,Ci , , C' η} 第五子单元，用于对处理后图像块集合C^i进行处理，取第1个图像小块中的第（1+nXd) 行，将取出的各行按顺序拼接在一起，得到处理后的第i个谱段的图像小块，1 = 1，2,…，n，d 为使(1+n X d)不大于a的任意整数； 第六子单元，用于将同谱段的处理后的图像小块按原顺序拼在一起，得到所有谱段的 处理后的完整图像。5. 根据权利要求4所述基于接触式图像传感器的多谱段图像重建系统，其特征在于:n 二 3 〇6. 根据权利要求4或5所述基于接触式图像传感器的多谱段图像重建系统，其特征在 于:第一子单元中，对各谱段图像每两行之间进行线性插值处理，其纵向分辨率增大为原图 的η倍，并保证第i个谱段的图像的原始数据位于放大图像的第(i+nXd)行，d=l，2,…，N，N 为原图行像素数。
【文档编号】G06T3/40GK105976316SQ201610270348
【公开日】2016年9月28日
【申请日】2016年4月27日
【发明人】周严, 鹿璇
【申请人】武汉卓目科技有限公司