CN111986320B - 一种面向智慧城市应用的dem与倾斜摄影模型空间贴合优化方法 - Google Patents
一种面向智慧城市应用的dem与倾斜摄影模型空间贴合优化方法 Download PDFInfo
- Publication number
- CN111986320B CN111986320B CN202010352400.4A CN202010352400A CN111986320B CN 111986320 B CN111986320 B CN 111986320B CN 202010352400 A CN202010352400 A CN 202010352400A CN 111986320 B CN111986320 B CN 111986320B
- Authority
- CN
- China
- Prior art keywords
- dem
- model
- oblique photography
- precision
- original
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 12
- 238000000034 method Methods 0.000 title claims description 43
- 238000011439 discrete element method Methods 0.000 claims abstract description 121
- 230000000694 effects Effects 0.000 claims abstract description 16
- 238000012937 correction Methods 0.000 claims abstract description 12
- 230000006872 improvement Effects 0.000 claims abstract description 9
- 230000004927 fusion Effects 0.000 claims description 16
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 102100022117 Abnormal spindle-like microcephaly-associated protein Human genes 0.000 claims 2
- 101000900939 Homo sapiens Abnormal spindle-like microcephaly-associated protein Proteins 0.000 claims 2
- ZPCCSZFPOXBNDL-ZSTSFXQOSA-N [(4r,5s,6s,7r,9r,10r,11e,13e,16r)-6-[(2s,3r,4r,5s,6r)-5-[(2s,4r,5s,6s)-4,5-dihydroxy-4,6-dimethyloxan-2-yl]oxy-4-(dimethylamino)-3-hydroxy-6-methyloxan-2-yl]oxy-10-[(2r,5s,6r)-5-(dimethylamino)-6-methyloxan-2-yl]oxy-5-methoxy-9,16-dimethyl-2-oxo-7-(2-oxoe Chemical compound O([C@H]1/C=C/C=C/C[C@@H](C)OC(=O)C[C@H]([C@@H]([C@H]([C@@H](CC=O)C[C@H]1C)O[C@H]1[C@@H]([C@H]([C@H](O[C@@H]2O[C@@H](C)[C@H](O)[C@](C)(O)C2)[C@@H](C)O1)N(C)C)O)OC)OC(C)=O)[C@H]1CC[C@H](N(C)C)[C@@H](C)O1 ZPCCSZFPOXBNDL-ZSTSFXQOSA-N 0.000 claims 2
- 238000012545 processing Methods 0.000 abstract description 10
- 238000010276 construction Methods 0.000 description 10
- 238000004519 manufacturing process Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 6
- 238000012876 topography Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000003475 lamination Methods 0.000 description 4
- 230000002093 peripheral effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000003780 insertion Methods 0.000 description 3
- 230000037431 insertion Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000013499 data model Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000009430 construction management Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000007499 fusion processing Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Graphics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Remote Sensing (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Processing Or Creating Images (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种面向智慧城市应用的DEM与倾斜摄影模型空间贴合优化算法,包括以下步骤:输入原始DEM、房屋矢量面、原始倾斜摄影模型;对原始倾斜摄影做纠偏处理得到位置正确的标准倾斜摄影高精模型;基于标准倾斜摄影高精模型与房屋矢量面,采用模型压平算子,形成近似地表倾斜摄影模型;基于近似地表倾斜摄影模型,采用DSM生成算子,形成倾斜摄影模型覆盖区域的近似高精度DEM;将原始DEM与近似高精度DEM镶嵌融合,形成融合后DEM;对比分析两种DEM的优劣,展示本发明的改进效果;结果输出与应用。本发明可以实现DEM和倾斜摄影模型空间匹配,服务于各项三维专题应用,实现智慧城市的空间数据在大尺度三维场景的表达,供智慧城市相关应用调取与使用。
Description
为实现DEM与倾斜摄影模型的贴合服务于三维场景的各专题应用,本发明提供了一种面向智慧城市应用的DEM与倾斜摄影模型空间贴合优化方法,包括如下:
技术领域
本发明涉及空间数据多源融合技术领域,具体涉及一种面向智慧城市应用的DEM与倾斜摄影空间贴合优化算法。
背景技术
随着“智慧城市”概念的快速推广,将城市规划建设管理从粗放模式向精细化、可视化、动态化转变是发展的必然趋势。在“智慧城市”的建设过程中,城市倾斜摄影模型、DSM、DEM等地理信息数据则是“智慧城市”建设的基础所在,没有精度高、效果好的空间数据支撑,“智慧城市”只能是空中楼阁。在实际应用中,需要将倾斜摄影模型叠加到具有DEM的场景中,然而受限于倾斜摄影模型与DEM生成方法、生产时间、采集精度不一致等原因,导致两种数据的叠加不匹配,出现了诸如DEM***倾斜摄影模型的现象。具体表现如下:
(1)生成方法不同:倾斜摄影数据模型主要通过飞行平台搭载多个相机,同时从垂直、倾斜等多个不同的角度采集影像并生成。DEM主要通过航空航天摄影测量和人工野外测量、现有的地形图数字化获取,也可直接通过遥感影像立体像对,根据视差模型选取同名点生成。受到观测技术、生成方法等因素的影响,DEM和倾斜摄影模型往往在空间分辨率、数据精度和覆盖范围等方面存在差异性,造成DEM生产和倾斜摄影模型生产存在平面和高程误差,高程精度误差一般受限于倾斜摄影模型控制点坐标。导致倾斜摄影模型和DEM贴合时存在交融、交叉、展示效果差的问题。
(2)生产时间不同:受限于倾斜摄影模型和DEM来源于多年度、多部门。如:DEM采集与生产于2009年,但倾斜摄影模型生产于2018年,而此段时间内,地貌现状已经发生了很大的变化,如:坑塘已经填平用于构建房屋、部分山坡已挖平用于修建公路。因此造成倾斜摄影模型与DEM不匹配。
(3)数据精度不同:受限于倾斜摄影模型和DEM分辨率精度不一致,如:城镇内部DEM空间分辨率精度一般在2m-10m,城镇外部DEM空间分辨率精度一般在5m-30m,而倾斜摄影模型的模型空间分辨率精度一般在2cm-10cm,因此造成倾斜摄影模型与DEM不匹配。
综上,针对现有DEM和倾斜摄影模型不匹配的问题,通过融合的手段,实现高精度倾斜摄影模型和DEM两者之间的准确贴合成为“智慧城市”建设中重点关注的问题,也是本专利的发明需要解决的重点问题。
发明内容
(一)解决的技术问题
倾斜摄影模型是一种包含了真实地形信息的数字地表模型,在实际应用中,常常需要将倾斜摄影模型叠加到已具有DEM的场景中,而由于倾斜摄影模型与原有DEM的生成方法、生产时间、采集精度不一致,导致两种数据的地形不匹配,出现了诸如原DEM***倾斜模型的现象。为了解决这类问题,本发明提供了一种面向智慧城市应用的DEM与倾斜摄影空间贴合优化算法,可以通过修改原有DEM达到目的。本发明采用选取模型控制点的方法按照1985国家高程标准对拿到的原始倾斜摄影模型进行纠偏处理,减小DEM生产和倾斜摄影模型生产存在数据平面误差和高程误差。采用模型压平算子和DSM生成算子,形成倾斜摄影模型覆盖区域的AH-DEM数据,采用O-DEM与AH-DEM镶嵌融合形成融合后DEM(AF-DEM),DEM数据融合能够克服单一数据的限制,提高数据的利用效率和原始DEM的数据精度,解决DEM生产和倾斜摄影模型生产时间和采集精度不同造成的匹配差异。本发明可提升倾斜摄影覆盖区域的O-DEM数据质量,使AF-DEM和ASPM可以有效贴合的展示在三维场景中,服务于智慧城市建设的各项三维专题应用。本专利实现了已有多种三维GIS技术的整合应用,属于现有技术的整合应用范畴,用于解决现实智慧城市项目建设过程中常见的DEM与倾斜摄影模型贴合效果不佳,从而影响用户体验的问题。
(二)技术方案
本发明所采用的技术方案是:
1.一种面向智慧城市应用的DEM与倾斜摄影模型空间贴合优化方法,包括以下步骤:
1)输入,输入数据包括:原始数字高程模型(后简称O-DEM)、房屋矢量面(后简称HVP)、原始倾斜摄影模型(后简称OPM)。
2)HMSP,获取标准倾斜摄影高精模型(后简称HMSP);
2.1)输入原始倾斜摄影模型(后简称OPM);
2.2)设置空间参考,根据1985国家高程基准;
2.3)选取控制点,采用OPM模型坐标控制点;
2.4)位置校正,对OPM做高度纠偏处理,获得HMSP。
3)AH-DEM,获取近似高精度DEM(AH-DEM);
3.1)模型压平算子,输入HMSP与HVP,采用模型压平算子,形成近似地表倾斜摄影模型(后简称ASPM);
3.2)DSM生成算子,基于ASPM,采用数字地表模型(后简称DSM)生成算子,形成OPM覆盖区域的近似高精度DEM(后简称AH-DEM)。
4)镶嵌融合,O-DEM与AH-DEM镶嵌融合,形成融合后DEM(后简称AF-DEM)。
4.1)栅格镶嵌算子,输入O-DEM与AH-DEM,采用栅格镶嵌算子,实现O-DEM与AH-DEM镶嵌融合;
4.2)梯度下降计算,O-DEM与AH-DEM的接边区域采用梯度下降模型;
4.3)形成AF-DEM。
5)对比,对比分析两种DEM的优劣,展示本发明的改进效果;
5.1)对比分析O-DEM与AF-DEM的优劣;
5.2)对比本发明专利改进效果。
6)输出,结果输出与应用;
6.1)输出AF-DEM,作为本专利方法的输出成果;
6.2)结果应用。
(三)有益效果
与现有技术相比,本发明的有益效果是:一种面向智慧城市应用的DEM与倾斜摄影空间贴合优化算法采用选取模型控制点的方法按照1985国家高程标准对拿到的原始倾斜摄影模型进行高度纠偏处理,减少由数据生成方法、机器性能、天气等因素造成的倾斜摄影模型和DEM模型的平面误差和高程偏差,得到位置正确的ASPM,并在此倾斜摄影模型基础上结合多种算法生成高度正确的倾斜摄影模型覆盖区域的近似高精度DEM(AH-DEM)。O-DEM与AH-DEM镶嵌融合获取到的AH-DEM。这种方法能够克服单一DEM数据的限制,提高数据的利用效率和原始DEM的数据精度。此外,AF-DEM可以与ASPM较好的贴合,解决了O-DEM和OPM无法准确贴合的问题。在智慧城市三维数据构建的实际应用中具有较好的应用价值。
附图说明
图1为本发明流程图。
图2为本发明OPM数据。
图3为本发明HVP数据。
图4为本发明O-DEM数据。
图5为本发明OPM数据存在高度偏差。
图6为本发明倾斜摄影模型生成过程。
图7为本发明倾斜摄影模型流程图。
图8为本发明航线平高点铺设方式。
图9为本发明矫正后的OPM。
图10为本发明ASPM模型数据。
图11为本发明压平算法流程。
图12为本发明AH-DEM地形数据。
图13为本发明O-DEM与HMSP贴合。
图14为本发明AH-DEM与HMSP贴合。
图15表格1为本发明反距离等分权重法内插权值分配表
图16表格2为本发明结果对比分析表
具体实施方式
为使本发明的实施目的、技术方案和优点更加清楚,下面将结合本发明实施方式中的附图,对本发明实施方式中的技术方案进行清楚、完整地描述,显然,所描述的实施方式是本发明一部分的实施方式,而不是全部的实施方式。基于本发明中的实施方式,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
结合图流程图1,一种面向智慧城市应用的DEM与倾斜摄影模型空间贴合优化方法,包括如下步骤:
1)输入,输入数据O-DEM、HVP、OPM,具体数据形式见图2、图3和图4。
2)HMSP,获取HMSP数据。
2.1)输入OPM模型控制点坐标;
原始倾斜摄影模型由于受到航测技术、生成方法、机器性能、天气等因素造成模型与实际地形存在高度偏差,如图5所示,为了得到高度正确的ASPM,采用选取模型控制点对高度进纠偏处理,具体做法如下:
倾斜摄影模型是指通过同一飞行平台搭载多个传感器,同时对从垂直、倾斜等不同角度采集的影像,并通过摄影测量原理,采用自动化建模软件对获得的倾斜影像数据进行几何处理、多视匹配、三角网(TIN)构建、自动赋予纹理输出城市实景三维模型成果,并采用人工修改完成模型局部效果纠正。其倾斜摄影模型生成过程如图6所示,流程图如图7所示。
通过自动化软件生成的倾斜摄影模型受到航测技术、模型生成方法、机器性能、天气等因素造成的倾斜摄影模型高度存在误差,本发明采用选取模型控制点的方法对产生的高度差进行纠偏处理。像控点应该选择在航摄像片上影像清晰、目标明显的像点,实地选点时,也应考虑侧视相机是否会被遮挡。对于弧形地物、阴影、狭窄沟头、水系、高程急剧变化的斜坡、圆山顶、跟地面有明显高差的房角、围墙角等以及航摄后有可能变迁的地方,均不应当做选择目标。控制点选取具体方法如下:
①确定测区所用的坐标***、投影方式、高程基准。
②为保障数据成果精度,对控制点的要求相对提高,需保证1平方公里内30个控制点。房屋顶部应适当增加控制点,从而使数据的精度有进一步的提高。
③采用航线两端及中间均隔一或两条航线布设平高点的方法。此方法既能保证成图精度,又能减少外业工作量。具体布点方式如图8所示。
2.2)设置空间参考,根据1985国家高程基准,采用OPM模型控制点坐标,对OPM做高度纠偏处理,获得HMSP;
中国测绘主管部门通过计算黄海平均海面,并以青岛验潮站1952年~1979年的潮汐观测资料为计算依据,称为“1985国家高程基准”。本发明根据1985国家高程标准。
2.3)选取控制点,采用OPM模型坐标控制点;
2.4)位置校正,对OPM做高度纠偏处理,获得HMSP。本发明专利基于控制点采用GIS桌面软件(本发明专利选取的是GIS软件是超图软件SuperMap iDesktop9D)对OPM做高度纠偏处理,得到纠偏后的OPM如图9所示
3)AH-DEM,获取AH-DEM数据;
3.1)模型压平算子,输入HMSP与HVP,采用模型压平算子,形成近似地表倾斜摄影模型ASPM;
由于倾斜摄影模型与原有地形的采集时间、采集精度不一致,导致两种数据的地形不匹配,出现了诸如原地形***倾斜模型的现象。为了解决这类问题,可以通过替换或修改原有地形数据达到目的。但在缺乏数据修改权限或能力时只能另寻办法。为此本发明提供了在小范围的倾斜摄影模型区域通过HVP,采用模型压平算子对上一步得到的HMSP进行计算得到近似地表倾斜摄影模型ASPM,如图10所示。之后可采用DSM生成算子,用ASPM生成倾斜摄影模型覆盖区域的近似高精度DEM。该方法通过将指定区域的原有倾斜摄影模型上的房屋等建筑按照高程压平,生成ASPM,解决地形***倾斜摄影模型不匹配的问题,且不会影响没有倾斜摄影模型区域的地形显示。倾斜摄影模型可以看成是一张表面覆盖了高分影像的连续的TIN三角网,通常获取到的倾斜摄影模型也是OSGB格式的三维缓存。因此可通过GIS二次开发组件(本发明专利选取的是GIS二次开发组件是超图软件SuperMap iObjects9D)对HMSP进行建筑面的压平计算,其算法流程图如图11所示。
3.2)DSM生成算子,基于ASPM,采用数字地表模型DSM生成算子,形成OPM覆盖区域的AH-DEM;
DSM(Digital Surface Model)即数字地表模型,是指包含了地表建筑物、桥梁、树木等要素的地表高程模型,与DEM的区别在于:DEM只包含了地形高程信息,并未包含其它地表信息,DSM在DEM的基础上,进一步涵盖了除地面以外的其它地表信息的高程,多应用于对建筑物高度有需求的领域。本发明采用改进的四点反距离等分权重内插算法生成DSM,算法主要原理是根据反距离权重内插方法,顾及倾斜摄影模型表面的复杂性,设定距离拟内插点最近的点予以一半的权重,其他点另一半权重,称之为“等分”。这种规则同样应用在远离内插点的三个点之间平衡其权重。可能搜索范围内只有不到四个点的特殊情况,因此需要把没有四个点的情况下缺少点的权值平均分配到各个点上,权值分配如表1所示,表1中以0.15m为距离单元。例如在8倍的距离单元处(1.2m)的点的权值只有0.1,而在0.15m之内的点具有100%的权值。四点反距离等分权重内插算法具体描述如下:
①根据最大搜索半径,建立空间索引网格,最大搜索半径由ASPM顶点间的平均距离确定,一般确定为平均距离的1.2倍。
②利用八邻域,建立DSM拟内插点备选集合。首先判断拟内插点在空间索引格网的编号(x,y,z),确定八邻域的编号,注意在研究区域边际的特殊情况。
③利用改进的反距离权重算法内插DSM格网点(Pxy)。规则如下:
a.如果周边点距离小于最小距离阈值,其值一般取搜索半径小一点的值,有大于0个点,则距离最小点的高程值直接赋给Pxy点。
b.使用这一算法,有搜索不到点的情况,一般是原始倾斜摄影模型存在数据“漏洞”,因此,在数据预处理时需要较长边的三角网,内插出DEM,用DEM点代替DSM。
c.对于个别少量点的“漏洞”区(一般少于3个点),需要搜索其周边的点,判断高程值的大小,如果高程值较大,需要内插周边的点,搜索所有周边的DSM格网点,高程值取算数平均值代替“漏洞”区的点值,以消除DSM的孔洞特征。本发明基于ASPM,采用数字地表模型DSM生成算子,形成OPM覆盖区域的AH-DEM,如图12所示。
4)镶嵌融合,O-DEM与AH-DEM镶嵌融合,形成AF-DEM;
4.1)栅格镶嵌算子,输入O-DEM与AH-DEM,采用栅格镶嵌算子,实现O-DEM与AH-DEM的镶嵌融合;
在DEM获取过程中,数据的精度和尺度不可避免地受到观测手段、处理方法和地表观测环境的影响,使单一来源的DEM在数据质量上往往不能完全满足应用的需求。面对DEM精度异常值、噪声、空洞等常见的质量问题,需要对原始DEM数据进行质量改善。本发明尝试将倾斜摄影模型生成的AH-DEM与O-DEM进行数据融合,挖掘两数据之间的互补信息,建立融合数据集的表达模型,以克服O-DEM数据质量的缺陷,提高数据精度。
DEM表达的是经过处理和解译得到的地表高程信息,由于地表高程起伏的稳定性,在中高分辨率的数据中很难区分短时间内的地形沉降、地物变化与数据的***和随机误差。并且,DEM的灰度变化反映的是地形的高低起伏,不同的DEM可能通过不同传感器在不同时间获取而造成的误差,观测手段和生成方法也非常多样化,因此DEM精度受到观测值本身和生成过程***误差和随机误差的共同影响。在数据融合中,需要结合数据本身的特性,并根据应用需求对多源数据的误差进行考虑。本发明采用将稀疏表示的框架引入与O-DEM和AH-DEM的融合过程,建立数据集的表达模型,形成融合后的AF-DEM数据模型。这种方法可提高O-DEM和AH-DEM重合区域的数据的分辨率。
4.2)梯度下降计算,O-DEM与AH-DEM的接边区域采用梯度下降模型;
大区域的O-DEM与小范围的高精度AH-DEM融合后形成的AF-DEM,在接边区域存在明显的高度差异。为了减小这种差异,使接边区域可以平缓过渡,本发明尝试在接边区域采用梯度下降模型(Gradient Descent,简称GD)对AF-DEM做优化。GD作为目前优化神经网络最基础的迭代算法,其最基本思想是:用在小批量的子集上计算出的梯度近似在整个数据集上计算出的真实的梯度。如公式(1-1),GD每一步用小批量样本迭代更新权重。
ωt+1=ωt+Δωt
Δωt=-ηΔωE(ωt) (1-1)
其中,η是算法的学习率,E(ωt)是关于第t次迭代权重ωt的损失函数,ΔωE(ωt)为权重w在t时刻关于损失函数的一阶梯度,简记为gt,ωt+1为t+1时刻的权重值ωt为t时刻的权重值,Δω为梯度算子,即每次迭代的权重更新部分。
4.3)形成AF-DEM;
本发明实现O-DEM与AH-DEM的镶嵌融合,并通过对接边区域采用梯度下降模型优化形成AF-DEM。
5)对比,对比分析两种DEM的优劣,展示本发明的改进效果;
5.1)对比分析O-DEM与AF-DEM的优劣;
由前面步骤可知,AF-DEM是AH-DEM与O-DEM融合的结果,与O-DEM相比,AF-DEM在倾斜摄影模型覆盖区域的地形拥有相对较好的现势性。为了做改进前后的贴合效果对比,分别将O-DEM与AF-DEM,同HMSP于三维场景中叠加展示,如图13、图14所示,对比指标具体包括:相同区域内DEM与HMSP交叉区域大小、DEM与HMSP交叉区域个数、道路面所在区域的模型贴合度,根据此三项指标,做本专利方法的用益效果对比,具体对比结果如表2所示。
5.2)对比本发明专利改进效果;
本发明的改进效果是:采用选取模型控制点的方法对拿到的原始倾斜摄影模型做纠偏处理,减少由数据生成方法、机器性能、天气等因素造成的倾斜摄影模型和DEM模型的平面误差和高程偏差,得到HMSP,并在此基础上结合多种算法生成位置正确的倾斜摄影模型覆盖区域的AH-DEM,其与O-DEM镶嵌融合获取到的AH-DEM。这种方法能够克服单一DEM数据的限制,提高数据的利用效率和原始DEM的数据精度,解决了O-DEM和OPM无法准确贴合的问题。在智慧城市三维数据构建的实际应用中具有较好的应用价值。
6)结果输出与应用;
6.1)输出AF-DEM,作为本专利方法的输出成果;
AF-DEM与HMSP空间叠加的展示效果如图14所示。可以明显看到倾斜摄影数据已经和地形数据有了较好的贴合。其成果可通过超图软件的Supermap iDesktop生成*.S3M格式三维缓存。
6.2)结果应用;
DEM数据能够描述地面的起伏形态特征,因此广泛应用在地形分析、地质、水文领域,在生态、农业、气象及军事应用中也扮演着重要的角色。倾斜摄影数据提供了现实世界建筑物的表面模型,地形数据提供了倾斜摄影数据在宏观环境中的位置信息及周边环境。将倾斜摄影数据与DEM数据进行融合,为用户提供更加真实、准确的可视化体验,以及三维空间分析数据支撑,服务于智慧城市建设的各项三维专题应用。
综上所述,本发明的一种面向智慧城市应用的DEM与倾斜摄影模型空间贴合优化方法,解决了小范围的倾斜摄影模型和大范围的DEM贴合效果存在交融、交叉、展示效果差的问题。实现了主城区倾斜摄影模型与地形数据色贴合,可用于智慧城市建设,在实际应用中具有较好的应用价值。
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施方式,凡属于本发明思路下的技术方案均属于本发明的保护范围。对于本技术领域的普通技术人员而言,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。
Claims (1)
1.一种面向智慧城市应用的DEM与倾斜摄影模型空间贴合优化方法,包括以下步骤:
1)输入原始数字高程模型(Original digital elevation model,简称O-DEM)、房屋矢量面(House vector polygon,简称HVP)、原始倾斜摄影模型(Original photographicmodel,简称OPM);
2)获取标准倾斜摄影高精模型(High precision model of standard photographic简称HMSP);
2.1)输入OPM模型坐标控制点;
2.2)根据1985国家高程基准,采用OPM模型坐标控制点,对OPM做高度纠偏处理,获得HMSP;
3)获取近似高精度DEM(AH-DEM);
3.1)输入HMSP与HVP,采用模型压平算子,形成近似地表倾斜摄影模型(Approximatesurface photographic model,简称ASPM);
3.2)基于ASPM,采用数字地表模型(Digital Surface Model,DSM)生成算子,形成OPM覆盖区域的近似高精度DEM(Approximate high-precision DEM of photographic modelcovered area,简称AH-DEM);
4)O-DEM与AH-DEM镶嵌融合,形成融合后DEM(After fusion DEM,简称AF-DEM);
4.1)输入O-DEM与AH-DEM,采用栅格镶嵌算子,实现O-DEM与AH-DEM镶嵌融合;
4.2)O-DEM与AH-DEM的接边区域采用梯度下降模型;
4.3)形成AF-DEM;
5)对比分析两种DEM的优劣,展示本方法的改进效果;
5.1)对比分析O-DEM与AF-DEM的优劣;
5.2)对比本方法改进效果;
6)结果输出与应用;
6.1)输出AF-DEM,作为本方法的输出成果;
6.2)结果应用。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010352400.4A CN111986320B (zh) | 2020-04-28 | 2020-04-28 | 一种面向智慧城市应用的dem与倾斜摄影模型空间贴合优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010352400.4A CN111986320B (zh) | 2020-04-28 | 2020-04-28 | 一种面向智慧城市应用的dem与倾斜摄影模型空间贴合优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111986320A CN111986320A (zh) | 2020-11-24 |
CN111986320B true CN111986320B (zh) | 2024-01-30 |
Family
ID=73441769
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010352400.4A Active CN111986320B (zh) | 2020-04-28 | 2020-04-28 | 一种面向智慧城市应用的dem与倾斜摄影模型空间贴合优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111986320B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117830539A (zh) * | 2022-09-27 | 2024-04-05 | 杭州海康威视***技术有限公司 | 地形数据更新方法、装置、设备、存储介质 |
CN117095151B (zh) * | 2023-10-19 | 2024-01-19 | 成都坤舆空间科技有限公司 | 基于共轭梯度的dem与倾斜摄影空间贴合初匹配优化方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102855810A (zh) * | 2012-09-04 | 2013-01-02 | 绍兴文理学院 | 一种基于利用卫星影像图进行数字线划地图的方法 |
CN107576311A (zh) * | 2017-08-23 | 2018-01-12 | 长江水利委员会长江科学院 | 一种基于三维gis的水库巡检实时监测方法 |
CN108052776A (zh) * | 2018-01-12 | 2018-05-18 | 河南省水利勘测设计研究有限公司 | 基于bim和三维实景模型的洪水模型构建与动态展示方法 |
CN108050995A (zh) * | 2017-12-14 | 2018-05-18 | 滁州学院 | 一种基于dem的倾斜摄影无像控点航摄测区归并方法 |
KR20180131932A (ko) * | 2017-06-01 | 2018-12-11 | 충남대학교산학협력단 | 드론과 공간정보를 이용한 하천지형정보 생성 방법 |
CN109840338A (zh) * | 2017-11-28 | 2019-06-04 | 南京国图信息产业有限公司 | 一种应用于三维不动产信息管理的三维楼盘模型构建方法 |
CN109934914A (zh) * | 2019-03-28 | 2019-06-25 | 东南大学 | 一种嵌入式城市设计场景仿真方法及*** |
CN110136259A (zh) * | 2019-05-24 | 2019-08-16 | 唐山工业职业技术学院 | 一种基于倾斜摄影辅助bim和gis的三维建模技术 |
CN110763205A (zh) * | 2019-11-05 | 2020-02-07 | 新疆维吾尔自治区测绘科学研究院 | 数字摄影测量***生成边境狭长区域正射影像图的方法 |
CN110955742A (zh) * | 2019-11-28 | 2020-04-03 | 中国电建集团成都勘测设计研究院有限公司 | 基于dem的升压站场坪雨水口设置区域分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2582971A1 (en) * | 2004-10-15 | 2006-04-20 | Ofek Aerial Photography International Ltd. | Computational solution of and building of three dimensional virtual models from aerial photographs |
-
2020
- 2020-04-28 CN CN202010352400.4A patent/CN111986320B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102855810A (zh) * | 2012-09-04 | 2013-01-02 | 绍兴文理学院 | 一种基于利用卫星影像图进行数字线划地图的方法 |
KR20180131932A (ko) * | 2017-06-01 | 2018-12-11 | 충남대학교산학협력단 | 드론과 공간정보를 이용한 하천지형정보 생성 방법 |
CN107576311A (zh) * | 2017-08-23 | 2018-01-12 | 长江水利委员会长江科学院 | 一种基于三维gis的水库巡检实时监测方法 |
CN109840338A (zh) * | 2017-11-28 | 2019-06-04 | 南京国图信息产业有限公司 | 一种应用于三维不动产信息管理的三维楼盘模型构建方法 |
CN108050995A (zh) * | 2017-12-14 | 2018-05-18 | 滁州学院 | 一种基于dem的倾斜摄影无像控点航摄测区归并方法 |
CN108052776A (zh) * | 2018-01-12 | 2018-05-18 | 河南省水利勘测设计研究有限公司 | 基于bim和三维实景模型的洪水模型构建与动态展示方法 |
CN109934914A (zh) * | 2019-03-28 | 2019-06-25 | 东南大学 | 一种嵌入式城市设计场景仿真方法及*** |
CN110136259A (zh) * | 2019-05-24 | 2019-08-16 | 唐山工业职业技术学院 | 一种基于倾斜摄影辅助bim和gis的三维建模技术 |
CN110763205A (zh) * | 2019-11-05 | 2020-02-07 | 新疆维吾尔自治区测绘科学研究院 | 数字摄影测量***生成边境狭长区域正射影像图的方法 |
CN110955742A (zh) * | 2019-11-28 | 2020-04-03 | 中国电建集团成都勘测设计研究院有限公司 | 基于dem的升压站场坪雨水口设置区域分析方法 |
Non-Patent Citations (4)
Title |
---|
Spatio-Temporal Landscape Analysis in MountainousTerrain by Means of Small Format Photography: AMethodological Approach;Janette Aschenwald;IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING;全文 * |
基于BIM的铁路建设管理平台及关键技术研究;智鹏;中国优秀博士学位论文全文数据库-工程科技II辑;全文 * |
无人机低空倾斜摄影测量技术在玛纳斯智慧城市建设中的应用;李杰;马强;;岩土工程技术(第04期);全文 * |
无人机倾斜摄影测量技术在桥梁施工现场中的应用研究;耿小平;王波;马钧霆;袁海天;闫亚妮;;现代测绘(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111986320A (zh) | 2020-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Balsa-Barreiro et al. | Generation of visually aesthetic and detailed 3D models of historical cities by using laser scanning and digital photogrammetry | |
US7944547B2 (en) | Method and system of generating 3D images with airborne oblique/vertical imagery, GPS/IMU data, and LIDAR elevation data | |
US8427505B2 (en) | Geospatial modeling system for images and related methods | |
CN102506824B (zh) | 一种城市低空无人机***生成数字正射影像图的方法 | |
CN100578525C (zh) | 生成物体三维电子模型的***和方法 | |
CN109238239B (zh) | 基于航空摄影的数字测量三维建模方法 | |
CN111724477A (zh) | 一种多源数据融合构建多层次三维地形模型方法 | |
CN105976426B (zh) | 一种快速的三维地物模型构建方法 | |
CN105046251B (zh) | 一种基于环境一号卫星遥感影像的自动正射校正方法 | |
WO2018061010A1 (en) | Point cloud transforming in large-scale urban modelling | |
Yan et al. | Integration of 3D objects and terrain for 3D modelling supporting the digital twin | |
CN101154313A (zh) | 三维仿真数字信息航道***及实现方法 | |
CN109472865B (zh) | 一种基于图像模型绘制的自由可量测全景再现方法 | |
US8395760B2 (en) | Unified spectral and geospatial information model and the method and system generating it | |
CN110660125B (zh) | 用于电力配网***的三维建模装置 | |
CN111986320B (zh) | 一种面向智慧城市应用的dem与倾斜摄影模型空间贴合优化方法 | |
US11922572B2 (en) | Method for 3D reconstruction from satellite imagery | |
CN103954970A (zh) | 一种地形要素采集方法 | |
Bolkas et al. | Creating a virtual reality environment with a fusion of sUAS and TLS point-clouds | |
CN116883604A (zh) | 一种基于天、空、地影像的三维建模技术方法 | |
CN111986074A (zh) | 一种真正射影像制作方法、装置、设备及存储介质 | |
CN111831778A (zh) | 一种三维地理信息***快速集成展现的方法 | |
Jiang et al. | Low–high orthoimage pairs-based 3D reconstruction for elevation determination using drone | |
KR100732915B1 (ko) | 디지털사진 측량기술 및 인공위성영상을 이용한 기본설계용도로노선의 3차원적 결정 방법 | |
CN111982077B (zh) | 电子地图绘制方法、***和电子设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |