区域土地利用时空变化图谱分析.docVIP

区域土地利用时空变化图谱分析 摘要：运用地学信息图谱理论和方法，在RS和GIS的支 持下，以海河流域40a 土地利用5期空间数据为基础，合成 了 40a变化过程图谱，并通过图谱重构，提取了海河流域40a 土地利用变化的“涨势”系列图谱、40a “涨势”变化过程 谱和40a 土地利用时空演变模式图谱。通过图谱分析，逐 一分析了不同时序单元、不同空间尺度土地利用格局的变化 规律，并总结了 40a 土地利用的分异规律。 关键词：土地利用图谱；图谱分析；海河流域 研究区概况 研究区位于东经1120?1190、北纬350?420之间，西 与山西高原与黄河流域接界，北与蒙古高原与内陆河流域接 界，南界黄河，东临渤海。包括河北、河南、山西、内蒙古、 北京和天津四省两市，流域面积23.2万平方公里，占全国 总面积的2. 4%。据1993年统计，海河流域内人口 1.176亿 人，约占全国总人口的10%;耕地1.63亿亩，占全国耕地的 11%o海河流域处于中国干旱和湿润气候的过渡地带，多年 年平均降水量为548mm。流域多年平均水资源总量为419亿 m3,产水系数（水资源总量与降水量之比）为24%。人均水 资源占有量仅有338m3,不足全国的1/6。流域西北高，东 南低，主要地貌类型为平原、山地、丘陵和台地。海拔最高 处3059m，最低地方仅为0m。 数据来源及处理 地学信息图谱单元的提取需要有具有一定时间间隔的 土地利用状况数据作为基础。本文以海河流域生态环境治理 工程项目，海河流域5期）土地利用数据为基础，参照2000 年矢量化底图，利用GIS和RS软件进行地图矢量化。对不 同期影像进行配准，保证误差不超过3个像元，参照不同期 遥感影像，对照统计数据，采用人机交互的方法对遥感影像 进行分类，最终得到5期土地利用数据。利用野外调研数据 进行检验，并采用误差矩阵进行评价，精度都达到90%以上。 在已有各期土地利用数据的基础上，对土地利用类型进 行重编码，以获取能最大反映土地利用变化的基本地物类 别。为了便于进行空间数据之间的地图代数运算和编码的提 取、综合以及图谱重构，以期能更好的在空间上反映各土地 利用类型的时序变化状况。 分析结果 3. 1 1970?2008年土地利用变化过程图谱分析 40a 土地利用变化过程图谱是由40a的5个采样时刻的 土地利用类型数据共同合成，图谱共产生4301类具有5位 编码的图谱单元。其中，461类图谱单元占到变化面积96%。 40年来，利用方式没有发生变化的耕地主要分布在海河流域 地势平坦的东部平原和地势较高的西部和南部平原。发生变 化的耕地主要集中在西北地势较高的平原、沿海以及大型城 市的周边地区。尤其是沿海地区耕地次生盐碱化，耕地和水 体湿地的频繁转换，以及城市建设，导致沿海地区的土地利 用类型变化剧烈。海河流域呈从西部、北部耕牧交错区到中 西部山区，经过草地过渡带到东部、南部的耕地、盐碱地、 水域直到海域的土地利用格局。在这种变化格局中，耕地占 据着优势，在北部山区的耕牧交错区，耕牧之间的转化最为 频繁。应该说，农牧交错区、林草交错区及沿海地区的土地 利用变化频繁，生态脆弱。 3. 2 1970?2008年土地利用变化“涨势”系列图谱分析 根据土地利用的转入状况，分析涨势系列图谱可知， 1970?1980年间，变化面积最大的图谱单元类型为新增林 地、园地，主要分布在海河流域中部太行山山脉一线，在流 域北部和南部山区也有零星分布。第二大变化图谱单元类型 为新增草地，主要分布在流域西北部到西南部山脉沿线。草 地增加大部分来源于林地逆向演替，也有一部分来源于耕地 转化。此外，新增耕地成为面积位居第三的土地利用变化类 型，集中分布在流域中南部平原向山区过度地带，在北部山 区及沿海也有增加，分布零散。1980?1990年间，土地利用 变化面积最大的是新增草地，草地增加主要集中在流域北部 张家口、赤诚和怀来一带，另外，沿太行山山脉一线也有零 星分布。与草地增长方式不同，林地增长多为成片增加，主 要分布在流域北部山区，平原上也有零星增加。1990?2000 年，草地增加趋势明显，集中分布在流域西北部山区，呈块 状分布。另外，林地也有明显增加，山区增加林地多成片状， 平原上多为零星散布。20⑻?2008年，林地增加占据主要地 位，主要来源于北部山区退耕还林以及平原区果园增加。第 二大变化图谱单元类型为新增草地，主要分布在流域西北部 山区，多为退耕还草，另外其它新增草地多为块状分散分布， 由林地采伐逆向演替为草地所致。 3. 3海河流域40年土地利用模式图谱特征分析 为了从宏观上了解海河流域土地利用转变方式时空特 征，根据图谱单元变化特点及其稳定性，归纳了 18类土地 利用时空演变模式特征。这样，同一图谱单元记录了相同空