一种基于PSR模型的稀土矿区环境压力量化评估方法及系统技术方案

本发明专利技术涉及稀土矿环境压力评估技术领域，特别是一种基于PSR模型的稀土矿区环境压力量化评估方法，包括以下步骤，步骤S101：数据获取，通过数据获取单元获取所需评估区域原始遥感影像、降雨量数据、地形数据、土壤数据、人口密度数据；步骤S102：数据预处理，将步骤S101获得的遥感影像通过几何校正模块校正影像，大气校正后再掩膜裁剪处理获取研究区遥感影像；步骤S103：遥感影像信息提取，通过遥感信息提取单元解析评估区域地物类别信息，并统计其面积；步骤S104：PSR模型分析。采用上述方法和系统后，本发明专利技术定量评估出稀土矿区的环境安全状况，为稀土矿区环境治理提供技术支持；相比于采用实地调查的方式，本发明专利技术具有省时、省力、高效的特点。

【技术实现步骤摘要】
一种基于PSR模型的稀土矿区环境压力量化评估方法及系统
本专利技术涉及稀土矿环境压力评估
，特别是一种基于PSR模型的稀土矿区环境压力量化评估方法及系统。
技术介绍
稀土作为中国特有的矿产资源，有“工业维生素”和新材料“金库”之称，是我国重要的战略资源。稀土的开采为国民经济发展做出了重要贡献，但同时也引发了严重的生态环境问题。稀土开采过程，也是生态环境破坏的过程，如植被破坏、土地荒漠化、水土流失、尾砂压占等生态环境问题，加剧了矿区生态系统的脆弱程度和退化速度，严重威胁矿区的生态安全。矿区一系列的生态环境问题，已引起人们对矿区资源开采模式的思考，只顾开采不管生态是不可持续的，要实现矿区的可持续开采和发展，必须对矿区的生态环境问题进行监测、评价和管理。我国矿区环境压力评估，更多的依赖于实地调查，对于矿区开采所带来的环境影响与破坏并不能完整估算。稀土矿大多位于偏远山区、山高林密、矿区分散、矿点众多，如果采用以往实地调查不但费时费力，工作效率也很低，对于定量评估矿区环境压力的难度也相对较大。
技术实现思路
本专利技术需要解决的技术问题是提供一种实时、快速对大面积稀土矿区生态环境压力的评估方法。为解决上述技术问题，本专利技术的一种基于PSR模型的稀土矿区环境压力量化评估方法，包括以下步骤，步骤S101：数据获取，通过数据获取单元获取所需评估区域原始遥感影像、降雨量数据、地形数据、土壤数据、人口密度数据；步骤S102：数据预处理，将步骤S101获得的遥感影像通过几何校正模块校正影像，大气校正后再掩膜裁剪处理获取研究区遥感影像；步骤S103：遥感影像信息提取，通过遥感信息提取单元解析评估区域地物类别信息，并统计其面积；步骤S104：PSR模型分析，确定PSR模型指标，并根据PSR模型指标构建PSR模型；通过构建的PSR模型计算稀土矿区环境压力并进行环境压力分级；最后根据稀土矿区环境压力空间分布情况进行分析。进一步的，所述步骤S104具体包括以下步骤，步骤S41：PSR模型确定，确定PSR模型的10个指标，然后分别计算出研究区该10个PSR模型的指标，并对其标准化处理；步骤S42：PSR模型指标加权，加权PSR模型指标，通过综合指数方法构建PSR模型；步骤S43：评估区域环境压力计算，通过构建的PSR模型计算出稀土矿区环境压力，并依据环境压力等级对研究区的环境压力分级；步骤S44：评估区域环境压力分析，依据步骤S43得到的稀土矿区环境压力空间分布情况对其进行分析。进一步的，步骤S01中所述原始遥感影像为Landsat8影像，其中包括波段2-7和波段10，所述地形数据为分辨率为30米的DEM数据。更进一步的，步骤S102中所述DEM数据需通过镶嵌地形数据，再通过裁剪处理获得。进一步的，所述步骤S102中几何校正模块选择地面控制点与影像同名点运用多项式模型几何校正，大气校正通过大气校正模块处理。进一步的，所述步骤S104中PSR模型包括压力层、状态层和响应层3个准则层；其中压力层包括人口密度、荒漠化指数和土壤侵蚀模数3个指标因子；状态层包括植被指数、生态弹性度、生物丰度指数3个指标因子；响应层包括景观破碎度、分维数、香农多样性、稀土矿区地表温度4个指标因子。本专利技术还包括一种基于PSR模型的稀土矿区环境压力量化评估系统，包括数据获取单元、数据处理单元、PSR模型构建单元、遥感信息提取单元和稀土矿区环境压力量化评估单元；所述数据获取单元与数据处理单元输入端相连接，所述数据处理单元输出端与PSR模型构建单元相连接；所述遥感信息提取单元分别与数据处理单元输入端、PSR模型构建单元输入端相连接；所述PSR模型构建单元输出端与稀土矿区环境压力量化评估单元输入端相连接。进一步的，所述遥感信息提取单元包括依次连接的定义训练样本模块、神经网络监督分类模块、评价分类结果模块和分类后处理模块；所述定义训练样本模块与数据处理单元相连接，所述分类后处理模块与PSR模型构建单元相连接。更进一步的，所述定义训练样本模块用于创建稀土矿区不同地物类型的感兴趣区，作为训练样本；所述神经网络监督分类模块是指利用神经网络算法对定义训练的样本进行图像分类；所述评价分类结果模块用于比较分析分类结果和地表真实信息，可采用ENVI提供的混淆矩阵包含的要素：总体分类精度、Kappa系数、混淆概率、错分误差、漏分误差、制图精度以及用户精度来反映；所述分类后处理模块是指利用Majority/Minority分析对分类结果中不可避免产生的一些面积很小的图斑进行剔除或重新分类处理。采用上述方法和系统后，本专利技术提出了结合层次分析法和PSR模型构建了一种矿区生态环境压力评估的方法，定量评估出稀土矿区的环境安全状况，为稀土矿区环境治理提供技术支持；相比于采用实地调查的方式，本专利技术具有省时、省力、高效的特点。附图说明下面结合附图和具体实施方式对本专利技术作进一步详细的说明。图1为本专利技术一种基于PSR模型的稀土矿区环境压力量化评估方法的流程图。图2为本专利技术一种基于PSR模型的稀土矿区环境压力量化评估系统的结构框图。图3为本专利技术PSR模型结构图。其中：1为数据获取单元，2为数据处理单元，3为PSR模型构建单元，4为稀土矿区环境压力量化评估单元，5为遥感信息提取单元具体实施方式如图1所示，本专利技术的一种基于PSR模型的稀土矿区环境压力量化评估方法，步骤S101：数据获取，通过数据获取单元获取所需评估区域原始遥感影像、降雨量数据、地形数据、土壤数据、人口密度数据。所述原始遥感影像为Landsat8影像，其中包括波段2-7和波段10，所述地形数据为分辨率为30米的DEM数据。所述评估区域土壤数据通过土壤数据模块直接导出；所述评估区域降雨量数据及人口密度数据分别通过查阅气象站及统计年鉴获得。步骤S102：数据预处理，将步骤S101获得的遥感影像通过几何校正模块校正影像，大气校正后再掩膜裁剪处理获取研究区遥感影像。所述DEM数据需通过镶嵌地形数据，再通过裁剪处理获得。几何校正模块，选择地面控制点与影像同名点运用多项式模型几何校正，大气校正通过大气校正模块处理，其采用ENVI软件中FLAASH大气校正模块进行，减少或消除大气对遥感影像的影响，然后通过研究区矢量图掩膜裁剪得到研究区影像；研究区地形数据位于3景相邻地形图中，因此需通过ArcCatalog软件将3景相邻地形图拼接成一个大范围、无缝的地形图，然后通过研究区矢量图掩膜裁剪得到研究区DEM数据。步骤S103：遥感影像信息提取，通过遥感信息提取单元解析评估区域地物类别信息，并统计其面积。步骤S104：PSR模型分析，确定PSR模型指标，并根据PSR模型指标构建PSR模型；通过构建的PSR模型计算稀土矿区环境压力并进行环境压力分级；最后根据稀土矿区环境压力空间分布情况进行分析。进一步的，所述步骤S104具体包括以下步骤，步骤S41：PSR模型确定，确定PSR模型的10个指标，然后分别计算出研究区该10个PSR模型的指标，并对其标准化处理。如图3所示，PSR模型包括压力层、状态层和响应层3个准则层；其中压力层包括人口密度、荒漠化指数和土壤侵蚀模数3个指标因子；状态层包括植被指数、生态弹性度、生物丰度指数3个指标因子；响应层包括景观破碎度、分维数、香农多样本文档来自技高网...

