基于GIS的河网分形研究.pdfVIP

西北大学学报(自然科学版) 2OO6年8月，第 36卷第4期，Aug．，2006，Vo1．36，No．4 Journal of Northwest University(Natural Science Edition) 基于 GIS的河网分形研究 张宏才 ，汤国安 (1．西北大学 新闻传播学院，陕西 西安 710069；2．南京师范大学 地理科学学院，江苏 南京 210097) 摘要：目的 探讨在 GIS支持下河网分维数的计算方法以及不同分维的意义。方法 应用数字高 程模型，在 GIS支持下提取陕北黄土高原韭园沟流域的河网，利用空间分析方法、分形理论以及统 计学知识进行计算、模拟。结果 基于霍顿水系定律得到韭园沟流域河网分维数是 1．662，基于分 形基本定义得出的是 1．589。结论 ①利用GIS计算河网分维数与传统方法相比，具有高效、简 捷、检验性好的优点；②对同一河网，使用不同方法得到的分维数不同，表明河网分形性质丰富多 彩；③分形维数对比时一般应满足计算方法一致，且在同一个无标度域内。 关 键 词：分形；DEM；GIS；河网；分维数 中图分类号：K343．1 文献标识码：A 文章编号：1000-274X(2006)04-0659-04 分形理论是人们认识 自然界的崭新工具，它为 人们提供了一条从有限到无限，从部分到整体，与系 统论截然相反的认识途径，使人们对自然界中许多 大尺度、大范围事物的认识在理论上成为可能。 目前，河网分形的研究基于两种思路：一种是将 水系作为整体，探讨河 网分形 ?，--13]；另一种是探 讨河网参数之间或水系参数的分形特征，如河道长 度与集水面积、流域周长与集水面积、河道纵剖面 等L3 J。其具 体 方法 一 是 基 于 霍 顿水 系定 理l ’．9．1卜 J，第 二 是 基 于 分 形 的 基 本 定 义u ’加川 。有关数据的获取，或实际测量，或在 一 定比例尺地形图上采用改变步长、手工量测获 得l ]，或依据别人有关的数据统计取得[3I6I9]。总 的来说 ，数据获取费时费力，计算周期长，且在同一 次量测过程中很难保证步长的一致性，同时也无法 满足测量标度趋于 0的条件，数据精度差，直接影响 计算结果的准确性，而且计算方法可检验性差，严重 影响了分形理论在水科学中的进一步研究与应用。 本文探讨在GIS支持下计算河网分维值的方法。 园沟流域，位于东经 110。15 00”～110。22 45”，北纬 37。32 00”～37。37 30”，总面积约 70．13 km 。该区 属典型的黄土丘陵沟壑区，地表形态复杂，沟谷纵 横，以发育充分的梁峁地貌为主，平均海拔高程 980．06 m，沟壑密度7．18 kin／kin ，平均地面粗糙度 为 1．18，水土流失严重。韭园沟流域是黄河水利委 员会水土流失重点治理实验区ll 。 1．2 实验数据 实验数据为 1：10 000 DEM数据。1：10 000 DEM数据是将 1：10 000比例尺地形图等高线扫描 数字化，再内插后生成。制作过程严格依照国家测 绘局关于建立数字高程模型的标准来进行。 1．3 应用软件 GeoScan，DEMix 2．0主要用于建立流域的 DEM 数据，Arc／INFO 8．0，ArcView 3．2用于 DEM 的建 立 ，以及后续分析，Excel是统计、模拟、计算的主要 软件。 2 实验原理与方法 1 实验基础 2．1 基于分形定义的河网分维数 1．1 实验样区 用一定尺度的方格网去覆盖分形图(曲线)时， 研究选择的实验样区是陕西榆林地区绥德县韭 非空方格数(含有图形)N(a)与方格的边长a具有 收稿日期：2004．10．15 基金项目：国家自然科学基金资助项 目40301027) 作者简介：张宏才(1970一)，男，陕西宝鸡人，西北大学讲师，硕士，从事地理信息系统方面研究。 一 660一 西北大学学报(自然科学版) 第 36卷 一 一 对应关系。若有 个不同的 a，(，=1，2，3， ? )，则有 肘 个不 同的 Ⅳ(口)值 ，即有：?r(a)～ a 一 。利用最小二乘法对 ln(Ⅳ)和In(口)作一元线性 回归便可得到分形维数D的估计值。利用此种方法 得到的分形维数叫盒维数(见图 1)。盒子可以取为 小区间、正方形、立方体、圆或球体等。 在 GIS中，矢量数据转换为栅格数据的过程，与 用某一尺寸的盒子去覆盖图形(曲线)是一致的。 通过 GIS的属性表查询栅格数据的数 目，即为非空 栅格的数 目。变换矢一栅转化时的 Output Grid Cel 图 1 盒维数计算原理图 Fig．1 Principle of calculating box dimension