留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

Classification of coal/rock based on Hyperspectral LiDAR calibration-free signals

He Zixin Shao Hui Guo Hang Chen Jie

何子辛, 邵慧, 郭航, 陈杰. 基于高光谱激光雷达信号强度免校准的煤岩分类[J]. 红外与激光工程, 2021, 50(10): 20200518. doi: 10.3788/IRLA20200518
引用本文: 何子辛, 邵慧, 郭航, 陈杰. 基于高光谱激光雷达信号强度免校准的煤岩分类[J]. 红外与激光工程, 2021, 50(10): 20200518. doi: 10.3788/IRLA20200518
He Zixin, Shao Hui, Guo Hang, Chen Jie. Classification of coal/rock based on Hyperspectral LiDAR calibration-free signals[J]. Infrared and Laser Engineering, 2021, 50(10): 20200518. doi: 10.3788/IRLA20200518
Citation: He Zixin, Shao Hui, Guo Hang, Chen Jie. Classification of coal/rock based on Hyperspectral LiDAR calibration-free signals[J]. Infrared and Laser Engineering, 2021, 50(10): 20200518. doi: 10.3788/IRLA20200518

基于高光谱激光雷达信号强度免校准的煤岩分类

doi: 10.3788/IRLA20200518
详细信息
  • 中图分类号: TN247

Classification of coal/rock based on Hyperspectral LiDAR calibration-free signals

