一、景觀生態風險指數計算
0 准備工作
數據:
1.不同時期的土地利用數據
2.研究區矢量邊界
插件:
格網裁剪影像插件tool(下載鏈接:https://pan.baidu.com/s/1r5Fjp8KeH5w1f-0xe3d5xQ 提取碼:7qy9)
軟件:arcgis、total commander、Fragstats4.2
1 生成研究區漁網
1.打開創建漁網工具,設置輸出路徑、模板范圍、像元高度(本文設置20km)、行列數(設置為0或者空值可進行自動計算)、幾何類型(面)
注:像元高度請結合研究區大小參考相關文獻,本文設置20km僅為方便計算。
2.按位置進行選擇,對研究區內的漁網進行選擇,導出。
2 分割漁網
1.打開屬性表,添加字段LID(文本型,長度20)如圖2-1;右鍵新建的表頭LID,打開字段計算器,進行如下字段計算如圖2-2;得到結果如圖2-3。
2-1 2-2 2-3
2.使用分割工具,進行一下設置,分割結果位於新建的數據庫中。
注意:目標工作空間是數據庫,需要在分割前新建一個文件地理數據庫。
3 影像分割
1.將上文中的插件添加到工具箱中,打開工具箱進行設置參數。計算得到若干tif相關文件(時間取決於分割小格的數量)。
注意:文件名一定設置為%Name%.tif
4 景觀格局指數計算
1.將文件中的0.tif(當然也可以是其他的.tif)導入到Fragstats4.2中,對其進行設置,點擊ok。
注意:左側窗口選擇GDAL Geo TIFF grid(.tif),右側背景值為999,文件的路徑必須是英文路徑,否則會出錯。(本文已對路徑進行了修改)
2.然后不進行任何操作進行,將其導出;點擊左側export batch,選擇路徑,生成2000.fbt文件(默認文件格式);將其后綴修改為2000.csv,雙擊打開;
3.使用total commander查找存儲的影像小格文件的路徑,點擊擴展名,擴展名為tif的全選,然后在標記中選中復制完整路徑及名稱。
4.復制在2000.csv中,並將其他列都修改
5.修改2000.csv的后綴名為2000.fbt,然后打開Fragstats4.2,點擊import batch,將2000.fbt導入,導入完成在左側會出現所有的影像斑塊,選擇合適的斑塊級別、景觀級別指標進行運算(Fragstats4.2軟件使用請在網上查閱資料),將計算結果另存為2000,系統會默認生成2000.class和2000.land文件
6.將生成的2000.class和2000.land文件后綴名分別修改為2000_class.csv和2000_land.csv,打開進行景觀指數運算,得到景觀生態風險指數並另存為2000ERI.csv。
(參考: [1]侯蕊,李紅波,高艷麗.基於景觀格局的武漢市江夏區土地利用生態風險評價
[2]侯蕊,李紅波,高艷麗.基於景觀格局的武漢市江夏區土地利用生態風險評價
[3陳心怡 謝跟蹤 張金萍.海口市海岸帶近30年土地利用變化的景觀生態風險評價
等文獻)
二、景觀生態風險制圖
0添加X、Y字段
1.打開屬性表,新建X、Y字段
1 計算質心
1.計算X、Y的質心坐標。分別右鍵X、Y字段,進行計算幾何。
2 將景觀生態風險指數與漁網連接
1.打開arcgis,將漁網和生態風險指數加載,連接2000ERI.csv(LID字段相同,以LID字段進行連接)
3 導出並以XY數據顯示,另存
將連接好的屬性表導出顯示,右鍵導出文件,選擇顯示XY數據。將生成的點文件另存為2000point.shp。
4 克里金插值
1.打開克里金插值的工具,將數據添加,設置參數,本文使用默認(根據需求進行修改)。
注意:需要將研究區20km加載,在處理范圍里的范圍和柵格分析里的掩膜進行選擇。
5出圖
1.加載邊界數據,進行掩膜,設置色帶顯示。
三、寫在最后
1.本文僅此為提供方法,很多細節沒有注意,望大家在學習借鑒時注意起來!!
2.這也是我第一次寫博客,有錯誤希望大家提出指正!!