期刊目录 | 《地球信息科学学报》第22卷第3期“数字地形分析”专栏

特邀主编介绍

汤国安教授

 奥地利理学博士，南京师范大学教授，从事GIS空间分析、基于DEM的数字地形分析等方面的科学研究工作，主持8项国家自然科学基金项目（含2项重点），获四项省部级以上科学研究奖励，以及“中国地理学会优秀地理科技工作者”、中共中央组织部“万人计划”领军人才（教学名师）及人力资源和社会保障部、教育部“全国模范教师”等荣誉称号。

熊礼阳副教授

       南京师范大学副教授，主要从事数字地形分析、黄土地貌演化建模与模拟等方面的科学研究工作，主持2项国家自然科学基金项目，发表论文50余篇，出版英文专著1部，获教育部自然科学奖二等奖、江苏省优秀博士论文、江苏省青年遥感与地理信息科技奖等奖励，担任国际数字地球学会数字山地专业委员会副主任委员。

0

黄土高原沟谷地貌发育演化研究

进展与展望

熊礼阳，汤国安*

摘要：黄土沟谷是黄土地貌中最有活力、最具变化、最富特色的对象单元，黄土高原千沟万壑的地貌形态以及触目惊心的侵蚀 状态也让区域内沟谷地貌的形成、发育及演化问题成为研究中焦点及前沿性科学问题。近年来，诸多学者采用地学测年法、 特征表达法、监测模拟法力图实现对黄土沟谷发育演化进程中“过去-现代-未来”的科学认知。这些研究在相当程度上丰富 了黄土沟谷发育过程的认知。本文梳理了黄土高原沟谷地貌演化相关研究的现状，并从黄土高原地貌演化、黄土沟谷发育、 基于DEM的沟谷信息提取与表达等研究进行了系统的回顾、梳理与分析。此外，本文提出“黄土沟道剖面群组”概念与方法， 试图从新的视角审视黄土沟谷地貌发育演化过程。沟道剖面在黄土沟谷发育演化进程中传递物质能量和累积地形动力，并 通过径流节点的串联实现剖面群的连接与组合，形成独特的剖面“群组”模式；该沟道剖面群组是集黄土沟谷地貌特征与过程 于一体的综合信息集成体，其三维空间结构是对黄土沟谷地貌发育演化的高度抽象与映射，并可望进一步丰富黄土高原数字 地形分析理论与方法体系，为黄土高原黄土地貌成因机理与空间分异格局带来创新的认识。

沟谷地貌演化研究不同层次

（此文被收录在陈述彭先生纪念专刊，于2020年第4期发表）

1

DEM分辨率对黄土侵蚀沟形态特征表达的不确定性分析

李思进，代文，熊礼阳*，汤国安

摘要：黄土侵蚀沟的地形表达是开展黄土沟谷侵蚀研究的基础工作,利用数字高程模型（DEM）定量描述侵蚀沟特征有助于研究侵蚀沟的形态变化和发育过程。基于DEM数据计算多种指标对黄土侵蚀沟特征进行描述是目前侵蚀沟研究中最为常用的方法。但是,受到格网DEM数据结构的限制,其计算结果会存在一定的不确定性。在侵蚀沟地形表达时,对形态特征的表达会受到DEM数据分辨率的影响,进而造成表达结果的不确定性。尤其在黄土高原地区,地形特征更为破碎,地形要素更为复杂,其表达结果受DEM分辨率的影响更为明显。本文以黄土高原典型样区为例,基于点云数据建立不同分辨率的DEM数据集,通过不同地形因子对侵蚀沟特征进行表达,分析DEM分辨率在黄土侵蚀沟形态特征表达时的不确定性。结果显示,分辨率的降低对主沟支沟比和纵比降等侵蚀沟形态特征因子产生了较大影响,且指标与分辨率多呈现线性变化关系。但是,随着侵蚀沟的横向扩张,DEM分辨率对其特征表达的影响逐渐被削弱。此外,在使用固定分析窗口进行侵蚀沟特征计算时,由于分辨率的降低,格网尺寸增大,其实际分析半径随之增大,使得计算范围内地表形态变化增加,导致沟谷切割深度随着分辨率的降低反而增加。同时,侵蚀沟主沟道区域受分辨率影响较小,沟头区域指标与分辨率的关系较弱。

