方位和品质对南疆果品空间特性光谱影响及反演方法研究
徐嘉翊1,3, 黄雪2,3, 罗华平1,3,*, 刘金秀1,3, 索玉婷1,3, 王长旭1,3
1.塔里木大学机械电气化工程学院, 新疆 阿拉尔 843300
2.塔里木大学植物科学学院, 新疆 阿拉尔 843300
3.新疆维吾尔自治区普通高等学校现代农业工程重点实验室, 新疆 阿拉尔 843300
*通讯作者 e-mail: luohuaping739@163.com

作者简介: 徐嘉翊, 1997年生, 塔里木大学机械电气化工程学院硕士研究生 e-mail: 1056939585@qq.com

摘要

高光谱无损检测技术在果品定量无损检测中应用广泛, 以冬枣、 红提、 香梨三种果品空间特性光谱为研究目标, 探索空间特性光谱的影响因素和反演方法, 为提高户外果品无损检测精度提供了一种新思路。 分别提取三种果品的光谱库并计算空间特性光谱, 依次使用马氏距离、 浓度残差等预处理方法以及竞争性自适应权重取样算法选取特征波长, 将处理后的三种果品空间特性光谱分别与品质(糖度、 水分)和方位(探测角、 方位角、 相位角)建模, 建模结果如下: 三种果品(按照冬枣、 红提、 香梨的顺序)与糖分模型的相关系数 r分别为: 0.853 3, 0.822 7和0.913 3; 水分模型的相关系数 r分别为: 0.741 3, 0.784 7和0.891 3; 探测角模型相关系数 r分别为: 0.985 6, 0.992 7和0.974 7; 方位角模型相关系数 r分别为: 0.941 8, 0.910 5和0.936 9; 相位角模型相关系数 r分别为: 0.960 9, 0.957 0和0.956 3。 可以看出, 不同果品方位模型相关性都明显高于品质模型相关性, 因此方位因素是影响空间特性光谱的主要原因。 使用Roujean模型和Walthall模型分别对不同方位的空间特性光谱进行反演, 反演结果如下: 使用Roujean模型反演三种果品(按照冬枣、 红提、 香梨的顺序)空间特性光谱时 R2分别为0.934 4, 0.928 1和0.830 6; r分别为0.990 2, 0.983 9和0.969 1; RMSEP分别为0.030 9, 0.048 7和0.062 7; 平均模型误差分别为7.27%, 11.02%和8.61%。 使用Walthall模型描述不同果品空间特性光谱时 R2分别为0.943 3, 0.859 7和0.839 0; r分别为0.991 8, 0.971 8和0.970 2; RMSEP分别为0.036 6, 0.066 1和0.068 7; 平均模型误差分别为6.19%, 15.40%和7.84%。 可以看出, Roujean模型可以很好的描述冬枣和红提的空间特性光谱, 也可以较好的描述香梨空间特性光谱; Walthall模型可以很好的描述冬枣空间特性光谱, 也可以较好的描述红提和香梨空间特性光谱。 综上所述, 在今后试验中可以使用Roujean模型反演红提和香梨的空间特性光谱, 使用Walthall模型反演冬枣的空间特性光谱, 进而提高户外果品户外果品无损检测精度。

关键词: 二向反射分布函数; 空间特性光谱; Roujean模型; Walthall模型
中图分类号:TS207.3 文献标志码:A
Effects of Orientation and Quality on Spatial Spectrum Characteristics of Fruits in Southern Xinjiang
XU Jia-yi1,3, HUANG Xue2,3, LUO Hua-ping1,3,*, LIU Jin-xiu1,3, SUO Yu-ting1,3, WANG Chang-xu1,3
1. College of Mechanical and Electrical Engineering, Tarim University, Alar 843300, China
2. College of Plant Science, Tarim University, Alar 843300, China
3. The Key Laboratory of Colleges and Universities under the Department of Education of Xinjiang Uygur Autonomous Region, Alar 843300, China
*Corresponding author
Abstract

