一种基于倾斜摄影城市三维模型的高大植被提取方法技术

技术编号:25398402 阅读:41 留言:0更新日期:2020-08-25 23:03
本发明专利技术涉及一种基于倾斜摄影城市三维模型的高大植被提取方法,具体包括以下步骤:S1:初始植被区域提取:先计算纹理空间的植被指数,再将三角网格中的每个三角面同绑定的UV值投影到对应纹理上,获取纹理上三角网覆盖区域的像素点集合,进一步的得到单个三角面的植被指数;S2:植被区域精化:将植被指数相似的空间邻接三角面聚集成区域,再通过最大流/最小割定理计算最小割路径,将区域分成对应一致性类别小区域,以实现对植被区域的提取;S3:高大植被提取:空间分割,单个植被对象的高程判断。本申请的一种基于倾斜摄影城市三维模型的高大植被提取方法,自动化程度高、无需人工干预,解决现有植被提取方法无法在三维场景中适用的技术问题。

【技术实现步骤摘要】
一种基于倾斜摄影城市三维模型的高大植被提取方法
本专利技术涉及测绘
,特别是涉及一种基于倾斜摄影城市三维模型的高大植被提取方法。
技术介绍
随着低空倾斜系统平台的快速发展与推广,出现了一种全新获取城市场景三维重建的方法:倾斜摄影。倾斜摄影是通过在同一飞行器平台上,搭载多台传感器,同时从一个垂直、四个倾斜等五个不同的角度采集影像,基于多视影像的三维建模方法,恢复目标对象的三维结构的同时获取三维模型表面的纹理信息。近几年,随着无人机技术的日益成熟,我们获取影像数据的成本也大大降低。因此,基于多视倾斜影像的城市三维场景的构建就变得更加的广泛。植被作为城市生态系统中的重要组成部分,其提取方法是城市规划、生态评估的基础。与此同时,基于多视影像的自动化重建结果由一整张三角网格以及相对应的纹理信息组成,虽然较为真实的表达了城市三维场景,但无对应植被属性信息。影像数据与实景三维数据示意图如图1所示,与常规影像数据相比,实景三维模型数据具有以下特点:1)由像素组成,MVS数据基本单元为三角网以及三角网对应一块纹理区域;2)与二维栅格影像不同,MVS数据具有在三维空间中的拓扑连接性,导致现有植被提取方法在此处无法适用。
技术实现思路
为了克服现有技术的上述不足,本专利技术提出了一种基于倾斜摄影城市三维模型的高大植被提取方法,解决现有三维模型中植被提取难度大、无法实现的技术问题。本专利技术是通过以下技术方案实现的:一种基于倾斜摄影城市三维模型的高大植被提取方法,具体包括以下步骤:S1:初始植被区域提取:将三维场景模型中的三角网格表示为M={fi},每张三角网格对应的纹理影像I={Ii},根据纹理空间的植被指数计算公式计算纹理空间的植被指数VI',再将三角网格中的每个三角面fi同绑定的UV值投影到对应纹理Ii上,获取纹理Ii上三角网fi覆盖区域的像素点集合进一步的得到单个三角面fi的植被指数VIi;S2:植被区域精化:将植被指数相似的空间邻接三角面聚集成区域,再通过最大流/最小割定理计算最小割路径,将聚集之后的区域分成对应一致性类别小区域,以实现对植被区域的提取;S3:高大植被提取:对提取出的植被对象进行空间分割,并将空间分割的单个对象的高程与设定的阈值进行比较,大于阈值,即为高大植被。进一步的,S1中纹理空间的植被指数VI'=(2*G'-R'-B')–(1.4*R'-G'),其中,G'=G/(R+G+B),R'=R/(R+G+B),B'=B/(R+G+B),R、G、B分别表示波段红、波段率、波段蓝。进一步的,S1中计算单个三角面fi的植被指数VIi通过面积规则化植被指数:其中,Si=Square(fi)表示三角面fi的平面面积。进一步的,S2中将植被指数相似的空间邻接三角面聚集成区域的具体方式是:迭代计算与当前区域相连接的三角面到区域中心的距离,判断其是否满足设定的判别条件,如果三角面满足条件,则加入到区域中,更新区域的植被指数参数,直到没有满足条件的相邻三角面加入为止。进一步的,判断三角面到当前区域中心距离是否满足判别条件具体操作为:以正态分布的概率密度函数f(μ,σ)=N(μ,σ)来描述区域的属性,用是否在μ±α·σ内作为三角面是否属于当前区域的判别准则,其中:μ为区域中心植被指数参数,σ为植被指数参数的变化程度,参数α与邻域的平滑程度相关。进一步的,S2中求解最小割路径的步骤:构造前景/背景二元标号的能量函数,能量函数由数据项和平滑项两部分组成,表示为:E(f)=Edata(f)+γ·Esmooth(f),数据项损失函数Edata(f)=∑p∈PDp(fp),平滑项损失函数Esmooth(f)=∑p∈PV{p,q}∈N(fp,fq),γ系数控制数据项与平滑项对总体损失函数的影响效力,计算能量函数的最小值,即为最小分隔路径。进一步的,以高斯函数f(v)=N(μ,σ)作为概率密度函数描述数据项损失,数据中心项为最理想情况下获取的植被。进一步的,平滑项损失V{p,q}∈N(fp,fq)=L(fp,fq)*G(fp,fq),其中,L(fp,fq)表示对象fp,fq之间公共边长度之和,G(fp,fq)表示fp,fq之间的距离。进一步的,S3中的阈值设为10m。与现有技术相比,本专利技术的有益效果在于:本专利技术提出的一种基于倾斜摄影城市三维模型的高大植被提取方法,实现了城市三维场景数据中高大植被的提取,为后续的城市规划、生态评估提供支持,而且无需人工干预、自动化程度高。说明书附图图1为
技术介绍
中所述影像数据与实景三维数据示意图;图2为本专利技术实施例所述纹理空间植被指数提取示意图;图3为本专利技术实施例所述三角面植被指数提取示意图;图4为本专利技术实施例所述网格空间聚类结果图;图5为本专利技术实施例所述植被提取结果示意图;图6为本专利技术实施例所述一种基于倾斜摄影城市三维模型的高大植被提取方法的流程图。具体实施方式展示一下实例来具体说明本专利技术的某些实施例,且不应解释为限制本专利技术的范围。对本专利技术公开的内容可以同时从材料、方法和反应条件进行改进,所有这些改进,均应落入本专利技术的精神和范围之内。如图6所示,本申请一种基于倾斜摄影城市三维模型的高大植被提取方法,具体包括以下步骤:S1:初始植被区域提取:将三维场景模型中的三角网格表示为M={fi},每张三角网格对应的纹理影像I={Ii};S11:参照文献VerificationofColorVegetationIndicesforAutomatedCropImagingApplications中体提及的纹理空间植被指数计算公式VI'=(2*G'-R'-B')–(1.4*R'-G'),其中,G'=G/(R+G+B),R'=R/(R+G+B),B'=B/(R+G+B),R、G、B分别表示波段红、波段率、波段蓝,计算纹理空间的植被指数VI′,空间植被指数提取结果如图2所示;S12:将三角网格中的每个三角面fi同绑定的UV值投影到对应纹理Ii上,获取纹理Ii上三角网fi覆盖区域的像素点集合再计算单个三角面fi的植被指数VIi通过面积规则化植被指数:其中,Si=Square(fi)表示三角面fi的平面面积;三角面植被指数提取结果如图3所示;S2:植被区域精化:S21:在三维网格中,每一个三角形对应的植被指数可以描述其局部区域性质,局部性质的变化可以通过比较当前三角形与周围三角形植被指数的差异来描述,将植被指数相似的空间邻接三角面聚集成区域;将植被指数相似的空间邻接三角面聚集成区域的具体方式是:迭代计算与当前区域相连接的三角面到区域中心的距离,判断其是否满足设定的判别条件,如果三角面满足条件,则加入到区域中,更新区域的植被指数参数,直到没有满足条件的相邻三角面加入为止;为统一判断方程,以正态分布的概率密度函数f(μ,σ)=N(μ,σ)来描述区域的属性,用是否在μ±α·σ内作为三角面是否属于当前本文档来自技高网
...