分辨率对沟谷切割深度的影响

2

基于ICESat/GLAS的山东省SRTM与ASTERGDEM高程精度评价与修正

秦臣臣，陈传法*，杨娜，高原，王梦樱

摘要：SRTM3和ASTER GDEM V2数据具有较高的空间分辨率和广泛的覆盖范围,对于地学研究具有重要意义;但在不同地形复杂度和地面覆盖物区域,两类数据的误差分布并不均匀。SRTM3和ASTER GDEM V2 数据自公布以来,其精度修正一直是研究热点。然而大范围区域精度验证缺乏有效手段,传统方法可靠性差且数据获取成本较高。自ICESat-1数据公开以来,它们已成为SRTM3和ASTER GDEM V2精度评定的主要检核点。为此,本文以山东省为研究区域,借助ICESat-1评估了SRTM3和ASTER GDEM V2的高程精度,并根据插值误差曲面对两种DEM进行了修正。分析表明,原始SRTM和ASTER高程中误差分别为5.57 m和7.20 m,均高于标称精度;随着坡度的增大,高程精度呈降低的趋势。通过分析土地覆盖类型与误差分布关系表明：农田、灌丛土地类型精度较高;森林、湿地精度较低。分别采用反距离加权、普通克里金、地形转栅格和自然邻域插值方法构建误差曲面。结果表明：不同的插值方法构建的误差曲面的特征和精度也不同。其中,反距离加权修正的效果最佳,其次是地形转栅格和自然邻域,而普通克里金修正的效果最差。

山东省原始SRTM和ASTER与修正后对比

3

规则格网DEM局地坡面凸凹性精度分析

王春*，徐燕，江岭，赵明伟

摘要：DEM地表形态精度分析理论与方法的建立,对DEM数据的生产和广泛应用具有重要意义。本文从局地坡面形态的凸凹性角度,剖析规则格网DEM格网点位置、格网分辨率对DEM局地坡面凸凹性的影响,以期进一步完善和发展DEM质量分析的理论与方法。论文首先阐述了DEM局地坡面凸凹性的基本概念,研究建立了规则格网DEM的局地坡面凸凹性量化分析方法,并以黄土丘陵5、10、15、25、……、155 m DEM为例,采用比较分析方法研究了局地坡面凸凹性随DEM格网点位置和格网分辨率的变化特征。研究表明：对于本研究中的1:5万DEM,10 m（跃变率≤ 0.3%）是其最佳的格网分辨率阈值,当DEM实际格网分辨率高于该阈值时,实际DEM与最佳格网分辨率DEM具有近乎相同的局地坡面凸凹性,主要在正地形与负地形的过渡区域会发生不同程度的坡面凸凹性变化;当DEM实际格网分辨率低于该阈值时,实际DEM的局地坡面凸凹性,会随着DEM格网点布设位置和DEM格网分辨率发生较大的不确定性变化。

DEM局地坡面凸凹性的误差空间分布

4

资源三号卫星DEM数据与全球开放DEM数据的误差对比

赵尚民*,程维明,蒋经天,沙文娟

摘要：全球开放DEM数据为数字地形分析提供了重要数据源。与已有的全球开放DEM数据相比,资源三号卫星具有更高的空间分辨率、更大的覆盖范围和更好的现势性。将资源三号卫星生成的DEM数据与全球开放DEM数据进行误差对比则为基于资源三号卫星的全球DEM数据研制提供科学依据。本文以山西省中部太原市为研究区,基于高精度激光点云数据生成DEM为参考数据,对资源三号卫星影像生成的DEM数据与全球典型的开放DEM数据（AW3D30、SRTM1和ASTER GDEM）的误差进行了对比分析,并获得了其在不同坡度等级下绝对误差与相对误差的平均值、平均绝对值、均方根值和标准偏差值。研究结果表明：① 4种DEM数据的误差分布均具有较好的对称性。同时,平均误差接近于0 m,SRTM1和ASTER GDEM数据更是如此。因此均方根误差值与标准偏差值近似一致;② 资源三号DEM具有最高的精度,误差最小（均方根误差4.6 m）。其次为AW3D30数据（均方根误差5.6 m）和SRTM1数据（均方根误差8.8 m）。ASTER GDEM数据误差最大（均方根误差12.6 m）,精度最差;③ 资源三号DEM、SRTM1和ASTER GDEM数据的误差均随坡度的变大而增大,而AW3D30数据误差随着坡度增加呈现先减小后增大的趋势。总体上,与其他3种DEM数据相比,资源三号DEM在所有坡度范围均具有最小的误差值。

