本发明涉及一种煤炭粉尘污染程度的量化方法，所述方法具体包括：获取待量化的各量化单元对应的原始数据；提取各量化单元所对应各量化指标的值，生成原始指标值矩阵；建立评语集：V＝{轻微，较轻，较重，严重}；设置各评语的分值，以及各评语对应的分级区间；构建隶属于各评语的隶属度函数；确定量化指标权重集；将原始指标值矩阵转化转换为量化指标分数矩阵；依据量化指标权重集、量化指标分数矩阵和隶属度函数确定各量化单元的煤炭粉尘污染程度。还涉及了装置，包括：数据获取单元、设置单元、第一计算单元、第二计算单元、第三计算单元。本发明的方法和装置，能够全面、准确地评估煤炭粉尘污染程度。
1.一种煤炭粉尘污染程度的量化方法，其特征在于，包括：获取待量化的各量化单元对应的原始数据；从获取的原始数据提取各量化单元所对应各量化指标的值，生成原始指标值矩阵；其中，所述量化指标有：煤炭粉尘污染强度算数平均值、煤炭粉尘污染强度中位数、年内持续影响范围、年内最大影响范围、静态产尘面积比、动态产尘面积比、静态万吨产尘量、动态万吨产尘量、矿山生产建设规模和所属生命周期阶段；
设置各评语的分值，以及各评语对应的分级区间；依照评语集中的顺序分别构建隶属于各评语的递增型梯形分布隶属度函数；
依据应用场景，采用模糊层次‑熵值组合赋权法确定量化指标权重集；
采用正向标准化公式将原始指标值矩阵转化为正向标准化指标值矩阵；确定各正向标准化指标值所属的分级区间，依据分级区间、评语、评语的分值的对应关系，将正向标准化指标值矩阵转换为量化指标分数矩阵；
对任一所述量化单元，分别将该量化单元的各量化指标分数代入各隶属度函数求解，得到该量化单元各量化指标的隶属度向量，各向量组合生成该量化单元的模糊矩阵；所述量化指标权重集与该模糊矩阵相乘，得到该量化单元的隶属度向量，该隶属度向量中最大的隶属度所对应的评语即为该量化单元的煤炭粉尘污染程度；其中，对任一所述量化单元，提取所述煤炭粉尘污染强度算数平均值、中位数包括：选择在该量化单元所含年份内拍摄的、所含地域的卫星遥感影像和/或无人机航拍影像，并预处理所选择的影像；
分别计算预处理后影像所含各像元的优化的增强型煤炭粉尘指数OECDI：其中，RNIR、RSWIR1、RSWIR2分别为NIR波段、SWIR1波段和SWIR2波段的反射率；
对各像元的OECDI进行异常值剔除与归一化处理；应用像元二分模型，逐像元计算各像元的煤炭粉尘污染强度；
统计所有所述像元的煤炭粉尘污染强度，根据得到的统计直方图，确定平均值和中位数，即为提取到的所述煤炭粉尘污染强度算数平均值、中位数；
对任一所述量化单元，提取所述年内持续影响范围、年内最大影响范围包括：选择在该量化单元所含年份内拍摄的、所含地域的卫星遥感影像和/或无人机航拍影像，并预处理所选择的影像；
若在预处理后的影像中，一个季节包含不只一期影像，将该季节的全部影像进行影像合成处理，得到该季节的单景影像；若一个季节只一期影像，将该期影像作为该季节的单景影像；
分别计算各季节单景影像所含像元的煤炭粉尘污染强度；
将各季节单景影像表达为数据立方体的形式，X维和Y维组成影像所覆盖的地域，Z维为时间维，与四季对应为：z＝1,2,3,4；
对于所覆盖的地域相同的像元P(x,y,1)、P(x,y,2)、P(x,y,3)、P(x,y,4)，i(x,y)表示四个像元中煤炭粉尘污染强度大于预先设定的阈值的像元的数量；取i(x,y)＝4对应的各像元，全部所取的像元所覆盖的地域即为年内持续影响范围；取i(x,y)＝1～4对应的像元，全部所取的像元所覆盖的地域年内最大影响范围；
对任一所述量化单元，提取所述静态产尘面积比、动态产尘面积比包括：通过以下公式计算所述静态产尘面积比、动态产尘面积比：其中，SContinuous为该量化单元的年内持续影响范围；SMaximum为该量化单元的年内最大影响范围；Sstope为该量化单元所含年份的露采场未充分复垦区域的面积；Souter‑dump为该量化单元所含年份的外排土场未充分复垦区域的面积；
对任一所述量化单元，提取所述静态万吨产尘量、动态万吨产尘量包括：通过以下公式计算静态万吨产尘量、动态万吨产尘量：
其中，SContinuous为该量化单元的年内持续影响范围；SMaximum为该量化单元的年内最大影响范围；P为该量化单元所含年份的煤炭产量值，单位为万吨。
2.根据权利要求1所述的方法，其特征在于，所述采用模糊层次‑熵值组合赋权法确定量化指标权重集包括：
步骤A1～F1采用模糊层次法获取各量化指标的权重：A1、将所述各量化指标划分为两层量化指标集；其中，一层量化指标集：U＝{煤炭粉尘污染强度，煤炭粉尘污染范围，产尘面积比，万吨产尘量，矿山开发利用状况}；二层量化指标集：U1＝{算数平均值，中位数}，U2＝{年内持续影响范围，年内最大影响范围}，U3＝{静态产尘面积比，动态产尘面积比}，U4＝{静态万吨产尘量，动态万吨产尘量}，U5＝{矿山生产建设规模，所属生命周期阶段}；
B1、依据应用场景，采用0.1‑0.9九级标度法，构造目标‑一层量化指标，以及一层量化指标‑二层量化指标的相应模糊互补判断矩阵F1～F6；
C1、模糊互补判断矩阵F1＝(fij)5×5做一致化处理，构造出矩阵F1相应的模糊一致判断矩阵R1＝(rij)5×5；
D1、对模糊一致判断矩阵R1进行行归一化处理，得到指标初始权向量W1：将R1转化为互反矩阵E1：
将W1作为初始向量V1 ，利用下式进行迭代：其中，
若|||V ||∞‑||V ||∞|≤0.0001，则迭代结束，最终得到最大特征向量V ；
E1、将V 进行归一化处理，得到目标‑一层量化指标的权重向量w1；
F1、分别将F2～F5代入步骤C1～E1，得到相应一层量化指标‑二层量化指标的权重向量T T T
w2＝[w21,w22] ，w3＝[w31,w32] ，…，w6＝[w61,w62] ；则采用模糊层次所得各量化指标的权重为：w11w21＝w′1，w11w22＝w′2，w12 w31＝w′3，w12 w32＝w′4，w13 w41＝w′5，w13 w42＝w′6，w14w51＝w′7，w14w52＝w′8，w15 w61＝w′9，w15w62＝w′10；
步骤A2～D2采用熵权法获取各量化指标的权重：A2、采用负向标准化公式将所述原始指标值矩阵中定量指标的各指标值，以及定性指标的量化指标分数，分别转化为负向标准化指标值aij；将aij按照公式 进行平移处理，得到负向标准化阵 其中，m代表待量化的量化单元数量，n代表量化指标数量；
B2、计算 所占比重，得到矩阵P＝(pij)m×n，其计算公式为：C2、计算量化指标的信息熵值，第j个量化指标的信息熵值计算公式为：D2、计算第j个量化指标的权重，第j个量化指标的权重计算公式为：采用相对熵原理对模糊层次分析法所得权重和熵权法所得权重进行修正，获取耦合权重最优解w″′i，计算公式为：
确定的量化指标权重集为：[w″′1 w″′2 … w″′n]。
3.一种煤炭粉尘污染程度的量化装置，其特征在于，包括：数据获取单元、设置单元、第一计算单元、第二计算单元、第三计算单元；其中，所述数据获取单元，用于获取待量化的各量化单元对应的原始数据；从获取的原始数据提取各量化单元所对应各量化指标的值，生成原始指标值矩阵；其中，所述量化指标有：煤炭粉尘污染强度算数平均值、煤炭粉尘污染强度中位数、年内持续影响范围、年内最大影响范围、静态产尘面积比、动态产尘面积比、静态万吨产尘量、动态万吨产尘量、矿山生产建设规模和所属生命周期阶段；
设置单元，用于建立评语集：V＝{轻微，较轻，较重，严重}；设置各评语的分值，以及各评语对应的分级区间；依照评语集中的顺序分别构建隶属于各评语的递增型梯形分布隶属度函数；
第一计算单元，用于依据应用场景，采用模糊层次‑熵值组合赋权法确定量化指标权重集；
第二计算单元，用于采用正向标准化公式将原始指标值矩阵转化为正向标准化指标值矩阵；确定各正向标准化指标值所属的分级区间，依据分级区间、评语、评语的分值的对应关系，将正向标准化指标值矩阵转换为量化指标分数矩阵；
第三计算单元，用于对任一所述量化单元，分别将该量化单元的各量化指标分数代入各隶属度函数求解，得到该量化单元各量化指标的隶属度向量，各向量组合生成该量化单元的模糊矩阵；所述量化指标权重集与该模糊矩阵相乘，得到该量化单元的隶属度向量，该隶属度向量中最大的隶属度所对应的评语即为该量化单元的煤炭粉尘污染程度；其中，对任一所述量化单元，提取所述煤炭粉尘污染强度算数平均值、中位数时，所述数据获取单元具体用于：
选择在该量化单元所含年份内拍摄的、所含地域的卫星遥感影像和/或无人机航拍影像，并预处理所选择的影像；
分别计算预处理后影像所含各像元的优化的增强型煤炭粉尘指数OECDI：其中，RNIR、RSWIR1、RSWIR2分别为NIR波段、SWIR1波段和SWIR2波段的反射率；
对各像元的OECDI进行异常值剔除与归一化处理；应用像元二分模型，逐像元计算各像元的煤炭粉尘污染强度；
统计所有所述像元的煤炭粉尘污染强度，根据得到的统计直方图，确定平均值和中位数，即为提取到的所述煤炭粉尘污染强度算数平均值、中位数；
对任一所述量化单元，提取所述年内持续影响范围、年内最大影响范围时，所述数据获取单元具体用于：
选择在该量化单元所含年份内拍摄的、所含地域的卫星遥感影像和/或无人机航拍影像，并预处理所选择的影像；
若在预处理后的影像中，一个季节包含不只一期影像，将该季节的全部影像进行影像合成处理，得到该季节的单景影像；若一个季节只一期影像，将该期影像作为该季节的单景影像；
分别计算各季节单景影像所含像元的煤炭粉尘污染强度；
将各季节单景影像表达为数据立方体的形式，X维和Y维组成影像所覆盖的地域，Z维为时间维，与四季对应为：z＝1,2,3,4；
对于所覆盖的地域相同的像元P(x,y,1)、P(x,y,2)、P(x,y,3)、P(x,y,4)，i(x,y)表示四个像元中煤炭粉尘污染强度大于预先设定的阈值的像元的数量；取i(x,y)＝4对应的各像元，全部所取的像元所覆盖的地域即为年内持续影响范围；取i(x,y)＝1～4对应的像元，全部所取的像元所覆盖的地域年内最大影响范围；
对任一所述量化单元，提取所述静态产尘面积比、动态产尘面积比时，所述数据获取单元具体用于：
通过以下公式计算所述静态产尘面积比、动态产尘面积比：其中，SContinuous为该量化单元的年内持续影响范围；SMaximum为该量化单元的年内最大影响范围；Sstope为该量化单元所含年份的露采场未充分复垦区域的面积；Souter‑dump为该量化单元所含年份的外排土场未充分复垦区域的面积；
对任一所述量化单元，提取所述静态万吨产尘量、动态万吨产尘量时，所述数据获取单元具体用于：
通过以下公式计算静态万吨产尘量、动态万吨产尘量：其中，SContinuous为该量化单元的年内持续影响范围；SMaximum为该量化单元的年内最大影响范围；P为该量化单元所含年份的煤炭产量值，单位为万吨。
4.根据权利要求3所述的装置，其特征在于，所述第一计算单元具体用于：步骤A1～F1采用模糊层次法获取各量化指标的权重：A1、将所述各量化指标划分为两层量化指标集；其中，一层量化指标集：U＝{煤炭粉尘污染强度，煤炭粉尘污染范围，产尘面积比，万吨产尘量，矿山开发利用状况}；二层量化指标集：U1＝{算数平均值，中位数}，U2＝{年内持续影响范围，年内最大影响范围}，U3＝{静态产尘面积比，动态产尘面积比}，U4＝{静态万吨产尘量，动态万吨产尘量}，U5＝{矿山生产建设规模，所属生命周期阶段}；
B1、依据应用场景，采用0.1‑0.9九级标度法，构造目标‑一层量化指标，以及一层量化指标‑二层量化指标的相应模糊互补判断矩阵F1～F6；
C1、模糊互补判断矩阵F1＝(fij)5×5做一致化处理，构造出矩阵F1相应的模糊一致判断矩阵R1＝(rij)5×5；
D1、对模糊一致判断矩阵R1进行行归一化处理，得到指标初始权向量W1：将R1转化为互反矩阵E1：
将W1作为初始向量V1 ，利用下式进行迭代：其中，
若|||V ||∞‑||V ||∞|≤0.0001，则迭代结束，最终得到最大特征向量V ；
E1、将V 进行归一化处理，得到目标‑一层量化指标的权重向量w1；
F1、分别将F2～F5代入步骤C1～E1，得到相应一层量化指标‑二层量化指标的权重向量T T T
w2＝[w21,w22] ，w3＝[w31,w32] ，…，w6＝[w61,w62] ；则采用模糊层次所得各量化指标的权重为：w11w21＝w′1，w11w22＝w′2，w12 w31＝w′3，w12 w32＝w′4，w13 w41＝w′5，w13 w42＝w′6，w14w51＝w′7，w14w52＝w′8，w15 w61＝w′9，w15 w62＝w′10；
步骤A2～D2采用熵权法获取各量化指标的权重：A2、采用负向标准化公式将所述原始指标值矩阵中定量指标的各指标值，以及定性指标的量化指标分数，分别转化为负向标准化指标值aij；将aij按照公式 进行平移处理，得到负向标准化阵 其中，m代表待量化的量化单元数量，n代表量化指标数量；
B2、计算 所占比重，得到矩阵P＝(pij)m×n，其计算公式为：C2、计算量化指标的信息熵值，第j个量化指标的信息熵值计算公式为：D2、计算第j个量化指标的权重，第j个量化指标的权重计算公式为：采用相对熵原理对模糊层次分析法所得权重和熵权法所得权重进行修正，获取耦合权重最优解w″′i，计算公式为：
确定的量化指标权重集为：[w″′1 w″′2 … w″′n]。
一种煤炭粉尘污染程度的量化方法及装置
技术领域
[0001] 本发明涉及环境监测技术领域，尤其是一种煤炭粉尘污染程度的量化方法及装置。
背景技术
[0002] 近年，我国工业生产水平不断提高，对煤炭资源开发和利用的力度也在不断加大，露天煤矿作为我国煤炭工业的重要组成部分，具有生产规模大、资源回收率高、安全条件高等优势，然而，露天开采势必会对原有的生态环境带来严重破坏，其中，对大气环境的破坏形式之一就是煤炭粉尘污染，该污染容易引起尘肺疾病、加速机器磨损、降低生产场所能见度等问题。而不同建设时期、不同矿山条件下的煤炭粉尘污染问题差异较大，目前，煤炭粉尘污染评估的方法有如下几种。
[0003] 使用采样器、测尘仪或粉尘监测站等设备对粉尘浓度进行测定。其中，主要测定对象包括总粉尘浓度、个体呼吸性粉尘浓度及粉尘中游离SiO2含量。粉尘浓度测定结果能够及时、准确地反映出相应测尘点的粉尘污染状况，但仅能测定静态、有限点位数据，难以全面、准确地评估矿山整体的煤炭粉尘污染情况。
[0004] 公告号为CN110414153B的发明专利提出了一种确定露天煤矿粉尘扩散范围的方法，首先利用经验公式获取露天煤矿生产过程中各个作业环节的产尘强度，其次依据大气
环境和尘源特点，确定煤炭粉尘的扩散模型和扩散参数，最后利用Matlab编制程序，模拟得到露天煤矿粉尘的扩散范围。即采用数值模拟方法，模拟研究不同前提条件下相关因素对
粉尘扩散的影响强度和影响机理，以预测露天煤矿开采过程中粉尘的扩散范围。但基于数
值模拟的不确定性无法评估，且由于粉尘污染存在多变性，难以建立出精确的粉尘传输扩
散模型，因而模型模拟结果并不能用作煤炭粉尘污染程度的实际量化指标。
[0005] 叶面降尘遥感监测方法，一方面基于实测叶面降尘量和航空航天遥感数据，建立叶面降尘量的反演与预测模型，另一方面通过构建新的指数(叶面降尘差值指数、比值指
数、归一化指数，尘源光谱指数等)或分析已有植被指数对叶面降尘量的敏感程度来间接评估叶面降尘的范围和强度，但叶面降尘量反演与预测模型的建立依赖于野外实测工作的进
行，费时费力，效率较低，且现有预测模型和指数大多针对自然降尘，应用到植被叶面稀疏的露天矿及其周边区域并不能准确反应煤炭粉尘污染状况。
[0006] 综上所述，目前的煤炭粉尘污染评估的方法并不能全面、准确地评估矿山整体的煤炭粉尘污染程度。
发明内容
[0007] 有鉴于此，本发明的主要目的在于提供一种煤炭粉尘污染程度的量化方法及装置，能够全面、准确地评估煤炭粉尘污染程度。
[0008] 为达到上述目的，一方面，本申请提供了一种煤炭粉尘污染程度的量化方法，包括：获取待量化的各量化单元对应的原始数据；从获取的原始数据提取各量化单元所对应
各量化指标的值，生成原始指标值矩阵；其中，所述量化指标有：煤炭粉尘污染强度算数平均值、煤炭粉尘污染强度中位数、年内持续影响范围、年内最大影响范围、静态产尘面积比、动态产尘面积比、静态万吨产尘量、动态万吨产尘量、矿山生产建设规模和所属生命周期阶段；
[0009] 建立评语集：V＝{轻微，较轻，较重，严重}；
[0010] 设置各评语的分值，以及各评语对应的分级区间；依照评语集中的顺序分别构建隶属于各评语的递增型梯形分布隶属度函数；
[0011] 依据应用场景，采用模糊层次‑熵值组合赋权法确定量化指标权重集；
[0012] 采用正向标准化公式将原始指标值矩阵转化为正向标准化指标值矩阵；确定各正向标准化指标值所属的分级区间，依据分级区间、评语、评语的分值的对应关系，将正向标准化指标值矩阵转换为量化指标分数矩阵；
[0013] 对任一所述量化单元，分别将该量化单元的各量化指标分数代入各隶属度函数求解，得到该量化单元各量化指标的隶属度向量，各向量组合生成该量化单元的模糊矩阵；所述量化指标权重集与该模糊矩阵相乘，得到该量化单元的隶属度向量，该隶属度向量中最
大的隶属度所对应的评语即为该量化单元的煤炭粉尘污染程度。
[0014] 在一个可能的实现中，对任一所述量化单元，提取所述煤炭粉尘污染强度算数平均值、中位数包括：
[0015] 选择在该量化单元所含年份内拍摄的、所含地域的卫星遥感影像和/或无人机航拍影像，并预处理所选择的影像；
[0016] 分别计算预处理后影像所含各像元的优化的增强型煤炭粉尘指数OECDI：
[0017]
[0018] 其中，RNIR、RSWIR1、RSWIR2分别为NIR波段、SWIR1波段和SWIR2波段的反射率；
[0019] 对各像元的OECDI进行异常值剔除与归一化处理；应用像元二分模型，逐像元计算各像元的煤炭粉尘污染强度；
[0020] 统计所有所述像元的煤炭粉尘污染强度，根据得到的统计直方图，确定平均值和中位数，即为提取到的所述煤炭粉尘污染强度算数平均值、中位数。
[0021] 在另一个可能的实现中，对任一所述量化单元，提取所述年内持续影响范围、年内最大影响范围包括：
[0022] 选择在该量化单元所含年份内拍摄的、所含地域的卫星遥感影像和/或无人机航拍影像，并预处理所选择的影像；
[0023] 若在预处理后的影像中，一个季节包含不只一期影像，将该季节的全部影像进行影像合成处理，得到该季节的单景影像；若一个季节只一期影像，将该期影像作为该季节的单景影像；
[0024] 分别计算各季节单景影像所含像元的煤炭粉尘污染强度；
[0025] 将各季节单景影像表达为数据立方体的形式，X维和Y维组成影像所覆盖的地域，Z维为时间维，与四季对应为：z＝1,2,3,4；
[0026] 对于所覆盖的地域相同的像元P(x,y,1)、P(x,y,2)、P(x,y,3)、P(x,y,4)，i(x,y)表示四个像元中煤炭粉尘污染强度大于预先设定的阈值的像元的数量；取i(x,y)＝4对应的各像元，全部所取的像元所覆盖的地域即为年内持续影响范围；取i(x,y)＝1～4对应的像元，全部所取的像元所覆盖的地域年内最大影响范围。
[0027] 在另一个可能的实现中，对任一所述量化单元，提取所述静态产尘面积比、动态产尘面积比包括：通过以下公式计算所述静态产尘面积比、动态产尘面积比：
[0028]
[0029]
[0030] 其中，SContinuous为该量化单元的年内持续影响范围；SMaximum为该量化单元的年内最大影响范围；Sstope为该量化单元所含年份的露采场未充分复垦区域的面积；Souter‑dump为该量化单元所含年份的外排土场未充分复垦区域的面积。
[0031] 在另一个可能的实现中，对任一所述量化单元，提取所述静态万吨产尘量、动态万吨产尘量包括：通过以下公式计算静态万吨产尘量、动态万吨产尘量：
[0032]
[0033]
[0034] 其中，SContinuous为该量化单元的年内持续影响范围；SMaximum为该量化单元的年内最大影响范围；P为该量化单元所含年份的煤炭产量值，单位为万吨。
[0035] 在另一个可能的实现中，所述采用模糊层次‑熵值组合赋权法确定量化指标权重集包括：
[0036] 步骤A1～F1采用模糊层次法获取各量化指标的权重：
[0037] A1、将所述各量化指标划分为两层量化指标集；其中，一层量化指标集：U＝{煤炭粉尘污染强度，煤炭粉尘污染范围，产尘面积比，万吨产尘量，矿山开发利用状况}；二层量化指标集：U1＝{算数平均值，中位数}，U2＝{年内持续影响范围，年内最大影响范围}，U3＝{静态产尘面积比，动态产尘面积比}，U4＝{静态万吨产尘量，动态万吨产尘量}，U5＝{矿山生产建设规模，所属生命周期阶段}；
[0038] B1、依据应用场景，采用0.1‑0.9九级标度法，构造目标‑一层量化指标，以及一层量化指标‑二层量化指标的相应模糊互补判断矩阵F1～F6；
[0039] C1、模糊互补判断矩阵F1＝(fij)5×5做一致化处理，构造出矩阵F1相应的模糊一致判断矩阵R1＝(rij)5×5；
[0040] D1、对模糊一致判断矩阵R1进行行归一化处理，得到指标初始权向量W1：
[0041]
[0042] 将R1转化为互反矩阵E1：
[0043]
[0044] 将W1作为初始向量V1(0)，利用下式进行迭代：
[0045]
[0046] 其中，
[0047] 若|||V(k+1)||∞‑||V(k)||∞|≤0.0001，则迭代结束，最终得到最大特征向量V(k+1)；
[0048] E1、将V(k+1)进行归一化处理，得到目标‑一层量化指标的权重向量w1；
[0049]
[0050] F1、分别将F2～F5代入步骤C1～E1，得到相应一层量化指标‑二层量化指标的权重T T T向量w2＝[w21,w22] ，w3＝[w31,w32] ，…，w6＝[w61,w62] ；则采用模糊层次所得各量化指标的权重为：w11w21＝w′1，w11w22＝w′2，w12 w31＝w′3，w12 w32＝w′4，w13 w41＝w′5，w13 w42＝w′6，w14w51＝w′7，w14w52＝w′8，w15 w61＝w′9，w15 w62＝w′10；
[0051] 步骤A2～D2采用熵权法获取各量化指标的权重：
[0052] A2、采用负向标准化公式将所述原始指标值矩阵中定量指标的各指标值，以及定性指标的量化指标分数，分别转化为负向标准化指标值aij；将aij按照公式 进
行平移处理，得到负向标准化阵 其中，m代表待量化的量化单元数量，n代表量
化指标数量；
[0053] B2、计算 所占比重，得到矩阵P＝(pij)m×n，其计算公式为：
[0054]
[0055] C2、计算量化指标的信息熵值，第j个量化指标的信息熵值计算公式为：
[0056]
[0057] D2、计算第j个量化指标的权重，第j个量化指标的权重计算公式为：
[0058]
[0059] 采用相对熵原理对模糊层次分析法所得权重和熵权法所得权重进行修正，获取耦合权重最优解w″′i，计算公式为：
[0060]
[0061] 确定的量化指标权重集为：[w″′1 w″′2 … w″′n]。
[0062] 另一方面，本申请还提供了一种煤炭粉尘污染程度的量化装置，包括：数据获取单元、设置单元、第一计算单元、第二计算单元、第三计算单元；其中，
[0063] 所述数据获取单元，用于获取待量化的各量化单元对应的原始数据；从获取的原始数据提取各量化单元所对应各量化指标的值，生成原始指标值矩阵；其中，所述量化指标有：煤炭粉尘污染强度算数平均值、煤炭粉尘污染强度中位数、年内持续影响范围、年内最大影响范围、静态产尘面积比、动态产尘面积比、静态万吨产尘量、动态万吨产尘量、矿山生产建设规模和所属生命周期阶段；
[0064] 设置单元，用于建立评语集：V＝{轻微，较轻，较重，严重}；设置各评语的分值，以及各评语对应的分级区间；依照评语集中的顺序分别构建隶属于各评语的递增型梯形分布隶属度函数；
[0065] 第一计算单元，用于依据应用场景，采用模糊层次‑熵值组合赋权法确定量化指标权重集；
[0066] 第二计算单元，用于采用正向标准化公式将原始指标值矩阵转化为正向标准化指标值矩阵；确定各正向标准化指标值所属的分级区间，依据分级区间、评语、评语的分值的对应关系，将正向标准化指标值矩阵转换为量化指标分数矩阵；
[0067] 第三计算单元，用于对任一所述量化单元，分别将该量化单元的各量化指标分数代入各隶属度函数求解，得到该量化单元各量化指标的隶属度向量，各向量组合生成该量
化单元的模糊矩阵；所述量化指标权重集与该模糊矩阵相乘，得到该量化单元的隶属度向
量，该隶属度向量中最大的隶属度所对应的评语即为该量化单元的煤炭粉尘污染程度。
[0068] 在一个可能的实现中，对任一所述量化单元，提取所述煤炭粉尘污染强度算数平均值、中位数时，所述数据获取单元具体用于：
[0069] 选择在该量化单元所含年份内拍摄的、所含地域的卫星遥感影像和/或无人机航拍影像，并预处理所选择的影像；
[0070] 分别计算预处理后影像所含各像元的优化的增强型煤炭粉尘指数OECDI：
[0071]
[0072] 其中，RNIR、RSWIR1、RSWIR2分别为NIR波段、SWIR1波段和SWIR2波段的反射率；
[0073] 对各像元的OECDI进行异常值剔除与归一化处理；应用像元二分模型，逐像元计算各像元的煤炭粉尘污染强度；
[0074] 统计所有所述像元的煤炭粉尘污染强度，根据得到的统计直方图，确定平均值和中位数，即为提取到的所述煤炭粉尘污染强度算数平均值、中位数。
[0075] 在另一个可能的实现中，对任一所述量化单元，提取所述年内持续影响范围、年内最大影响范围时，所述数据获取单元具体用于：
[0076] 选择在该量化单元所含年份内拍摄的、所含地域的卫星遥感影像和/或无人机航拍影像，并预处理所选择的影像；
[0077] 若在预处理后的影像中，一个季节包含不只一期影像，将该季节的全部影像进行影像合成处理，得到该季节的单景影像；若一个季节只一期影像，将该期影像作为该季节的单景影像；
[0078] 分别计算各季节单景影像所含像元的煤炭粉尘污染强度；
[0079] 将各季节单景影像表达为数据立方体的形式，X维和Y维组成影像所覆盖的地域，Z维为时间维，与四季对应为：z＝1,2,3,4；
[0080] 对于所覆盖的地域相同的像元P(x,y,1)、P(x,y,2)、P(x,y,3)、P(x,y,4)，i(x,y)表示四个像元中煤炭粉尘污染强度大于预先设定的阈值的像元的数量；取i(x,y)＝4对应的各像元，全部所取的像元所覆盖的地域即为年内持续影响范围；取i(x,y)＝1～4对应的像元，全部所取的像元所覆盖的地域年内最大影响范围。
[0081] 在另一个可能的实现中，对任一所述量化单元，提取所述静态产尘面积比、动态产尘面积比时，所述数据获取单元具体用于：
[0082] 通过以下公式计算所述静态产尘面积比、动态产尘面积比：
[0083]
[0084]
[0085] 其中，SContinuous为该量化单元的年内持续影响范围；SMaximum为该量化单元的年内最大影响范围；Sstope为该量化单元所含年份的露采场未充分复垦区域的面积；Souter‑dump为该量化单元所含年份的外排土场未充分复垦区域的面积。
[0086] 在另一个可能的实现中，对任一所述量化单元，提取所述静态万吨产尘量、动态万吨产尘量时，所述数据获取单元具体用于：
[0087] 通过以下公式计算静态万吨产尘量、动态万吨产尘量：
[0088]
[0089]
[0090] 其中，SContinuous为该量化单元的年内持续影响范围；SMaximum为该量化单元的年内最大影响范围；P为该量化单元所含年份的煤炭产量值，单位为万吨。
[0091] 在另一个可能的实现中，所述第一计算单元具体用于：
[0092] 步骤A1～F1采用模糊层次法获取各量化指标的权重：
[0093] A1、将所述各量化指标划分为两层量化指标集；其中，一层量化指标集：U＝{煤炭粉尘污染强度，煤炭粉尘污染范围，产尘面积比，万吨产尘量，矿山开发利用状况}；二层量化指标集：U1＝{算数平均值，中位数}，U2＝{年内持续影响范围，年内最大影响范围}，U3＝{静态产尘面积比，动态产尘面积比}，U4＝{静态万吨产尘量，动态万吨产尘量}，U5＝{矿山生产建设规模，所属生命周期阶段}；
[0094] B1、依据应用场景，采用0.1‑0.9九级标度法，构造目标‑一层量化指标，以及一层量化指标‑二层量化指标的相应模糊互补判断矩阵F1～F6；
[0095] C1、模糊互补判断矩阵F1＝(fij)5×5做一致化处理，构造出矩阵F1相应的模糊一致判断矩阵R1＝(rij)5×5；
[0096] D1、对模糊一致判断矩阵R1进行行归一化处理，得到指标初始权向量W1：
[0097]
[0098] 将R1转化为互反矩阵E1：
[0099]
[0100] 将W1作为初始向量V1(0)，利用下式进行迭代：
[0101]
[0102] 其中，
[0103] 若|||V(k+1)||∞‑||V(k)||∞|≤0.0001，则迭代结束，最终得到最大特征向量V(k+1)；
[0104] E1、将V(k+1)进行归一化处理，得到目标‑一层量化指标的权重向量w1；
[0105]
[0106] F1、分别将F2～F5代入步骤C1～E1，得到相应一层量化指标‑二层量化指标的权重T T T向量w2＝[w21,w22] ，w3＝[w31,w32] ，…，w6＝[w61,w62] ；则采用模糊层次所得各量化指标的权重为：w11 w21＝w′1，w11 w22＝w′2，w12 w31＝w′3，w12 w32＝w′4，w13 w41＝w′5，w13 w42＝w′6，w14 w51＝w′7，w14 w52＝w′8，w15 w61＝w′9，w15 w62＝w′10；
[0107] 步骤A2～D2采用熵权法获取各量化指标的权重：
[0108] A2、采用负向标准化公式将所述原始指标值矩阵中定量指标的各指标值，以及定性指标的量化指标分数，分别转化为负向标准化指标值aij；将aij按照公式 进
行平移处理，得到负向标准化阵 其中，m代表待量化的量化单元数量，n代表量化
指标数量；
[0109] B2、计算 所占比重，得到矩阵P＝(pij)m×n，其计算公式为：
[0110]
[0111] C2、计算量化指标的信息熵值，第j个量化指标的信息熵值计算公式为：
[0112]
[0113] D2、计算第j个量化指标的权重，第j个量化指标的权重计算公式为：
[0114]
[0115] 采用相对熵原理对模糊层次分析法所得权重和熵权法所得权重进行修正，获取耦合权重最优解w″′i，计算公式为：
[0116]
[0117] 确定的量化指标权重集为：[w″′1 w″′2 … w″′n]。
[0118] 基于上述，本发明提供的一种煤炭粉尘污染程度的量化方法及装置，具有以下优点和特点：
[0119] 1、能够用于绿色矿山评估、环境损害鉴定、矿山长时序监测等，为其实现的基础和依据。
[0120] 2、更加清晰地界定了煤炭粉尘污染强度和范围，有助于为环境损害鉴定相关工作的针对性、规范性和科学性提供参考。
[0121] 3、不仅能够实现单个矿山长时序大面域尺度条件下，煤炭粉尘污染程度的多指标综合量化评估，同时能够实现多个矿山煤炭粉尘污染程度的差异评估，确保矿山运营的可
持续发展，为绿色矿山的评比和遴选提供量化数据支撑。
附图说明
[0122] 图1为本发明实施例一种煤炭粉尘污染程度的量化方法的流程示意图。
具体实施方式
[0123] 具体的，本发明实施例一种煤炭粉尘污染程度的量化方法的流程如图1所示，包括步骤101～106。
[0124] 步骤101：获取待量化的各量化单元对应的原始数据；从获取的原始数据提取各量化单元所对应各量化指标的值，生成原始指标值矩阵。
[0125] 其中，所述量化指标有：煤炭粉尘污染强度算数平均值、煤炭粉尘污染强度中位数、年内持续影响范围、年内最大影响范围、静态产尘面积比、动态产尘面积比、静态万吨产尘量、动态万吨产尘量、矿山生产建设规模和所属生命周期阶段。
[0126] 步骤102：建立评语集：V＝{轻微，较轻，较重，严重}。
[0127] 步骤103：设置各评语的分值，以及各评语对应的分级区间；依照评语集中的顺序分别构建隶属于各评语的递增型梯形分布隶属度函数。
[0128] 步骤104：依据应用场景，采用模糊层次‑熵值组合赋权法确定量化指标权重集。
[0129] 步骤105：采用正向标准化公式将原始指标值矩阵转化为正向标准化指标值矩阵；确定各正向标准化指标值所属的分级区间，依据分级区间、评语、评语的分值的对应关系，将正向标准化指标值矩阵转换为量化指标分数矩阵。
[0130] 步骤106：对任一所述量化单元，分别将该量化单元的各量化指标分数代入各隶属度函数求解，得到该量化单元各量化指标的隶属度向量，各向量组合生成该量化单元的模
糊矩阵；所述量化指标权重集与该模糊矩阵相乘，得到该量化单元的隶属度向量，该隶属度向量中最大的隶属度所对应的评语即为该量化单元的煤炭粉尘污染程度。
[0131] 这里，所述量化单元包含年份和地域，即对该地域在该年份的煤炭粉尘污染程度进行量化。
[0132] 步骤103中，所述各评语的分值为：轻微1分、较轻2分、较重3分、严重4分。
[0133] 对于定量量化指标煤炭粉尘污染强度算数平均值、煤炭粉尘污染强度中位数、年内持续影响范围、年内最大影响范围、静态产尘面积比、动态产尘面积比、静态万吨产尘量、动态万吨产尘量，所述评语轻微、较轻、较重、严重与分级区间[0，0.25]、(0.25，0.5]、(0.5，
0.75]、(0.75，1]顺序对应。
[0134] 对于定性量化指标，正向标准化指标值即为分值区间；对于定性量化指标矿山生产建设规模，正向标准化指标值—、小型、中型、大型与所述评语轻微、较轻、较重、严重顺序对应；对于定性量化指标所属生命周期阶段，正向标准化指标值衰退期、稳产期、达产期+丰产期、投产期与所述评语轻微、较轻、较重、严重顺序对应。
[0135] 所述隶属度函数分别为：
[0136] 隶属于轻微的隶属度函数：
[0137]
[0138] 隶属于较轻的隶属度函数：
[0139]
[0140] 隶属于较重隶属度函数：
[0141]
[0142] 隶属于严重的隶属度函数：
[0143]
[0144] 其中，bij为量化单元i的量化指标j的量化指标分数。
[0145] 步骤105中，所述正向标准化公式为：
[0146]
[0147] 其中，x为原始指标值，xnew为正向标准化指标值，xmin为全部量化单元中相应量化指标的各原始指标值中的最小值，xmax为全部量化单元中相应量化指标的各原始指标值中的最大值。
[0148] 对于定性量化指标的原始指标值，不需要进行正向标准化运算，即，将定性量化指标的原始指标值作为正向标准化指标值。
[0149] 第i个量化单元的所述模糊矩阵表示为：
[0150]
[0151] 其中，n为量化指标的数量。
[0152] 所述量化指标权重集与第i个量化单元的模糊矩阵相乘，得到第i个量化单元的隶属度向量Gi表示为：
[0153]
[0154] 其中，n为量化指标的数量。
[0155] 在一个可能的实现中，对任一所述量化单元，提取所述煤炭粉尘污染强度算数平均值、中位数包括：
[0156] 选择在该量化单元所含年份内拍摄的、所含地域的卫星遥感影像和/或无人机航拍影像，并预处理所选择的影像；
[0157] 分别计算预处理后影像所含各像元的优化的增强型煤炭粉尘指数OECDI：
[0158]
[0159] 其中，RNIR、RSWIR1、RSWIR2分别为NIR波段、SWIR1波段和SWIR2波段的反射率；
[0160] 对各像元的OECDI进行异常值剔除与归一化处理；应用像元二分模型，逐像元计算各像元的煤炭粉尘污染强度；
[0161] 统计所有所述像元的煤炭粉尘污染强度，根据得到的统计直方图，确定平均值和中位数，即为提取到的所述煤炭粉尘污染强度算数平均值、中位数。
[0162] 这里，所述各量化单元对应的原始数据包括：卫星遥感影像和/或无人机航拍影像。
[0163] 卫星遥感影像可选择Landsat TM/ETM+/OLI、Sentinel‑2MSI等；可在以下下载平台下载卫星遥感影像 ：G oogle Earth Eng ine平台 (https ://
code.earthengine.google.com/)、美国地质调查局(https://www.usgs.gov/)、欧空局
(https://scihub.copernicus.eu/)、地理空间数据云(http://www.gscloud.cn/)等。
[0164] 所述选择影像为：依据实际需要，选择各量化单元所含年份内各季节所有可用的卫星遥感影像和/或无人机航拍影像；或者，各量化单元所含地域的风向风级指标，在各季节所有可用的卫星遥感影像和/或无人机航拍影像中，选择单期或少量的无风天气或具有
季节代表性风向的不同风级天气下的影像。
[0165] 所述预处理具体为：预处理卫星遥感影像和/或无人机航拍影像；其中，
[0166] 预处理卫星遥感影像，包括：时间和云量筛选、辐射校正、几何校正、波段合成、云及阴影掩膜和影像裁剪等。这里，可选择预处理平台Google Earth Engine平台、ENVI软件、Sen2Cor软件或SNAP软件等，实现卫星遥感影像预处理。
[0167] 预处理无人机航拍影像，包括：影像配准、辐射校正、几何校正、影像拼接和裁剪等。
[0168] 这里，对增强型煤炭粉尘指数(Enhanced Coal Dust Index，ECDI)进行改进，得到优化的增强型煤炭粉尘指数(Optimized Enhanced Coal Dust Index,OECDI)，能够较好地区分煤炭粉尘与周边其他地物，且解决了ECDI在裸煤地区容易出现负异常值，以致影响煤
炭粉尘污染状况评估准确性和可靠性的问题。
[0169] 所述对各像元的OECDI进行异常值剔除采用拉依达3σ法则实现。
[0170] 所述对各像元的OECDI进行归一化处理的计算公式为：
[0171]
[0172] 其中，OECDImin和OECDImax分别为各像元OECDI中的最小值和最大值。
[0173] 所述煤炭粉尘污染强度的计算公式为：
[0174]
[0175] 其中， 和 分别为全粉尘覆盖和无粉尘覆盖像元的*
分别取累积频率为2％和98％的OECDI值作为 和
[0176] 在另一个可能的实现中，采用时间加权平均总粉尘浓度来印证所述煤炭粉尘污染强度的可用性，具体为：
[0177] 步骤a：测定时间加权平均总粉尘浓度；包括：布设采样点：以待测区域各露天矿采场为中心，从采场中心向外沿西、西北、北、东北、东、东南、南、西南8个方向设置辐射状的调查样线，采场内采样间隔为100m，距采场边界0～2km范围内采样间隔为200m，2～4km范围内采样间隔为400m；沿待评价单元煤炭运输道路布设采样点，采样间隔为100m～200m。
[0178] 步骤b：根据《GBZ/T 192.1工业场所空气中粉尘测定第1部分：总粉尘浓度》标准，采用一个工作班连续滤膜采样称重法对各采样点的时间加权平均总粉尘浓度进行测定，各采样点分别测定3～5次，取算数平均值结果作为该采样点时间加权平均总粉尘浓度；
[0179] 步骤c：基于采样点实测的时间加权平均总粉尘浓度，利用克里金插值法预测待测区域时间加权平均总粉尘浓度的空间分布图；对比分析同一时间煤炭粉尘污染强度与时间
加权平均总粉尘浓度的空间分布图；若在空间分布上保持近似一致，则表明获取的煤炭粉
尘污染强度指标具有可用性，可继续进行下一步的操作；若在空间分布上具有很大差异，则表明获取的煤炭粉尘污染强度不具有可用性，需要重新计算煤炭粉尘污染强度。
[0180] 在另一个可能的实现中，所述时间加权平均总粉尘浓度也可以作为量化指标之一。
[0181] 在另一个可能的实现中，对任一所述量化单元，提取所述年内持续影响范围、年内最大影响范围包括：
[0182] 选择在该量化单元所含年份内拍摄的、所含地域的卫星遥感影像和/或无人机航拍影像，并预处理所选择的影像；
[0183] 若在预处理后的影像中，一个季节包含不只一期影像，将该季节的全部影像进行影像合成处理，得到该季节的单景影像；若一个季节只一期影像，将该期影像作为该季节的单景影像；
[0184] 分别计算各季节单景影像所含像元的煤炭粉尘污染强度；
[0185] 将各季节单景影像表达为数据立方体的形式，X维和Y维组成影像所覆盖的地域，Z维为时间维，与四季对应为：z＝1,2,3,4；
[0186] 对于所覆盖的地域相同的像元P(x,y,1)、P(x,y,2)、P(x,y,3)、P(x,y,4)，i(x,y)表示四个像元中煤炭粉尘污染强度大于预先设定的阈值的像元的数量；取i(x,y)＝4对应的各像元，全部所取的像元所覆盖的地域即为年内持续影响范围；取i(x,y)＝1～4对应的像元，全部所取的像元所覆盖的地域年内最大影响范围。
[0187] 这里，所述阈值为μ+σ，μ和σ分别为各季节单景影像所含的全部像元的煤炭粉尘污染强度的算数平均值和标准差。
[0188] 所述影像合成处理采用均值合成算法，可通过Google Earth Engine平台实现。
[0189] 在另一个可能的实现中，对任一所述量化单元，提取所述静态产尘面积比、动态产尘面积比包括：通过以下公式计算所述静态产尘面积比、动态产尘面积比：
[0190]
[0191]
[0192] 其中，SContinuous为该量化单元的年内持续影响范围；SMaximum为该量化单元的年内最大影响范围；Sstope为该量化单元所含年份的露采场未充分复垦区域的面积；Souter‑dump为该量化单元所含年份的外排土场未充分复垦区域的面积。
[0193] 在另一个可能的实现中，对任一所述量化单元，提取所述静态万吨产尘量、动态万吨产尘量包括：通过以下公式计算静态万吨产尘量、动态万吨产尘量：
[0194]
[0195]
[0196] 其中，SContinuous为该量化单元的年内持续影响范围；SMaximum为该量化单元的年内最大影响范围；P为该量化单元所含年份的煤炭产量值，单位为万吨。
[0197] 在另一个可能的实现中，所述采用模糊层次‑熵值组合赋权法确定量化指标权重集包括：
[0198] 步骤A1～F1采用模糊层次法获取各量化指标的权重：
[0199] A1、将所述各量化指标划分为两层量化指标集；其中，一层量化指标集：U＝{煤炭粉尘污染强度，煤炭粉尘污染范围，产尘面积比，万吨产尘量，矿山开发利用状况}；二层量化指标集：U1＝{算数平均值，中位数}，U2＝{年内持续影响范围，年内最大影响范围}，U3＝{静态产尘面积比，动态产尘面积比}，U4＝{静态万吨产尘量，动态万吨产尘量}，U5＝{矿山生产建设规模，所属生命周期阶段}；
[0200] B1、依据应用场景，采用0.1‑0.9九级标度法，构造目标‑一层量化指标，以及一层量化指标‑二层量化指标的相应模糊互补判断矩阵F1～F6；
[0201] C1、模糊互补判断矩阵F1＝(fij)5×5做一致化处理，构造出矩阵F1相应的模糊一致判断矩阵R1＝(rij)5×5；
[0202] D1、对模糊一致判断矩阵R1进行行归一化处理，得到指标初始权向量W1：
[0203]
[0204] 将R1转化为互反矩阵E1：
[0205]
[0206] 将W1作为初始向量V1(0)，利用下式进行迭代：
[0207]
[0208] 其中，
[0209] 若|||V(k+1)||∞‑||V(k)||∞|≤0.0001，则迭代结束，最终得到最大特征向量V(k+1)；
[0210] E1、将V(k+1)进行归一化处理，得到目标‑一层量化指标的权重向量w1；
[0211]
[0212] F1、分别将F2～F5代入步骤C1～E1，得到相应一层量化指标‑二层量化指标的权重T T T向量w2＝[w21,w22] ，w3＝[w31,w32] ，…，w6＝[w61,w62] ；则采用模糊层次所得各量化指标的权重为：w11 w21＝w′1，w11 w22＝w′2，w12 w31＝w′3，w12 w32＝w′4，w13 w41＝w′5，w13 w42＝w′6，w14 w51＝w′7，w14 w52＝w′8，w15 w61＝w′9，w15w62＝w′10；
[0213] 步骤A2～D2采用熵权法获取各量化指标的权重：
[0214] A2、采用负向标准化公式将所述原始指标值矩阵中定量指标的各指标值，以及定性指标的量化指标分数，分别转化为负向标准化指标值aij；将aij按照公式 进
行平移处理，得到负向标准化阵 其中，m代表待量化的量化单元数量，n代表量
化指标数量；
[0215] B2、计算 所占比重，得到矩阵P＝(pij)m×n，其计算公式为：
[0216]
[0217] C2、计算量化指标的信息熵值，第j个量化指标的信息熵值计算公式为：
[0218]
[0219] D2、计算第j个量化指标的权重，第j个量化指标的权重计算公式为：
[0220]
[0221] 采用相对熵原理对模糊层次分析法所得权重和熵权法所得权重进行修正，获取耦合权重最优解w″′i，计算公式为：
[0222]
[0223] 确定的量化指标权重集为：[w″′1 w″′2 … w″′n]。
[0224] 这里，0.1‑0.9九级标度如下表：
[0225] 表1
[0226]
[0227] 所述应用场景可以是绿色矿山评估、环境损害鉴定、矿山时序监测等。其中，应用场景为绿色矿山评估时，主要量化指标为矿山开发利用状况、煤炭粉尘污染强度，其次是产尘面积比、万吨产尘量和煤炭粉尘污染范围；应用场景为环境损害鉴定时，煤炭粉尘污染范围、煤炭粉尘污染强度作为首要量化指标，其余指标为次要量化指标；应用场景为矿山时序监测时，若量化单元为像元，首要考虑煤炭粉尘污染强度指标，若量化单元为任意损毁单元或其他任意地域范围，量化指标重要性从高到低依次为：煤炭粉尘污染强度、煤炭粉尘污染范围、产尘面积比、万吨产尘量、矿山开发利用状况。
[0228] 主要目标比其余目标重要，构造模糊互补判断矩阵F1～F6时，参与调查的专家依据应用场景以及表1确定重要程度fij。
[0229] 本发明的应用实例：
[0230] 以实现绿色矿山评估为应用场景，以胜利矿区为各量化单元的地域，参与评价的年份包括2006年、2008年、2009年、2011年、2014年～2019年，分别作为各量化单元的年份。
[0231] 量化指标有：煤炭粉尘污染强度算数平均值、煤炭粉尘污染强度中位数、年内持续影响范围、年内最大影响范围、静态产尘面积比、动态产尘面积比、静态万吨产尘量、动态万吨产尘量、矿山生产建设规模和所属生命周期阶段。
[0232] 具体评价步骤如下：
[0233] 1、数据获取。原始遥感影像数据均来源于美国地质调查局，通过GEE平台获取胜利+矿区参评年份共275景Landsat TM/ETM/OLI地表反射率影像数据。
[0234] 2、预处理。通过GEE平台对原始遥感影像数据进行预处理，包括去云处理(SimpleCloudScore算法)、地形校正(基于DEM的朗伯体反射率模型的SCS+C校正法)和影像
裁剪(Quality Mosaic算法)。
[0235] 3、提取各量化指标，得到原始指标值矩阵，见表2。
[0236] 4、计算权重。构造模糊互补判断矩阵F1～F6，见表3～表8，模糊层次法所得权重见表9。
[0237] 计算熵权法所得权重。对原始指标值进行负向标准化处理，得到负向标准化指标值矩阵，见表10。熵权法所得各量化指标权重见表11。
[0238] 5、计算各量化指标权重集，见表12.
[0239] 6、综合评价。基于原始指标值矩阵，得到量化指标分数矩阵，见表13，煤炭粉尘污染程度评价结果见表14。
[0240] 表2原始指标值矩阵
[0241]
[0242] 表3目标层(A)‑一层(B)模糊互补判断矩阵
[0243]
[0244] 表4一层(B2)‑二层(C)模糊互补判断矩阵
[0245]
[0246] 表5一层(B3)‑二层(C)模糊互补判断矩阵
[0247]
[0248] 表6一层(B4)‑二层(C)模糊互补判断矩阵
[0249]
[0250] 表7一层(B5)‑二层(C)模糊互补判断矩阵
[0251]
[0252] 表8一层(B6)‑二层(C)模糊互补判断矩阵
[0253]
[0254] 表9模糊层次法所得量化指标权重
[0255]
[0256] 表10负向标准化指标值矩阵
[0257]
[0258] 表11熵权法所得量化指标权重
[0259]
[0260] 表12量化指标权重集
[0261]
[0262] 表13量化指标分数矩阵
[0263]
[0264] 表14综合评价结果
[0265]
[0266] 以上所述，仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。