Hyperspectral nondestructive testing technology is widely used in quantitative nondestructive testing of fruit. In this paper, the spatial characteristic spectra of jujube, grape and pear are taken as the research objectives, and the influencing factors and inversion methods of spatial characteristic spectra are explored, which provides a new idea for improving the accuracy of outdoor fruit nondestructive testing. The spectral library of three kinds of fruits was extracted, and the spatial characteristic spectra were calculated. The characteristic wavelengths were selected by Mahalanobis distance, concentration residual and competitive adaptive weight sampling algorithm. Model characteristic spectra of three kinds of fruits after pretreatment with quality parameters and positional parameters respectively.The modeling results are as follows: In the sugar model, the R of jujube, grape and pear were 0.853 3, 0.822 7 and 0.913 3 respectively; In the water model, the R were 0.741 3, 0.784 7 and 0.891 3 respectively; In the detection angle model, the R were 0.985 6, 0.992 7 and 0.974 7 respectively; In the azimuth angle model, the R were 0.941 8, 0.910 5 and 0.936 9 respectively; In the phase angle model, the R were 0.960 9, 0.957 0 and 0.956 3 respectively. In summary, the correlation of different fruit positional models was significantly higher than quality models. Therefore, the positional factor is the main reason affecting the characteristic spectrum. Therefore, the roujean model and waltall model are used to invert the characteristic spatial spectrum of different directions. The inversion results are as follows: roujean model is used when retrieving the spatial characteristic spectra of three kinds of fruits (in the order of jujube, grape and pear), R2 is 0.934 4, 0.928 1 and 0.830 6 respectively; R is 0.990 2, 0.983 9 and 0.969 1 respectively; RMSEP is 0.030 9, 0.048 7 and 0.062 7 respectively; the average model error is 7.27%, 11.02% and 8.61% respectively. The results showed that R2 was 0.943 3, 0.859 7, 0.839 0; R was 0.991 8, 0.971 8, 0.970 2; RMSEP was 0.036 6, 0.066 1, 0.068 7; the average model error was 6.19%, 15.40%, 7.84%. It can be seen that roujean model can well describe the characteristic spatial spectrum of jujube and grape, and also can better describe the characteristic spatial spectrum of pear; waltall model can well describe the characteristic spatial spectrum of jujube, and also can better describe the characteristic spatial spectrum of grape and pear. In conclusion, roujean model can be used to invert the characteristic spatial spectrum of grape and pear, and waltall model can be used to invert the characteristic spatial spectrum of jujube to improve the accuracy of outdoor fruit nondestructive testing.

Keyword: Bidirectional reflectivity distribution function; Background spectrum; Roujean model; Walthall model
引言

新疆是久负盛名的“ 瓜果之乡” , 水果品种繁多, 质地优良, 盛产石榴、 葡萄、 无花果、 杏、 梨、 哈密瓜等水果, 2019年水果产量1 118.72万吨, 同比增长5.6%, 位于全国前列, 其中红枣、 葡萄、 香梨的种植面积和产量均排全国第一, 特色林果业已经成为新疆农民致富的新的支柱产业。 高光谱成像技术近年来日益兴起, 相比较于传统的光谱技术, 高光谱成像技术在快速获得待测样品光谱信息的同时, 还可以获取样品的空间分布信息以及图像特征[1], 在食品、 农产品、 烟草等方面无损检测中应用广泛, 而如何提高检测速度和精度是当前的热点问题。 在户外检测时, 探测器接收到的能量主要是果品反射和散射后的光, 也有背景反射光和杂散光的影响, 而空间特性光谱可以描述这种影响。 本文以三种南疆特色果品的空间特性光谱为研究目标, 探索空间特性光谱的影响因素和反演方法, 从而达到提高检测精度的目的。