4种DEM数据的相对误差计算结果

5

一种用于DEM数据认证与篡改定位的感知哈希算法

张鑫港，闫浩文*，张黎明

摘要：DEM数据作为重要的基础地理信息数据,其数据完整性问题不容忽视。基于DEM数据完整性认证的要求,以及相关认证算法的欠缺,本文运用感知哈希技术设计了一种DEM数据认证算法,并可实现篡改定位。因DEM数据具有数据量大、细节丰富的特点,首先对其进行规则格网划分,将其划分为互不重叠的格网单元;然后对格网单元数据进行DCT分解,提取数据的特征信息以生成特征向量矩阵,并对特征向量矩阵进行摘要化处理;随后,使用Logistic混沌系统对简化后的特征向量矩阵进行置乱;对置乱矩阵进行量化、编码后,便可生成感知哈希序列。在数据认证时,首先计算原始数据与待验证数据的高程相对中误差,再将二者的感知哈希序列进行归一化汉明距离度量,结合判定阈值,即可对DEM数据进行数据认证与篡改定位。该算法对DEM数据的格式转换、水印嵌入等攻击有较强的鲁棒性,对各类改变内容的操作具有敏感性,并可实现DEM数据微小篡改的识别与定位。与已有的DEM完整性认证方法相比,将DEM数据的“内容”作为完整性度量的重要标准,在具体应用中更具有实用价值。

篡改定位测试

6

多模型协同下的城郊地区DEM构建

方法研究

赵明伟*，金永林，江岭，王春，杨灿灿，徐燕

摘要：人类活动对地表的改造使其呈现形态多样化、不连续等特征,此时传统的DEM构建方法难以满足这些区域DEM精度要求。为此,本文提出了一种多模型协同构建DEM的思路,首先按照形态特征和语义信息对地形进行分类,然后对不同类型的地形区域,选择、设计适宜的方法分别进行DEM构建,最后将不同区域构建DEM结果融合拼接形成区域完整的DEM结果。本文选择江苏省南京市城市郊区某区域为实验区,以1:500比例尺地形图为基本数据源进行DEM构建实验。实验结果表明,与传统经典DEM构建方法相比,本文提出的多模型协同的DEM构建方法能够有效表达实验区域不同的地形特征,特别是对于人工改造的地形（如道路、边坡等区域）,本文方法构建的DEM其形态精度优势显著;同时,基于验证点法的高程精度分析结果表明,本文方法构建DEM的高程精度亦优于传统DEM构建方法,特别是对于一些形态规则而高程信息相对稀少的区域,以边坡区域为例,经典DEM构建法平均误差均超过5 m,而本文构建结果平均误差为0.26 m,精度优势非常明显。研究表明本文提出的多模型协同的DEM构建方法适用于人类活动改造或显著影响的区域的DEM构建。

不同方法构建DEM对应的坡度

7

基于傅里叶变换的谷间距特征信息提取及其影响因素研究

蔡顺*，耿豪鹏*，郑炜珊，潘保田

摘要：谷间距（Valley Spacing）是描述相邻沟谷距离的特征参数,它能较好地反映沟谷的次序等级和空间分布特征。本研究以黑河正义峡和兰州大岭岘2个样区为例,利用ArcGIS软件,将无人机测绘获取的0.12 m分辨率地表高程数据通过重采样生成不同分辨率的数字地表模型。通过MATLAB软件,将不同分辨率、不同空间域的数字地表模型作为二维空间域信号进行傅里叶变换。通过傅里叶变换及频谱分析研究地形的频谱特征与地表谷间距之间的转换关系。分析结果显示：① 当区域内只有一级沟谷时,频谱中谷间距特征信号的有效提取要求地形分辨率至少优于1/5谷间距,分辨率的粗略化则直接影响着地形频谱中谷间距特征信号的识别,但是分析空间域对频谱谷间距特征信号的捕获影响较小;② 当区域内有多级沟谷时,分辨率优于1/3谷间距时即可有效提取到该级沟谷的谷间距特征信号,分辨率的粗略化和空间域的增大都会使得频谱中较低序次等级沟谷的谷间距特征信号减弱,而较高序次等级沟谷的谷间距特征信号增强。

