【ERDAS模型】基于DEM、Slope和AOI的监督分类后处理

1.处理对象

        厦门市行政区域监督分类后的数据。

2.分类（Value）

1：水体（蓝色） 2：林地（绿色） 3：耕地（黄色） 

4：建设用地（紫色） 5：滩涂（灰色） 6：未利用地（白色）

注：若value不是这样需要改正。

图 1Value改正步骤

3.处理原理

（1）高程较高、坡度较陡的耕地其实为林地。

（2）高程较高、坡度较陡的建筑用地其实为林地。

图 2例子

（3）原始图像的滩涂可能会被归于建筑用地或未利用地。

图 3例子

4.处理思路

（1）耕地中（坡度>12°，高程>120）–> 林地

（2）建筑用地中（坡度>12°，高程>140）–> 林地

（3）建筑物且原始图像滩涂部分 –> 滩涂

        在原始图像中圈出滩涂部分，保存为AOI文件，用该文件对建筑物进行裁剪，找出应是滩涂的建筑物。

5.最终结果

（1）水体最终结果 = 原始水体

（2）未利用地最终结果 = 原始未利用地 – 滩涂中的未利用地

（3）林地最终结果 = 原始林地 + 耕地中的林地 + 建筑用地中的林地

（4）耕地最终结果 = 原始耕地 – 耕地中的林地

（5）滩涂最终结果 = AOI中的滩涂

（6）建筑用地最终结果 = 原始建筑用地 – 建筑用地中的林地 + AOI以外的滩涂

6.数据准备

（1）厦门行政区域DEM数据

（2）厦门行政区域坡度数据

（3）各年份的滩涂AOI文件（厦门行政区域的原始图像（未监督分类时的影像）滩涂）

AOI文件制作步骤：

①打开原始图像

②打开AOI工具

图 4打开AOI工具

③圈出滩涂区域

图 5圈出滩涂区域

④保存AOI文件

图 6保存AOI文件

7.校正公式

model01:

（1）分离各类

原始图像 == n (其中n=1,2,3,4,5,6)

（2）耕地中的林地

$n10_gengdi AND $n18_dem30_xiamen > 120  AND $n19_slope > 12

（3）建设用地中的林地

$n11_jianshe AND $n18_dem30_xiamen > 140  AND $n19_slope > 12

（4）最终林地

$n9_lindi  OR  $n16_linj  OR  $n17_ling

（5）最终耕地和建设用地

$n10_gengdi  -  $n16_linj

$n11_jianshe  -  $n17_ling

（6）各类合成

CONDITIONAL { ($n8_shuiti) 1 , ($n21_f_lindi) 2 , ($n24_f_gengdi) 3 , ($n25_f_jianshe) 4 , ($n12_tantu) 5 , ($n13_weiliyong) 6}

model02:

1. 滩涂外区域

$n27_PROMPT_USER - $n28_PROMPT_USER

（2）滩涂外区域的分类

原始图像 == n (其中n=1,2,3,4,5,6)

（3）各类合成

        将原始区域外的滩涂赋值为建设用地，AOI区域的滩涂赋值为滩涂。

CONDITIONAL { ($n37_e_shuiti) 1 , ($n38_e_lindi) 2 , ($n39_e_gengdi) 3 , ($n40_e_jianshe) 4 , ($n41_e_tantu) 3 , ($n42_e_weiliyong) 6 , ($n28_PROMPT_USER) 5 }

8.模型使用步骤

（1）打开model01模型

（2）输入监督分类后的图像

图 7输入监督分类后的图像

（3）输入DEM和slope数据

图8输入DEM和slope数据

（4）更改路径

        更改成自己的路径，数据命名可以改变。

图 9更改路径

（5）运行模型，输出数据

        数据类型为4bit，若不是，要改正。

图 10数据类型为4bit

（6）打开model02模型

（7）输入model01生成的数据

图 11输入model01生成的数据

（8）输入滩涂的AOI文件

图 12输入model01生成的数据

图 13输入滩涂的AOI文件

（9）更改路径

        更改成自己的路径，数据命名可以改变。

（10）输出图像，输出最终数据

        数据类型为4bit，若不是，要改正。