Roujean模型是一种表面双向反射率模型, 该模型遵循半经验方法, 被广泛应用于描述植被材料BRDF特征, 也适用于各类异质新表面。 Yang等[2]将其应用在研究聚四氟乙烯薄板稀释剂的双向性能上; Inoue等[3]将其用于评估纸张表面光泽度上; Luongo等[4]将其用于评估打印机控制表面反射率能力上; Zhang等[5]将其用于水下激光探测系统的设计中。 Yeom等[6]将其用于地面薄云探测以提高薄云像素探测效率; Farhad等[7]将其用于为卫星传感器开发一种新的交叉校准技术提高数据协调性。 Aval等[8]将其用于提高水下环境的目标识别。 Zhou等[9]将其应用在沙漠的研究中。

Walthall模型是一种具有互异性的经验模型, 由Walthall等提出, 以大量野外实验数据为基础形成的模型, 被广泛应用于模拟土壤的BRDF, 后也作为普通的BRDF模型用于其他地物。 Wu等[10]将其应用于雪地地表反射率研究; 有报道将其应用于MODIS资料进行BRDF订正的气溶胶光学厚度反演试验研究; Gong等[11]将其应用于多光谱传感器采集坡地栽培葡萄角度依赖性研究; Damian等[12]将其应用于无人机近红外图像序列双向反射分布函数的质量评价; Roujean模型和Walthall模型都属于线性模型, 其特点是模型参数的反演可以通过矩阵变换来实现, 因此反演速度较快, 并且这种反演还适用于混合像元, 因而使用两种模型对不同角度的空间特性光谱反演比较。

1 实验部分
1.1 材料

实验使用的冬枣、 红提、 香梨全部采集自新疆维吾尔自治区农一师阿拉尔市十团团场(北纬40° 35'10″, 东经81° 18'12″), 采集大小不同表面无损伤的冬枣、 红提、 香梨各120个, 分别进行擦拭、 标号、 称重等。

1.2 仪器

户外光谱采集使用四川双利合谱科技有限公司的Image-λ -N17E-N3 型增强型近红外高光谱相机, 其光谱范围900~1 750 nm、 256个波段、 光谱分辨率5 nm; 暗箱光谱采集使用北京卓立汉光公司的Hyperspspectral Sorting System推扫式高光谱分选系统; 以及HSIA-CT-150× 150型标准白板、 A2003型电子天平、 GZX-9140MBE电热鼓风干燥箱、 GMK-701R型手持式糖度计等。

1.3 高光谱图像采集

首先进行户外图像采集, 将冬枣、 红提、 香梨按从左到右的顺序摆放在试验架上, 在不同位置摆放两块标准白板进行定标, 如图1所示。 将高光谱相机调到与试验架同一高度并放置在试验架前1.5 m处接线调试, 测量户外温度21 ℃、 空气湿度55%、 东风2级。 实验开始于北京时间2020年10月1日13:00, 使用高光谱相机对试验架上的果品和白板进行扫描, 扫描时测量太阳高度角和方位角。 扫描之后使用激光测距仪测量相机高度、 样品高度、 到白板的距离角度、 到样品的距离角度, 以便计算入射角、 探测角、 方位角、 相位角。 之后进行暗箱图像采集, 将样品送入高光谱分选仪中进行扫描, 调节相机高度42 cm、 曝光时间11.5 ms、 传送带移动速度2.0 m· min-1, 冬枣和红提每次扫描20个样本, 香梨每次扫描6个样本。

图1 试验样品摆放图Fig.1 Layout of test samples

1.4 理化值测定

将样品扫描完光谱的一面切片后一分为二, 一半测量水分, 一半测量糖度。 切片时避开果核部分。

水分测定使用国标(GB/T 5009.3— 2016)的烘干减质量方法, 每个样品切取扫描光谱位置的果肉约3 g放入电热鼓风干燥箱中进行烘干, 温度102 ℃。 每隔两小时将样品称重一次, 直至重量变化小于0.001 g。