兰州大岭岘地区1414 m×2071 m空间域不同分辨率地形特征谷间距提取

8

基于等高线空间关系的鞍部点提取方法

黄楠，杨昕*，刘海龙

摘要：鞍部点是反映地表形态起伏变化的重要地形特征点之一,准确地提取鞍部点有利于地形的空间关系和结构特征分析。现有的鞍部点提取方法通常是直接基于规则格网DEM数据,无法顾及鞍部点与周围地形的空间拓扑关系和复杂地形对其的影响,不仅产生大量的伪鞍部点,而且忽略一些关键地区的鞍部点。本文根据鞍部点的地形形态特征,设计了一种基于等高线数据的鞍部点提取算法。该算法利用等高线闭合的特征,将等高线按照一定规则转成等高面数据,再利用等高面之间的相邻拓扑关系实现递归查找并自动提取鞍部点。实验结果显示：① 鞍部点的数量和位置与等高距的大小显著相关,在一定尺度范围内,等高距越小,提取出鞍部点越多,位置精度也逐渐提高;② 与基于规则格网DEM数据提取方法相比,该方法能更有效的过滤大量伪鞍部点,提高了鞍部点的提取精度,同时也降低了鞍部点提取算法的复杂度;与基于等高线的增量缓冲方法（Incremental Buffering Algorithm）相比,本文的方法能有效提高鞍部点提取的完整性,更适用于本文DEM的尺度即5 m DEM数据。

实验区不同方法提取结果对比

9

基于DEM的一体化山地特征要素提取

胡金龙，唐梦鸽，罗明良*，魏兰，

晏自红，秦子晗

摘要：山顶点和山脊线等特征地形要素是构成地表地形及其起伏变化的基本框架，对地形在地表的空间分布具有控制作用。基于DEM研究山顶点、山脊线及其空间组合关系，是DEM地表形态特征研究的重要内容，也是衔接从地形特征分析向山峰等地貌学本源语言的途径之一。本文以四川盆地西南缘与青藏高原过渡地带的川西凉山山原为例，基于山峰—山脊线—控制范围一体化构建的算法策略，识别了山峰和山脊线及其等级、主山脊及其范围。结果表明，研究区内有主峰9座，次峰53 座，平均高程2540m；山脊线230条，其中主山脊9条，平均长度60km；9大山系，近南北走向，平均控制面积1017km2。研究用模糊隶属度方法对算法所提取的主峰、主脉进行精度验证，隶属度介于0.98~1.00 和0.37~0.57 时提取的主峰、主脉基本吻合算法提取的结果。研究采用一体化山地特征要素提取方法，实现了各山地要素间紧密联系、总体结构与区域地貌特征相对吻合的目标；完成了由栅格单元向地理对象的转变；可以应用于协助地貌类型划分，协同区域地理规划等。

凉山山原区隶属度提取与一体化算法提取主峰对比

10

基于DEM的黄土地貌逼近度因子构建及特征分析

雷雪，周毅*，李阳，王泽涛

摘要：沟谷源点作为沟沿线上最为活跃的部位，其分别到上游分水线、下游沟谷线的流线空间比对关系，是表征三线空间结构的重要突破口。它在空间上向流域分水线逼近的程度，是量化黄土流域地貌系统发育程度的重要切入点。为探究黄土地貌区流域沟谷源点向流域分水线逼近的程度，揭示黄土流域地貌发育进程及该进程所表现的主要侵蚀方式，本文从水平和垂直2 个维度，构建量化三线空间结构关系的核心因子——逼近度（PI），其中包括水平逼近度（HPI）和垂直逼近度（VPI），基于5 m分辨率数字高程模型，在陕北黄土高原遴选了包含16 种地貌类型的42 个样区，利用数字地形分析方法，探讨其平均值（MHPI, MVPI）的空间分异规律。另外，选择其中南北序列分别代表黄土塬、残塬、梁状丘陵沟壑和峁状丘陵沟壑区的淳化、宜君、甘泉和绥德4个地区，完备包含一至五级别沟谷的典型流域为重点实验样区，探讨陕北黄土高原流域尺度平均逼近度变异指数（MPIV）序列性特征。实验结果表明：①陕北黄土高原平均逼近度（MPI）存在强烈的空间自相关，MHPI在南北序列上先增大后减小，东西方向上逐渐减小，在黄河沿岸峡谷丘陵区MHPI达到最大；MVPI由西南—东北先减小后增加，由西北—东南逐渐减小，在渭北黄土台塬区达到最小；② 在流域尺度上，MPIV值的正负与黄土塬区、丘陵沟壑区敏感相关；③4个重点样区MHPI、MVPI与其它地形因子存在南北序列上的一致性。104个外部汇流区平均水平逼近度与平均坡度相关性较好（P=0.43, a<0.001），平均垂直逼近度与面积高程积分强烈相关（P=0.75, a<0.001）。平均逼近度指标综合考量了黄土高原地区最典型的三条具有结构控制意义特征线的空间关系，对黄土地貌的发育程度有明显的指示性意义。

