ERDAS+ArcGis論文技術(shù)匯總

1、碩士畢業(yè)論文過程中所遇到的技術(shù)問題及解答匯總2016.02.231：遙感影像的下載（1）原始數(shù)據(jù)下載網(wǎng)站：地理空間數(shù)據(jù)云（2）打開地理空間數(shù)據(jù)云網(wǎng)站，登陸后高級(jí)檢索選擇行政區(qū)（安徽省六安市）、時(shí)間（2005,2015年）、數(shù)據(jù)集（2005對(duì)應(yīng)landsat7 ETM SLC-off, 2015對(duì)應(yīng)landsat8 OLI-TIRS）搜索。（3）確定所需區(qū)域的地理位置，具體可以在操作更多顯示地理位置。選擇云量較小的數(shù)據(jù)進(jìn)行下載。如果下載數(shù)據(jù)不能包含所有的研究區(qū)域，則需要下載多幅影像在面進(jìn)行拼接。（4）對(duì)于landsat7 ETM SLC-off的遙感圖像需要進(jìn)行帶條修復(fù)模型，按照下圖設(shè)置好相關(guān)的

2、信息，進(jìn)行條帶修復(fù)。 進(jìn)行土地利用類型分析，選擇345波段，其他時(shí)像數(shù)據(jù)中選擇云量較小的數(shù)據(jù)最后在模型服務(wù)我的計(jì)算任務(wù)，查看任務(wù)進(jìn)度。并將數(shù)據(jù)下載下來整理文件夾。2：遙感影像波段融合（1） 打開Erdas，新建一個(gè)viewer. （2）開始對(duì)之前下載的某一時(shí)段的遙感影像中的543波段進(jìn)行融合。具體操作如下主菜單MainImage InterpreterUtilitiesLayerstark（2） 添加需要融合的波段，每添加一個(gè)波段需要點(diǎn)擊Add添加按鈕，指定輸出文件夾。合并完成之后點(diǎn)擊OK退出。（3） 在Viewer中打開合并好的文件，使用放大縮小圖標(biāo)進(jìn)行查看。 3：遙感影像的條帶處理（1）在

3、Erdas菜單欄上選擇MainImage InterpreterRadiometric EnhancementDestripe TM ,按照下圖設(shè)置參數(shù)。（2）完成之后點(diǎn)擊OK，在Viewer中查看去條帶的影像。4：在GIS中創(chuàng)建適應(yīng)于裁剪遙感圖像的行政區(qū)邊界及面域。（1） 所需資料：大范圍的shp文件（六安市鄉(xiāng)鎮(zhèn)邊界線圖層矢量圖）。（2） 打開ARCMAP，創(chuàng)建文件夾鏈接，以便打開儲(chǔ)存文件。具體方法：主菜單窗口目錄文件夾鏈接右鍵選擇相應(yīng)文件夾 （3） 將所需文件拖入窗口中加載，將相應(yīng)圖層可以編輯模式。具體做法：右鍵圖層編輯要素開始編輯出現(xiàn)編輯窗口選擇不必要的線，刪除。（4） 將范圍外面的線打

4、斷，具體做法：先選中一條待分割的線，點(diǎn)擊編輯器的分割工具，在端點(diǎn)處打斷。將多余的線刪除（5） 同時(shí)檢查斷口是否連接閉合。雙擊需要連接的線調(diào)整，將確定連接上的兩條線并。（6） 線轉(zhuǎn)面：確定所有線閉合之后，地理處理Arctoolbox數(shù)據(jù)管理工具要素要素轉(zhuǎn)面形成面之后，單擊面域，檢查是否形成單獨(dú)面域。形成單獨(dú)面域5：坐標(biāo)系及投影變換：查看坐標(biāo)系，并將gis坐標(biāo)系轉(zhuǎn)化成為ERDAS可用的坐標(biāo)系。從GCS-WGS-1984轉(zhuǎn)變?yōu)閁TM/1984 具體做法：查看shp文件的坐標(biāo)系，右鍵點(diǎn)擊目錄中的shp文件查看屬性XY坐標(biāo) 確定需要轉(zhuǎn)換的坐標(biāo)系及區(qū)間，具體做法：在工具欄中單擊ArcToolbox數(shù)據(jù)管理