含水量=(m1-m2)m1×100%(1)

式(1)中: m1为烘干前样品质量; m2为烘干后样品质量。 糖度测定使用手持式数显糖度计。 测量开始前先使用校正液进行基准校正。

1.5 空间特性光谱提取

由于受到系统探测响应、 电子信号噪声、 环境温度变化等因素影响, 探测器采集到的图像不准确, 因此使用ENVI5.1软件对采集到的高光谱图像进行平场校正和黑白校正。 由于高光谱分选仪处于暗室, 并且有四个均匀分布的高亮度卤素灯照射, 因此将其测得的光谱近似看作样品真实光谱, 通过校正后的户外高光谱图像和高光谱分选仪图像计算每个样品的空间特性光谱并建立空间特性光谱库, 空间特性光谱的计算公式为

Sback=S1S2(2)

式(2)中: Sback为样品的空间特性光谱; S1为样品的户外光谱; S2为样品的暗室光谱。

对建立的空间特性光谱库进行包络线去除(continuum removal)以突出空间特性光谱的特征并减小计算量。 部分样品空间特性光谱如图2所示, 其中红色为冬枣空间特性光谱, 绿色为红提空间特性光谱, 蓝色为香梨空间特性光谱。

图2 部分样品空间特性光谱图Fig.2 Background spectra of some samples

1.6 数据预处理

在建模之前, 剔除异常样本, 因为其中可能包含异常光谱和异常理化值, 与其他样本存在显著差异, 如果异常样本参与建模, 必然会影响建模的精度和可靠性。 异常光谱使用马氏距离剔除, 异常理化值使用浓度残差剔除。 图3和图4分别是马氏距离和浓度残差方法剔除异常红枣光谱和糖度值, 分别设定阈值2.5和2, 共剔除14个异常波段和6个异常样本。

图3 马氏距离剔除异常冬枣光谱Fig.3 Eliminating abnormal spectra of winter Jujube by Mahalanobis distance

图4 浓度残差剔除异常冬枣糖度样本Fig.4 Eliminating abnormal sugar content samples of winter Jujube by concentration residual

由于高光谱数据的波长间往往会包含一些其他的多余信息, 增加大量的计算工作, 研究表明无关信息的存在会极大地降低模型的稳健性, 因此对全波波段的波长变量进行优选和精简十分必要。 研究中使用竞争性自适应权重取样法(CARS)选取特征波长, 图5为CARS方法选取冬枣糖分特征波长, 选取后剩余21个特征波长。

图5 CARS方法选取冬枣糖分特征波长
(a): 迭代次数与波长变量关系; (b): 不同迭代次数的RMSECV值; (c): 变量PLS回归系数值
Fig.5 CARS method selects the characteristic wavelength of winter Jujube sugar
(a): Relationship between interation number and wavelength variable; (b): RMSECV value for different interations; (c): Variable PLS regression coefficient value

2 结果与讨论
2.1 冬枣糖分模型建立

使用TQ Analyst软件进行建模, 采用偏最小二乘法(partial least squares, PLS), 光程分别使用恒定(constant)、 多元散射校正(multiplicative singal correction, MSC)、 标准正态变换(standard normal variate, SNV), 处理方法选择原始光谱(spectrum)、 一阶导数(frist derivative, 1st)、 二阶导数(second derivative, 2st)。 评价参数分别为: 相关系数(related coefficient, r)、 校正标准偏差(standard error of calibration, RMSEC)、 预测标准偏差(standard error of prediction, RMSEP), 计算公式为

r=i=1nxi-x-yi-y-i=1n(xi-x-)2i=1n(yi-y-)2(3)

式(3)中, xi, yi(i=1, 2, …, n)为两个变量xy的样本值; x-y-为两个变量的样本平均值; n为两个变量的样本值的个数。