水平与垂直方向逼近度指标构建原理

11

基于地形剖面的黄土塬梁峁地形各向异性奇异强度研究

曹建军*，汤国安，方炫，刘永娟，那嘉明，朱莹，李金莲

摘要：黄土高原不同地貌类型对应的地形起伏变化特征各不相同，地形起伏在不同方向上的变化特征具有很大差异性，即地形各向异性。本文引入多重分形谱参数描述不同地貌类型地形各向异性的变化规律，以反映不同黄土地貌地形各向异性的局部和整体变化格局。本研究以5 m分辨率的DEM数据为基础，选取淳化、甘泉、绥德3 个典型黄土塬梁峁地貌类型研究样区，对其地形各向异性的变化特征进行分析。研究表明：①3 种黄土地貌类型地形各向异性变化具有显著的多重分形特征，地形各向异性变化的奇异强度由强到弱依次为淳化、绥德、甘泉；②淳化样区地形各向异性奇异强度在1.4 的概率最大，在奇异强度为2.2 附近又有一个小的峰值，总体以奇异强度1.4 为主；③甘泉样区的奇异强度在1.8 有一个小的峰值，总体以奇异强度0.7 为主；④绥德样区奇异强度为0.8 的概率最大，且较为均匀地分布在0.8 附近，地形各向异性的概率分布呈现为对称分布。研究结果为黄土塬梁峁地形各向异性变化奇异强度为黄土塬最大，黄土峁次之，黄土梁最小。该研究可为黄土地貌的精确分类提供定量信息支持，为区域水土侵蚀评价提供地形指标。

研究样区地形各向异性奇异强度概率分布

12

基于随机森林的黄土地貌分类研究

曹泽涛，方子东，姚瑾，熊礼阳*

摘要：地貌分类在指导人类建设活动的规模与布局中有着重要的意义。然而，传统的基于数字高程模型（DEM）的地貌分类方法使用的地形因子和考虑到的地貌特征往往比较单一。本文提出了一种基于流域单元的地貌分类方法，该方法考虑了流域单元的多方面特征，包括基本地形因子统计量、地形特征点线统计量、小流域特征和纹理特征。本研究首先基于DEM进行水文分析将研究区域划分成不同的小流域。然后利用数字地形分析提取29个不同方面的特征来表征流域的形态，并基于随机森林（RF）算法进行了特征选择和参数标定。RF是一种基于决策树算法的集成分类器，能有效地处理高维数据，分类精度高。最后选择训练集小流域对RF分类器进行训练，使用训练完成的分类器对整个研究区域的地貌进行分类，研究地貌分异的规律。该实验在我国陕北黄土高原典型黄土地貌区域的地貌分类中取得了较好的结果，结果表明不同的地貌之间存在明显的区域界线，特定的地貌类型在空间上表现出明显的聚集性。通过人工判读进行验证的分类精度达到了85%，Kappa系数为0.83。

基于随机森林的小流域单元分类结果

13

基于BEMD分解的地貌分类研究

顾文亚，孟祥瑞，朱晓晨，邱新法*

