技术领域
[0001] 本
发明涉及
图像处理技术领域，具体而言，涉及一种高光谱图像预处理方法、装置及电子设备。
背景技术
[0002] 目前常规的对无人机航测图像的处理技术大多存在不
稳定性，使得处理之后的图像存在
辐射畸变，光谱定标之后的光谱曲线严重失真。
发明内容
[0003] 有鉴于此，本发明
实施例的目的在于提供一种高光谱图像预处理方法、装置及电子设备，以解决
现有技术对图像处理之后存在辐射畸变和导致光谱曲线失真的问题。
[0004] 本发明实施例提供了一种高光谱图像预处理方法，所述方法包括：
[0005] 获得多个待处理高光谱图像，各所述待处理高光谱图像包括反映图像尺寸的行数；
[0006] 查找出所述多个待处理高光谱图像中行数最小的待处理高光谱图像，将该待处理高光谱图像作为基准图像，根据该基准图像的最小行数对所述多个待处理高光谱图像中的其余待处理高光谱图像进行修正，将修正之后的待处理高光谱图像和所述基准图像作为修正版高光谱图像；
[0007] 针对每个修正版高光谱图像，获得该修正版高光谱图像的
亮度和调整系数，根据所述调整系数对该修正版高光谱图像的亮度进行调整，得到调整版高光谱图像；
[0008] 获得定标公式，针对每个调整版高光谱图像，采用所述定标公式对该调整版高光谱图像进行光谱定标处理，得到光谱定标图像；
[0009] 针对每个光谱定标图像，获取该光谱定标图像的光谱反射率，根据所述光谱反射率对该光谱定标图像进行平滑处理，得到平滑高光谱图像。
[0010] 可选地，根据该基准图像的最小行数对所述多个待处理高光谱图像进行修正的步骤，包括：
[0011] 针对所述多个待处理高光谱图像中的其余待处理高光谱图像中的每个待处理高光谱图像，根据最小行数对该待处理高光谱图像的行数进行剔除处理，使经过剔除处理的该待处理高光谱图像的行数与所述最小行数一致。
[0012] 可选地，各所述待处理高光谱图像还包括多个波段数，所述各波段数包括多个像元值，所述各像元值包括对应的亮度，针对每个修正版高光谱图像，获得该修正版高光谱图像的亮度和调整系数，根据调整系数对该修正版高光谱图像的亮度进行调整，得到调整版高光谱图像的步骤，包括：
[0013] 获得该修正版高光谱图像的所有波段数中每个波段数的亮度，计算所有波段数的亮度的平均值；
[0014] 根据所述平均值计算得到该修正版影像的调整系数；
[0015] 根据所述调整系数对该修正版高光谱图像的每个波段中的每个像元值所对应的亮度进行调整，得到调整版高光谱图像。
[0016] 可选地，根据所述调整系数对该修正版高光谱图像的每个波段中的每个像元值所对应的亮度进行调整，得到调整版高光谱图像的步骤，包括：
[0017] 采用调整公式对该修正版高光谱图像的每个波段中的每个像元值所对应的亮度进行调整，得到调整版高光谱图像，所述调整公式为：
[0018] ak(i,j)＝bk(i,j)×t(i,j)
[0019] 其中，
[0020] i为该修正版高光谱图像的列数，i的取值范围为1到列数设定值；
[0021] j为该修正版高光谱图像的行数，j的取值范围为1到行数设定值；
[0022] k为波段数，k的取值范围为1到波段数设定值；
[0023] ak(i,j)为调整过亮度的像元值；
[0024] bk(i,j)为未调整过亮度的像元值；
[0025] t(i,j)为调整系数，t(i,j)通过所述平均值计算得到。
[0026] 可选地，所述定标公式通过以下步骤获得：
[0027] 获得白板定标值和黑板定标值；
[0028] 根据所述白板定标值和所述黑板定标值计算得到定标公式；
[0029] 所述定标公式为：
[0030] R k(i,j)＝(ak(i,j)-dbk)/(dwk-dbk)
[0031] 其中，
[0032] dwk为第k波段对应的白板定标值；
[0033] dbk为第k波段对应的黑板定标值。
[0034] 可选地，针对每个光谱定标图像，获取该光谱定标图像的光谱反射率，根据所述光谱反射率对该光谱定标图像进行平滑处理，得到平滑高光谱图像的步骤，包括：
[0035] 获得该光谱定标图像的每个像元值在每个波段数中的光谱反射率；根据所述多个光谱反射率和所述多个波段数获得该像元值的光谱曲线；
[0036] 根据所述光谱曲线对该像元值进行滤波处理；
[0037] 计算滤波处理后的像元值和滤波处理前的像元值的差值，若该差值不超过预设值，将滤波处理前的像元值作为平滑像元值，将该平滑像元值对应的光谱定标图像作为平滑高光谱图像；若该差值超过预设值，将滤波处理后的像元值作为平滑像元值，将该平滑像元值对应的光谱定标图像作为平滑高光谱图像。
[0038] 可选地，所述方法还包括：
[0039] 针对所述光谱曲线，查找出该光谱曲线中
信噪比低于预设标准值的波段范围，对该波段范围内的每个像元值进行滤波处理；
[0040] 计算滤波处理后的像元值和滤波处理前的像元值的差值，若该差值不超过预设值，将滤波处理前的像元值作为平滑像元值，将该平滑像元值对应的波段范围作为平滑处理后的波段范围，将该平滑处理后的波段范围对应的光谱定标图像作为平滑高光谱图像；若该差值超过预设值，将滤波处理后的像元值作为平滑像元值，将该平滑像元值对应的波段范围作为平滑处理后的波段范围，将该平滑处理后的波段范围对应的光谱定标图像作为平滑高光谱图像。
[0041] 本发明实施例还提供了一种高光谱图像预处理装置，所述高光谱图像预处理装置包括：
[0042] 获取模
块，用于获得多个待处理高光谱图像，各所述待处理高光谱图像包括行数；
[0043] 行数修正模块，用于查找出所述多个待处理高光谱图像中行数最小的待处理高光谱图像，将该待处理高光谱图像作为基准图像，根据该基准图像的最小行数对所述多个待处理高光谱图像中的其余待处理高光谱图像进行修正，将修正之后的待处理高光谱图像和所述基准图像作为修正版高光谱图像；
[0044] 亮度调整模块，用于针对每个修正版高光谱图像，获得该修正版高光谱图像的亮度和调整系数，根据所述调整系数对该修正版高光谱图像的亮度进行调整，得到调整版高光谱图像；
[0045] 光谱定标模块，用于获得定标公式，针对每个调整版高光谱图像，采用所述定标公式对该调整版高光谱图像进行光谱定标处理，得到光谱定标图像；
[0046] 平滑处理模块，用于针对每个光谱定标图像，获取该光谱定标图像的光谱反射率，根据所述光谱反射率对该光谱定标图像进行平滑处理，得到平滑高光谱图像。
[0047] 本发明实施例还提供了一种电子设备，包括
存储器、处理器及存储在存储器上并可在处理器上运行的
计算机程序，所述处理器执行所述程序时实现上述高光谱图像预处理方法。
[0048] 本发明实施例还提供了一种计算机可读存储介质，所述可读存储介质包括计算机程序，所述计算机程序运行时控制所述可读存储介质所在电子设备执行上述高光谱图像预处理方法。
[0049] 本发明实施例提供的高光谱图像预处理方法、装置及电子设备能根据调整系数对图像的亮度进行调整，保证了每张图的总体亮度一致。对光谱定标后的图像进行平滑处理，能对光谱曲线失真的部分进行修正，进而减少失真和辐射畸变。
[0050] 为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附
附图，作详细说明如下。
附图说明
[0051] 为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
[0052] 图1为本发明较佳实施例提供的一种电子设备10的方框示意图。
[0053] 图2为本发明较佳实施例提供的一种高光谱图像预处理方法的
流程图。
[0054] 图3为一实施方式中图2所示步骤S22包括的另一子步骤的示意图。
[0055] 图4为一实施方式中图4所示步骤S24包括的另一子步骤的示意图。
[0056] 图5为本发明实施例所提供的一种高光谱图像预处理装置20的模块
框图。
[0057] 图标：10-电子设备；11-存储器；12-处理器；13-网络模块；20-高光谱图像预处理装置；21-获取模块；22-行数修正模块；23-亮度调整模块；24-光谱定标模块；25-平滑处理模块。
具体实施方式
[0058] 经研究发现，目前常规的对无人机航测获得的图像进行处理的
软件，如ENVI、PCI、ERDAS等，在图像的辐射校正方面的功能模块，只针对参见的国内外陆地卫星影像，没有针对国产新型旋翼无人机载高光谱
传感器的选项，所以不能解决此类传感器获取的光谱图像上出现的辐射畸变问题，在进行图像处理后，光谱曲线的形态严重偏离实际光谱。
[0059] 基于上述研究，本发明实施例提供了一种能够对无人机航测获得的高光谱图像进行畸变处理和平滑处理的方法。
[0060] 下面将结合本发明实施例中附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的
选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0061] 应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
[0062] 如图1所示，是本发明较佳实施例提供的电子设备10的方框示意图。本发明实施例中的电子设备10可以为具有数据存储、传输、处理功能的服务端。如图1所示，电子设备10包括：存储器11、处理器12、网络模块13及高光谱图像预处理装置20。
[0063] 所述存储器11、处理器12以及网络模块13相互之间直接或间接地电性连接，以实现数据的传输或交互。例如，这些元件相互之间可通过一条或多条通讯总线或
信号线实现电性连接。存储器11中存储有高光谱图像预处理装置20，所述高光谱图像预处理装置20包括至少一个可以软件或
固件(firmware)的形式存储于所述存储器11中的软件功能模块，所述处理器12通过运行存储在存储器11内的软件程序以及模块，如本发明实施例中的高光谱图像预处理装置20，从而执行各种功能应用以及
数据处理，即实现本发明实施例中的
视频点播方法。
[0064] 其中，所述存储器11可以是，但不限于，
随机存取存储器(Random Access Memory，RAM)，
只读存储器(Read Only Memory，ROM)，可编程只读存储器(Programmable Read-Only Memory，PROM)，可擦除只读存储器(Erasable Programmable Read-Only Memory，EPROM)，电可擦除只读存储器(Electric Erasable Programmable Read-Only Memory，EEPROM)等。其中，存储器11用于存储程序，所述处理器12在接收到执行指令后，执行所述程序。
[0065] 所述处理器12可能是一种集成
电路芯片，具有数据的处理能
力。上述的处理器12可以是通用处理器，包括
中央处理器(Central Processing Unit，CPU)、网络处理器(Network Processor，NP)等。可以实现或者执行本发明实施例中公开的各方法、步骤及逻辑框图。通用处理器可以是
微处理器或者该处理器也可以是任何常规的处理器等。
[0066] 网络模块13用于通过网络建立电子设备10与其他通信终端之间的通信连接，实现网络信号及数据的收发操作。上述网络信号可包括无线信号或者有线信号。
[0067] 可以理解，图1所示的结构仅为示意，电子设备10还可包括比图1中所示更多或者更少的组件，或者具有与图1所示不同的配置。图1中所示的各组件可以采用
硬件、软件或其组合实现。
[0068] 本发明实施例还提供一种计算机可读存储介质，所述可读存储介质包括计算机程序。所述计算机程序运行时控制所述可读存储介质所在电子设备10执行下面的高光谱图像预处理方法。
[0069] 请参阅图2，是本发明较佳实施例提供的一种高光谱图像预处理方法的流程图。所述方法有关的流程所定义的方法步骤应用于电子设备10，可以由所述处理器12实现。下面将对图2所示的具体流程进行详细阐述。
[0070] 步骤S21，获得多个待处理高光谱图像，对各待处理高光谱图像的统一性进行处理，得到多个修正版高光谱图像。
[0071] 在本实施例中，由于人工控制的旋翼无人机光谱成像，在
图像采集过程中，图像行数不一致，因此需要对每次航测获取的多个待处理高光谱图像的行数进行调整，便于后续的统计分析，为后续的图像亮度的统计和调整做准备。
[0072] 可选地，针对n个待处理高光谱图像，查找出n各待处理高光谱图像中行数最小的待处理高光谱图像，将该最小行数作为标准，例如，最小行数为1023，将1023作为标准。将其他待处理高光谱图像的上下边缘各剔除等量的行数，使各个待处理高光谱图像的行数都为1023。由于图像分幅之间存在较大的重叠度，此过程实际只是去除了极少数的重叠区影像条带，不会造成图像的缺失。
[0073] 步骤S22，针对每个修正版高光谱图像，获得该修正版高光谱图像的亮度和调整系数，根据调整系数对该修正版高光谱图像的亮度进行调整，得到调整版高光谱图像。
[0074] 请结合参阅图3，本实施例中通过步骤S221、步骤S222和步骤S223列举了步骤S22的其中一种实现方式。
[0075] 步骤S221，获得该修正版高光谱图像的所有波段数中的每个波段数的亮度，计算所有波段数的亮度的平均值。
[0076] 在本实施例中，各待处理高光谱图像包括多个波段数，各波段数包括多个像元值，各像元值包括对应的亮度。经过步骤S21的处理后，各修正版高光谱图像的波段数、像元值和亮度没有改变。
[0077] 例如，每个修正版高光谱图像可以包括176个波段，每个波段可以包括960*1023个像元值，每个像元值包括对应的亮度。
[0078] 计算176个波段的亮度的平均值M，计算176个波段的总体亮度图像p(i,j)。
[0079] 步骤S222，根据所述平均值计算得到该修正版影像的调整系数。
[0080] 根据平均值M，计算调整系数t(i,j)，计算公式如下：
[0081] t(i,j)＝M/p(i,j)
[0082] 其中，
[0083] t(i,j)为调整系数；
[0084] i为列数，在本实施例中，i的取值范围可以是[0,959]；
[0085] j为行数，在本实施例中，j的取值范围可以是[0,1022]。
[0086] 步骤S223，根据所述调整系数对该修正版高光谱图像的每个波段中的每个像元值所对应的亮度进行调整，得到调整版高光谱图像。
[0087] 根据调整系数，采用调整公式对该修正版高光谱图像的每个波段中的每个像元值所对应的亮度进行调整，得到调整版高光谱图像，调整公式为：
[0088] ak(i,j)＝bk(i,j)×t(i,j)
[0089] 其中，
[0090] k为波段数，在本实施例中，k的取值范围可以是[0,175]；
[0091] ak(i,j)为调整过的像元值；
[0092] bk(i,j)为未调整过的像元值。
[0093] 通过上述公式，可以获得每个波段下的每个像元值经过亮度调整过后的像元值ak(i,j)。
[0094] 步骤S23，获得定标公式，针对每个调整版高光谱图像，采用所述定标公式对该调整版高光谱图像进行光谱定标处理，得到光谱定标图像。
[0095] 遥感成像时，需要测定白板的反射辐射强度，这种辐射来自于被白板反射的
太阳辐射，还需要测定黑板辐射强度，以便进行光谱定标处理。
[0096] 在本实施例中，可以采用双利合谱的V10成像光谱仪在定标测试过程中产生两个文件，*.wcor和*.dcor，这两个文件分别对应白板和黑板的定标记录，根据定标记录可以得到定标公式如下所示：
[0097] Rk(i,j)＝(ak(i,j)-dbk)/(dWk-dbk)
[0098] 其中，
[0099] dWk为第k波段对应的白板定标值；
[0100] dbk为第k波段对应的黑板定标值；
[0101] Rk(i,j)为光谱定标之后的像元值。
[0102] 通过上述公式，可以对经过亮度调整的像元值进行光谱定标处理，得到Rk(i,j)。
[0103] 在本实施例中，Rk(i,j)的结果是浮点型数据，理论值介于0.0～1.0之间，为便于分析，可以将Rk(i,j)放大10000倍，以整形方式记录。
[0104] 步骤S24，针对每个光谱定标图像，获取该光谱定标图像的光谱反射率，根据所述光谱反射率对该光谱定标图像进行平滑处理，得到平滑高光谱图像。
[0105] 请结合参阅图4，本实施例中通过步骤S241和步骤S242列举了步骤S24的其中一种实现方式。
[0106] 步骤S241，获得该光谱定标图像的每个像元值在每个波段数中的光谱反射率，根据所述多个光谱反射率和所述多个波段数获得该像元值的光谱曲线。
[0107] 在本实施例中，每个像元值在每个波段数中的光谱反射率可以作为纵坐标，波段数可以作为横坐标，根据横、纵坐标可以得到该像元值的光谱曲线。
[0108] 步骤S242，根据所述光谱曲线对该像元值进行滤波处理。
[0109] 在本实施例中，对像元值进行滤波处理的公式如下：
[0110] MRk(i,j)＝mean(Rk(i,j)，N)
[0111] 其中，
[0112] MRk(i,j)为滤波之后的像元值；
[0113] Rk(i,j)为光谱定标之后的像元值，可以理解为滤波之前的像元值；
[0114] N为均值滤波的核大小，在本实施例中N的取值可以为7；
[0115] mean(,)为中值滤波计算函数。
[0116] 通过上述公式，可以对光谱定标后的像元值Rk(i,j)进行滤波处理，得到滤波之后的像元值MRk(i,j)。
[0117] 进一步地，计算各像点在每个波段的像元值滤波前后的差值Δk(i,j)，若差值Δk(i,j)不超过预设值D，则保持滤波前的像元值不变，若差值Δk(i,j)超过预设值D时，将滤波前的像元值替换为滤波后的像元值，如此判别的意义在于保留某些细微的变化。
[0118] 例如，判定过程如下：
[0119] Δk(i,j)＝Rk(i,j)-MRk(i,j)
[0120] if(Δk(i,j)≥D)，then Rk(i,j)＝MRk(i,j)
[0121] 可选地，可针对某一光谱曲线，查找出该光谱曲线中信噪比低于预设标准值的波段范围，对该波段范围内的每个像元值进行上述的滤波处理。
[0122] 例如，分析某一光谱曲线，查找出
波长介于400～500nm、700～1000nm之间的光谱曲线信噪比低于预设标准值，则对波长介于400～500nm、700～1000nm之间的光谱曲线对应的像元值进行滤波处理，又例如，判定过程如下：
[0123] if(Δk(i,j)≥D)and(k≤k1or k≥k2)
[0124] then Rk(i,j)＝MRk(i,j)
[0125] 其中，k1和k2为波长小于500nm和大于700nm处对应的波段序号，预设标准值根据实际需求设置。
[0126] 如此设置，根据波长范围和变化幅度有限制的调整，既能有效地平滑信噪比低的区间的曲线，又能保留信噪比高的区间的光谱细节特征。
[0127] 在上述
基础上，如图5所示，本发明实施例提供了一种高光谱图像预处理装置20，所述高光谱图像预处理装置20包括获取模块21、行数修正模块22、亮度调整模块23、光谱定标模块24和平滑处理模块25。
[0128] 获取模块21，用于获得多个待处理高光谱图像，各所述待处理高光谱图像包括行数。
[0129] 由于获取21和图2中步骤S21的实现原理类似，因此在此不作更多说明。
[0130] 行数修正模块22，用于查找出所述多个待处理高光谱图像中行数最小的待处理高光谱图像，将该待处理高光谱图像作为基准图像，根据该基准图像的最小行数对所述多个待处理高光谱图像进行修正，将修正之后的待处理高光谱图像和所述基准图像作为修正版高光谱图像。
[0131] 由于行数修正模块22和图2中步骤S21的实现原理类似，因此在此不作更多说明。
[0132] 亮度调整模块23，用于针对每个修正版高光谱图像，获得该修正版高光谱图像的亮度和调整系数，根据所述调整系数对该修正版高光谱图像的亮度进行调整，得到调整版高光谱图像。
[0133] 由于亮度调整模块23和图2中步骤S22的实现原理类似，因此在此不作更多说明。
[0134] 光谱定标模块24，用于获得定标公式，针对每个调整版高光谱图像，采用所述定标公式对该调整版高光谱图像进行光谱定标处理，得到光谱定标图像。
[0135] 由于光谱定标模块24和图2中步骤S23的实现原理类似，因此在此不作更多说明。
[0136] 平滑处理模块25，用于针对每个光谱定标图像，获取该光谱定标图像的光谱反射率，根据所述光谱反射率对该光谱定标图像进行平滑处理，得到平滑高光谱图像。
[0137] 由于平滑处理模块25和图2中步骤S24的实现原理类似，因此在此不作更多说明。
[0138] 本发明实施例中的高光谱图像预处理方法、装置及电子设备，会对每张待处理高光谱图像在每个波段下的像元值进行亮度调整，基于大量高光谱图像的统计，分析传感器在推扫过程中出现的调整系数的变化趋势，作出对应的调整，保障的各个波段的总体亮度一致。对光谱定标后的图像进行平滑处理，根据波长范围和变化幅度进行灵活调整，既能有效的平滑信噪比低的区间的噪音，又能保留信噪比高的区间的光谱细节特征，解决了辐射畸变和光谱定标之后的光谱曲线失真的问题。
[0139] 在本发明实施例所提供的几个实施例中，应该理解到，所揭露的装置和方法，也可以通过其它的方式实现。以上所描述的装置和方法实施例仅仅是示意性的，例如，附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
[0140] 另外，在本发明各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。
[0141] 所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，电子设备10，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动
硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。
在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
[0142] 以上所述仅为本发明的可选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何
修改、等同替换、改进等，均应包含在本发明的保护范围之内。
