1、第 卷 第 期 年 月 北京测绘 引文格式：刘梁美子基于 的番禺区地表覆盖时空变化分析北京测绘，（）：，（）：收稿日期 作者简介 刘梁美子（），女，湖北黄冈人，硕士，研究方向为 和 应用技术研究。：基于 的番禺区地表覆盖时空变化分析刘梁美子（广州市城市规划勘测设计研究院，广东 广州）摘 要 本文基于全球 地表覆盖数据（），采用地表覆盖转移矩阵、重心迁移方法，分析了番禺区 年不同地表覆盖类型的相互转化及其空间分布。结果表明：过去 年间，人造地表持续占用耕地、森林和草地，其中耕地最为明显。同时，年人造地表扩张速度较 年更为迅速。从相互转化的空间分布上看，番禺区人造地表在 年呈点状增加，之后 快速蔓
2、延至全区。人造地表和水体的重心向东迁移，耕地的重心转移则逐步西移，呈“一进一退”的趋势。关键词；地表覆盖变化；重心迁移；番禺区中图分类号 文献标识码 文章编号（）引言土地利用 覆盖变化（，）是地表综合集成研究的一个学科基础和前沿领域。遥感技术已成为对地表覆盖变化动态监测的重要手段，各类 基础数据集层出不穷，典型的主要包括美国地质调查局（，）的 数据集（，）、欧盟联合研究中心（，）的 土地覆盖数据集（，）、美国波士顿大学（，）的 数据集（，）等。全球 地表覆盖数据产品（）是我国自主研制的首套 分辨率、全球尺度的地表覆盖产品，在生态、环境、气候等研究领域均得到广泛应用，数据精度亦得到验证。近年来，
3、土地利用 覆盖变化相关研究中，国内外均关注运用遥感和地理信息系统技术进行 动态监测、对生态系统服务影响等方面，时空变化、驱动力分析、模型模拟，以及 对径流的影响也一直是国内关注的焦点。相应地，国内研究热点区域亦多集中于生态脆弱区、流域和人类活动频繁地区（如珠三角、长三角等）。自 年广州市政府撤销番禺市设立番禺区后，广州城市空间逐渐向南发展。与此同时，人口集聚以及快速城市化带来的土地资源矛盾日益凸显，建设用地供需失衡。然而，目前涉及番禺区土地 地表覆盖相关的研究大多聚焦于较大空间尺度，如粤港澳大湾区，珠三角城市群，广东省，广州市，但总体上看，基于小尺度的、针对番禺区的研究很少，更鲜有长时间序列的
4、时空尺度的分析，这对于区（县）域的政府了解本地生态、资源等问题，制定相关政策，实现区域可持续发展是不利的。鉴于此，本文基于 、年三期地表覆盖数据，在番禺区开展地表覆第 卷 第 期刘梁美子基于 的番禺区地表覆盖时空变化分析盖时空演变特征的研究，以期为未来番禺区国土空间规划和城市治理提供科学依据。研究区域和研究方法 研究区概况广州市番禺区地处广东省中南部的珠三角中心地带（、），东邻东莞，西邻佛山。境内地势平坦，河网众多，东江自东、北部入境，上游来水及本区水系均归依珠江三大口门即虎门、蕉门、洪奇门出海，见图。市莲路以北为市桥台地，南部是冲积平原，最高峰大夫山海拔.。作为广州市“南拓”发展战略重点区域
5、，近年来社会经济发展迅速。石楼镇化龙镇石碁镇南村镇沙湾街道石壁街道大龙街道洛浦街道钟村街道大石街道新造镇沙头街道桥南街道小谷围街道东环街道市桥街道北江干流水道江珠珠江陈村水道大夫山番禺区镇街河流水系图 研究区地理位置 研究方法 数据获取本文使用的番禺区、年地表覆盖数据来自 官方网站（：？）。数据共包括 个类型，分别是耕地、森林、草地、灌木地、湿地、水体、苔原、人造地表、裸地、冰川和永久积雪，数据空间分辨率为，官方公布总体精度在 以上。本文所有数据的统一坐标为 坐标系，投影方式为 投影。地表覆盖转移矩阵地表覆盖转移矩阵来源于土地利用转移矩阵，能够刻画区域地表覆盖变化的结构特征，反映人类活动引导的
6、地表覆盖变化方向。该方法来源于系统分析中对系统状态与状态转移的定量描述，反映在一定时间间隔下，一个亚稳定系统从 时刻向 时刻状态转化的过程，数学表达形式为|（）式中，为研究期初与研究期末的地表覆盖状态；为地表覆盖的类型数。地表覆盖转移矩阵中向量可以是地表覆盖类型面积，也可以是期初地表覆盖类型向期末地表覆盖类型转化的概率。本文采用地表覆盖类型面积进行分析。利用 空间分析功能逐栅格计算并统计不同地表覆盖类型相互转换的面积，得到研究区的地表覆盖转移矩阵。地表覆盖重心迁移地表覆盖重心迁移借鉴人口重心分布变化研究。通过 软件的空间统计工具箱，可以获得各年份不同地表覆盖类型的重心，通过分析各覆盖类型的重心
7、迁移规律，可以揭示土地和资源空间格局变化，反映时空变化过程和速率。计算公式如式（）、式（）所示。k （kk）k（）Yk （kYk）k（）式中，k为第 k 类地表覆盖类型的重心 坐标；Yk为在第 k 类地表覆盖类型的重心 Y 坐标；k为第k 类地表覆盖类型中第 个栅格面积；k为第 k类地表覆盖类型中第 个栅格 坐标的平均值；Yk为第 k 类地表覆盖类型中第 个栅格的 Y 坐标的平均值；k为第 k 类地表覆盖类型总面积。结果分析 地表覆盖变化分析不同时间地表覆盖类型能够指示区域地表覆盖基 本 情 况 及 其 数 量 变 化 趋 势。番 禺 区 年不同地表覆盖类型情况如表 所示。年，番禺区主要地表覆
8、盖类型以耕地（）、人造地表和水体（均）为主，北京测绘第 卷 第 期这 年间，各地表覆盖类型变化较为平稳，耕地减少相对明显，年均降幅为.，森林、水体和人造地表均小幅增加，年均变化率均未超过。至 年，人造地表占比由 年的.增至.，耕地、森林和草地占比则分别从.降 至.，.降 至.，.降至.，年均变化率均大于，其中草地降幅.，人造地表增幅达.。这表明 年 年间，番禺区城市扩张速度较 年开始加快，表现为人造地表面积急剧增加，耕地、森林和草地不同程度被挤占。年以来，水体面积持续平稳增加，这一方面与农业结构调整和居民生活生产方式改变有关，另一方面，从 年开始，番禺区大力开展水系的综合整治工作，陆续编制了番
9、禺区生态水系整体规划蓝图番禺区水环境治理进行“大兵团作战”工作方案等政策文件，极大地改善了 世纪八九十年代至 年初水体恶化的情况。表 番禺区 年不同地表覆盖类型情况土地利用类型 年 年 年 年 年 年面积 占比 面积 占比 变化量 年均变化率 面积 占比 面积 占比 变化量 年均变化率 耕地.森林.草地.水体.人造地表.地表覆盖转移分析地表覆盖的变化一定程度可反映人类活动对土地的改造情况。表、表 中的对角线表示未变化的地表覆盖类型，横向反映某一地表覆盖类型转换为其他类型的面积，纵向反映其他地表覆盖类型转变为某一类型的面积。通过分析 年番禺区地表覆盖转移矩阵，结果显示：（）年间，减少的耕地面积为
10、.，其中.转化为了水体；增加的耕地面积为.，.来自水体；森林转化为其他覆盖类型的面积为.，主要转化为人造地表.，其他地表覆盖类型转化为森林的面积为.，来自耕地（.）、水体（.）和人造地表（.）。草地面积略有减少，表现为与人造地表和水体间的相互转化。水体与耕地间的转化较为频繁，减少的面积（.）中.转为耕地，增加面积为.，.来源于耕地。人造地表有.转出，转化为少部分森林（.）和草地（.），同时有.转入，主要来自耕地（.）和草地（.）。（）年间，耕地转化为其他地表覆盖类型的面积为.，其中.转化为人造地表，.转化为水体，从其他地表覆盖类型转化为耕地的面积仅为.，主要来源于水体（.）。与耕地大量转为人造
11、地表相似，森林和草地减少面积分别为.、.，其中分别.、.转向人造地表，两者增加面积仅为.和.。水体面积转移情况与 年大体一致，主要体现为与耕地间的转化，不同的是，人造地表对水体的占用更为明显，水体减少面积中.转为人造地表，但由于耕地等转化为水体面积较多，使得水体总面积持续增加。人造地表的转入远大于转出，仅.转出，转化为森林、耕地和草地，相应地，转入面积则达.，来自耕地（.）、森林（.）和草地（.），约一半人造地表由耕地转化而来。结果表明，番禺第 卷 第 期刘梁美子基于 的番禺区地表覆盖时空变化分析区城市扩张对耕地的占用十分明显，这也是 年地表覆盖相互转化较 年差异最为明显的部分，人造地表净增加
12、面积为上一阶段（年）的 倍，说明过去 间，人类建设占用耕地，森林和草地，农用地流失的速度加快，城市化的进程持续加速。表 年番禺区不同地表覆盖类型转移矩阵单位：土地类型耕地森林草地水体人造地表总计耕地.森林.草地.水体.人造地表.总计.表 年番禺区不同地表覆盖类型转移矩阵单位：土地类型耕地森林草地水体人造地表总计耕地.森林.草地.水体.人造地表.总计.地表覆盖变化空间格局图 和图 描述了不同地表覆盖类型相互转化的空间分布，结果显示，番禺区人造地表由 年的点状增加，之后快速蔓延，至 年基本填充全区。耕地向人造地表的转化集中在东北（化龙镇、新造镇），西北（石壁街道、洛浦街道）和南部（石碁镇、桥南街道
13、），森林、草地转向人造地表则发生在中部、西部和东部，主要包括南村镇、石楼镇、沙湾街道等。人造地表的扩张方向以沿江、河涌延伸模式为主。区域内各地表覆盖类型的相互转化导致各自重心的变化（图）。从不同地表覆盖类型重心的转移方向上看，年间，森林重心迁移不明显，人造地表和水体的重心向东迁移，与此同时，耕地的重心转移方向与之相反，逐步西移，形成“一进一退”的趋势。根据番禺区城乡发展规划（年），番禺区以“两轴两带四中心”作为引导城市空间发展的基本框架。其中，化龙镇汽车城以及位于石楼镇的广州亚运城分别属于“四中心”的国际创新服务中心和亚运城服务中心，由此可能导致番禺区人造地表重心东迁。同时，退耕、围海等人为活
14、动的进行也使得水体重心同步向东迁移。相应地，由于人造地表和水体均占用耕地，尤其是城镇用地的扩张，造成耕地重心呈现向西移动的趋势。北京测绘第 卷 第 期石楼镇化龙镇石碁镇南村镇沙湾街道石壁街道大龙街道洛浦街道钟村街道大石街道新造镇沙头街道桥南街道小谷围街道东环街道市桥街道石楼镇化龙镇石碁镇南村镇沙湾街道石壁街道大龙街道洛浦街道钟村街道大石街道新造镇沙头街道桥南街道小谷围街道东环街道市桥街道石楼镇化龙镇石碁镇南村镇沙湾街道石壁街道大龙街道洛浦街道钟村街道大石街道新造镇沙头街道桥南街道小谷围街道东环街道市桥街道石楼镇化龙镇石碁镇南村镇沙湾街道石壁街道大龙街道洛浦街道钟村街道大石街道新造镇沙头街道桥南
15、街道小谷围街道东环街道市桥街道2330N22570N石楼镇化龙镇石碁镇南村镇沙湾街道石壁街道大龙街道洛浦街道钟村街道大石街道新造镇沙头街道桥南街道小谷围街道东环街道市桥街道2330N22570N森林不变森林水体森林草地森林人造地表森林耕地草地不变草地森林草地耕地草地人造地表草地水体水体不变水体森林水体耕地水体人造地表人造地表不变人造地表水体人造地表草地人造地表森林人造地表耕地耕地不变耕地森林耕地水体耕地人造地表(e)人造地表(d)水体(c)草地(a)耕地(b)森林图 年不同地表覆盖类型相互转化空间分布石楼镇化龙镇石碁镇南村镇沙湾街道石壁街道大龙街道洛浦街道钟村街道大石街道新造镇沙头街道桥南街道
16、小谷围街道东环街道市桥街道石楼镇化龙镇石碁镇南村镇沙湾街道石壁街道大龙街道洛浦街道钟村街道大石街道新造镇沙头街道桥南街道小谷围街道东环街道市桥街道石楼镇化龙镇石碁镇南村镇沙湾街道石壁街道大龙街道洛浦街道钟村街道大石街道新造镇沙头街道桥南街道小谷围街道东环街道市桥街道石楼镇化龙镇石碁镇南村镇沙湾街道石壁街道大龙街道洛浦街道钟村街道大石街道新造镇沙头街道桥南街道小谷围街道东环街道市桥街道石楼镇化龙镇石碁镇南村镇沙湾街道石壁街道大龙街道洛浦街道钟村街道大石街道新造镇沙头街道桥南街道小谷围街道东环街道市桥街道森林不变森林水体森林草地森林人造地表森林耕地草地不变草地森林草地耕地草地人造地表草地水体水体不变水体森林水体耕地水体人造地表人造地表不变人造地表水体人造地表草地人造地表森林人造地表耕地耕地不变耕地森林耕地水体耕地人造地表(a)耕地(b)森林(c)草地(d)水体(e)人造地表图 年不同地表覆盖类型相互转化空间分布第 卷 第 期刘梁美子基于 的番禺区地表覆盖时空变化分析石楼镇化龙镇石碁镇南村镇沙湾街道石壁街道大龙街道洛浦街道钟村街道大石街道新造镇沙头街道桥南街道小谷围街道东环街道市桥街道石楼镇