5、工具投影和變換要素投影，打開投影變化框輸出坐標(biāo)系單擊右側(cè)的瀏覽按鈕，進(jìn)入到“空間參考屬性”對(duì)話框，在對(duì)話框中單擊“選擇”，打開“瀏覽坐標(biāo)系”對(duì)話框，依次找到Projected Coordinate SystemsUTMWGS1984Northern/ Southern HemisphereWGS 1984 UTM Zone *.prj，其中的UTM Zone *代表不同的投影分度帶，根據(jù)城市的經(jīng)緯度坐標(biāo)選擇，具體選擇方法如下：   北半球地區(qū)，選擇最后字母為“N”的

6、帶（在Northern Hemisphere文件夾中），南半球地區(qū)選擇最后字母為“S”的帶（在Southern Hemisphere文件夾中）；帶數(shù)=經(jīng)度整數(shù)位/6的整數(shù)部分+31（東經(jīng)為正值，西經(jīng)為負(fù)值）如：北京約在東經(jīng)116°24'27.09"， 帶數(shù)=116/6+31=50，50N，即WGS 1984 UTM ZONE 50N ，選中后依次單擊確定返回。六安市也是WGS 1984 UTM ZONE 50N。 單擊“確定”執(zhí)行投影變換，自動(dòng)加載

7、到列表內(nèi)容中。6：圖像增強(qiáng)光譜增強(qiáng) 具體做法：點(diǎn)擊interpreterSpectral Enhancementprincipal compent，彈出對(duì)話框，按照對(duì)話框提示輸入數(shù)據(jù)。單擊OK確定。7：裁剪圖像在ERDAS中裁剪遙感圖像。利用行政面域矢量圖shp文件，（采用增強(qiáng)后的影像作為待裁剪影像）（額好吧，增強(qiáng)之后還不如不增強(qiáng)的，換回原來的）打開待裁剪的遙感圖像，再打開六安市規(guī)劃區(qū)矢量圖，發(fā)現(xiàn)覆蓋在遙感圖像上。選中六安市矢量圖，圖像變?yōu)榱咙S色，點(diǎn)擊菜單欄的AOICopy selelacation to AOI，再次選中六安市矢量圖，出現(xiàn)矩形邊框。 將AOI保存到指定文件夾。Filesave

8、AOI layeras，點(diǎn)擊主菜單中的Data reparationsubset Image，在AOI處點(diǎn)擊選擇指定文件AOI file。確定裁剪圖像。8：監(jiān)督分類（1） 打開需要進(jìn)行分類的圖像（2） 打開分類模板編輯器，具體步驟： 點(diǎn)擊主菜單上的Classifier圖標(biāo)ClassificationSignature Editor窗口主面板AOItools（3） 點(diǎn)擊主菜單上的ClassiferSignature Editor，跳出Signature Editor面板。點(diǎn)擊主窗口菜單欄上的AOITools，用多邊形工具畫出一類用地的訓(xùn)練區(qū)，點(diǎn)擊加入特征文件編輯面板。（4）畫

9、一塊，添加一塊，每類用地在重點(diǎn)地域和全區(qū)域選5-10塊。點(diǎn)擊，可將選擇的幾塊地歸為一類。歸并后，點(diǎn)擊Editdelete將原有的刪除。并重命名地類。 （5）將全區(qū)域土地分為幾種地類。保存分類特征信息：在Signature Editor面板上，點(diǎn)擊FileSave As，在彈出的對(duì)話框中選擇路徑并為文件命名，點(diǎn)擊OK。（7） 評(píng)價(jià)分類模板：可能性矩陣評(píng)價(jià)工具（8） 執(zhí)行監(jiān)督分類，具體做法：ERDAS圖標(biāo)面板工具條:點(diǎn)擊ClassifierClassification菜單 Supervised Classification菜單項(xiàng)（9） 定義閾值，確定哪些像元

10、最沒可能正確分類，進(jìn)而對(duì)分類結(jié)果進(jìn)一步優(yōu)化。ClassificationThreshold9：建模（1）查看各類用地的屬性值：點(diǎn)擊主窗口的RasterProfile Tools，在跳出的對(duì)話框選擇Spatial，跳出SPATIAL PROFILE 面板，用在要查詢的用地上拉條線，即可顯示用地屬性值。建筑用地 林地 農(nóng)地 水域 （2）構(gòu)建模型，具體做法：點(diǎn)擊主菜單ModelerModel MakerNew Model窗口，用New Model窗口工具條構(gòu)建運(yùn)算模型，頭和尾分別輸入文件和輸出文件，中間部分雙擊中間過程，跳出Funciton Definitio