【技术保护点】
1.一种基于倾斜摄影城市三维模型的高大植被提取方法,其特征在于,具体包括以下步骤:/nS1:初始植被区域提取:将三维场景模型中的三角网格表示为M={f

【技术特征摘要】
1.一种基于倾斜摄影城市三维模型的高大植被提取方法,其特征在于,具体包括以下步骤:
S1:初始植被区域提取:将三维场景模型中的三角网格表示为M={fi},每张三角网格对应的纹理影像I={Ii},根据纹理空间的植被指数计算公式计算纹理空间的植被指数VI′,再将三角网格中的每个三角面fi同绑定的UV值投影到对应纹理Ii上,获取纹理Ii上三角网fi覆盖区域的像素点集合进一步的得到单个三角面fi的植被指数VIi;
S2:植被区域精化:将植被指数相似的空间邻接三角面聚集成区域,再通过最大流/最小割定理计算最小割路径,将聚集之后的区域分成对应一致性类别小区域,以实现对植被区域的提取;
S3:高大植被提取:对提取出的植被对象进行空间分割,并将空间分割的单个对象的高程与设定的阈值进行比较,大于阈值,即为高大植被。


2.根据权利要求1所述的一种基于倾斜摄影城市三维模型的高大植被提取方法,其特征在于,S1中纹理空间的植被指数VI′=(2*G′-R′-B′)-(1.4*R′-G′),其中,G′=G/(R+G+B),R′=R/(R+G+B),B′=B/(R+G+B),R、G、B分别表示波段红、波段率、波段蓝。


3.根据权利要求1所述的一种基于倾斜摄影城市三维模型的高大植被提取方法,其特征在于,S1中计算单个三角面fi的植被指数VIi通过面积规则化植被指数:其中,Si=Square(fi)表示三角面fi的平面面积。


4.根据权利要求1所述的一种基于倾斜摄影城市三维模型的高大植被提取方法,其特征在于,S2中将植被指数相似的空间邻接三角面聚集成区域的具体方式是:迭代计算与当前区域相连接的三角面到区域中心的距离,判断其是否满足设定的判别条件,如果三角面满足条件,则加入到区域中,...

【专利技术属性】
技术研发人员:高云龙
申请(专利权)人:武汉大势智慧科技有限公司
类型:发明
国别省市:湖北;42

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1