摘要：地貌是指地势高低起伏的变化，即地表的形态。地貌划分对气温、降水、太阳辐照等诸多应用领域都有着重要作用。本文选择空间分辨率为90 m的福建省的数字高程模型（DEM）数据作为地理信号，运用二维经验模态分解（BEMD）进行分解处理，得到多个具有不同尺度、不同物理意义的本征模函数（BIMF1—BIMF3）以及对应余量ORIG。这些BIMF分量对应不同尺度的微观地形，ORIG余量表现为该研究区的地貌分布趋势，体现了平原、丘陵与山地的大致分布区域。运用变点分析法确定最佳计算单元，利用地形起伏度对各个信号区域进行一级分类，依据绝对高度进行二级分类，最后将一级分类与二级分类相结合，实现对地形的分类，这一分类过程体现了研究地区地貌组合复杂的特征。结果表明：① 叠加BIMF分量，提取出分量和大于74 m的区域为高频信号区域。该区域以小起伏度的低山为主，并伴随有丘陵和小起伏度中山。②ORIG余量中余量高度小于等于340 m的区域，去除其中包含的高频信号区域后为低频信号区域，该区域以平原、丘陵为主。③剩余区域定义为中频信号区域，该区域的地貌以平丘陵和小起伏的山地为主。研究成果表明福建地貌可分为7种主要类型：低频平原，低频丘陵，中频丘陵，高频丘陵，中频小起伏低山，高频小起伏低山，高频小起伏中山。

基于2000 年福建省DEM数据地形分解的高频词提取

14

面向复杂地形的坡位K-means聚类

划分研究

高海峰，葛莹*，张杰，肖胜昌，陈科

摘要：土壤植被研究建立在精准坡位划分的基础上。但现有的坡位大多采用手工划分的方式，存在着自动化程度低、划分精度不高且耗时较长等问题。本文提出一种顾及复杂地形的坡位自动划分算法，尝试采用机器学习K-means 方法解决高海拔山区坡位划分的问题，并在山峰区域提取、聚类数确定、以及初始聚类中心选取等关键技术进行了算法的优化。为了验证算法的有效性，以云南省姚安县为研究区，运用提出的算法对研究区坡位进行自动划分，再采用Calinski-Harabasz 聚类评价指标、调整兰德系数ARI和误差平方和SSE等一系列方法对坡位K-means 聚类划分实验进行分析和评价。研究结果表明，利用该算法所生成的复杂地形坡位与研究区实测等高线相匹配。其次，再从姚安县规划风电场任选4 个场址，比较13×13、25×25、37×37三种适宜窗口下坡位自动划分结果，结果表明选取25×25适宜窗口进行坡位划分可靠性最强。再者，计算的规划风电场内山脊、坡肩及背坡比例高达57.13%，也从一个侧面证实了利用该算法划分的坡位结果良好。

坡位划分

15

神农架林区植被分布与地形的关系研究

张兴航，张百平*，王晶，姚永慧，余付勤

摘要：神农架林区是我国物种多样性最为丰富的地区之一，地形地貌复杂，对植被分布影响巨大。本文利用该地区2007 年数字高程数据、2007 年植被分布图以及2017 年野外实地调查数据，基于最大熵模型和空间分析理论，从植被类型和种群两个角度研究该地区不同尺度植被空间分布的地形特征，分别量化植被类型和种群空间分布的地形范围，得到植被类型与地形因子关系模型、植被种群与地形因子关系模型。结果表明：① 神农架林区影响植被空间分布的地形因子不同，其中影响针叶林分布的最重要的地形因子是高程和高程变异系数，影响阔叶林分布的是高程和坡向，影响灌丛分布的是坡向变率和坡向，影响草丛分布的较为分散；② 典型植被种群分布的地形范围和植被类型的基本一致，其中90%针叶林分布在高程1600~2600 m间，典型种群巴山冷杉和华山松主要分布在高程1700~3200 m和1700~2200 m；85%的阔叶林分布在高程1000~2000 m间，典型种群青冈类和鹅耳枥主要分布在高程1200~2200 m间；95%的灌丛分布在坡向变率0~40°间，典型种群杜鹃和蔷薇主要分布在坡向变率小于40°的范围，但相应的关系模型存在差异，植被类型与地形因子为高斯模型，典型种群与地形因子关系模型相对复杂，不同种群的分布模式不同；③ 虽然坡度常作为数字地形的重要因子，但本文研究发现该地区坡度对植被类型和种群分布的影响不明显。研究结果可为神农架林区植被保护和恢复，以及植被规划和管理提供基础参考。