RMSEC=i=1n(yi, actual-yi, qredicted)2n-1(4)

RMSEP=i=1n(yi, actual-yi, predicted)2m-1(5)

式(4)和式(5)中, yi, actual为第i样品参考方法的测定值, yi, predicted为用所建模型对校正集中第i样品的预测值, n为校正集的样品数, m为验证集的样本数。 RMSEP越小, 表明所建模型的预测能力越强, 通常RMSEP大于RMSEC。

将剩余的114个冬枣按糖度从大到小排列, 按照3∶ 1的比例依次划分校正集和预测集, 使校正集糖分梯度包含预测集糖分梯度, 以便有更好的预测效果。 分别使用Con, MSC, SNV, Con+1st, MSC+1st, SNV+1st, Con+2st, MSC+2st和SNV+2st处理方法对冬枣空间特性光谱与糖分建模, 建模结果如表1所示, 可以看出, 不经过处理的冬枣空间特性光谱建模相关系数r只有0.549 93, 而使用Con+1st的相关系数提升至0.853 31, 结果表明冬枣空间特性光谱与糖分有较好的相关性, 图6为空间特性光谱与糖分Con+1st建模的结果。

表1 不同处理方法下冬枣空间特性光谱与糖分建模结果 Table 1 Results of background spectrum and sugar modeling of winter Jujube under different treatment methods

图6 冬枣空间特性光谱与糖分con+1st建模结果
(a): 模型样本回归; (b): 模型样本误差
Fig.6 Background spectra of winter jujube and results of con+1st modeling of sugar content
(a): Model sample regression; (b): Model sample error

2.2 冬枣理化值模型建立

依次使用马氏距离、 浓度残差、 CARS算法对水分、 探测角、 方位角、 相位角进行预处理, 按理化值梯度排序并划分校正集后使用以上九种处理方法建模, 分别选取建模效果最好的三组处理方法列于表2。 可以看出, 冬枣的内部品质(糖度、 水分)r可以达到0.853 3和0.741 3, 而探测方位(探测角、 方位角、 相位角)r可以达到0.985 6, 0.941 8和0.960 9。 可以看出, 冬枣的内部品质与空间特性光谱建模效果较好, 而探测方位与空间特性光谱有更好的建模效果。

表2 冬枣不同理化值建模结果 Table 2 Modeling results of different winter jujube physicochemical values
2.3 红提、 香梨理化值模型建立

使用相同的方法对红提和香梨空间特性光谱进行预处理。 马氏距离阈值设定为2.5, 剔除了红提11个波段和香梨11个波段。 浓度残差阈值设定为2, 分别剔除了8, 6, 5, 6和7个红提异常理化值样本和8, 6, 5, 6和7个香梨异常理化值样本。 使用CARS选取特征波长后导入TQ Analyst划分校正集和预测集, 进行PLS建模, 建模结果如表3表4所示。

表3 红提不同理化值建模结果 Table 3 Modeling results of different Red grapes physicochemical values
表4 香梨不同理化值建模结果 Table 4 Modeling results of different Fragrant Pear physicochemical values

对比冬枣、 红提、 香梨三种果品的建模结果可以看出, 红提和香梨的内部品质r分别为0.822 7, 0.784 7和0.913 3, 0.891 3, 而方位r分别为0.992 7, 0.910 5, 0.957 0和0.974 7, 0.936 9, 0.956 3。 对比发现, 果品空间特性光谱的方位模型相关性明显好于品质模型。 在品质模型中, 糖度模型相关性都高于水分模型。 而在方位模型中, 探测角模型相关性高于相位角模型高于方位角模型。

2.4 Roujean模型反演空间特性光谱

由于空间特性光谱与探测角、 方位角、 相位角建模都有很好的相关性, Roujean模型和Walthall模型都与BRDF角度相关, 因此使用两种模型对空间特性光谱进行反演。 模型反演评价参数主要有相关系数r、 决定系数R2、 预测标准差RMSEP、 模型误差, 计算公式如式(6)