More Information
    Author Bio:

    何子辛,男,硕士生,主要从事激光雷达数据处理方面的研究

    Corresponding author: 邵慧,女,副教授,博士,主要从事激光雷达与数字图像处理方面的研究。
  • 摘要: 安全性对于深井开采至关重要。基于激光雷达的扫描探测技术可以有效监测巷道和深部煤矿现场的围岩安全状况。新兴的高光谱激光雷达不仅可以提供空间几何信息还可以提供丰富的光谱数据,在深井煤矿安全检测和精细结构分析方面具有良好的应用前景,而精确的煤岩分类是监测分析的基础。在实际应用中,雷达强度信号易受仪器属性和环境因素的影响,需校准才能使用。由于深井煤矿粉尘污染严重,常规校准方法难以达到理想效果。针对这个问题,提出一种信号强度免校准的方法,从激光雷达回波信号中提取新特征实现煤岩精确分类。首先,使用高光谱激光雷达获得煤/岩石样本的回波强度信息,并计算出波形熵(WE)和联合偏斜度-峰度系数(JSKF)作为新分类特征参数。其次,采用随机森林(RF)与支持向量机(SVM)分类器实现煤/岩石分类。最后,笔者进行了光谱分段测试,对特征分类性能进行优化。结果表明,所提的免校准方法,提高高光谱激光雷达直接应用能力的同时能够保持良好的分类性能。
  • Figure  1.  Schematic diagram of HSL

    Figure  2.  (a) Waveform distribution diagram of WE; (b) Waveform distribution diagram of JSKF

    Table  1.   Specification of AOTF-HSL

    ParameterValue
    Spectral range/nm650-1100
    Spectral resolution/nm5
    Output efficiency>40%
    PolarizationLine polarization
    Beam divergence/mrad0.4
    下载: 导出CSV

    Table  2.   Full channels classification accuracy

    DataRFSVM
    AccuracyAccuracy
    Intensity100%100%
    WE100%100%
    JSKF100%100%
    Calibrated100%100%
    下载: 导出CSV

    Table  3.   Minimum number of random channels

    DataRFSVM
    Number of channelsNumber of channels
    Intensity2120
    WE1515
    JSKF66
    Calibrated64
    下载: 导出CSV

    Table  4.   5-channel classification accuracy

    DataRFSVM
    AccuracyAccuracy
    WE88%90%
    JSKF96%100%
    Calibrated98%100%
    下载: 导出CSV

    Table  5.   720-820 nm classification results

    ClassifierRFSVM
    ChannelsAccuracyChannelsAccuracy
    WE2095.02%2095.34%
    JSKF3100%4100%
    Calibrated2100%3100%
    下载: 导出CSV

    Table  6.   1000-1100 nm classification results

    ClassifierRFSVM
    ChannelsAccuracyChannelsAccuracy
    WE11100%7100%
    JSKF19100%7100%
    Calibrated10100%6100%
    下载: 导出CSV
  • [1] Qun Feng, Hong Chen. The safety-level gap between China and the US in view of the interaction between coal production and safety management [J]. Safety Science, 2013, 54: 80-86. doi:  10.1016/j.ssci.2012.12.001
    [2] Qi Qingjie, Zhao Xiaoliang, Song Baichao. Pre-evaluation method of coal mine safety based on continental distance model with varying weight [J]. Procedia Earth and Planetary Science, 2009, 1(1): 180-185. doi:  10.1016/j.proeps.2009.09.030
    [3] Song Weihua, Zhang Hongwei. Regional prediction of coal and gas outburst hazard based on multi-factor pattern recognition [J]. Procedia Earth and Planetary Science, 2009, 1(1): 347-353. doi:  10.1016/j.proeps.2009.09.055
    [4] Khatik V M, Nandi A K. A generic method for rock mass classification [J]. Journal of Rock Mechanics and Geotechnical Engineering, 2018, 10(1): 106-120.
    [5] Weidner L, Walton G, Kromer R A. Classification methods for point clouds in rock slope monitoring: A novel machine learning approach and comparative analysis [J]. Engineering Geology, 2019, 263: 105326. doi:  10.1016/j.enggeo.2019.105326
    [6] Ryan A K, D Jean H, Lato M J, et al. Identifying rock slope failure precursors using LiDAR for transportation corridor hazard management [J]. Engineering Geology, 2015, 195: 93-103. doi:  10.1016/j.enggeo.2015.05.012
    [7] Hartzell P, Glennie C, Biber K, et al. Application of multispectral LiDAR to automated virtual outcrop geology [J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2014, 88: 147-155.
    [8] Chen Yuwei, Jiang Changhui, Hyyppa Juha, et al. Feasibility study of ore classification using active hyperspectral LiDAR [J]. IEEE Geoscience and Remote Sensing Letters, 2018, 15(11): 1785-1789. doi:  10.1109/LGRS.2018.2854358
    [9] Liu Bingyi, Zhuang Quanfeng, Qin Shengguang, et al. Aerosol classification method based on high spectral resolution lidar [J]. Infrared and Laser Engineering, 2017, 46(4): 0411001. (in Chinese) doi:  10.3788/IRLA201746.0411001
    [10] Dong Junfa, Liu Jiqiao, Zhu Xiaolei, et al. Splitting ratio optimization of spaceborne high spectral resolution lidar [J]. Infrared and Laser Engineering, 2019, 48(S2): S205001. (in Chinese) doi:  10.3788/IRLA201948.S205001
    [11] Xu Junjie, Bu Lingbing, Liu Jiqiao, et al. Airborne high-spectral-resolution lidar for atmospheric aerosol detection [J]. Chinese Journal of Lasers, 2020, 47(7): 0710003. doi:  10.3788/CJL202047.0710003
    [12] Morsdorf, Nichol, Malthus, et al. Assessing forest structural and physiological information content of multi-spectral LiDAR waveforms by radiative transfer modelling [J]. Remote Sensing of Environment, 2009, 113(10): 2152-2163. doi:  10.1016/j.rse.2009.05.019
    [13] Chen Yuwei, Räikkönen E, Kaasalainen S, et al. Two-channel hyperspectral LiDAR with a supercontinuum laser source [J]. Sensors, 2010, 10(7): 7057-7066. doi:  10.3390/s100707057
    [14] Sanna K, Anttoni J, Mikko K, et al. Analysis of incidence angle and distance effects on terrestrial laser scanner intensity: search for correction methods [J]. Remote Sensing, 2011, 3(10): 2207-2221. doi:  10.3390/rs3102207
    [15] Balduzzi M A F, Dimitry V D Z, Stuckens J, et al. The properties of terrestrial laser system intensity for measuring leaf geometries: A case study with conference pear trees (Pyrus Communis) [J]. Sensors, 2011, 11(2): 1657-1681. doi:  10.3390/s110201657
    [16] Hoefle B, Pfeifer N. Correction of laser scanning intensity data: Data and model-driven approaches [J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2007, 62(6): 415-433. doi:  10.1016/j.isprsjprs.2007.05.008
    [17] Yan W Y, Shaker A, Habib A, et al. Improving classification accuracy of airborne LiDAR intensity data by geometric calibration and radiometric correction [J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2012, 67: 35-44.
    [18] Shao Hui, Chen Yuwei, Yang Zhirong, et al. A 91-channel hyperspectral LiDAR for coal/rock classification [J]. IEEE Geoscience and Remote Sensing Letters, 2020, 17(6): 1052-1056. doi:  10.1109/LGRS.2019.2937720
    [19] Chen Yuwei, Li Wei, Juha H, et al. A 10-nm spectral resolution hyperspectral LiDAR system based on an acousto-optic tunable filter [J]. Sensors, 2019, 19(7): 1620. doi:  10.3390/s19071620
    [20] CoifmanR R, Wickerhauser M V. Entropy-based algorithms for best basis selection [J]. IEEE Transactions on Information Theory, 1992, 38(2): 713-718. doi:  10.1109/18.119732
    [21] Diner D J, Xu F, GarayM J, et al. The Airborne Multiangle Spectro Polarimetric Imager (AirMSPI): A new tool for aerosol and cloud remote sensing [J]. Atmospheric Measurement Techniques, 2013, 6(8): 2007-2025. doi:  10.5194/amt-6-2007-2013
    [22] Chen Jikai, Li Guoqing. Tsallis wavelet entropy and its application in power signal analysis [J]. Entropy, 2014, 16(6): 3009-3025. doi:  10.3390/e16063009
    [23] Hovi A, Korhonen L, Vauhkonen J, et al. LiDAR waveform features for tree species classification and their sensitivity to tree- and acquisition related parameters [J]. Remote Sensing of Environment, 2016, 173: 224-237. doi:  10.1016/j.rse.2015.08.019
    [24] Su Dianpeng, Yang Fanlin, Ma Yue, et al. Classification of coral reefs in the south China sea by combining airborne LiDAR bathymetry bottom waveforms and bathymetric features [J]. IEEE Transactions on Geoscience and Remote Sensing, 2019, 57(2): 815-828. doi:  10.1109/TGRS.2018.2860931
    [25] Gu Zhujun, Cao Sen, Sanchez-Azofeifa G A. Using LiDAR waveform metrics to describe and identify successional stages of tropical dry forests [J]. ITC Journal, 2018, 73: 482-492.
    [26] Guo Lei, Chang Weiwei, Fu Chaoyang. Band selection of optimal for hyperspectral image fusion [J]. Journal of Astronautics, 2011, 32(2): 374-379.
    [27] Chan C W, Desiré P. Evaluation of random forest and adaboost tree-based ensemble classification and spectral band selection for ecotope mapping using airborne hyperspectral imagery [J]. Remote Sensing of Environment, 2008, 112(6): 2999-3011. doi:  10.1016/j.rse.2008.02.011
    [28] Ghasemian N, Akhoondzadeh M. Introducing two Random Forest based methods for cloud detection in remote sensing images [J]. Advances in Space Research, 2018, 62(2): 288-303. doi:  10.1016/j.asr.2018.04.030
    [29] Zhao Chunhui, Gao Bing, Zhang Lejun, et al. Classification of hyperspectral imagery based on spectral gradient, SVM and spatial random forest [J]. Infrared Physics and Technology, 2018, 95: 61-69.
    [30] Emma I V, Raúl Z M. An evaluation of guided regularized random forest for classification and regression tasks in remote sensing [J]. International Journal of Applied Earth Observations and Geoinformation, 2020, 88: 102051. doi:  10.1016/j.jag.2020.102051
    [31] Pedregosa F, Varoquaux G, Gramfort A, et al. Scikit-learn: Machine learning in python [J]. Journal of Machine Learning Research, 2011, 12: 2825-2830.
    [32] Shao Hui, Chen Yuwei, Li Wei, et al. An investigation of spectral band selection for hyperspectral LiDAR technique [J]. Electronics, 2020, 9(9): 148.
    [33] Yan Shouxun, Zhang Bing, Zhao Yongchao, et al. Summarizing the VIS-NIR spectra of minerals and rocks [J]. Remote Sensing Technology and Application, 2003, 18(4): 191-201.
  • [1] 田汶鑫, 陈育伟, 唐伶俐, 李子扬, 邱实, 吴昊昊, 张慧静, 陈林生, 蒋长辉, 胡佩纶, 贾建鑫, 孙海彬, 王一程, 胡以华.  高光谱激光雷达后向散射强度的粗糙表面二向反射模型 . 红外与激光工程, 2023, 52(10): 20230108-1-20230108-14. doi: 10.3788/IRLA20230108
    [2] 邵慧, 撒贝宁, 李伟, 陈育伟, 刘璐, 陈杰, 孙龙, 胡玉霞.  一种适用于古建筑建模的全波形高光谱激光雷达设计与实现 . 红外与激光工程, 2022, 51(8): 20210786-1-20210786-10. doi: 10.3788/IRLA20210786
    [3] Fu Yanjun, Cai Xiaoqi, Zhong Kejun, Ma Baiheng, Yan Zhanjun.  Method for phase-height mapping calibration based on fringe projection profilometry . 红外与激光工程, 2022, 51(4): 20210403-1-20210403-9. doi: 10.3788/IRLA20210403
    [4] Yang Sen.  Hybrid PSO-AMLS-based method for data fitting in the calibration of the infrared radiometer . 红外与激光工程, 2021, 50(8): 20200471-1-20200471-9. doi: 10.3788/IRLA20200471
    [5] 刘舒扬, 张晨, 赵安娜, 李奇峰, 贾晓东.  基于光谱特性的高判别准确度激光近场探测传感器研究 . 红外与激光工程, 2020, 49(4): 0403009-0403009-5. doi: 10.3788/IRLA202049.0403009
    [6] 胡善江, 贺岩, 陶邦一, 俞家勇, 陈卫标.  基于深度学习的机载激光海洋测深海陆波形分类 . 红外与激光工程, 2019, 48(11): 1113004-1113004(8). doi: 10.3788/IRLA201948.1113004
    [7] 张东阁, 傅雨田.  基于一类支持向量机的盲元检测方法 . 红外与激光工程, 2018, 47(4): 404001-0404001(7). doi: 10.3788/IRLA201847.0404001
    [8] 马跃, 张文豪, 张智宇, 马昕, 李松.  基于半解析模型的激光测高回波海水海冰波形分类方法 . 红外与激光工程, 2018, 47(5): 506005-0506005(7). doi: 10.3788/IRLA201847.0506005
    [9] 刘志青, 李鹏程, 郭海涛, 张保明, 陈小卫, 丁磊, 赵传.  基于相关向量机的机械LiDAR点云数据分类 . 红外与激光工程, 2016, 45(S1): 98-104. doi: 10.3788/IRLA201645.S130006
    [10] 李峰, 梁汉东, 米晓楠, 卫爱霞.  Landsat8卫星影像的多子区决策树土地覆被分类方法 . 红外与激光工程, 2015, 44(7): 2224-2230.
    [11] 翟文超, 徐骏, 郑小兵, 丁蕾, 陆俊桦.  应用于超光谱定标的光谱可调光源光学设计 . 红外与激光工程, 2014, 43(3): 950-955.
    [12] 石轶, 刘常杰, 郭寅, 叶声华, 石松.  基于双目视觉的接触网几何参数测量系统 . 红外与激光工程, 2014, 43(6): 1936-1942.
    [13] 徐骏, 孟炳寰, 翟文超, 丁蕾, 郑小兵.  基于热红外标准辐亮度计的常温黑体定标技术 . 红外与激光工程, 2014, 43(3): 716-721.
    [14] 李斌, 吴海英, 王文涛, 刘鹏.  一种简易红外标定系统的设计 . 红外与激光工程, 2014, 43(2): 458-463.
    [15] 普晗晔, 王斌, 张立明.  基于流形学习的新高光谱图像降维算法 . 红外与激光工程, 2014, 43(1): 232-237.
    [16] 杨国洪, 李军, 韦敏习, 侯立飞, 易涛, 刘慎业.  X射线分光晶体高阶衍射效率标定方法 . 红外与激光工程, 2013, 42(12): 3325-3329.
    [17] 孙高飞, 张国玉, 高玉军, 王凌云, 苏拾, 付芸, 王向东.  星敏感器地面标定设备的设计 . 红外与激光工程, 2013, 42(10): 2822-2827.
    [18] 张晓娟, 张国玉, 孙高飞, 刘石, 郑茹, 高玉军.  大视场LCOS拼接星模拟器标定方法 . 红外与激光工程, 2013, 42(11): 2996-3001.
    [19] 陈应航, 隋左宁, 孙红胜, 魏建强, 李敬峰, 李世伟, 孙广尉.  紫外-真空紫外光谱辐照度校准系统 . 红外与激光工程, 2013, 42(4): 998-1002.
    [20] 徐骏, 孟炳寰, 郑小兵, 翟文超, 李健军.  红外标准辐亮度计的研制及定标 . 红外与激光工程, 2013, 42(4): 863-868.
  • 加载中
图(2) / 表(6)
计量
  • 文章访问数:  247
  • HTML全文浏览量:  98
  • PDF下载量:  28
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-12-09
  • 修回日期:  2021-03-10
  • 刊出日期:  2021-10-20

Classification of coal/rock based on Hyperspectral LiDAR calibration-free signals

doi: 10.3788/IRLA20200518
    作者简介:

    何子辛,男,硕士生,主要从事激光雷达数据处理方面的研究

    通讯作者: 邵慧,女,副教授,博士,主要从事激光雷达与数字图像处理方面的研究。
  • 中图分类号: TN247

摘要: 安全性对于深井开采至关重要。基于激光雷达的扫描探测技术可以有效监测巷道和深部煤矿现场的围岩安全状况。新兴的高光谱激光雷达不仅可以提供空间几何信息还可以提供丰富的光谱数据,在深井煤矿安全检测和精细结构分析方面具有良好的应用前景,而精确的煤岩分类是监测分析的基础。在实际应用中,雷达强度信号易受仪器属性和环境因素的影响,需校准才能使用。由于深井煤矿粉尘污染严重,常规校准方法难以达到理想效果。针对这个问题,提出一种信号强度免校准的方法,从激光雷达回波信号中提取新特征实现煤岩精确分类。首先,使用高光谱激光雷达获得煤/岩石样本的回波强度信息,并计算出波形熵(WE)和联合偏斜度-峰度系数(JSKF)作为新分类特征参数。其次,采用随机森林(RF)与支持向量机(SVM)分类器实现煤/岩石分类。最后,笔者进行了光谱分段测试,对特征分类性能进行优化。结果表明,所提的免校准方法,提高高光谱激光雷达直接应用能力的同时能够保持良好的分类性能。

English Abstract

    • Coal is an important energy in the world and it plays an important role in the development of the economy. Although a lot of new energy sources are applied in production and living, it will still occupy the dominant position for a long time[1]. As the demand for coal increases and surface resources decrease, deep mines are concerned and carried out. However, disasters caused by deep mining threaten safety operation seriously. Many methods had been implemented to ensure the operation safety, such as safety pre-evaluation[2], regional pre-diction[3], rock mass classification[4] and so on. Whereas, traditional geological monitoring methods usually require extensive field work, which is time-consuming and labor-intensive. At present, remote sensing has become one of the most effective tools for monitoring the condition of underground spaces. Rock classification based on light detection and ranging (LiDAR) technology attracted wide attention due to its high classification accuracy and good recognition precision[5].

      LiDAR is an effective remote sensing technique that it can not only assess stability of rock conditions accurately[6], but also collect surface geometry information of rock and distinguish different rock properties by acquiring 3D point cloud data[7]. Whereas, traditional LiDAR sensors operate at single or several wavelengths that can only obtain limited spectral information, which restrict the performance of application[8]. With development of remote sensing, the emergence of hyperspectral LiDAR (HSL), fusion of hyperspectral data and LiDAR, significantly enhanced the quantitative and qualitative analysis capabilities of LiDAR, which has been widely used in atmospheric detection, forest protection and artificial target detection[9-13].

      However, the performance of HSL depends on the quality of intensity signal received by LiDAR receiver. Many aspects impact the accuracy of LiDAR intensity signals in on-site applications, such as, instrument properties[14] and environmental factors[15]. The common solution is to calibrate the intensity signal to eliminate data deviation. In 2007, Hoefle et al. calibrated the intensity signal using data-driven and model-driven correction, which provided support for surface classification and multi-temporal analysis[16]. In 2011, Kaasalainen et al. studied the effect of range and incidence angle on intensity signals, conducted instrument calibration and target surface features calibration researches[14]. In the same year, Yan et al. evaluated the influence of geometric calibration and radiometric calibration of LiDAR signal on land surface classification, improved the accuracy of LiDAR data by eliminating parameter deviations and correcting the scanning angle[17]. Most calibration methods need a whiteboard with a calculated reflectance as reference to correct LiDAR signals. While the environment of signal collection site is seriously polluted and the reference whiteboard is easily covered by dust, which leads to inaccurate calculations and is hard to ensure the accuracy of calibration.

      To address this issue, we propose a new method to classify coal/rock samples without calibration. The instrument we employed is previously proposed an HSL based on acousto-optic tunable filter (AOTF)[18], which offers a quicker tuning speed and broader wavelength ranges. First, waveform entropy (WE) and joint skewness-kurtosis figure (JSKF) are calculated based on HSL intensity signals, as new extracted features. Then, coal/rock samples are classified with WE and JSKF by random forest (RF) and support vector machine (SVM) classifiers and the classification results are compared with the intensity signals. Finally, with spectral segmentation test, the classification performances are optimized by selecting the fewer channels and compared results with calibrated intensity signal.

    • The instrument we employed is the previously designed AOTF-HSL with a spectral 5 nm resolution, using a supercontinuum laser source (YSL® SC-OEM)[18], as Fig.1 shown. A computer triggers laser source to transmit supercontinuum laser pulses and filtered by the AOTF device to select specific wavelength. The specifications of AOTF module can be found in Tab.1. With beam expander, the laser pulses are collimated to transmit to the target. The echo signals are collected by the optical receiving module and focused them on avalanche photodiode (APD) model, which converts the laser echo signals into electronic signals and amplifies them. The output signals of APD are sampled and converted to digital signals by a high-speed A/D converter and transmitted to computer for data processing.

      Figure 1.  Schematic diagram of HSL

      Table 1.  Specification of AOTF-HSL

      ParameterValue
      Spectral range/nm650-1100
      Spectral resolution/nm5
      Output efficiency>40%
      PolarizationLine polarization
      Beam divergence/mrad0.4
    • The measurements were conducted in a controlled laboratory environment to obtain hyperspectral information through AOTF-HSL. Due to the low sensitivity of APD and the low transmitted power intensity of supercontinuum lasers below 650 nm[19], the measurement of the spectrum channels was selected in the range of 650 nm to 1100 nm. The waveforms of pulses reflected from coal/rock samples were collected by an oscilloscope sampling at 20 GHz.

      Besides coal and gangue-rock, we also added rock from roof layer and rock from floor layer as samples, which come from deep mines[18]. All samples were scanned at the same distance from the instrument that range differences between samples can be neglectable.

    • The precision of the intensity calibration of HSL system is a prerequisite for many quantitative applications, and it has become an important research object[14]. The equation (1) describes parameters related to the signal received by the LiDAR sensor:

      $$ {P}_{\rm r}=\frac{{P}_{\rm t}{D}_{\rm r}^{2}}{4\pi {R}^{4}{\beta }_{\rm t}^{2}}\sigma $$ (1)

      Where, Pr is the power receiver, Pt is the transmitted power, Dr is the receiver optical aperture diameter, R is the distance, βt denotes transmitter beam width, and σ is the backscatter cross section.

      The backscatter cross section can be expressed equation (2):

      $$ \sigma =\frac{4\pi }{\varOmega }\rho {A}_{S} $$ (2)

      Where, Ω is a cone of a solid angle into which the incoming radiation is scattered uniformly, ρ is the surface reflectance, and AS denotes the scattering receiving area.

      The equations indicate that the accuracy of intensity signals not only depends on distance, scattering cross-section, and atmospheric parameters, but also relates to reflectance of targets and directionality of the incident angle[15]. Therefore, it is necessary to calibrate intensity signals. An efficient calibration method is to set a whiteboard with a constant reflectance for correction[18]. Whereas, in deep coal mines, the conventional calibration method is hard to valid because of dust pollution.

    • The entropy concept is introduced to measure uncertainty of random variables[20]. A metric of the dispersion degree of echo signals with respect to their variables is selected as a new feature, namely waveform entropy.

      While the different target is present, the entropy value corresponds different, further, WE parameter of the same target is different under different spectral laser scanning[21]. Therefore, we can distinguish different type sequences based on WE, which is vary with scatter distribution, wavelength and other features[21-22].

      In order to adopt the concept of entropy, we suppose digital echo signals as Xjt= xjt(i), i= n,···,N, t=1,···4, j=1,···m. i denotes the echo signal from initial point to the ending point, t is the different type of coal sample, and j is the number of channels, setting,

      $$ {X}_{j}^{t}=\sum\limits _{i=n}^{N}\left|{x}_{j}^{t}\left(i\right)\right| $$ (3)
      $$ {p}_{j}^{t}\left(i\right)=\frac{{\left|{x}_{j}^{t}\left(i\right)\right|}^{2}}{{\left({X}_{j}^{t}\right)}^{2}} $$ (4)

      Where, pjt(i) is an energy distribution of the vector at individual components, WE is defined as,

      $$ {E}_{j}^{t}\left(x\right)=-\sum\limits _{i=n}^{N}{p}_{j}^{t}\left(i\right)\mathit{\rm ln}{p}_{j}^{t}\left(i\right) $$ (5)

      Evidently, WE is closely related with amplitude, and vary with signal amplitude fluctuation. The more uniform probability pjt(i) is, the larger WE value is. Consequently, WE of a digital waveform depends on the waveform itself, we will apply the distinguished values as feature to classify the different coal samples.

    • Digital waveform provides more specific information potentially which easily available for operation[23]. Statistical variables of waveform energy were extracted from LiDAR echo signal waveform, which has been widely applied in the discriminated classes of coastal habitats and species[24], forest species[25]. In 2011, Guo et al.[26] proposed a JSKF model by combining skewness and kurtosis for band selection and achieved good results. We employ JSKF as a new feature to classify coal/rock.

      Skewness is a characteristic number that characterizes the degree of asymmetry of the probability distribution density curve relative to the average, which is expressed by the third-order standardized moments. The equation is:

      $$ S=\frac{E{[X-\mu ]}^{3}}{{\sigma }^{3}} $$ (6)

      Where, E(X) is the expectation of vector X, μ is the mean value of vector X, and σ is the standard deviation of vector X. The larger the skewness, the more asymmetrical the distribution of random variables.

      Kurtosis is a characteristic number that characterizes the peak of the probability density distribution curve at average value, which is expressed by the ratio of the fourth-order central moment of the random variable to the square of the variance. Kurtosis reflects the sharpness of the peak of the probability density distribution curve. The larger the kurtosis value, the sharper the probability density distribution curve. The equation is:

      $$K=\frac{E{[X-\mu ]}^{4}}{{\sigma }^{4}}-3 $$ (7)

      Skewness and kurtosis represent the asymmetry of random distribution, which could not only measure the difference between different bands of target characteristics, but also evaluate non-Gaussianity of data samples. JSKF uses the product of skewness and kurtosis as an indicator to measure the amount of information deviating from the normal distribution of the size. It is defined as

      $$ {F}_{\rm JSKF}=S\cdot K $$ (8)

      That is:

      $$ {F}_{\rm JSKF}=\frac{1}{{\sigma }^{7}}\left[E{\left(X-\mu \right)}^{4}-3{\sigma }^{4}\right]\cdot E{(X-\mu )}^{3} $$ (9)

      It is obvious that JSKF is related to the expectation of the waveform. The expectation is related to mean and standard deviation. The mean and standard deviation depend on the waveform itself. Therefore, the JSKF of the digital waveform also depends on the waveform itself that has no relation to external conditions. We employ the waveform of JSKF as another feature to classify different coal and rock samples in the next section.

    • In order to explore coal/rock classification based on AOTF-HSL, RF and SVM are employed as classifiers in our experiments. The main advantage of RF is only a few parameters and manual intervention needed to require high stability[27], which is usually applied for remote sensing image analysis[28]. SVM can provide high classification accuracy and it has become a very popular kernel-based classification algorithm in hyperspectral image classification[29]. Both of them have been widely applied in remote sensing researches and their efficiencies have proven in remote sensing data classification[30]. Multi-label classification is implemented by the scikit-learn Python package[31].

    • The experiments were conducted under a controlled laboratory environment. We employed HSL of 91 channels to collect different coal/rock sample data and calculated WE and JSKF based on echo signals. Then, intensity signals, WE, JSKF, and calibrated intensity signals were classified by RF and SVM classifiers to observe the performance.

    • With the number of channels in HSL increases, the instrument complexity and the amount of collected data also increase at the same time. Complex large-scale instruments and huge data volumes are hard for on-site operation and data processing. Therefore, we hope that a miniaturization HSL with less channels can also achieve accurate classification results. In order to simplify equipment complexity to save equipment resources and improve efficiency, we have improved the experiment: select a part of channels from 91 channels for classi-fication by random selection method and test accuracy. The number of channels increases from 1 in turn until the precision reaches 100%. We evaluate the consequences by comparing the number of channels required.

    • To further compare the capacity of WE, JSKF and calibrated intensity signals, we attempt to randomly extract 5 channels from 91 channels and assess the property of them respectively. At the same time, our previous research confirmed that the data of different bands have different properties[32], so we conduct a spectral segmentation test based on WE and JSKF to find the optimal band for classification.

    • The classification results of all channels are shown in the Tab.2. The accuracy of all features can reach 100% by using the full channels spectral information. The results prove that with sufficient spectral information, the classification of coal and rock can achieve ideal results.

      Table 2.  Full channels classification accuracy

      DataRFSVM
      AccuracyAccuracy
      Intensity100%100%
      WE100%100%
      JSKF100%100%
      Calibrated100%100%
    • With sufficient experiment tests, the channels of random selected classification results are shown in the Tab.3. For the RF classifier, while accuracy reaches 100%, the least number of channels needed for intensity signal is 21, which is maximum compared with others. WE needs 15 channels, the channels of JSKF and calibrated intensity signal are both 6. Compared with the intensity signal, WE, JSKF and calibrated intensity signal reduce 6, 15, 15 channels respectively. For the SVM classifier, when accuracy reaches 100%, the intensity signal needs 20 channels at least. WE needs 15 channels, JSKF needs 6 channels, and the calibrated intensity signal needs 4 channels. Compared with the intensity signal, WE, JSKF and calibrated intensity signal reduce 5, 14, 16 channels respectively. With comparison, we can draw a conclusion: the four groups of data can still achieve the desired performances by using fewer spectrum channels. Meanwhile, the number of channels needed for WE and JSKF classification is significantly less than that of intensity signal. Although WE and JSKF do not achieve the classification ability as same as calibrated intensity signal, they can provide enhancement of classification effect significantly than the intensity signal, which indicate that the HSL intensity signal calibration-free method we proposed is feasible for improving the classification performance.

      Table 3.  Minimum number of random channels

      DataRFSVM
      Number of channelsNumber of channels
      Intensity2120
      WE1515
      JSKF66
      Calibrated64
    • 5-channel random extraction classification results are shown in Tab.4. The accuracies are calculated by averaging multiple experiments.

      Table 4.  5-channel classification accuracy

      DataRFSVM
      AccuracyAccuracy
      WE88%90%
      JSKF96%100%
      Calibrated98%100%

      From the results, we can see that with 5 channels, the accuracy of calibrated intensity signal and JSKF is 98% and 96% respectively by RF, and they are both 100% by SVM. For WE, although more channels are needed to achieve accurate classification, the precision still reaches 90% by extracting 5 channels randomly. The reason is that the 91-channel HSL covers a wide spectrum that different spectral bands have different classification properties[32]. Therefore, we conduct a spectrum segmentation test to find the optimal interval of classification.

      Figure 2 show the waveform distribution diagram of WE and JSKF.

      Figure 2.  (a) Waveform distribution diagram of WE; (b) Waveform distribution diagram of JSKF

      We can see that WE and JSKF present different properties in different bands. For WE, the entropy changes slowly at 650-1000 nm and increases significantly at 1000-1100 nm in Fig.2(a). JSKF is used as an indicator to measure the amount of information that deviates from the normal distribution. The data are approximately at the lowest value between 720 nm and 820 nm in Fig.2(b). That means, the data deviation is relatively low; and when it reaches the maximum at 1000-1100 nm, the data deviation is relatively high. Thus, we test the 720-820 nm band and the 1000-1100 nm band respectively and compared them with the calibrated intensity signal. There are 20 channels in each band. The results are shown in Tab.5 and Tab.6.

      Table 5.  720-820 nm classification results

      ClassifierRFSVM
      ChannelsAccuracyChannelsAccuracy
      WE2095.02%2095.34%
      JSKF3100%4100%
      Calibrated2100%3100%

      Table 6.  1000-1100 nm classification results

      ClassifierRFSVM
      ChannelsAccuracyChannelsAccuracy
      WE11100%7100%
      JSKF19100%7100%
      Calibrated10100%6100%

      From Tab.5, we can see that in the 720-820 nm spectrum band, for the WE, the precision is only reach 95% using all of 20 channels, whereas the performance of JSKF is improved further compared to random channels selection method. The 100% accuracy only needs 3 channels with RF classifier, and needs 6 channels with random channels selection. And the 100% accuracy only needs 4 channels with SVM classifier, where efficiency increased 33% than random channels selection. In this band, the number of channels that accuracy reach 100% needed for JSKF is basically as same as the number of channels required for calibrated intensity signal.

      From Tab.6, we can see that in the 1 000 - 1 100 nm spectrum band, the performance of WE has been improved further than random channel selection method. The 100% accuracy only needs 11 channels with RF classifier, and increase 27% than channel selection randomly, and it needs 7 channels with SVM classifier, which improve 53% than channel selection randomly. For JSKF, the performance is slightly reduced. The number of channels of RF and SVM classifiers is reduced by 217% and 40% compared with random channel selection. In this band, the number of channels needs to reach 100% accuracy with WE is basically as same as the number of channels required with calibrated intensity signal.

      Based on the above experiments, we can see that WE and JSKF have different classification results in different spectral band. The reason is as following: WE characterizes the degree of dispersion for variables. Due to the different characteristic wavelength of object, the absorption characteristic of the spectrum is different. The characteristic wavelength of coal/rock is in the near-infrared spectrum[33], where the absorption characteristic is obvious, the degree of dispersion for variables is large and WE changes significantly (Fig.2(a)). Thus, the classification performance of WE is better in the near-infrared band. JSKF measures the deviation of the variable from the normal distribution. The smaller the data deviation, the better the data quality. The data collection results show that the deviation of the data is small in 720-820 nm band (Fig.2(b)). Thus, the classification performance of JSKF is better in 720-820 nm band. Meanwhile, both of them can basically reach the classification performance of calibrated intensity signal in their corresponding ideal classification band. While realizing the intensity signal calibration-free, it maintains excellent classification performances. The results show that our calibration-free classification method is feasible.

    • In this paper, we proposed a new method based on the classification of coal/rock in deep coal mines. Without any reference, we employed WE and JSKF to classify four different types of coal/rock, explored the classification performances of different features in different spectral bands and evaluated the classification performance by comparing with calibrated intensity signal. The following conclusions can be drawn from the results:

      (1) When the spectral information is sufficient, all features can achieve accurate classification.

      (2) WE and JSKF require less spectral channels than intensity signals to achieve accurate classification.

      (3) For different spectral bands, WE and JSKF have different classification performance, and they can achieve great classification performance as calibrated intensity signal in their ideal band.

      This research not only solves the problem that intensity signals cannot be calibrated in the deep coal mines, but also simplify the complexity of the equipment, laying the foundation for miniaturized LiDAR. In the future, we will do the further research for instrument miniaturization package and deep mine field application.

参考文献 (33)

目录

    /

    返回文章
    返回