神农架南北坡植被优势种群在高程上的分布

16

Jason-2 测高卫星对湖泊水位的监测精度评价

何飞，刘兆飞*，姚治君

摘要：湖泊水位是评估湖泊水量变化的重要指标。本文以洪泽湖、高邮湖及洞庭湖为研究对象，利用集中度的概率密度函数方法（CPDF）来提高Jason-2 测高数据精度，分析了降水量与各个湖泊水位变化的相关性，并基于实测水位数据对比评价了Jason-2 测高卫星原始GDR数据和CPDF方法处理后的卫星数据的精度。结果表明：① Jason-2 原始GDR数据点的分布存在疏密之分，大部分数据分布相对集中，且有一定的周期变化，但评价结果显示精度较差，故原始GDR数据不能直接用于湖泊水位监测；② CPDF方法可以极大提高测高卫星的水位数据精度，洪泽湖与高邮湖的均方根误差分别由1.92 m与1.74 m减少到了0.32 m和0.36 m，相关系数由0.28 和0.04 提高到了0.85 和0.72。对于南北宽度较窄且日水位变化较大的湖泊（如洞庭湖），CPDF方法提高原始GDR结果的精度有限；③ 洞庭湖降水与水位相关性最强，高邮湖次之，而洪泽湖降水与水位成不显著的负相关，是洪泽湖水利工程对于水位的调节导致了这一结果。本研究对于利用测高卫星获得湖泊水位值，进而对湖泊进行动态监控，特别是在填补资料匮乏地区湖泊水位数据方面具有重要意义。

CPDF方法流程

17

基于CUDA的地表水动态模拟并行方法

吴钱娇，陈玉敏*

摘要：如何高效地动态监测、模拟和预测地表水过程是防灾减灾中亟待解决的问题，也是科学化的国土整治、区域规划、环境保护和水资源管理的基础。因此，本文利用统一计算设备架构（CUDA）对基于不规则三角网（TIN）的地表水动态模拟算法进行并行化改进，提出了一种基于CUDA的地表水动态模拟并行方法，旨在对任意时刻的地表水进行快速、高精度的动态模拟，从而满足实际的应用需求。该算法从高精度的数字高程模型（DEM）中提取地形特征点和流域线，生成受流域线约束的TIN。在此基础上，根据TIN表面的三角面坐标数据获取水流方向，再结合任意位置的降雨源点追踪得到流水线网络。基于曼宁公式，利用流水线流速计算核函数得到每条流水线上雨滴的流速，结合预设的时间，利用汇流量统计核函数得到该时刻的地表汇流量。具体的并行化过程包括数据的传输，CUDA中线程的划分和流水线流速计算核函数和汇流量统计核函数的实现。本文选取位于加拿大新布伦瑞克西北部的BBW（Black Brook Watershed）流域作为实验样区，将该算法与改进前的方法进行对比表明，该模型在保持精度的同时，大幅提高了模拟效率，加速比达到11.2；另外，通过与SWAT（Soil and Water Assessment Tool）模型进行对比，结果表明，该并行方法的Nash系数提高了88%，相关系数提高了5%。

不同降雨源点尺度下的流水线网络（阈值为0.5 m）

（原文请到中国知网下载阅读）

《地球信息科学学报》是由中国科学院地理科学与资源研究所和中国地理学会联合主办，由中国科学院资源与环境信息系统国家重点实验室、虚拟地理环境教育部重点实验室（南京师范大学）、三维信息获取与应用教育部重点实验室（首都师范大学）共同协办的学术期刊。学报为月刊，是中国科学引文数据库（CSCD）核心期刊、中国科技核心期刊、全国中文核心期刊、北大《中文核心期刊要目总览》来源期刊。《2018年版中国科技期刊引证报告（核心版）》显示，《地球信息科学学报》影响因子位列测绘科学技术期刊第3名，综合评价总分位列第4名。

转载自《地球信息科学学报》

经作者授权转载

文章仅代表作者观点，与本公众号无关，版权归原作者所有

原文标题：数字地形分析：理论与方法创新

图文排版：罗湘阳

责任编辑：田可

审编：王波涛 李倩

终审：顾伟男 田巍 梁龙武

猜你喜欢

都是成年人了，要记得告诉地小联你也“在看”哦~👇