基于IDL和Mapinfo的子流域划分与信息提取方法

基于IDL和Mapinfo的子流域划分与信息提取方法
【专利摘要】本发明公开一种基于IDL和Mapinfo的子流域划分与信息提取方法，属于地球物理科学门类的水文科学类。本发明是一种通过IDL程序语言实现利用高精度的DEM资料对流域进行升尺度的子流域划分方法，很好地利用了IDL面向矩阵，运算效率高的特点。通过对DEM的分析，获得子流域划分结构、子流域间拓扑关系、子流域控制点、控制范围、集水面积、子流域间河道长度、坡度等地质、地貌特征，为流域径流模拟提供支撑。同时，利用强大的地理信息系统软件MapInfo对子流域划分结果显示。基于IDL生成了子流域拓扑关系、子流域控制点位置，利用MapInfo的编程语言MapBasic展示子流域的拓扑关系。
【专利说明】基于IDL和Mapinfo的子流域划分与信息提取方法
[0001]
【技术领域】
[0002]本发明涉及一种利用编程语言IDL和MapInfo进行子流域划分、分析、显示的方法，属于地球物理下水文分支【技术领域】。
【背景技术】
[0003]分布式水文模型是预测流域内水文情势的空间分布，评价人类活动(水资源开发利用、土地利用、面源污染等)的影响，为流域规划与治理、水资源可持续开发利用科学指导的有效方法。基于EffiM的分布式水文模型是水文模型研究的热点问题。基于EffiM的流域内子单元的划分对模拟精度有重要的影响，根据需要进行有效的流域划分具有重要意义。传统的子单元划分多是基于格网的划分，精度低，并不能满足现实需要。日本学者Yamazaki等于2009年提出基于格网的自然子流域划分，自然子流域划分具有单元内和单元间水文过程清晰的特点，子流域作为独立单元，容易进行产汇流计算，很多模型选择自然子流域划分，因此研究子流域划分及地形、地貌信息对流域模拟很重要。子流域划分最基本的方法是利用ArcGIS进行子流域划分，但随着人们对子流域划分结果的要求越来越高，单独依据ArcGIS进行的子流域划分并不能满足应用要求。本发明利用IDL和Mapinfo相结合的技术结合日本学者提出的方法进行子流域划分。

【发明内容】

