CN111833445B - 一种基于多源数据的区域地形分割及数字高程模型获取方法 - Google Patents
一种基于多源数据的区域地形分割及数字高程模型获取方法 Download PDFInfo
- Publication number
- CN111833445B CN111833445B CN202010590427.7A CN202010590427A CN111833445B CN 111833445 B CN111833445 B CN 111833445B CN 202010590427 A CN202010590427 A CN 202010590427A CN 111833445 B CN111833445 B CN 111833445B
- Authority
- CN
- China
- Prior art keywords
- digital
- elevation model
- digital elevation
- terrain
- regional
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000011218 segmentation Effects 0.000 title claims abstract description 17
- 238000004519 manufacturing process Methods 0.000 claims abstract description 12
- 238000001914 filtration Methods 0.000 claims abstract description 6
- 238000000605 extraction Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000012876 topography Methods 0.000 claims description 7
- 230000000903 blocking effect Effects 0.000 claims description 3
- 238000013075 data extraction Methods 0.000 claims description 3
- 230000001172 regenerating effect Effects 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 2
- 238000007621 cluster analysis Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
Abstract
本发明公开了一种基于多源数据的区域地形分割及数字高程模型获取方法,包括以下步骤:对地形进行分割,并形成相应的等级;对不同等级的区域采取自适应的方式获取与之相匹配的数字地表模型,并进行滤波;对不同级别的成果进行融合,获取可用于数字正影像图生产的初始数字高程模型。本发明采用地形分割的方法获取分级数字高程模型,在此基础上对数字高程模型进行处理,可以极大程度上减轻后期数字高程模型的人工处理工作量,提高数字高程模型生产及数字正射影像生产的自动化程度。
Description
技术领域
本发明属于数字正射影像图生成技术领域,具体涉及一种基于多源数据的区域地形分割及数字高程模型获取方法。
背景技术
数字高程模型生产是数字正射影像生产的关键技术流程之一,也是数字正射影像生产中人工干预工作量占比重的环节。目前数字正射影像生产中用于数字微分纠正的数字高程模型一般是在历史数字高程模型或者按统一标准参数提取滤波出的数字高程模型基础上进行人工处理获得的,该方式自动化程度低,人工处理工作量大。
发明内容
本发明提供了一种基于多源数据的区域地形分割及数字高程模型获取方法,提高生产的自动化及信息化程度。
本发明的技术方案为:一种基于多源数据的区域地形分割及数字高程模型获取方法,包括以下步骤:
S1、进行区域地形的分割,并形成相应的等级;
S2、对不同等级的区域采取自适应的方式获取与之相匹配的数字地表模型,并进行滤波;
S3、对不同类型的成果进行融合,获取可用于数字正射影响图生产的初始数字高程模型。
作为优选,所述步骤S1包括:根据区域地形特点,将地形区分成三种等级:L0,L1,L2,其中,L0级代表所在区域划分成平地区域,L1级代表所在区域划分成山地区域,L2级代表所在区域划分成平山结合处。
作为优选,所述步骤S1中,L0级通过坡度改正计算式(1)计算,
i=2*h/l*100% (1);
若即将相应区域归到L0;
其中h为地形高差,l为水平距离,i为坡度计算值。
作为优选,所述步骤S1中,L1级通过最新地理国情地表覆盖要素代码的聚类分析进行归算,建立聚类函数,见式(2),
L=x1+x2+x3+.....xn (2);
其中x1,x2,x3.....xn为与山地有关的代码要素标识。
本发明中将与山地有关的代码要素都合并到相应的层里,并将相应区域归算到L1级。
作为优选,所述步骤S2-1中,将L0级及L1级提取出的区块与最新时相的历史数字正射影像、数字线划图进行对比纠正,最终确定相应层级的成果,利用式(3),将剩余地形区域归算到L2级,
L2=L总-L0-L1 (3)。
作为优选,所述步骤S2包括:
S2-1、数字表面模型获取:利用L0层级的成果范围面作为约束范围面获取数字表面模型,提取规则设置为山地,为便于后期数据的加载和处理,在点云数据提取过程中作分块处理;
S2-2、数字高程模型获取:利用L1层级的成果范围作为约束范围面并利用时相最接近的数字高程模型成果作为约束条件获取山地区域初始的数字高程模型,综合考虑后期对L2成果的处理,提取规则设置为丘陵地。
作为优选,所述步骤S3包括:
S3-1、对步骤S2-1中获取的分块点云数据进行滤波处理;
S3-2、对步骤S2-2中初始的数字高程模型进行植被高程改正;
S3-3、获取L2层级的数字高程模型:对步骤S3-1中滤波后的点云数据进行高程内插,生成平地区域的栅格数据,用于内插平山分界线,由此获取三维特征线,从而处理平山结合处的数字高程模型格网点;
S3-4、将步骤S3-1中获取的点云数据、步骤S3-2中的数字高程模型成果、步骤S3-3中的三维特征线联合处理并进行高程内插,重新生成经初始处理后的数字高程粗成果。
作为优选,步骤S3-1中涉及到人工构筑物阈值的确定以及地形坡度变化阈值的设置,人工构筑物阈值的确定通过提取时相最接近的数字线划图中的居民地矢量数据层的房屋要素为样本进行设定,地形坡度变化阈值以L0层的初值为基础适当进行调整。
作为优选,步骤S3-2中,进行植被高程改正时,通过式(4)计算高程值,
其中,h为成果高程值,h0为提取出的高程值,h1为对步骤S2-2中作为约束条件的初始值。
作为优选,所述数字表面模型的格网间距为1米。
与现有技术相比,本发明的有益效果体现在:
本发明中利用与现有数据源时相最接近的地表覆盖数据、数字正射影像图数据、数字线划图等数据结合代码聚类分析、特征提取及坡度改正计算完成区域地形的分割并形成相应的等级,对不同等级的区域采取自适应的方式获取与之相匹配的数字地表模型并进行滤波,对不同类型的成果进行融合获取可用于后续数字正射影像图生产的初始数字高程模型,可以极大程度上减轻后期数字高程模型的人工处理工作量,提高数字正射影像生产的自动化程度。
附图说明
图1为本发明的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
本实施例为一种基于多源数据的区域地形分割及数字高程模型获取方法,如图1所示,主要包括以下步骤:
I:地形分割
根据区域地形特点,将地形区分成三种等级:L0,L1,L2。L0级代表所在区域划分成平地区域,L1级代表所在区域划分成山地区域,L2级代表所在区域划分成平山结合处。
(1)L0级通过坡度改正计算式计算,
i=2*h/l*100% (1)
若即将相应区域归到L0;
其中h为地形高差,l为水平距离,i为坡度计算值。
(2)L1级通过最新的国情地表要素代码的聚类分析进行归算,建立聚类函数,见式(2),将阔叶林等与山地有关的代码要素都合并到相应的层里,并将相应区域归算到L1级,
L=x1+x2+x3+.....xn (2);
其中x1,x2,x3.....xn为与山地有关的代码要素标识。
(3)将L0级及L1级提取出的区块与最新时相的历史数字正射影像、数字线划图进行对比纠正,最终确定相应层级的成果。
(4)利用式(3),将剩余地形区域归算到L2级。
L2=L总-L0-L1 (3)。
II:点云提取
(1)数字表面模型获取:利用L0层级的成果范围面作为约束范围面获取格网间距为1米的数字表面模型,提取规则设置为山地,为便于后期数据的加载和处理,在点云数据提取过程中作分块处理。
(2)数字高程模型获取:利用L1层级的成果范围作为约束范围面并利用时相最接近的数字高程模型成果作为约束条件获取山地区域初始的数字高程模型,综合考虑后期对L2成果的处理,提取规则设置为丘陵地。
(3)对所有区域进行点云数据获取,提取规则设置为平地。
III:点云融合及数字高程模型获取
(1)对II(1)中获取的分块点云数据进行滤波处理:该步骤涉及到人工构筑物阈值的确定及地形坡度变化阈值的设置,人工构筑物阈值的确定通过提取时相最接近的数字线划图中的居民地矢量数据层的房屋要素为样本进行设定,地形坡度变化阈值以L0层的初值为基础适当进行调整。
(2)对II(2)中初始的数字高程模型进行植被高程改正:通过式(4)计算高程值,
其中,h为成果高程值,h0为提取出的高程值,h1为对II(2)中作为约束条件的初始值。
(3)获取L2层级的数字高程模型:对III(1)中滤波后的点云数据进行高程内插,生成平地区域的栅格数据,用于内插平山分界线,由此获取三维特征线,从而处理平山结合处的数字高程模型格网点。
(4)将III(1)中获取的点云数据、III(2)中的数字高程模型成果、III(3)中的三维特征线联合处理并进行高程内插,重新生成经初始处理后的数字高程粗成果。
Claims (7)
1.一种基于多源数据的区域地形分割及数字高程模型获取方法,其特征在于,包括以下步骤:
S1、进行区域地形的分割,并形成相应的等级;
S2、对不同等级的区域采取自适应的方式获取与之相匹配的数字地表模型,并进行滤波;
S3、对不同类型的成果进行融合,获取可用于数字正射影像图生产的初始数字高程模型
其中,步骤S1包括:根据区域地形特点,将地形区分成三种等级:L0,L1,L2,其中,L0级代表所在区域划分成平地区域,L1级代表所在区域划分成山地区域,L2级代表所在区域划分成平山结合处;
步骤S2包括:
S2-1、数字表面模型获取:利用L0层级的成果范围面作为约束范围面获取数字表面模型,提取规则设置为山地,为便于后期数据的加载和处理,在点云数据提取过程中作分块处理;
S2-2、数字高程模型获取:利用L1层级的成果范围作为约束范围面并利用时相最接近的数字高程模型成果作为约束条件获取山地区域初始的数字高程模型,综合考虑后期对L2成果的处理,提取规则设置为丘陵地;
步骤S3包括:
S3-1、对步骤S2-1中获取的分块点云数据进行滤波处理;
S3-2、对步骤S2-2中初始的数字高程模型进行植被高程改正;
S3-3、获取L2层级的数字高程模型:对步骤S3-1中滤波后的点云数据进行高程内插,生成平地区域的栅格数据,用于内插平山分界线,由此获取三维特征线,从而处理平山结合处的数字高程模型格网点;
S3-4、将步骤S3-1中获取的点云数据、步骤S3-2中的数字高程模型成果、步骤S3-3中的三维特征线联合处理并进行高程内插,重新生成经初始处理后的数字高程粗成果。
2.如权利要求1所述的基于多源数据的区域地形分割及数字高程模型获取方法,其特征在于,所述步骤S1中,L0级通过坡度改正计算式(1)计算,
i=2*h/l*100% (1);
若即将相应区域归到L0;
其中h为地形高差,l为水平距离,i为坡度计算值。
3.如权利要求2所述的基于多源数据的区域地形分割及数字高程模型获取方法,其特征在于,所述步骤S1中,L1级通过最新地理国情地表覆盖要素代码的聚类分析进行归算,建立聚类函数,见式(2),
L=x1+x2+x3+.....xn (2);
其中x1,x2,x3.....xn为与山地有关的代码要素标识。
4.如权利要求3所述的基于多源数据的区域地形分割及数字高程模型获取方法,其特征在于,所述步骤S2-1中,将L0级及L1级提取出的区块与最新时相的历史数字正射影像、数字线划图进行对比纠正,最终确定相应层级的成果,利用式(3),将剩余地形区域归算到L2级,
L2=L总-L0-L1 (3)。
5.如权利要求4所述的基于多源数据的区域地形分割及数字高程模型获取方法,其特征在于,步骤S3-1中涉及到人工构筑物阈值的确定以及地形坡度变化阈值的设置,人工构筑物阈值的确定通过提取时相最接近的数字线划图中的居民地矢量数据层的房屋要素为样本进行设定,地形坡度变化阈值以L0层的初值为基础适当进行调整。
6.如权利要求5所述的基于多源数据的区域地形分割及数字高程模型获取方法,其特征在于,步骤S3-2中,进行植被高程改正时,通过式(4)计算高程值,
其中,h为成果高程值,h0为提取出的高程值,h1为对步骤S2-2中作为约束条件的初始值。
7.如权利要求1所述的基于多源数据的区域地形分割及数字高程模型获取方法,其特征在于,所述数字表面模型的格网间距为1米。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010590427.7A CN111833445B (zh) | 2020-06-24 | 一种基于多源数据的区域地形分割及数字高程模型获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010590427.7A CN111833445B (zh) | 2020-06-24 | 一种基于多源数据的区域地形分割及数字高程模型获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111833445A CN111833445A (zh) | 2020-10-27 |
CN111833445B true CN111833445B (zh) | 2024-05-31 |
Family
ID=
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101006729B1 (ko) * | 2010-07-23 | 2011-01-10 | (주)동광지엔티 | 수치 표고모델 제작 방법 및 시스템 |
JP2012243226A (ja) * | 2011-05-23 | 2012-12-10 | Nec System Technologies Ltd | 地盤高データ生成装置、地盤高データ生成方法、及びプログラム |
CN103970837A (zh) * | 2014-04-11 | 2014-08-06 | 北京工业大学 | 基于城市用地和竖向规划的不连续dem分类制作方法 |
CN105205328A (zh) * | 2015-09-28 | 2015-12-30 | 辽宁工程技术大学 | 一种基于广义分维和特征参数的地质动力环境评价方法 |
CN107305701A (zh) * | 2017-05-14 | 2017-10-31 | 杭州师范大学 | 一种基于数字高程模型的城市洼地提取方法 |
CN108961286A (zh) * | 2018-06-27 | 2018-12-07 | 武汉大学 | 一种顾及建筑物三维及边缘形状特征的无人机影像分割方法 |
CN109579784A (zh) * | 2018-11-26 | 2019-04-05 | 青岛国测海遥信息技术有限公司 | 基于数字表面模型的城市地区建筑物高度的自动获取方法 |
CN110427857A (zh) * | 2019-07-26 | 2019-11-08 | 国网湖北省电力有限公司检修公司 | 一种基于遥感数据融合的输电线路地质灾害分析方法 |
CN110763205A (zh) * | 2019-11-05 | 2020-02-07 | 新疆维吾尔自治区测绘科学研究院 | 数字摄影测量系统生成边境狭长区域正射影像图的方法 |
CN110941689A (zh) * | 2019-11-18 | 2020-03-31 | 云南瀚哲科技有限公司 | 一种基于ArcGIS的地貌类型划分方法 |
CN111103595A (zh) * | 2020-01-02 | 2020-05-05 | 广州建通测绘地理信息技术股份有限公司 | 一种数字线划图的生成方法和装置 |
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101006729B1 (ko) * | 2010-07-23 | 2011-01-10 | (주)동광지엔티 | 수치 표고모델 제작 방법 및 시스템 |
JP2012243226A (ja) * | 2011-05-23 | 2012-12-10 | Nec System Technologies Ltd | 地盤高データ生成装置、地盤高データ生成方法、及びプログラム |
CN103970837A (zh) * | 2014-04-11 | 2014-08-06 | 北京工业大学 | 基于城市用地和竖向规划的不连续dem分类制作方法 |
CN105205328A (zh) * | 2015-09-28 | 2015-12-30 | 辽宁工程技术大学 | 一种基于广义分维和特征参数的地质动力环境评价方法 |
CN107305701A (zh) * | 2017-05-14 | 2017-10-31 | 杭州师范大学 | 一种基于数字高程模型的城市洼地提取方法 |
CN108961286A (zh) * | 2018-06-27 | 2018-12-07 | 武汉大学 | 一种顾及建筑物三维及边缘形状特征的无人机影像分割方法 |
CN109579784A (zh) * | 2018-11-26 | 2019-04-05 | 青岛国测海遥信息技术有限公司 | 基于数字表面模型的城市地区建筑物高度的自动获取方法 |
CN110427857A (zh) * | 2019-07-26 | 2019-11-08 | 国网湖北省电力有限公司检修公司 | 一种基于遥感数据融合的输电线路地质灾害分析方法 |
CN110763205A (zh) * | 2019-11-05 | 2020-02-07 | 新疆维吾尔自治区测绘科学研究院 | 数字摄影测量系统生成边境狭长区域正射影像图的方法 |
CN110941689A (zh) * | 2019-11-18 | 2020-03-31 | 云南瀚哲科技有限公司 | 一种基于ArcGIS的地貌类型划分方法 |
CN111103595A (zh) * | 2020-01-02 | 2020-05-05 | 广州建通测绘地理信息技术股份有限公司 | 一种数字线划图的生成方法和装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Metz et al. | Efficient extraction of drainage networks from massive, radar-based elevation models with least cost path search | |
CN110675392B (zh) | 一种基于路面连续激光点云的车辙精细三维特征提取方法 | |
CN108257142B (zh) | Dem中斜坡单元提取方法 | |
CN111047704B (zh) | 一种改进区域生长算法的多波束测深数据粗差自动清除方法 | |
CN110335352B (zh) | 一种机载激光雷达点云的双基元多分辨率层次滤波方法 | |
WO2013121340A1 (en) | Digital elevation model | |
WO2011085435A1 (en) | Classification process for an extracted object or terrain feature | |
CN111833445B (zh) | 一种基于多源数据的区域地形分割及数字高程模型获取方法 | |
CN109741451B (zh) | 一种基于地形图构建三维地表bim模型的系统 | |
CN111382715A (zh) | 一种地形图中地物要素的提取方法 | |
KR101547943B1 (ko) | 폴리곤과 dem을 이용한 유역추출 장치 및 방법 | |
Al-Muqdadi et al. | Automated watershed evaluation of flat terrain | |
CN110660027A (zh) | 一种针对复杂地形的激光点云连续剖面地面滤波方法 | |
CN115828054B (zh) | 一种通过改进拉普拉斯自动识别南支槽的方法 | |
CN116612245B (zh) | 基于视频图像的海滩地形地貌构建方法、系统、存储介质 | |
CN111833445A (zh) | 一种基于多源数据的区域地形分割及数字高程模型获取方法 | |
CN113838059B (zh) | 一种基于要素级别的数字正射影像图生成方法 | |
Ghazal et al. | Change Detection Study Of Al Razaza Lake Region Utilizing RemoteSensing And GIS Technique | |
CN111833446A (zh) | 一种基于特征线提取的架空地物快速纠正方法 | |
CN115937676A (zh) | 一种基于dem的大尺度沙丘地貌精细分类方法 | |
WO2011085437A1 (en) | Extraction processes | |
Mohammed et al. | Using Remote Sensing and GIS to Study Morphological Analysis of Kirkuk Province | |
Li et al. | Building boundaries extaction from point clouds using dual-threshold Alpha Shapes | |
CN111178372B (zh) | 基于遥感影像和地形数据的大区域尺度黄土塬面提取方法 | |
Novák et al. | The Potential and Implications of Automated Pre-Processing of LiDAR-Based Digital Elevation Models for Large-Scale Archaeological Landscape Analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |