基于ERDAS模型校正分类结果

1. 基本原理

1.1 校正依据

分类优先级：
1：水体 2：林地 3：滩涂 4：耕地 5：未利用地 6：建筑用地

规定:
耕地中【坡度>12，高程>120】=> 林地
建筑用地中【坡度>12，高程>140】=> 林地

最终结果校正公式：
水体最终结果 = 原始水体
未利用地最终结果 = 原始未利用地
林地最终结果 = 原始林地 + 建筑用地中的林地 + 耕地中的林地
耕地最终结果 = 原始耕地 - 耕地中的林地
滩涂最终结果 = AOI中的滩涂
建筑用地最终结果 = 原始建筑用地–建筑用地中的林地 + AOI以外的滩涂

1.2 分离六类地物算法

算法公式：$n1_PROMPT_USER == n (其中n=1,2,3,4,5,6)

1）分离水体：$n1_PROMPT_USER == 1

2）分离滩涂：$n1_PROMPT_USER == 3

输出文件命名：
1：水（water）2：林地（lindi） 3：滩涂（tantu） 4：耕地（gengdi） 5：未利用地（weiliyong） 6：建筑用地（building）

1.3 利用坡度和高程分离地物算法

算法原始公式： EITHER IF ( ) OR OTHERWISE
1）提取耕地中的林地：
EITHER $n11_gengdi IF ( $n14_slope > 12 or $n16_dem30 > 120 ) OR 0 OTHERWISE
2）提取建筑用地中的林地：
EITHER $n13_building IF ( $n14_slope > 12 or $n16_dem30 > 140 ) OR 0 OTHERWISE

1.4 图像合并或去除算法

1）去除耕地中的林地：$n11_gengdi - $n17_liningeng
2）合并原始林地、建筑用地中的林地和耕地中的林地：
$n9_lindi OR $n23_lininbuild OR $n17_liningeng
3）原始林地、去除林地后建筑用地与滩涂中的建筑用地合并：
$n13_building - $n23_lininbuild + $n32_buildintan

1.5 利用AOI分类

分离滩涂与建设用地：

使用AOI裁出原始分类结果滩涂部分:


二值化滩涂最终结果（最后合并六类分类结果时需要）

1.6 校正后分类结果合并

算法原始公式：
CONDITIONAL { (<>) <> , (<>) <> , … }
1）合并六类结果：公式如下图
说明：（<>）中填入分级后影像名称，<>中填入优先级数字，需与原始图像相同

1.7 模型示意图

2. 模型制作过程

2.1 分离六类地物模型

说明： 勾选【Prompt User for File at Run Time】模型运行时，可跳出对话框选择文件
Data Type为Unsigned 8-bit


说明： 文件名类别与分类一致输出结果只有一类地物，原始图像只需判断是否为该类地物，因此Data Type为Unsigned 1-bit，File Type为Thematic。
（文件名与算法公式根据地物类型不同进行更改）

2.2 利用坡度和高程分离地物

2.3 利用AOI分类（详见1.5）

2.4 校正剩余类别

林地校正：


耕地校正：


建筑用地校正：

2.5 合并最终分类结果

说明： 输入的六个分类结果都为校正后


说明： Data Type为Unsigned 4-bit，File Type为Thematic

3. 模型使用

点击运行后需输入两次原始分类结果图（两次输入图像应为同一张影像）：


运行过程中若出现以下提示可点击【OK to ALL】->【Dismiss】