[0004]发明目的:考虑到以往子流域划分不灵活，提取信息不完善等缺陷，本发明提供一种基于IDL和Mapinfo的子流域划分与信息提取方法，利用面向矩阵的IDL编程语言进行子流域划分，并通过MapInfo对划分结果进行展示，同时，利用MapInfo辅助语言MapBasic显示子流域拓扑结构。
[0005]技术方案:一种基于IDL和Mapinfo的子流域划分与信息提取方法利用面向矩阵运算的程序语言IDL对高精度DEM进行提取分析，获得流域的子流域划分、子流域拓扑关系等指标，利用MapInfo的MapBasic编程语言读取子流域拓扑信息、出口点信息，对子流域划分结果进行显示。具体包括如下步骤:
第一步:从中国科学院数据应用环境网站上下载研究区域的高精度DEM数据((30m或90m))，利用ArcGIS对研究区数据进行合并处理，基于研究流域边界提取研究流域内的DEM。利用Hydrology 模块的填挂功能生成研究区无挂地DEM,利用Hydrology模块的生成流向、汇流累积量功能生成研究区的流向、汇流累积量。利用ArcGIS中栅格转ASCII工具，将流域DEM、流向、汇流累积量分别导出为ASCII码，供IDL处理。
[0006]第二步:编写生成子流域控制点及拓扑关系的程序，按照网格划分，利用IDL的where函数获得网格内汇流累积量最大的点的位置，作为子流域控制点，最终生成每个大网格的子流域控制点位置文件。利用控制点沿流向向下游搜索，并记录搜索路径，获得该网格控制点的下游控制点，记录控制点之间的拓扑关系，即为子流域间拓扑关系，生成拓扑关系文件。利用控制点位置文件及DEM信息获得控制点高程，作为子流域的代表高程。利用控制点的搜索路径，计算上游子流域到下游子流域间的河道长度。利用控制点间河道长度和控制点高程差，得到河道坡度信息。
[0007]第三步:编写子流域生成程序。首先，利用子流域拓扑关系生成子流域的上游子流域，将子流域控制点的汇流累积量将去上游子流域的汇流累积量，获得子流域的汇流面积。其次，对子流域集水面积进行由小到大排序，利用IDL编写基于矩阵的流向追踪程序，依次对排序单元进行流向追踪，创建初始矩阵，判断流向流入控制点的网格，将子流域控制区内的网格赋值为子流域编号。最后，得到每个子流域的控制区域，生成子流域控制区域文件。
[0008]第四步:利用第三步中生成的子流域控制矩阵文件，按照各子流域控制区的边界，生成子流域的mif文件格式,方便MapInfo导入。
[0009]第五步:打开MapInfo,利用工具----表----1mport工具将mif文件导入到
MapInfo中，生成流域的子流域划分tab文件。
[0010]第六步:编写MapBasic程序，读取子流域拓扑关系文件、子流域控制点位置文件，生成子流域的拓扑关系图。与第六步的子流域划分图叠加在一起，就生成了流域的拓扑关系图。
[0011]有益效果:与现有技术相比，本发明提供的基于IDL和Mapinfo的子流域划分与信息提取方法，利用面向矩阵的IDL编程语言进行子流域划分，并通过MapInfo对划分结果进行展示；编程语言IDL简单、方便、计算效率高，容易被研究者所掌握，研究者可根据自己的需要进行不同的子流域划分。同时，本发明利用高精度DEM提取了子流域的有效信息，对径流模拟提供有效支撑。利用Ma·pInfo对子流域化法结果进行展示、出图，方便研究者使用。
【专利附图】

【附图说明】
[0012]图1为本发明实施例的研究区王家坝站以上流域边界；
图2为本发明实施例的王家坝站以上流域子流域划分；
图3为本发明实施例的王家坝站以上流域子流域拓扑关系展示；
图4为本发明实施例的研究区外洲站以上流域边界；
图5为本发明实施例的外洲站以上流域子流域划分；
图6为本发明实施例的外洲站以上流域子流域拓扑关系展示。
【具体实施方式】
[0013]下面结合具体实施例，进一步阐明本发明，应理解这些实施例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
[0014]I软件安装 Cl) IDL 安装:
第一步:下载IDL7.0，解压缩到本地文件。
[0015]第二步:双击解压缩文件里面的安装文件idl70winx86_32.exe,点击下一步。
[0016]第三步:同意IDL协议，点击下一步。[0017]第四步:选择安装路径或默认安装路径，点击下一步。
[0018]第五步:选择IDL扩展模块，根据需要选择安装。点击下一步，开始安装。
[0019]第六步:提示是否要导入License,选择是，进入License导入界面,选择第三项，从本机导入License。点击next,选择本机上IDL安装包里面的License,导入License。注意:Li cense里面的计算机名要改为本机名。
[0020]第七步:从开始菜单里面，启动IDL License管理器，选择Start/Stop/Reread,先点击 Stop 停止 License,接着 Start License,点击 Reread License,在下方出现 ServerStart Successful。表不启动成果。
[0021]第八步:点击IDL启动界面，启动IDL。
[0022](2) Maplnfo安装:第一步:下载正版MapInfol0.0版本软件,解压缩到本地文件。
[0023]第二步:双击解压缩文件里面的安装文件setup, exe,点击下一步。
[0024]第三步:输入安装包里面的序列号和注册码，点击下一步。
[0025]第四步:选择安装路径，点击下一步，完成安装。
[0026]第五步:在完成安装后，将安装包里面的micore.dll覆盖到安装目录，完成安装。
[0027]第六步:启动Maplnfo。
[0028](3) MapBasic 安装:
第一步:下载MapBasic软件，将其解压到本地文件夹。
[0029]第二步:双击SETUP.EXE，进入安装页面，点击下一步。
[0030]第三步:进入软件设置界面，输入MapBasic的注册码，点击下一步。
[0031]第四步:设置软件安装路径，最后完成安装。
[0032]2具体实施步骤:
第一步:从中国科学院数据应用环境网站上下载研究区域的高精度DEM数据(90m)，利用ArcGIS对研究区数据进行合并处理，基于王家坝流域边界提取研究流域内的DEM。利用Hydrology模块的填挂功能生成研究区无挂地DEM,利用Hydrology模块的生成流向、汇流累积量功能生成研究区的流向、汇流累积量。利用ArcGIS中栅格转ASCII工具，将流域DEM、流向、汇流累积量分别导出为ASCII 5马，供IDL处理。
[0033]第二步:编写生成子流域控制点及拓扑关系的程序，按照网格划分，利用IDL的where函数获得网格内汇流累积量最大的点的位置，作为子流域控制点，最终生成每个大网格的子流域控制点位置文件。利用控制点沿流向向下游搜索，并记录搜索路径，获得该网格控制点的下游控制点，记录控制点之间的拓扑关系，即为子流域间拓扑关系，生成拓扑关系文件。利用控制点位置文件及DEM信息获得控制点高程，作为子流域的代表高程。利用控制点的搜索路径，计算上游子流域到下游子流域间的河道长度。利用控制点间河道长度和控制点高程差，得到河道坡度信息。
[0034]第三步:编写子流域生成程序。首先，利用子流域拓扑关系生成子流域的上游子流域，将子流域控制点的汇流累积量将去上游子流域的汇流累积量，获得子流域的汇流面积。其次，对子流域集水面积进行由小到大排序，利用IDL编写基于矩阵的流向追踪程序，依次对排序单元进行流向追踪，创建初始矩阵，判断流向流入控制点的网格，将子流域控制区内的网格赋值为子流域编号。最后，得到每个子流域的控制区域，生成子流域控制区域文件。
[0035]第四步:利用第三步中生成的子流域控制矩阵文件，按照各子流域控制区的边界，生成子流域的mif文件格式,方便Maplnfo导入。
[0036]第五步:打开Maplnfo,利用工具----表----1mport工具将mif文件导入到
Maplnfo中，生成流域的子流域划分tab文件。
[0037]第六步:编写MapBasic程序，读取子流域拓扑关系文件、子流域控制点位置文件，生成子流域的拓扑关系图。与第六步的子流域划分图叠加在一起，就生成了流域的拓扑关系图。
[0038]3实列介绍:
实例1:选择王家坝水文站以上流域(如图1所示)进行划分，王家坝水文站位于淮河中游豫皖两省交界处，集水面积为3.6万km2，是淮河干流第一大站，它是淮河的晴雨表，在历次淮河流域抗洪抢险工作中发挥了巨大的作用。
[0039]淮河地处我国南北气候过渡带，淮河以北属暖温带区，淮河以南属北亚热带区，气候温和，多年平均降水量约为830mm，年内降水时空分布不均，其中最大年降水量约为最小年降水量的两倍，汛期(6、月)降水约占全年降水的60%~ 80%，集水面积为27万km2。
[0040]第一步:从国际科学数据服务平台(http://datamirror.csdb.cn/dem/f iles/ys.jsp)下载DEM数据，利用ArcGIS合并功能将数据合并，并利用王家坝站的边界提取出王家坝站流域内DEM。利用ArcGIS的填洼、流向、汇流累积量功能生成无洼地DEM，流域内流向、汇流累积量，并导出为相应ASCII文件。
[0041]第二步:利用生成子流域控制点及拓扑关系的程序，生成王家坝站子流域控制点及拓扑关系。生成子流域的高程、子流域间河道长度、河道坡度等信息，结果分别如表1、2、3所示，其中数据所在行列即为图2大网格化法子流域所在位置的数据。下面的表同理。
[0042]第三步:利用子流域分割程序，对流域进行分割，获得子流域控制面积(如表3)，子流域控制区矩阵。
[0043]第四步:利用子流域控制区矩阵，生成子流域的边界文件mif文件。
[0044]第五步:将mif文件导入Maplnfo,生成子流域的tab文件,如图2所示。
[0045]第六步:利用编写的MapBasic程序，读取生成的子流域控制点文件及子流域拓扑关系文件，生成子流域拓扑关系图，如图3所示。叠加在子流域图上，即为最终图形。
[0046]表1王家坝以上流域提取的子流域面积信息(km2)
【权利要求】
1.一种基于IDL和Mapinfo的子流域划分与信息提取方法，其特征在于，具体步骤如下: 步骤一:获取研究区域DEM ； 步骤二:利用ArcGIS的数据合并功能将获取的分块数据合并为一个大的DEM ;基于研究区的流域边界提取研究区的内的DEM ;利用ArcGIS的hydrology模块的填挂、生成流向、汇流累积量功能获得全流域无洼地DEM、流向、汇流累积量； 步骤三:利用ArcGIS的栅格转ASCII工具，将生成的高精度DEM、流向、汇流累积量结果导为ASCII码，以方便IDL读取； 步骤四:利用IDL编写程序，读取汇流累积量，对流域进行网格划分，利用IDL提取各网格范围内的汇流累积量信息，利用where函数确定汇流累积量最大的点，即子流域的出口点，称之为控制点；利用各网格的子流域划分控制点，生成流域子流域出口控制点文件；利用控制点沿流向向下游搜索，并记录搜索路径，获得该网格控制点的下游控制点，记录控制点之间的拓扑关系，即为子流域间拓扑关系，生成拓扑关系文件；利用控制点位置文件及DEM信息获得控制点高程，作为子流域的代表高程；利用控制点的搜索路径，计算上游子流域到下游子流域间的河道长度；利用控制点间河道长度和控制点高程差，得到河道坡度信息； 步骤五:子流域边界生成，利用子流域拓扑关系生成子流域的上游子流域，将子流域控制点的汇流累积量减去上游子流域的汇流累积量，获得子流域的汇流面积；对子流域汇流面积进行由小到大排序，生成的子流域出口控制点的汇流累积量和；利用IDL编写基于矩阵的流向追踪程序，依次对排序单元进行流向追踪，创建初始矩阵，判断流向流入控制点的网格，将子流域控制区内的网格赋值为子流域编号；最后，得到每个子流域的控制区域，生成子流域控制区域文件； 步骤六:利用步骤四中生成的子流域控制矩阵文件，生成流域的mif文件； 步骤七:将步骤六中生成的mif文件导入到MapInfo中，生成流域的子流域划分文件；步骤八:编写MapBasic程序，读取子流域拓扑关系文件、子流域控制点位置文件，生成子流域的拓扑关系图。
2.如权利要求1所述的基于IDL和Mapinfo的子流域划分与信息提取方法，其特征在于，获取研究区域DEM，是从中国科学院数据应用环境网站上下载研究区域DEM。
3.如权利要求1所述的基于IDL和Mapinfo的子流域划分与信息提取方法，其特征在于，编写MapBasic程序，读取子流域拓扑关系文件、子流域控制点位置文件，生成子流域的拓扑关系图；与第六步的子流域划分图叠加在一起，就生成了流域的拓扑关系图。
【文档编号】G06F9/44GK103677826SQ201310662799
【公开日】2014年3月26日 申请日期:2013年12月9日 优先权日:2013年12月9日
【发明者】刘京京, 吴志勇, 何海, 鞠琴, 董亮, 吴晓韬 申请人:河海大学