R2=1-i=1n(yi, actual-yi, predicted)2i=1n(yi, actual-y-i, actual)2(6)

式(6)中, yi, actual为第i样品参考方法的测定值, y-actual为校正集或验证集所有样品参考方法测定值的平均值, yi, predicted为用所建模型对校正集中第i样品的预测值, n为校正集的样品数。 在浓度范围相同的前提下, r越接近1, 回归或预测效果应越好。

由于户外实验光源是太阳光, 太阳光是平行光, 因此入射角都相同, 经计算均为51.52° 。 由于试验架限制, 探测角和方位角梯度变化不是很明显, 因此每种果品均选择摆放位置较远的8个样品来反演。 不同果品1号样品反演结果如图7所示, 可以看出每种果品反演都有较好地精度, 但在1 000~1 350 nm波段反演效果并不理想, 全部反演结果如表5所示。

图7 Roujean模型反演结果对比
(a): 冬枣; (b): 红提; (c): 香梨
Fig.7 Comparison of inversion results of Roujean model
(a): Winter jujube; (b): Red grapes; (c): Fragrant pear

表5 Roujean模型反演结果对比 Table 5 Comparison of inversion results of roujean model

经计算, Roujean模型描述冬枣、 红提、 香梨空间特性光谱时平均R2分别为0.934 4, 0.928 1和0.830 6; 平均r分别为0.990 2, 0.983 9和0.969 1; 平均RMSEP分别为0.030 9, 0.048 7和0.062 7; 平均模型误差分别为7.27%, 11.02%和8.61%; 单波段最大模型误差分别为28.7%, 42.92%和48.21%; 单波段最小模型误差分别为0, 0和0.01%。可以看出, Roujean模型反演空间特性光谱精度从高到低依次为冬枣、 红提、 香梨。 综上所述, Roujean可以很好地描述冬枣和红提的空间特性光谱, 可以较好地描述香梨空间特性光谱。

2.5 Walthall模型反演空间特性光谱

使用Walthall模型对三种果品反演, 不同果品1号样品反演结果如图8(a, b, c)所示, 全部反演结果如表6所示。

图8 Walthall模型反演结果对比
(a): 冬枣; (b): 红提; (c): 香梨
Fig.8 Comparison of inversion results of Walthall model
(a): Winter jujube; (b): Red grapes; (c): Fragrant pear

表6 Walthall模型反演结果对比 Table 6 Comparison of inversion results of Walthall model

经计算, Walthall模型描述冬枣、 红提、 香梨空间特性光谱时平均R2分别为0.943 3, 0.859 7和0.839 0; 平均r分别为0.991 8, 0.971 8和0.970 2; 平均RMSEP分别为0.036 6, 0.066 1和0.068 7; 平均模型误差分别为6.19%, 15.40%和7.84%; 单波段最大模型误差分别为28.73%, 44.91%和20.45%; 单波段最小模型误差分别为0, 0.01%和0.01%。 可以看出, Walthall模型反演空间特性光谱精度从高到低依次为冬枣、 红提、 香梨。 综上所述, Walthall可以很好地描述冬枣空间特性光谱, 可以较好地描述红提和香梨空间特性光谱。

2.6 Roujean模型与Walthall模型反演结果对比

Roujean模型与Walthall模型反演结果对比见图9(a, b, c)和图10(a, b, c)。

图9 Roujean模型反演不同果品误差对比
(a): 冬枣; (b): 红提; (c): 香梨
Fig.9 Error comparison of different fruits retrieved by Roujean model
(a): Winter jujube; (b): Red grapes; (c): Fragrant pear

图10 Walthall模型反演不同果品误差对比
(a): 冬枣; (b): 红提; (c): 香梨
Fig.10 Error comparison of different fruits retrieved by Walthall model
(a): Winter jujube; (b): Red grapes; (c): Fragrant pear

