毕业实习报告-矿山地质环境评价

PAGEPAGE12龙岩学院毕业实习报告（矿山地质环境评价）院系：班级：10测绘工程姓名：学实习单实习起止日期：2013.7—2014.4指导老师：目录1.序言 21.1实习项目名称 21.2项目概述 21.3实习时间 21.4实习地点 21.5实习任务 21.6范围及组织情况 22.测区概况 22.1测区范围与概述 22.2已有资料 23.技术依据 24.影像校正 34.1显示图像文件 34.2启动几何纠正模型 34.3启动控制点工具 34.4采集地面控制点 34.5计算转换模型 34.6图像重采样 34.7保存几何纠正模式 34.8检验校正结果 35.数据准备 45.1居民地密度分析 45.2NDVI分析 45.3坡度分析 45.4主要开采方式、主要开采矿种、地质灾害隐患的提取 55.4.1主要开采方式 55.4.2主要开采矿种 55.4.3.地质灾害隐患 55.5降雨量分析 56.综合因子计算 56.1收集评价因子 56.2计算综合因子 57.实习中发生、发现的问题及处理情况 68.实习体会 69.成果附表 7摘要：在矿山地质环境调查的基础上，依据相关标准，采用定量或半定量的方法，评定矿山地质环境问题的危害程度，预测矿山地质环境问题的变化趋势。本文选择安溪县潘田重点工作区为研究区,先收集、整理矿山地质环境调查的资料，然后对相关评价因子进行综合处理，对该地区的地质环境做出等级分类，统计各等级面积及分布情况。1.序言1.1实习项目名称矿山地质环境评价1.2项目概述矿产资源的开发利用,促进我国国民经济的发展和人民生活水平的提高。但随着矿产资源开发规模、强度、深度的不断扩大,也带来了一系列环境问题,如各类矿山地质灾害层出不穷,矿山及矿山周边的生活环境受到不同程度的影响等。为了了解福建省内各地矿产资源开发对环境带来的影响程度，对破坏严重地区采取治理措施，对全省各重点工作区进行矿山地质环境评价，为矿山地质环境保护、恢复治理及监测提供地质依据。1.3实习时间2013年7月-2014年4月1.4实习地点福建省地质测绘院遥感中心1.5实习任务根据相关的规定和要求，完成重点工作区的矿山地质环境评价并统计出各等级面积及分布情况。1.6范围及组织情况测区位于东经117°36′－118°17′，北纬24°50′－25°26′。由5个正式员工及1个实习生每人承担2个工作区共同完成。2.测区概况2.1测区范围与概述本次研究是安溪县重点工作区位于安溪县北部，山峦起伏，山峰林立。安溪县属南、中亚热带海洋性季风气候，境内地表水资源主要来自降雨量。此地植被属于中亚热带常绿阔叶林植被带，气候温和，四季分明，雨量充沛，农业自然资源丰富。境内层峦叠嶂，村落分散，交通不便；梯田层迭，田畴狭隘，耕作条件差。2.2已有资料安溪县2011年Rapideye遥感卫星影像、坡度、降雨量、居民地密度、岩性组合、主要开采方式、主要开采矿种以及地质灾害隐患。3.技术依据3.1《矿山地质环境保护与恢复治理方案编制规范》（DZ/T0223-2011）；3.2《崩塌、滑坡、泥石流监测规范》（DZ/T0221-2006）；3.3《地质灾害分类分级》（DZ0238-2004）；3.4《区域环境地质勘查遥感技术规程》（DZ/T0190-1997）。4.影像校正4.1显示图像文件首先，在ERDAS图标面板中点击Viewer图表两次，打开两个视窗(Viewer1/Viewer2)。在Viewer1中打开需要校正的RE图像,在Viewer2中打开参考底图。然后，将两个视窗平铺放置。4.2启动几何纠正模型选择多项式几何校正模型：Polynomial,同时打开GeoCorrectionTools对话框和PolynomialModelProperties对话框；在PolynomialModelProperties对话框中，定义多项式模型参数以及投影参数。4.3启动控制点工具打开GCPToolReferenseSetup对话框，然后在GCPToolReferenseSetup对话框中选择视窗采点模式；打开ViewerSelectionInstructions指示器。在显示参考底图的Viewer2中点击左键，打开referenceMapInformation提示。4.4采集地面控制点在图像几何校正过程中，采集控制点是一项非常重要和繁重的工作，具体过程如下：首先，在GCP工具对话框中，点击SelectGCP图表，进入GCP选择状态；在Viewer1中移动关联方框位置，寻找明显的地物特征点，作为输入GCP。在Viewer2中，移动关联方框位置，寻找对应的地物特征点，作为参考GCP。然后，在GCP工具对话框中，点击CreateGCP图标，系统将自动将参考点的坐标（X、Y）显示在GCP数据表中。完成后，在GCP对话框中，点击SelectGCP图标，重新进入GCP选择状态，并将光标移回到Viewer1中，准备采集另一个输入控制点。不断重复以上步骤，采集若干控制点GCP，直到满足所选定的几何模型为止。此后，每采集一个InputGCP，系统就自动产生一个Ref.GCP，通过移动Ref.GCP可以优化校正模型。4.5计算转换模型在控制点采集过程中，一般是设置为自动转换计算模型。所以随着控制点采集过程的完成，转换模型就自动计算生成。在Geo-CorrectionTools对话框中，点击DisplayModelProperties图表，可以查阅模型。4.6图像重采样重采样过程就是依据未校正图像的像元值，计算生成一幅校正图像的过程。原图像中所有栅格数据层都在ImageResample对话框中，定义重采样参数。4.7保存几何纠正模式在Geo-CorrectionTools对话框中点击Exit按钮，退出几何校正过程，按照系统提示，可以选择保存图像几何校正模式，并定义模式文件，以便下一次直接利用。4.8检验校正结果同时在两个视窗中打开两幅图像，一幅是矫正以后的图像，一幅是当时的参考图像，通过视窗地理连接功能，及查询光标功能进行目视定性检验。或者用两幅图像叠加显示，通过“拉窗帘”来目视检验。5.数据准备5.1居民地密度分析在ARCMAP中，从数据库中导出地类图斑（DLTB），将工作区涉及到的地类图斑复制到一个图层，再将图斑做一个merge，按工作区范围裁切clip地类图斑，在属性表中添加字段，建制镇、城市为1，村庄为2，剩下的为3；再对添加的字段做一个dissolve，即得到居民地密度。5.2NDVI分析按工作区范围建矢量裁切影像，按裁切好的影像进行NDVI的提取：在PCI里将img转为pix文件；增加一个波段，按re数据源的波段加权方法（调入NDVI_re.EAS）进行计算NDVI；将NDVI按植被覆盖度0-30%、30%-80%以及大于80%进行重分类：对应多少百分比的植被覆盖度查询相应的NDVI植被指数，并将利用对应的NDVI植被指数在arcmap里进行重分类，具体步骤如下：植被覆盖度0-30%、30%-80%以及大于80%对应生成的新影像的像元值分别为1、2、3。将新栅格数据直接转为矢量数据，对生成的矢量数据重新计算面积，将1万面积<5000平方米的图斑去掉。对工作区内分块计算的合并图层（merge）,并最后按1、2、3的属性字段做dissolve5.3坡度分析按工作区范围将各自的DEM数据裁切下来，在arcmap里将裁切好的dem数据进行坡度计算，对生成的坡度文件再进行重分类，坡度按3个等级进行重分类：0-20⁰、20⁰-35⁰以及35⁰以上，生成新的像元值分别对应3、2、1。最后按1、2、3的属性字段做dissolve。5.4主要开采方式、主要开采矿种、地质灾害隐患的提取5.4.1主要开采方式在矢量成果文件里，找到文件名带KM的文件，在KM文件的属性表里，字段KMKCFS里有相应的开采方式代码，如下表：代码值开采方式说明26001露天开采126002地下开采226003联合开采0按该KM文件重新输出一个开采方式的新文件,并与该工作区范围进行union，在属性字段里增加一个新字段“kcfs_dm”,将联合开采、露天、硐采的属性字段分别设置为0、1、2。而其余空的值设置为3。5.4.2主要开采矿种在矢量成果文件里，找到文件名带KM的文件，在KM文件的属性表里，字段KMLX里有相应的主要开采矿种代码，具体开采矿种代码参考“地质矿产术语分类代码”；按该KM文件重新输出一个开采矿种的新文件,并与该工作区范围进行union，在属性字段里增加一个新字段“kckz_dm”,将主要开采矿种分为3类：能源一类、金属和非金属一类，以及无开采一类，分别将这三类的值设置为1、2、3。5.4.3.地质灾害隐患在矢量成果文件里，找到文件名带HP的文件，在HP文件的属性表里，字段HPPG里有相应的滑坡等级描述。按该HP文件重新输出一个地质灾害的新文件,并与该工作区范围进行union，在属性字段里增加一个新字段“dzzh_dm”,将地质灾害隐患分为3类：严重、轻微和无三类，分别将这三类的值设置为1、2、3。（注：无图斑的地方设置为无。）5.5降雨量分析降雨量数据的收集：将该工作区重新转为一个新的矢量数据，并添加一个新的属性字段”jyl_dm”,并将该值赋值为1。6.综合因子计算6.1收集评价因子收集好该工作区的所有因子（坡度、降雨量、NDVI、居民地密度、岩性组合、主要开采方式、主要开采矿种以及地质灾害隐患）。6.2计算综合因子将每个因子的矢量数据直接转为IMG文件：按每个矢量数据的赋值的0、1、2、3代码转为像元值将第一步生成的每个因子的栅格数据调入ERDAS模型中：首先使用gmd模型（有地质灾害隐患的工作区使用含有地质灾害隐患的模型），打开gmd模型。然后分别将每个因子的栅格文件调入该模块中，并在“成果输出”里双击设置输出的文件，最后点击运行。如果生成的栅格数据在arcmap里没有显示正确的投影信息，需让该栅格数据先在erdas中去投影信息。然后在arcmap里转为矢量数据，并将该矢量数据按像元值显示出来：将上述的像元值结合所在工作区的遥感影像重新分为四个区间，这四个区间分别对应严重区、较严重区、一般区、较好区。其中综合值越低，该处环境评价等级越差。每个评价等级区间的值需要参考遥感影像及相关工作区资料。注：一般地情况下，严重区——地质灾害隐患处，矿山开采严重且无植被覆盖较严重区——矿山开采较严重地方且植被覆盖相对较少处一般区——居民地及矿山建筑等区域较好区——农田、水系和植被覆盖较好、坡度较缓及硬质岩区域（一般情况下，每个工作区的这4个等级的面积比例是逐渐递增的。）找出4个区间的阈值之后，在该矢量数据中增加一个text属性字段“评价等级”，长度设置为50；并将4个区间连同确定的阈值一起赋值到对应的属性字段中，如一块图斑的属性为严重区，其阈值区间为43-50，则该属性填“严重区”，其余的对应为“较严重区”、“一般区”、“较好区”。最后再增加一个字段“pjdj_dm”，并将严重区、较严重区、一般区和较好区分别赋值为4、3、2、1。将该矢量重新按“pjdj_dm”字段在arcmap里重新显示，结合遥感影像进行检查修改。最后对矢量进行圆滑图斑处理，如果工作区内有空的地方，需要结合该工作区整个范围进行一次union操作，再对union操作的成果对原来空的地方进行eliminate的操作。然后对矢量进行dissolve处理。最后的成果包含“评价等级”及“pjdj_dm”字段。7.实习中发生、发现的问题及处理情况在处理NDVI的过程中，由于数据量过大ARCMAP无法运行出来。必须把栅格数据给裁切成几个部分分别运行，然后再将各个成果merge起来，得到所需的数据。8.实习体会时光飞逝，大四为期9个月的实习已经结束了。在实习的过程中，我参与了许多遥感项目，也从中得到了许多收获，是我人生中一笔宝贵的财富。但是在收获的同时，我也发现了很多自己的不足之处。之前在学校里面