【技术保护点】
一种基于PSR模型的稀土矿区环境压力量化评估方法，其特征在于，包括以下步骤，步骤S101：数据获取，通过数据获取单元获取所需评估区域原始遥感影像、降雨量数据、地形数据、土壤数据、人口密度数据；步骤S102：数据预处理，将步骤S101获得的遥感影像通过几何校正模块校正影像，大气校正后再掩膜裁剪处理获取研究区遥感影像；步骤S103：遥感影像信息提取，通过遥感信息提取单元解析评估区域地物类别信息，并统计其面积；步骤S104：PSR模型分析，确定PSR模型指标，并根据PSR模型指标构建PSR模型；通过构建的PSR模型计算稀土矿区环境压力并进行环境压力分级；最后根据稀土矿区环境压力空间分布情况进行分析。

【技术特征摘要】
1.一种基于PSR模型的稀土矿区环境压力量化评估方法，其特征在于，包括以下步骤，步骤S101：数据获取，通过数据获取单元获取所需评估区域原始遥感影像、降雨量数据、地形数据、土壤数据、人口密度数据；步骤S102：数据预处理，将步骤S101获得的遥感影像通过几何校正模块校正影像，大气校正后再掩膜裁剪处理获取研究区遥感影像；步骤S103：遥感影像信息提取，通过遥感信息提取单元解析评估区域地物类别信息，并统计其面积；步骤S104：PSR模型分析，确定PSR模型指标，并根据PSR模型指标构建PSR模型；通过构建的PSR模型计算稀土矿区环境压力并进行环境压力分级；最后根据稀土矿区环境压力空间分布情况进行分析。2.按照权利要求1所述的一种基于PSR模型的稀土矿区环境压力量化评估方法，其特征在于，所述步骤S104具体包括以下步骤，步骤S41：PSR模型确定，确定PSR模型的10个指标，然后分别计算出研究区该10个PSR模型的指标，并对其标准化处理；步骤S42：PSR模型指标加权，加权PSR模型指标，通过综合指数方法构建PSR模型；步骤S43：评估区域环境压力计算，通过构建的PSR模型计算出稀土矿区环境压力，并依据环境压力等级对研究区的环境压力分级；步骤S44：评估区域环境压力分析，依据步骤S43得到的稀土矿区环境压力空间分布情况对其进行分析。3.按照权利要求1所述的一种基于PSR模型的稀土矿区环境压力量化评估方法，其特征在于：步骤S01中所述原始遥感影像为Landsat8影像，其中包括波段2-7和波段10，所述地形数据为分辨率为30米的DEM数据。4.按照权利要求3所述的一种基于PSR模型的稀土矿区环境压力量化评估方法，其特征在于：步骤S102中所述DEM数据需通过镶嵌地形数据，再通过裁剪处理获得。5.按照权利要求1所述的一种基于PSR模型的稀土矿区环境压力量化评估方法，其特征在于：所述步骤S102中几何校正模块选择地面控制点与影像...

【专利技术属性】
技术研发人员：李恒凯，杨柳，雷军，王秀丽，
申请(专利权)人：江西理工大学，
类型：发明
国别省市：江西,36