基于gis的土地利用动态变化空间融合方法.docxVIP

基于gis的土地利用动态变化空间融合方法 土地使用和屋顶变化（luc）一直是全球环境变化研究的热点和前沿，研究单元大小。李晓文等研究的单元是研究区所辖的较高级别行政区, 通过各行政区的土地利用综合动态度分析不同行政区土地利用数量的变化;而刘纪远等[10,11,12,13,14,15,16,17,18,19,20,21,22,23,24]研究的单元是研究区所辖的整个行政区, 通过综合得到的土地利用综合动态度表征该区整体的土地利用变化。这两种方式均以行政区为单位进行分析, 虽能清楚描述研究区域的土地利用变化, 但在一定程度上还是损失掉一些细节信息 (即更精细尺度上的变化信息) , 如研究区内较低级别行政区的土地利用变化。本文拟以2001年、2004年的土地利用数据为基础, 研究房山区最低级别的行政区 (即村镇) 的综合动态度, 最大限度保留精细尺度上的土地利用变化信息, 然后按精细尺度上对应的动态度空间分布, 在粗略尺度上突破行政区划的界限, 综合得到宏观层次上的土地利用变化程度空间分布。这种跨行政区综合的方法在尽可能少丢失信息的情况下, 能更好地反映土地利用变化程度的空间分布。 1 基层社会空间分布 北京市房山区位于39°30′~39°55′N、115°25′~116°15′E, 处于华北平原与太行山交界地段, 西部和北部是山地、丘陵, 约占全区总面积的2/3, 南部为沃野平原。房山区总面积2 019km2, 辖8个街道、11个镇、4个地区、6个乡。根据已有的2001年、2004年房山区土地利用矢量数据, 其属性字段只有相关的二级地类编码, 为了便于后续的研究, 根据《地类编码参照表》将其分为耕地、草地、林地、建设用地、水域、未利用地六大类, 从而得到2001年、2004年房山区六大地类的空间分布图。房山区的东部区域以林地和未利用地为主, 西部区域以耕地和建设用地为主 (图1, 见封3) 。 2 土地使用与行政区域相关的空间整合方法 2.1 土地利用动的数量动态分析方法 本文选用的土地数量变化模型是土地利用变化分析中的基本模型之一, 土地资源数量变化可以用土地利用动态度表示。本文的土地利用动态度即综合土地利用动态度, 是指某一地区某一时段内综合的土地利用类型的数量变化情况, 其表达式为: 式中:LUi为监测起始时刻第i类土地利用类型的面积;ΔLUi-j为监测时段第i类土地利用类型转为第j类土地利用类型面积的绝对值;T为监测时段长度。 由上述模型可知:土地利用动态度的描述对象是某一区域。传统的方法一般是以行政或自定义的网格为单位开展研究, 而本文研究的区域有473个村镇, 数量较多, 且网格化不利于表现土地利用空间变化的自然分布规律。村镇作为最精细的统计单位, 能够最大限度地保留土地利用变化程度的细节信息, 在这个基础上能更好地进行后续的研究, 因此, 本文以村镇为对象计算综合动态度。 2.2 基于综合动态度下的多态性分析 为便于土地利用动态度在空间分布上进行融合, 需要对综合动态度进行分类。村镇总数为473个, 为避免因组数过多加大分析的难度和复杂性以及因组数过少造成融合后得到的粗略尺度上空间格局不够客观, 本文按照各村综合动态度值分布区间的特点 (图2) , 并保证不同类型间的地块数大体相同, 将综合动态度分为六类, 并对每一类赋值 (表1) 。 2.3 土地利用变化程度的指标综合动态度 每个斑块都有一个代表该区域综合动态度所属类型的值, 综合考虑其空间临近性及属性临近性, 跨越行政区的限制, 将类型值相对比较接近且空间临近的图斑融合, 尽可能使变化程度相近的斑块融合到同一区块, 融合出一个粗略尺度上的跨行政区的空间变化格局分布 (图3) 。具体的操作方法是:斑块的类型值为1、2、3、4、5、6, 相邻类型值间距为1, 类型值2和4的间距为2, 表示属性值较为接近, 类型值2和6间距为4, 表示属性值相差很大。在后期的操作中, 对于每个斑块, 将其类型值与其一阶临近斑块的类型值相比较, 并尽量将该斑块融合到与其公共边长最大的类型值间距较小的图斑中。如果处理斑块的类型值很大 (如6) , 而其周围斑块的类型值都很小 (如1) , 基于提取粗略尺度上土地利用变化程度分布格局的目的, 直接将该处理斑块融合到周围的斑块中。以此类推, 得到综合动态度跨行政区的格局分布, 这样就能够在粗略尺度上反映土地利用变化程度的分布趋势。 提取出的各区块, 其组成斑块基本保持属性相近的特征, 能够反映该区块土地利用程度变化的信息。如果某一区块对应的原图斑的综合动态度多数小于0.5%, 则说明该区块的土地利用非常稳定, 如果大多数介于0.5%~1%, 则说明该区块的土地利用不稳定;如果介于1%~4%, 则说明该区块的土地利用变化较为剧烈;如果大于4%