3 结论

(1)对比冬枣、 红提、 香梨三种果品的Walthall建模结果可以看出, 果品空间特性光谱与探测方位有很强的相关性, 而与果品品质相关性较弱。 在户外果品近红外无损检测时, 探测器接收的能量主要有三部分: ① 果品反射和散射后到达探测器的太阳光, 包括果品表面镜面反射以及内部体散射的能量; ② 背景反射和散射后到达探测器的太阳光; ③ 大气散射后直接进入探测器的太阳光。 以上三种能量均与入射、 探测方位有关, 而只有果品镜面散射和体散射与果品品质有关, 说明角度模型的相关性好于品质模型的相关性。 因此, 在户外果品无损检测时, 要重点考虑角度对光谱的影响。

(2)Roujean可以很好地描述冬枣和红提的空间特性光谱, 可以较好地描述香梨空间特性光谱; Walthall可以很好地描述冬枣空间特性光谱, 可以较好地描述红提和香梨空间特性光谱。 两种模型均在980~1 120, 1 180~1 330和1 500~1 670 nm波长范围拟合效果最好, 1 340~1 440 nm范围效果最差。 由于在1 400 nm附近有水份吸收峰, 因此导致该范围拟合效果较差。

(3)由于单次试验太阳入射角单一, 样品摆放时探测角梯度不明显等原因会影响试验精度。 后期可通过多次试验、 增加试验样本量、 增加试验角度分布、 使用不同方位模型(如Rahman模型、 Ross-Li模型、 指数模型等)等方法对多角度下的空间特性光谱进行反演分析, 以确定适合每种果品的模型, 分别建立每种果品不同方位下的空间特性光谱库, 通过空间特性光谱来反演果品光谱, 为提高户外果品高光谱无损检测精度提供了一种新思路。

参考文献
[1] WANG Dong, WANG Kun, WU Jing-zhu, et al(王冬, 王坤, 吴静珠, ). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2021, 41(1): 52. [本文引用:1]
[2] Yang Peiyan, Zhuo min, Zhang M. International Journal of Heat and Mass Transfer, 2020, 148: 118991. [本文引用:1]
[3] Shinichi Inoue, Masanori Maki, Norimichi Tsumura. JAPAN TAPPI JOURNAL, 2020, 73(10): 1022. [本文引用:1]
[4] Luongo A, Falster V, Doest M B, et al. Computer Graphics Forum, 2020, 39(1): 347. [本文引用:1]
[5] ZHANG Ying-luo, WANG Ying-min, HUANG Ai-ping, et al(张英珞, 王英民, 黄爱平, ). Optoelectronics Letters(光电子快报), 2020, 16(2): 137. [本文引用:1]
[6] Jong-Min Yeom, Jean-Louis Roujean, Kyung-Soo Han, et al. Remote Sensing of Environment, 2020, 239: 111610. [本文引用:1]
[7] Farhad M M, Morakot Kaewmanee, Larry Leigh, et al. Remote Sensing, 2020, 12: 806. [本文引用:1]
[8] Aval J, Alfalou A, Brosseau C. Optical Society of America, 2019, 11(2): 10. [本文引用:1]
[9] Zhou Qu, Tian Liqiao, Li Jian, et al. ISPRS Journal of Photogrammetry and Remote Sensing, 2020, 160(12): 733. [本文引用:1]
[10] WU Hong-yi, TONG Ling(吴宏伊, 童玲). Chinese Optics Letters, 2011, 9(10): 83. [本文引用:1]
[11] Gong Chizhang, Henning Budden baum, Rebecca Retzlaff, et al. Remote Sensing, 2019, 11(21): 2561. [本文引用:1]
[12] Wierzbicki D, Kedzierski M, Fryskowska A, et al. Remote Sensing, 2018, 10(9): 1348. [本文引用:1]