畢業(yè)實(shí)習(xí)報(bào)告-礦山地質(zhì)環(huán)境評(píng)價(jià)

PAGEPAGE12龍巖學(xué)院畢業(yè)實(shí)習(xí)報(bào)告（礦山地質(zhì)環(huán)境評(píng)價(jià)）院系：班級(jí)：10測(cè)繪工程姓名：學(xué)實(shí)習(xí)單實(shí)習(xí)起止日期：2013.7—2014.4指導(dǎo)老師：目錄1.序言 21.1實(shí)習(xí)項(xiàng)目名稱 21.2項(xiàng)目概述 21.3實(shí)習(xí)時(shí)間 21.4實(shí)習(xí)地點(diǎn) 21.5實(shí)習(xí)任務(wù) 21.6范圍及組織情況 22.測(cè)區(qū)概況 22.1測(cè)區(qū)范圍與概述 22.2已有資料 23.技術(shù)依據(jù) 24.影像校正 34.1顯示圖像文件 34.2啟動(dòng)幾何糾正模型 34.3啟動(dòng)控制點(diǎn)工具 34.4采集地面控制點(diǎn) 34.5計(jì)算轉(zhuǎn)換模型 34.6圖像重采樣 34.7保存幾何糾正模式 34.8檢驗(yàn)校正結(jié)果 35.數(shù)據(jù)準(zhǔn)備 45.1居民地密度分析 45.2NDVI分析 45.3坡度分析 45.4主要開采方式、主要開采礦種、地質(zhì)災(zāi)害隱患的提取 55.4.1主要開采方式 55.4.2主要開采礦種 55.4.3.地質(zhì)災(zāi)害隱患 55.5降雨量分析 56.綜合因子計(jì)算 56.1收集評(píng)價(jià)因子 56.2計(jì)算綜合因子 57.實(shí)習(xí)中發(fā)生、發(fā)現(xiàn)的問題及處理情況 68.實(shí)習(xí)體會(huì) 69.成果附表 7摘要：在礦山地質(zhì)環(huán)境調(diào)查的基礎(chǔ)上，依據(jù)相關(guān)標(biāo)準(zhǔn)，采用定量或半定量的方法，評(píng)定礦山地質(zhì)環(huán)境問題的危害程度，預(yù)測(cè)礦山地質(zhì)環(huán)境問題的變化趨勢(shì)。本文選擇安溪縣潘田重點(diǎn)工作區(qū)為研究區(qū),先收集、整理礦山地質(zhì)環(huán)境調(diào)查的資料，然后對(duì)相關(guān)評(píng)價(jià)因子進(jìn)行綜合處理，對(duì)該地區(qū)的地質(zhì)環(huán)境做出等級(jí)分類，統(tǒng)計(jì)各等級(jí)面積及分布情況。1.序言1.1實(shí)習(xí)項(xiàng)目名稱礦山地質(zhì)環(huán)境評(píng)價(jià)1.2項(xiàng)目概述礦產(chǎn)資源的開發(fā)利用,促進(jìn)我國(guó)國(guó)民經(jīng)濟(jì)的發(fā)展和人民生活水平的提高。但隨著礦產(chǎn)資源開發(fā)規(guī)模、強(qiáng)度、深度的不斷擴(kuò)大,也帶來了一系列環(huán)境問題,如各類礦山地質(zhì)災(zāi)害層出不窮,礦山及礦山周邊的生活環(huán)境受到不同程度的影響等。為了了解福建省內(nèi)各地礦產(chǎn)資源開發(fā)對(duì)環(huán)境帶來的影響程度，對(duì)破壞嚴(yán)重地區(qū)采取治理措施，對(duì)全省各重點(diǎn)工作區(qū)進(jìn)行礦山地質(zhì)環(huán)境評(píng)價(jià)，為礦山地質(zhì)環(huán)境保護(hù)、恢復(fù)治理及監(jiān)測(cè)提供地質(zhì)依據(jù)。1.3實(shí)習(xí)時(shí)間2013年7月-2014年4月1.4實(shí)習(xí)地點(diǎn)福建省地質(zhì)測(cè)繪院遙感中心1.5實(shí)習(xí)任務(wù)根據(jù)相關(guān)的規(guī)定和要求，完成重點(diǎn)工作區(qū)的礦山地質(zhì)環(huán)境評(píng)價(jià)并統(tǒng)計(jì)出各等級(jí)面積及分布情況。1.6范圍及組織情況測(cè)區(qū)位于東經(jīng)117°36′－118°17′，北緯24°50′－25°26′。由5個(gè)正式員工及1個(gè)實(shí)習(xí)生每人承擔(dān)2個(gè)工作區(qū)共同完成。2.測(cè)區(qū)概況2.1測(cè)區(qū)范圍與概述本次研究是安溪縣重點(diǎn)工作區(qū)位于安溪縣北部，山巒起伏，山峰林立。安溪縣屬南、中亞熱帶海洋性季風(fēng)氣候，境內(nèi)地表水資源主要來自降雨量。此地植被屬于中亞熱帶常綠闊葉林植被帶，氣候溫和，四季分明，雨量充沛，農(nóng)業(yè)自然資源豐富。境內(nèi)層巒疊嶂，村落分散，交通不便；梯田層迭，田疇狹隘，耕作條件差。2.2已有資料安溪縣2011年Rapideye遙感衛(wèi)星影像、坡度、降雨量、居民地密度、巖性組合、主要開采方式、主要開采礦種以及地質(zhì)災(zāi)害隱患。3.技術(shù)依據(jù)3.1《礦山地質(zhì)環(huán)境保護(hù)與恢復(fù)治理方案編制規(guī)范》（DZ/T0223-2011）；3.2《崩塌、滑坡、泥石流監(jiān)測(cè)規(guī)范》（DZ/T0221-2006）；3.3《地質(zhì)災(zāi)害分類分級(jí)》（DZ0238-2004）；3.4《區(qū)域環(huán)境地質(zhì)勘查遙感技術(shù)規(guī)程》（DZ/T0190-1997）。4.影像校正4.1顯示圖像文件首先，在ERDAS圖標(biāo)面板中點(diǎn)擊Viewer圖表兩次，打開兩個(gè)視窗(Viewer1/Viewer2)。在Viewer1中打開需要校正的RE圖像,在Viewer2中打開參考底圖。然后，將兩個(gè)視窗平鋪放置。4.2啟動(dòng)幾何糾正模型選擇多項(xiàng)式幾何校正模型：Polynomial,同時(shí)打開GeoCorrectionTools對(duì)話框和PolynomialModelProperties對(duì)話框；在PolynomialModelProperties對(duì)話框中，定義多項(xiàng)式模型參數(shù)以及投影參數(shù)。4.3啟動(dòng)控制點(diǎn)工具打開GCPToolReferenseSetup對(duì)話框，然后在GCPToolReferenseSetup對(duì)話框中選擇視窗采點(diǎn)模式；打開ViewerSelectionInstructions指示器。在顯示參考底圖的Viewer2中點(diǎn)擊左鍵，打開referenceMapInformation提示。4.4采集地面控制點(diǎn)在圖像幾何校正過程中，采集控制點(diǎn)是一項(xiàng)非常重要和繁重的工作，具體過程如下：首先，在GCP工具對(duì)話框中，點(diǎn)擊SelectGCP圖表，進(jìn)入GCP選擇狀態(tài)；在Viewer1中移動(dòng)關(guān)聯(lián)方框位置，尋找明顯的地物特征點(diǎn)，作為輸入GCP。在Viewer2中，移動(dòng)關(guān)聯(lián)方框位置，尋找對(duì)應(yīng)的地物特征點(diǎn)，作為參考GCP。然后，在GCP工具對(duì)話框中，點(diǎn)擊CreateGCP圖標(biāo)，系統(tǒng)將自動(dòng)將參考點(diǎn)的坐標(biāo)（X、Y）顯示在GCP數(shù)據(jù)表中。完成后，在GCP對(duì)話框中，點(diǎn)擊SelectGCP圖標(biāo)，重新進(jìn)入GCP選擇狀態(tài)，并將光標(biāo)移回到Viewer1中，準(zhǔn)備采集另一個(gè)輸入控制點(diǎn)。不斷重復(fù)以上步驟，采集若干控制點(diǎn)GCP，直到滿足所選定的幾何模型為止。此后，每采集一個(gè)InputGCP，系統(tǒng)就自動(dòng)產(chǎn)生一個(gè)Ref.GCP，通過移動(dòng)Ref.GCP可以優(yōu)化校正模型。4.5計(jì)算轉(zhuǎn)換模型在控制點(diǎn)采集過程中，一般是設(shè)置為自動(dòng)轉(zhuǎn)換計(jì)算模型。所以隨著控制點(diǎn)采集過程的完成，轉(zhuǎn)換模型就自動(dòng)計(jì)算生成。在Geo-CorrectionTools對(duì)話框中，點(diǎn)擊DisplayModelProperties圖表，可以查閱模型。4.6圖像重采樣重采樣過程就是依據(jù)未校正圖像的像元值，計(jì)算生成一幅校正圖像的過程。原圖像中所有柵格數(shù)據(jù)層都在ImageResample對(duì)話框中，定義重采樣參數(shù)。4.7保存幾何糾正模式在Geo-CorrectionTools對(duì)話框中點(diǎn)擊Exit按鈕，退出幾何校正過程，按照系統(tǒng)提示，可以選擇保存圖像幾何校正模式，并定義模式文件，以便下一次直接利用。4.8檢驗(yàn)校正結(jié)果同時(shí)在兩個(gè)視窗中打開兩幅圖像，一幅是矯正以后的圖像，一幅是當(dāng)時(shí)的參考圖像，通過視窗地理連接功能，及查詢光標(biāo)功能進(jìn)行目視定性檢驗(yàn)?；蛘哂脙煞鶊D像疊加顯示，通過“拉窗簾”來目視檢驗(yàn)。5.數(shù)據(jù)準(zhǔn)備5.1居民地密度分析在ARCMAP中，從數(shù)據(jù)庫中導(dǎo)出地類圖斑（DLTB），將工作區(qū)涉及到的地類圖斑復(fù)制到一個(gè)圖層，再將圖斑做一個(gè)merge，按工作區(qū)范圍裁切clip地類圖斑，在屬性表中添加字段，建制鎮(zhèn)、城市為1，村莊為2，剩下的為3；再對(duì)添加的字段做一個(gè)dissolve，即得到居民地密度。5.2NDVI分析按工作區(qū)范圍建矢量裁切影像，按裁切好的影像進(jìn)行NDVI的提?。涸赑CI里將img轉(zhuǎn)為pix文件；增加一個(gè)波段，按re數(shù)據(jù)源的波段加權(quán)方法（調(diào)入NDVI_re.EAS）進(jìn)行計(jì)算NDVI；將NDVI按植被覆蓋度0-30%、30%-80%以及大于80%進(jìn)行重分類：對(duì)應(yīng)多少百分比的植被覆蓋度查詢相應(yīng)的NDVI植被指數(shù)，并將利用對(duì)應(yīng)的NDVI植被指數(shù)在arcmap里進(jìn)行重分類，具體步驟如下：植被覆蓋度0-30%、30%-80%以及大于80%對(duì)應(yīng)生成的新影像的像元值分別為1、2、3。將新柵格數(shù)據(jù)直接轉(zhuǎn)為矢量數(shù)據(jù)，對(duì)生成的矢量數(shù)據(jù)重新計(jì)算面積，將1萬面積<5000平方米的圖斑去掉。對(duì)工作區(qū)內(nèi)分塊計(jì)算的合并圖層（merge）,并最后按1、2、3的屬性字段做dissolve5.3坡度分析按工作區(qū)范圍將各自的DEM數(shù)據(jù)裁切下來，在arcmap里將裁切好的dem數(shù)據(jù)進(jìn)行坡度計(jì)算，對(duì)生成的坡度文件再進(jìn)行重分類，坡度按3個(gè)等級(jí)進(jìn)行重分類：0-20?、20?-35?以及35?以上，生成新的像元值分別對(duì)應(yīng)3、2、1。最后按1、2、3的屬性字段做dissolve。5.4主要開采方式、主要開采礦種、地質(zhì)災(zāi)害隱患的提取5.4.1主要開采方式在矢量成果文件里，找到文件名帶KM的文件，在KM文件的屬性表里，字段KMKCFS里有相應(yīng)的開采方式代碼，如下表：代碼值開采方式說明26001露天開采126002地下開采226003聯(lián)合開采0按該KM文件重新輸出一個(gè)開采方式的新文件,并與該工作區(qū)范圍進(jìn)行union，在屬性字段里增加一個(gè)新字段“kcfs_dm”,將聯(lián)合開采、露天、硐采的屬性字段分別設(shè)置為0、1、2。而其余空的值設(shè)置為3。5.4.2主要開采礦種在矢量成果文件里，找到文件名帶KM的文件，在KM文件的屬性表里，字段KMLX里有相應(yīng)的主要開采礦種代碼，具體開采礦種代碼參考“地質(zhì)礦產(chǎn)術(shù)語分類代碼”；按該KM文件重新輸出一個(gè)開采礦種的新文件,并與該工作區(qū)范圍進(jìn)行union，在屬性字段里增加一個(gè)新字段“kckz_dm”,將主要開采礦種分為3類：能源一類、金屬和非金屬一類，以及無開采一類，分別將這三類的值設(shè)置為1、2、3。5.4.3.地質(zhì)災(zāi)害隱患在矢量成果文件里，找到文件名帶HP的文件，在HP文件的屬性表里，字段HPPG里有相應(yīng)的滑坡等級(jí)描述。按該HP文件重新輸出一個(gè)地質(zhì)災(zāi)害的新文件,并與該工作區(qū)范圍進(jìn)行union，在屬性字段里增加一個(gè)新字段“dzzh_dm”,將地質(zhì)災(zāi)害隱患分為3類：嚴(yán)重、輕微和無三類，分別將這三類的值設(shè)置為1、2、3。（注：無圖斑的地方設(shè)置為無。）5.5降雨量分析降雨量數(shù)據(jù)的收集：將該工作區(qū)重新轉(zhuǎn)為一個(gè)新的矢量數(shù)據(jù)，并添加一個(gè)新的屬性字段”jyl_dm”,并將該值賦值為1。6.綜合因子計(jì)算6.1收集評(píng)價(jià)因子收集好該工作區(qū)的所有因子（坡度、降雨量、NDVI、居民地密度、巖性組合、主要開采方式、主要開采礦種以及地質(zhì)災(zāi)害隱患）。6.2計(jì)算綜合因子將每個(gè)因子的矢量數(shù)據(jù)直接轉(zhuǎn)為IMG文件：按每個(gè)矢量數(shù)據(jù)的賦值的0、1、2、3代碼轉(zhuǎn)為像元值將第一步生成的每個(gè)因子的柵格數(shù)據(jù)調(diào)入ERDAS模型中：首先使用gmd模型（有地質(zhì)災(zāi)害隱患的工作區(qū)使用含有地質(zhì)災(zāi)害隱患的模型），打開gmd模型。然后分別將每個(gè)因子的柵格文件調(diào)入該模塊中，并在“成果輸出”里雙擊設(shè)置輸出的文件，最后點(diǎn)擊運(yùn)行。如果生成的柵格數(shù)據(jù)在arcmap里沒有顯示正確的投影信息，需讓該柵格數(shù)據(jù)先在erdas中去投影信息。然后在arcmap里轉(zhuǎn)為矢量數(shù)據(jù)，并將該矢量數(shù)據(jù)按像元值顯示出來：將上述的像元值結(jié)合所在工作區(qū)的遙感影像重新分為四個(gè)區(qū)間，這四個(gè)區(qū)間分別對(duì)應(yīng)嚴(yán)重區(qū)、較嚴(yán)重區(qū)、一般區(qū)、較好區(qū)。其中綜合值越低，該處環(huán)境評(píng)價(jià)等級(jí)越差。每個(gè)評(píng)價(jià)等級(jí)區(qū)間的值需要參考遙感影像及相關(guān)工作區(qū)資料。注：一般地情況下，嚴(yán)重區(qū)——地質(zhì)災(zāi)害隱患處，礦山開采嚴(yán)重且無植被覆蓋較嚴(yán)重區(qū)——礦山開采較嚴(yán)重地方且植被覆蓋相對(duì)較少處一般區(qū)——居民地及礦山建筑等區(qū)域較好區(qū)——農(nóng)田、水系和植被覆蓋較好、坡度較緩及硬質(zhì)巖區(qū)域（一般情況下，每個(gè)工作區(qū)的這4個(gè)等級(jí)的面積比例是逐漸遞增的。）找出4個(gè)區(qū)間的閾值之后，在該矢量數(shù)據(jù)中增加一個(gè)text屬性字段“評(píng)價(jià)等級(jí)”，長(zhǎng)度設(shè)置為50；并將4個(gè)區(qū)間連同確定的閾值一起賦值到對(duì)應(yīng)的屬性字段中，如一塊圖斑的屬性為嚴(yán)重區(qū)，其閾值區(qū)間為43-50，則該屬性填“嚴(yán)重區(qū)”，其余的對(duì)應(yīng)為“較嚴(yán)重區(qū)”、“一般區(qū)”、“較好區(qū)”。最后再增加一個(gè)字段“pjdj_dm”，并將嚴(yán)重區(qū)、較嚴(yán)重區(qū)、一般區(qū)和較好區(qū)分別賦值為4、3、2、1。將該矢量重新按“pjdj_dm”字段在arcmap里重新顯示，結(jié)合遙感影像進(jìn)行檢查修改。最后對(duì)矢量進(jìn)行圓滑圖斑處理，如果工作區(qū)內(nèi)有空的地方，需要結(jié)合該工作區(qū)整個(gè)范圍進(jìn)行一次union操作，再對(duì)union操作的成果對(duì)原來空的地方進(jìn)行eliminate的操作。然后對(duì)矢量進(jìn)行dissolve處理。最后的成果包含“評(píng)價(jià)等級(jí)”及“pjdj_dm”字段。7.實(shí)習(xí)中發(fā)生、發(fā)現(xiàn)的問題及處理情況在處理NDVI的過程中，由于數(shù)據(jù)量過大ARCMAP無法運(yùn)行出來。必須把柵格數(shù)據(jù)給裁切成幾個(gè)部分分別運(yùn)行，然后再將各個(gè)成果merge起來，得到所需的數(shù)據(jù)。8.實(shí)習(xí)體會(huì)時(shí)光飛逝，大四為期9個(gè)月的實(shí)習(xí)已經(jīng)結(jié)束了。在實(shí)習(xí)的過程中，我參與了許多遙感項(xiàng)目，也從中得到了許多收獲，是我人生中一筆寶貴的財(cái)富。但是在收獲的同時(shí)，我也發(fā)現(xiàn)了很多自己的不足之處。之前在學(xué)校里面