采用中分辨率遥感影像提取土地利用变化信息

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
Leabharlann Baidu
0



土地利用动态遥感监测是国土资源部新一轮国土资 源大调查的重要组成部分,它的实施对于加强土地资源 管理,确保耕地总量动态平衡和占补平衡目标的实现有 着重要意义。随着航天遥感技术的快速发展和应用领域 的拓宽,利用多种信息源、不同分辨率的卫星影像进行 土地利用遥感动态监测也正在逐步推广应用。目前我国 的土地利用遥感动态监测,主要是对耕地和建设用地等 土地利用变化情况进行及时、 直接、 客观的定期监测[1,14], 在遥感监测方法上目前主要采用的是目视解译法在遥感 影像上对土地利用变化信息进行发现并提取,这就对影 像分辨率要求较高,但由于各种原因目前大部分都是以 中分辨率影像进行监测为主 [3-7,15-19],在中分辨率影像用 目视解译法比较困难,而且费时,人为因素的干扰造成 数据不准确。 针对以上问题本文以在不同时相、不同空间分辨率 的 TM、SPOT2 的影像上分析岳池县花园镇的基本农田 变化为例,探索出一种变化信息自动发现提取的方法, 即在分析比较的基础上采用图像差值法来发现变化信 息,然后运用基于光谱知识法(域值法)[20]来提取变化 信息的方法,从而确定变化范围及面积,以解决目视解 译法在中分辨率影像中费时且不易提取等问题。
图 2 SPOT1-TM2 中提取出的土地变化信息 Fig.2 Change information extracted from the difference of band 1 of SPOT and band 2 of TM
利用 ERDAS IMAGINE 将提取出的变化信息首先二 值化,然后对二值图像进行矢量化,如图 3,其中绿色信 息为自动矢量化后的变化信息。
3.1 影像预处理 3.1.1 影像校正 本文所获得的影像已经经过几何校正和辐射校正, 但是两影像同名地物重合得不够好,有移位现象。文中 选用 TM 作为标准图像,分别在两影像上选取控制点, 用三次卷积内插运算对 SPOT 影像进行重采样,将 SPOT
86
农业工程学报
2008 年
影像与其配准, 并重采样成同空间分辨率。 然后将 TM 和 SPOT 影像组合为一个影像文件, 便于下面的土地利用变 化信息识别工作。 3.1.2 相对辐射校正 相对辐射校正能够减小多时相遥感图像之间由大 气、照度、物候和传感器标度等差异造成的影响,有利 于提高动态监测精度。一般以一个时相的图像作为参考 图像,即主图像,对另一时相的图像进行辐射校正,即 从图像,把这两个时相的图像称为主-从像对。即利用图 像像元灰度值,建立多时相遥感图像各个波段之间的校 正方程,对遥感图像进行归一化处理。 由于最终以 2001 年的 TM 图像作为确定变化信息的 基年数据,因此,以 TM 图像作为参考图像,校正 2006 年的 SPOT 图像. 由于 SOPT 的 1、2、3 波段分别类似对 应于 TM 的 2、3、4 波段,因此,分别以 TM 的 2、3、4 波段来校正 SOPT 的 1、2、3 波段,分别建立其线形回 归方程后,其线形回归方程分别如下: m = x× n + y 式中 m从图像某波段的归一化像元灰度值,n 主图像某波段的像元灰度值, x , y 分别 归一化常 数[13]。经运算,SOPT 的 1、2、3 波段与 TM 的 2、3、4 波段的线形回归方程分别为: SPOT1=1.02× TM2 + 43.56 SPOT2=0.70× TM3 + 69.70 SPOT3=0.80× TM4 + 84.39 经过运算后得到校正后的 SOPT 的 1、2、3 波段。 3.2 土地利用变化信息的发现 本文采用图像差值法来发现土地利用变化信息,将 两个时相的遥感影像按波段进行逐像元相减,从而生成 一幅新的、代表两个时相间光谱变化的差值影像。假设 辐射值的显著变化代表了土地覆盖变化, 在差值影像中 接近于零的像元视为未变化, 大于或小于零的像元表示 其覆盖状况发生了变化。 由于 SPOT 的 1 波段与 TM 2 波 段对应,因此,将 SPOT1-TM2,在差值中接近于零的像 元视为未变化,图 1 中红色部分的像元的表示土地利用 变化信息。
2
技术方法及技术路线
1
1.1
研究区及资料收集概况
研究区概况 四 川 省 岳 池 县 花 园 镇 位 于 北 纬 29° 57'43" ~
修订日期:2008-07-28
技术方法:基于遥感信息自动提取为主的技术方法。 即在中低分辨率影像中采用基于光谱知识的自动提取法 对研究区的变化信息进行提取。 技术路线:由于采用的是不同时相、不同遥感平台、 不同分辨率的 TM 和 SPOT 数据,因此首先进行图像几 何纠正、图像间的匹配、重采样、相对辐射校正等处 理[13,16,17];然后识别两期影像中土地利用变化信息;接下 来进行提取;最后进行精度评价和总结。本文以花园镇 的 TM 和 SPOT 影像来识别并提取从 2001 年到 2006 年 该区基本农田变化信息,并进一步讨论中低分辨率影像 中土地利用变化信息的提取方法。
3
中空间分辨率影像中土地利用变化信息的提取
收稿日期:2008-03-27
作者简介:吴 玺(1975-) ,博士生,土地资源管理专业,主要从事土地 整理方面工作。北京 中国农业大学资源与环境学院,100094。 Email:028-85539903@163.com ※通讯作者:王小燕(1981-) ,硕士,地理信息系统专业,主要从事遥感 与 GIS 的应用方面工作。 成都 四川省农业科学院遥感应用研究所, 610066。 Email:wxyswallow@tom.com
本文针对以上发现土地利用变化信息,采用域值法 即基于光谱知识的方法进行自动提取法。 经过上面变化信息的识别,基本确定了其变化的位 置及范围,采用图像差值法后的某个波段影像进行自动 提取。因为 SPOT 的 1、2、3 波段分别类似对应于 TM2、 3、 4 波段, 因此, 分别将相对辐射校正后的 SPOT1-TM2、 SPOT2-TM3、SPOT3-TM4,经分析,只有 SPOT1-TM2 后的波段影像能较好的反映出变化信息。 再经过反复试验,在 SPOT1-TM2 波段影像上,当灰 度值在 0 到 5 之间就是变化了土地利用信息。因此,利 用域值法将灰度值在 0 到 5 的像元信息提取出来,如下 图 2,为提取出的变化信息的黑白二值图像,其中白色为 变化信息,黑色为背景信息。与图 1 相比,图 1 主要是 发现和确定变化信息的大致范围和位置,图 2 是对变化 了的信息进行了准确提取。
30° 02'47",东经 104° 56'53"~105° 01'04";全镇幅员面积 2 44.05 km ;典型中亚热带季风气候,气候温暖,雨量丰 沛;位于新华夏构造体系第三沉降带内,境内地质构造 基底差异微弱,较稳定的地块呈整体地貌类型;地貌类 型以丘陵为主,北部深丘山区地带,蕴藏丰富页岩资源。 1.2 资料收集 影像数据为 2001 年 5 月 22 日的 TM(行列号为 128/39) 影像及 2006 年 8 月的 SOPT2 (行列号为 266/289) 多光谱影像,空间分辨率分别为 30 m 和 20 m。其中, SOPT2 的三个单色波段 1 波段(0.5~0.59 μ m) 、2 波段 (0.61~0.68 μ m) 、3 波段(0.79~0.89 μ m)分别类似对 应于 TM2、3、4 波段[1];除影像资料外,还收集了研究 区相关的其他基础资料。
Fig.1
图 1 SPOT 1 波段与 TM 2 之差 Difference of band 1 of SPOT and band 2 of TM
图 3 自动矢量化后的变化信息 Fig.3 Outo-vectorised change information
3.3 土地利用变化信息的自动提取 3.3.1 基于光谱知识的变化信息的自动提取
于本人时间有限,此方法在具体运用中还存在一些问题, 需在以后的工作中进一步完善。
[参 考 文 献] 蒲体信,陈龙华,冯德俊,等.遥感技术在土地利用调查 与监测中的应用[M].成都:电子科技大学出版社,2005: 96-146. 任国业,湛洪举,周华茂,等.利用快鸟卫星影像开展县 级土地利用更新调查方法探讨[J].西南农业学报,2005, 18(1):76. 梅安新,彭望琭,秦其明,等.遥感导论[M].北京:高 等教育出版社,2001:128. freckle. 土 地 覆 盖 / 土 地 利 用 遥 感 监 测 与 变 化 信 息 提 取 . http : //www.gissky.net/blog/user1/1230/archives/2005/7185. html. 尤淑撑,刘顺喜.土地资源遥感监测产业化应用方法探讨 [J].农业工程学报.2004,20(2):264-266. 张继贤, 程 烨. 3S 技术支持的土地利用现状图更新[J]. 中 国土地科学,2002,16(1):20-23. 史培君,宫 鹏,李晓兵,等.土地利用/土地覆盖变化研 究的方法与实践[M].北京:科学出版社,2000:47-51. 莫源富,周立新.TM 数据在土地利用动态监测中的应用 [J].国土资源遥感,2000,44(2):13-15. 杨清华,齐建伟,孙永军.高分辨率卫星妖怪数据在土地 利用动态监测中的应用研究[J].国土资源遥感, 2001 , 50(4):20-27. 王海晖,彭嘉雄,吴 巍.一种多传感器遥感影像的配准 方法[J].2002,30(8):1-3. 王 强,彭嘉雄.基于轮廓的多源图像的配准[J].计算机 与数字工程,2002,30(5):2-6. 陈志伟, 郭志民, 陈永宝. 遥感图像的配准、 校正方法[J]. 福 建水土保持,2001,13(1):54-56. 丁丽霞,周 斌,王人潮.遥感监测中 5 种相对辐射校正 方法研究[J].浙江大学学报,2005,31(3):269-276. 冯德俊,李永树,蔡国林,等.多方法集成的土地利用变 化信息提取[J].成都理工大学学报(自然科学版) ,2005, 32(3):50-52. 张继贤.论土地利用与覆盖变化遥感信息提取技术框架 [J].中国土地科学,2003,17(4):105-109. 张鹏强,余旭初,刘 智,等.多时相遥感图像相对辐射 校正[J].遥感学报,2006,5:339-344. 钟耀武,刘良云,王纪华,等.基于矩匹配算法的山区影像 地形辐射校正方法研究[J].地理与地理信息学.2006,1:31 -34. 黄 永,何 立,王君华.卫星遥感图像的增强处理[J].广 西气象,2006,13:61-63. Sunar F, Musaoglu N. Merging multire solution SPOT Pan LANDSAT TM data: The Effects and Advantages[J]. INT J Remote Sensing,1998,19(2):219-224. 子工业出版社,2003.
第 24 卷 增刊 2 2008 年 9 月
农 业 工 程 学 报 Transactions of the CSAE
Vol. 24 Supp. 2 Sep. 2008
85
采用中分辨率遥感影像提取土地利用变化信息
吴 玺 1,2,王小燕 2※,魏 来 2,李章成 2,卿明福 3,任国业 2
(1.中国农业大学资源与环境学院,北京 100094; 2.四川省农业科学院遥感应用研究所,成都 6100663; 3.四川省土地统征整理事务中心,成都 610041) 摘 要:为了更好的利用中分辨率遥感影像监测土地利用变化信息,为今后土地利用遥感动态监测提供更准确快捷的方 法,该文介绍了一种基于光谱知识自动提取中分辨率影像中土地利用变化的方法,以岳池县花园镇为例,采用 2001 年 TM 和 2006 年 SPOT2 影像进行对比分析、识别、运算,然后提取。经分析提取精度较高,达 98.2%,并得出该方法的 适用范围。 关键词:土地利用;变化信息;图像差值;光谱知识;遥感影像 中图分类号:TP753 文献标识码:A 文章编号:1002-6819(2008)-(Supp.2)-0085-04 吴 玺, 王小燕, 魏 来, 等. 采用中分辨率遥感影像提取土地利用变化信息[J]. 农业工程学报, 2008, 24(Supp.2):85-88. Wu Xi, Wang Xiaoyan, Wei Lai, et al. Extraction of land use changing information using middle resolution remote sensing images[J]. Transactions of the CSAE, 2008,24(Supp.2):85-88.(in Chinese with English abstract)
由于同谱异物、同物异谱,以及两时相影像的月份 不一致等问题,即使经过相对辐射校正也会会导致变化
增刊

玺等:采用中分辨率遥感影像提取土地利用变化信息
87
了的信息未被提取、没变化的信息被提取、变化信息的 边界不准确等问题,这需要在室内采取必要措施尽量减 少这些问题的出现,提高数据精度,减少野外工作量。 将矢量图与用其他自动识别方法得到的变化信息图 以及影像在遥感影像处理软件中进行叠加显示。用手工 删除未发生变化的图斑,用绘图工具勾绘发生了变化而 未被提取出边界的图斑,同时,还需对一些变化图斑边 界不一致的地方进行编辑。如图 4 是手工去除干扰信息 后的变化信息与影像的叠加显示。
相关文档
最新文档