r
划分土地AOI编辑利用类型,制定分II
52中地葩S3
遥感处理组
.4:
y.倒蹈
2、过程解析
结合流程图以及大赛遥感平台提供的培训视频,我们对所有的操作步
骤进行了整块的分类综合。
首先是:
影像预处理模块,主要包括数据转换、影像融合、影像增强等操作。
其次是:
影像分类模块,主要包括AOI区编
辑、土地利用类型分类、分类后处理、矢量化编辑等相关操作;最后是:
结果输出模块。
以下是对各模块的详细介绍。
2.1、影像预处理模块
2.1.1、数据转换
因从网站下载获取的原数据均为tif格式,又因大赛指定软件为Mapgis,为方便使用和操作,将tif格式全部转换为软件自带msi格式。
具体操作如下:
打开MapGISK9软件遥感处理系统,在菜单栏中单击“数据转换”,选择“影像数据转换”。
在弹出的对话框中单击“添加文件”选择需要转换的原始数据,在“转换设置”中选择结果类型,单击“转换”,即可实现影像数据转换。
转换设置如下图所示:
2.1.2、波段合成
打开MapGISK9软件遥感处理系统,在菜单栏中单击“影响变换”,选择“波段合成”。
在弹出的对话框中单击“输入影像”,选择需要合成的原始数据,“波段选择”中选择“所有波段”单击“添加”,选择“输出影像”地址,单击确定。
2.1.3、创建数据库
打开MapGISK9GB吐业管理器,在MapGISLocal上右击,
创建数据库:
MapGISK9烛翟数捋阵实装器
点开数据库,找到“简单要素类”右击“导入其他数据”,导入已经下载好的“地区界”和“县界”。
p;***備甘
耳.4:
■.>■.倒蹈
42中tts^iaa
找到“栅格数据集”,右击“导入影像”,导入之前“波段合成”
的影像,点击转换:
打开MapGISK9地图编辑器。
从MapGlS数据目录中找到刚才导
入的“地区界”和“县界”上拉,并调整使得“县界”位于“地区界”
下方:
2.1.2、影像截图
鉴于影像数据覆盖面积较大,并且本组所选地区为赣州市和南康市两个面积超过2000千米的市级行政单位。
因此本组将在大赛已提供的数据库中提取并下载全国县级城市行政图,具体见数据显示:
打开MapGISK9地图编辑器。
从MapGlS数据目录中找到刚才导入的“地区界”和“县界”上拉,并调整使得“县界”位于“地区界”下方:
在图上找到南康市和赣州市的界限:
中地麵酉
赣州市南康市及其周边地区
中地匪OH
2CHMi:
3vc¥©eR
选中“地区界”和“县界”,右击选中“编辑”,点击菜单栏中
的“线编辑选择线”,选出南康市和赣州市的线边界:
在“新地图”上右击“新建图层”选择简单要素类,“确定”选择“线
类型”:
f3棗槌叭哉
.4:
-71倒蹈
数据源类别:
圍层类别:
伐件
MapGISS^的图屋期!
1
£简单要青囹层
囲连记囹层
#几何网络组图层
ECAD类图层
曰栅格數橢集层
lira据目录层
匕1对象类图匡
诸要需图层
Q三錐要臺类囹层
简单要秦圏层
名称:
南康市和解州市边畀
feS:
GDBP:
//MapGISLocal/J20137018遥底组01”测览…
I确產二弓取消
选择菜单栏上的“通用编辑”选择“拷贝图元”:
“:
-'"
右击“南康市和赣州市边界”图层选择“当前编辑”,菜单栏上的“通用编辑”选择“粘贴图元”。
保持“当前编辑状态”选中所有线,点击菜单栏的“线编辑”
“线
拓扑处理线拓扑差错”,没有报错,单击“线编辑”“线转弧造区
JLUE-r-^
线柘扑处8(E)
线拓扑查诸(Q)
线转弧直産
濟關徵短线啟)
文件名:
南康市和鵝州市边畀崔区
文件类型:
简单要耒类「乞区文件石二P)
保存
取消
示卖<-j
当前编辑“南康市和赣州市边界造区”影像,点击菜单栏“区编辑合并区”:
K用线分割区(刃画线分割区(U合并区01)
投影转换:
菜单栏“工具”“投影变换类投影”:
LA_祸Xr*KW
i
投影变撫(E)卜
单点投8US}
生成梯形图框(工〕►
厂类投孰E)
生成矩形圉框伍}卜
剛投囲遍L
源简单要素类,选择“南康市和赣州市边界
造区”,“目的简
单要素类”中“目的空间参照系”选择波段合成
.msi自定义,点击
中地匪SB
rrcrwc-i^
确定:
类授影皓按
将波段合成影像和投影转换的影像上拉至图层,调整位置:
使波段合成影像位于上方:
中
影像裁剪:
打开MapGISK9遥感处理系统,菜单栏影像预处理分块裁剪。
输入文件为波段合成的影像,设定裁剪区域中的裁剪模式设置为“矢量区裁剪”,单击下方的“简单要素区”选择之前“投影转换的区域”:
汁——…”;-”FW
毘僅的裁更与取子集
文件:
E:
谨感技能比赛0\L2212'.,2001E7M+...i疲段:
当前影稼第以玄碱段'.J
铀
文件:
(.J
丿设畫0为结果无效值「输岀多幅尉橡
设定载團区械
栽粵模式;矢量区我舅*
坐标类型:
地理”
显示
整景
正方羽
对开
简单要素£
关诃
裁剪完后的效果图:
赣州市+南康市4,3,2假彩色影像截图
以相同方法裁剪该区域第八波段的图像。
第8波段截图
2.1.3、影像融合
影像融合是将同一目标或场景的用不同传感器获得的,或用同种
传感器以不同成像方式,或在不同成像时间获得的不同影像,融合成
一幅影像,在保持多光谱影像辐射信息的同时提高了影像的空间分辨率的遥感影像处理方法。
根据遥感平台提供的几种融合方法,我们将上一节的多光谱加彩色截图与第8波段的全色影像截图进行融合并一一做了操作并最终发现“加权法”得到的结果影像在采用4,3,2波段假彩色显示后与实际地物最相似并且最容易辨认。
本小组将进行2001年与2009年的影像的分别融合并提取细节以供分析参考。
遥感处理组
细节1(主要地物为林地植被)
细节2(主要地物为城市地区,水体)
(以上两幅是2001年数据融合后的细节影像)
中
24
遥感处理组
族夫钊無
细节展示1
(主要地物为城市地区,水体)
细节展示2
(主要地物为耕地,农村居民点,
细节展示3
(主要地物为森林,农村居民点,耕地)
(以上两幅是2009年数据融合后的细节影像)
2.2、影像分类模块
2.2.1、AOI区编辑
在做AOI区编辑之前需要对已经预处理好的影像进行目视解译,对所选地区的土地利用做一个初步的判读,再根据影像数据的4,3,2加彩色显示,并根据4,3,2波段的红,绿,蓝三色的分布与实际地物的对比,对赣州市和南康市作AOI
编辑制作方案,最后在实际情况的基础上作出AOI区域调整
分类,准确选择训练区。
综合以上所有基础准备工作,本组制作出如下2001年和2009年的两个年份的分类对照方案:
26
万案
土地类型
组成
AOI区编辑
植被
耕地、林地、草地
水体
河流、湖泊、池塘、水库
城镇用地
住宅用地、厂房用地、城市建筑
交通运输用地
公路、铁路、乡间道路、飞机场
农村居民点用地
设施农用地
其他
裸地、空闲用地
遥感处理组
族丸444r
(2001年分类对照表)
27
遥感处理组
方案土地类型组成
宦中地葩祖
28
遥感处理组
(2009年分类对照表)
注:
因做A0I编辑的数据源是4,3,2,假彩色显示,且元数据的光
谱分辨率和空间分辨率不高导致多种地物的像元值相近或相同,造成
分类出现很大误差,精度评价不高,因此本小组将结合Googlearth
软件,利用其高空间分辨率类观察所选区域的实际地物,便于较精确
的选取训练区。
打开MapGISK9遥感处理系统,菜单栏影像分类AOI编辑,进行分类编辑:
2.2.2、土地利用类型分类
鉴于遥感的影像分类方法非常多,且各有特点,但总体上来说还是分为监督分类和非监督分来两大类,本小组将对影像作非监督的ISODATA分类和监督分类的最大似然法分类和LVQ2神经网络分类,并对其他分类方法进行尝试。
「is■卖劲舫
29
遥感处理组
ISADATA分类(分类效果不好,几乎看不到地物特征)
广义距离分类
30
-■JS-*创舫
遥感处理组
(2001年最大似然法分类结果)
2009年最大似然法分类(效果最好)
中地麵SB
*劲舫
31
遥感处理组
保存退出AOI编辑。
点击菜单栏影像分类进行监督分类,我们经过比较发现最大似然法效果较好。
监督分类方法,首先需要从研究区域选取有代表性的训练场地作为样本。
根据样本进行参数设置,对样本进行分类和归类。
监督分类最为关键的操作过程是选取训练区,其选取的训练区的准确与否将深刻影响精度评价结果和最终的分类结果。
32
223分类后处理
分类后处理主要包括分类结果编辑、小区处理、精度评价、面积统计、自动生成图例等。
分类结果编辑主要是类型区域合并。
在操作中本小组最终将土地利用类型定为五类,分别是:
耕地、林地、草地、水域、城镇用地、农村居民点用地、其他用地。
采用小区合并进行碎斑的处理,选择像元为30,不宜过大要基本
保持原基础上的分类结果,处理后的效果还是比较好的
分类后结果
打开MapGISK9软件,菜单栏专题产品分类后比较法:
输入i殳置
变化前分冀影像E:
陆感技能比塞數擱122斗2\2D01ETM4Afi:
固
殳化后分类影像曰谨感技能比寿数据U2242-副本课彳像分&
输出设養
输出检测后影像已曜感技能比赛数据卩22伯2D01bTM+llN固输出椅测后振告巳逼感技能比器数据订22舱凶01ETW*M凸监测报告设蛊
标题麥化检测报告
变化后转移矩阵标题变化后
穽化前转移拒阵标题变化前
原始影像单位|0鮭J报吿输出单位口椽素二
确走6厂取消(亍
tJ
分类后处理结果图
赣州市和南康市2001年--2009年土地利用变化
34
-■JS-*劲树
遥感处理组
3、结果分析
52中地曲13
35
1变化转移矩阵
(单位:
像素)
变化后
水体
林地
草地
蜒績用地
农村居民点
交通用丸
林
0857.
23P
灿.
2617.
328D
22S5.
13W.
仍.
lim
112&
371444&
10644M
14141L
】⑹聊
11701,
20357.
5155101
1906.
63902.
69J330.
81336.
110M?
.
349344.
463(7.
2MO5.
137M79.
1241
22KZ7L
652348.
WS4.
724S5.
190707,
157ZL
1282314.
交通用地
1371
韵1
処
5D3.
19813,
3331
3B11
654.
31916.
城績居民
M
6452.
1859.
W2&
M51
8608-.
1W.
T425.
1165.
】刃如一
农牖民
施
備
25522.
2W27.
34701
1U560.
143153
M26Z.
652130.
其它
2337.
183221.
8511?
?
.
S3711.
12570.
B5RW.
41fi47.
27649.
1451401,
合计
HtiSOl
4219433,
3591002,
4«4Z1
应5351
961971
1B3445:
1Q45J3,
mis&w.
2变化统计报表
(单位:
像素)
36
类别
代码
类别面积
变化数
变化率(%)
(A-B)/B*100
变化前但)
变化后(A)
水体
1
131836.
84654.
-47182.
-35.788:
林地
2
5720045.
4221055.H
-1498990.
-26.206
耕地
3
1520940.
3591337.
2070397.
136.13
草地
4
1420236.
446440.
-973796.
-68.566
交通用地
5
35333.
626028.:
590695.
1671.8「
城镇居民地
6
169828.
963049.
793221.
467.07
农村居民地
7
722570.
183454.
-539116.
-74.611
其它
8
1617123.
104545.
-1512578.
-93.535
遥感处理组
表2
赣州市和南康市1999年-2009年的土地利用类型在数量上的变化
有以下规律:
(1)从2001年到2009年,赣州市和南康市土地利用类型中,建筑用地和林地面积是增加的,而水域、草地,未利用地在减少。
这一规律我们可以从表1的变化量一栏中得出。
赣州市和南康市个土地利用类型的变化量中建筑用地和林地是正值,而水域、草地,未利用地均为负值。
(2)从2011年到2009年,赣州市和南康市土地利用类型中的建筑
用地增幅最大。
表2中,建筑用地的变化率在这六种土地利用
类型中达到最大。
而其他五种土地利用类型变化率均在50%以
下。