11、n對(duì)話框，在Funciton選擇Conditional，再選擇第二個(gè)條件EITHER <arg1> IF (<test> ) OR <arg2> OTHERWISE語句。其中text部分輸入文件名=前面計(jì)算的值，如EITHER 1 IF ($n1_2005fl =9) OR 0 OTHERWISE，點(diǎn)擊OK。再點(diǎn)擊紅色閃電圖標(biāo)運(yùn)行。 相同步驟輸入建筑。林地、農(nóng)地、水體的值，運(yùn)行完成， 建模后的成果10：將柵格數(shù)據(jù)轉(zhuǎn)換成矢量數(shù)據(jù)點(diǎn)擊主菜單上的VectorRaster to 

12、;Vector，在跳出的對(duì)話框中選擇輸入和輸出的文件，點(diǎn)擊OK，跳出Raster to ARC/INFO Coverage對(duì)話框，選擇Map和多邊形Polygon，點(diǎn)擊OK.生成info和四種類型的五個(gè)文件的矢量文件。11：使用gis統(tǒng)計(jì)數(shù)據(jù)并制作專題圖（1）導(dǎo)入剛剛轉(zhuǎn)換好的矢量數(shù)據(jù)（2）打開圖層屬性表，選出屬性為1的數(shù)據(jù)，實(shí)際上就是選擇實(shí)際的用地，去掉邊框底圖。（3）在屬性表主菜單中中打開Options，按屬性選擇，選擇GRID-CODE為1的數(shù)據(jù)。并顯示所選數(shù)據(jù)，剔除不需要的數(shù)據(jù)。選中AREAStatistic統(tǒng)計(jì)，彈出對(duì)話框以此類推各年份各種數(shù)據(jù)的屬性表20

13、15建筑2015林地2015農(nóng)地2015水域（4） 將矢量圖轉(zhuǎn)換成為shp格式輸出成果加載之后可調(diào)整顏色以及邊框樣式（5）修改Girdcode值，便于不同地類的區(qū)分 選中圖層打開屬性表，選中GIRD-CODE 列，右擊選擇Field calculator,（不同圖層設(shè)置不同的值）。（5） 合并圖層，具體步驟地理處理Arctoolbox數(shù)據(jù)管理工具常規(guī)合并改變合成圖像的顏色，具體步驟：右鍵圖層屬性，調(diào)整好之后，點(diǎn)擊應(yīng)用確定 12：FRAGSTATS景觀格局知識(shí)分析將矢量數(shù)據(jù)轉(zhuǎn)換成為柵格數(shù)據(jù)選擇arctoolbox轉(zhuǎn)換柵格要素轉(zhuǎn)柵格選取18個(gè)景觀格局指數(shù)注意：文件的讀

14、取路徑中不能含有中文字符13：Idrisi景觀格局預(yù)測(cè)（1） 數(shù)據(jù)格式的轉(zhuǎn)換：矢量文件轉(zhuǎn)化為柵格文件：10-11步驟將shp文件設(shè)置不同的ID字段，土地利用類型分別編號(hào)為1,2,3.。在arcgis中合并shp，將帶有屬性的合并矢量文件轉(zhuǎn)化為柵格文件。具體做法：柵格文件轉(zhuǎn)為ASCII文件：Arcgis的raster文件轉(zhuǎn)為IDRISI的raster文件，需要用到Arcgis的ASCII文件進(jìn)行轉(zhuǎn)換,arctoolbox轉(zhuǎn)換由柵格轉(zhuǎn)出柵格轉(zhuǎn)ASCII在IDRISI中將ASCII文件轉(zhuǎn)化為IDRISI的raster文件：（2） 重分類轉(zhuǎn)換后的柵格具有的屬性值是土地利用類型的編碼，沒有屬性的值即背景值是默認(rèn)的-9999，因此驚喜土地利用重分類，將-9999轉(zhuǎn)換成為0，即去除背景值。具體步驟為GIS AnalysisDatabase query reclass第一列是賦予的新值，第二列和第三列是舊地類編碼的范圍。由于柵格中的用地代碼是1,2,3,4，于是寫成如上圖所示，背景值為-99990，重新賦予新值0（3） 獲取馬爾科夫鏈矩陣方法:ModelingEnvironmental/Simulation mo