首页> 中国专利> 基于格式塔心理学准则和多标签图割最小化的城市建筑可视化的方法

基于格式塔心理学准则和多标签图割最小化的城市建筑可视化的方法

摘要

为了使大规模城市建筑的可视化结果尽可能符合人们空间认知,本发明提出基于格式塔心理学准则和多标签图割最小化的城市建筑可视化的方法。首先,基于格式塔心理学邻近性、相似性和规则性等准则,实现城市建筑轮廓图的预聚类,为了解决建筑轮廓图聚类间的冲突,利用多标签图割建立能量优化函数,最优化求解计算城市建筑优化分组;通过不同层次的聚类阈值,实现了每一组内建筑轮廓的多细节层次综合;利用索引树SceneTree实现了多层次细节模型的快速索引,利用并行计算技术实现大规模城市场景的快速可视化;与已有方法相比,本发明对大规模复杂城市建筑模型具有很高的综合精度和可视化效率,在城市导航、空间位置服务等方面具有广阔的应用前景。

著录项

  • 公开/公告号CN104063893A

    专利类型发明专利

  • 公开/公告日2014-09-24

    原文格式PDF

  • 申请/专利权人 北京师范大学;

    申请/专利号CN201410146260.X

  • 发明设计人 张立强;王跃宾;张良;

    申请日2014-04-14

  • 分类号G06T17/00(20060101);G06T15/00(20110101);

  • 代理机构

  • 代理人

  • 地址 100875 北京市海淀区新街口外大街19号科技处

  • 入库时间 2023-12-17 01:34:31

法律信息

  • 法律状态公告日

    法律状态信息

    法律状态

  • 2019-04-05

    未缴年费专利权终止 IPC(主分类):G06T17/00 授权公告日:20170308 终止日期:20180414 申请日:20140414

    专利权的终止

  • 2017-03-08

    授权

    授权

  • 2014-10-22

    实质审查的生效 IPC(主分类):G06T17/00 申请日:20140414

    实质审查的生效

  • 2014-09-24

    公开

    公开

说明书

一、技术领域

本发明涉及基于格式塔心理学准则和多标签图割最小化的城市建筑可视化的方法,属于 空间信息技术领域。

二、背景技术

城市建筑模型的快速三维交互式可视化作为“数字城市”重要的组成部分,受到越来越 多的关注。可视化表达大场景的建筑实体,使之满足人们的视觉习惯,符合人们对城市的印 象,帮助人们准确而快速的获取城市空间信息,是地理空间信息认知和计算机视觉领域的重 要研究内容。为了加快大场景、高密度城市建筑的可视化速度,实现实时交互操作,需要对 建筑模型进行综合并建立细节层次模型(LOD)。相对于矢量数据,建筑物模型呈现出不同的 几何特性,每栋建筑由一个或几个多边形组成,而每个多边形包含很少,或不包含几何信息 冗余。同时城市是一个高度人性化的系统,建筑的综合和可视化表达不仅要考虑几何和拓扑 上的关系,更要考虑人的认知规律。

为了使城市建筑可视化的结果符合人们空间认知的原理,本发明提出了一种基于格式塔 心理学与多标签图割最小化实现大规模城市建筑综合与可视化的方法。首先,基于城市建筑 格式塔心理学准则的邻近性、相似性和规则性,实现城市建筑模型的预分组;为了解决建筑 轮廓格式塔分组的冲突,通过多标签图割建立能量函数,并使用标签代价进行最优求解,获 得了城市建筑最优化分组;通过不同层次的聚类阈值,利用索引树SceneTree实现多层次细 节模型的的快速合并与可视化。与已有方法相比,本发明对大规模复杂城市建筑模型具有较 高的综合精度和较高的可视化效率,并能保证标志性建筑在综合前后结构形态不会发生改变。 本发明在城市导航、空间位置服务等方面具有广阔的应用前景。

三、发明内容

1、目的:以往综合方法生成的三维城市场景难以兼顾建筑仿真和城市意象的可视化应用,本 发明提出了一种基于格式塔心理学与多标签图割最小化实现大规模城市建筑综合的方法。与 已有方法相比,本发明对大规模复杂城市建筑模型具有较高的综合精度和较高的可视化效率, 并能保证标志性建筑在综合前后结构形态不会发生改变。

2、技术方案:

基于格式塔心理学准则和多标签图割最小化的城市建筑可视化的方法,其特征在于,包 括如下步骤:

步骤一:建筑潜在格式塔聚类的提取

对城市建筑分类时,依据格式塔心理学邻近性、相似性和规则性等准则提取城市建筑空 间分布模式,以实现满足格式塔心理学规则的建筑聚类。对街区内的建筑集合 G={B0,B1,...,Bn}依据上述格式塔心理学准则进行聚类分析,把属于同一类型的建筑划分为一 类。

邻近性:用建筑轮廓多边形之间的邻近距离衡量建筑间的邻近性程度,把距离小于阈值 tp的建筑归为一类GP

GP={(Bi,Bj)|d(Bi,Bj)<tp}(1)

式(1)中,d(Bi,Bj)为建筑Bi,Bj之间的欧氏距离,tp为距离阈值。

相似性:结合建筑轮廓的两个特点提取相似性特征:第一,两栋建筑轮廓多边形面积大 小比值应在阈值和之内,得到聚类GS1;如果它们最小外包矩形长和宽有较大差异,比 较两栋建筑轮廓多边形的最小外包矩形,得到聚类GS2。最后得到相似性的一类GS

GS1={(Bi,Bj)|ts2<SBi/SBj<ts1}---(2)

式(2)中,SBi,SBj分别为建筑Bi,Bj平面轮廓图的面积。

GS2={(Bi,Bj)|R(Hi,Hj)+R(Wi,Wj)2>ts}---(3)

GS=GS1∩GS2(4)

其中,W和H表示建筑轮廓图的最小外包矩形的宽度和高度,并且

R(a,b)=a/bifa<bb/aotherwise---(5)

对于规则性特征,侧重于建筑线状方向一致性。建筑主朝向一致,主朝向的夹角应小于 阈值to,得到聚类GR1;另外,计算聚类中初始建筑B0的中心到当前建筑Bi中心的方向di, 与建筑Bi的中心到聚类中上一幢建筑Bi-1重心的方向di-1,对于di和di-1间夹角θdir,只要 θdirdir,则认为Bi满足线状方向一致性,得到聚类GR2

GR1={(Bi,Bj)|θ(OBi,OBj)<to}---(6)

GR=GR1∩GR2(7)

以B0为起始建筑,初始时G={B0},将满足上述所定义的格式塔心理学邻近性、相似性 和规则性等准则的建筑B1,...,Bn加入G中,将集合G中元素个数不小于3的保留下来。

步骤二:基于多标签图割最小化的建筑格式塔聚类的优化

使用上述格式塔心理学准则之后,每个街区内的建筑被划分成不同的类别,有的建筑同 时满足不同的准则,这就产生了冲突。基于多标签图割对格式塔心理学准则的上述特征集成, 通过下面的能量函数最小化解决上述冲突。

E(f)=ΣpPD(p,f)+Σp,qNVp,q+ΣlLhl·δl(f)---(8)

在上式(8)中,δl(f)=1p:fp=l0otherwise,p表示某一栋建筑,fp表示一种格式塔心理学准 则,Dp(fp)表示数据代价,Vp,q(fp,fq)表示平滑代价,用来表示标注代价。

数据代价定义了一个标签赋给一个集合元素的代价值。一个元素被赋予的标签代价值越 大,属于这类标签所定义的元素可能性就越小。为了便于衡量不同格式塔心理学准则所描述 的代价值,把每个元素所对应的代价值统一映射到0-1之间。

格式塔心理学邻近性准则的数据代价:一个元素p到所属集合Gp元素的最小距离。通过 邻近距离反映建筑物在二维平面上的空间关系,因此,赋给一元素邻近性的代价值可表示为:

Dp(fp)=mind(p,q)p,q∈Gp(9)

格式塔心理学相似性准则的数据代价:获得相似性分组Gs中元素的平均面积计算 每个元素p或q的面积与平均面积差值的绝对值,数据代价可以表示为:

Dp(fp)=|SP-Sarea|pGs---(10)

其中,表示集合元素的平均面积。

格式塔心理学规则性准则的数据代价:计算规则性分组GR中所有建筑元素轮廓多边形的 重心,通过最小二乘法把这些坐标线性拟合为一条直线,求得重心坐标到拟合直线的距离, 作为对应建筑规则性格式塔分组的代价值,具体可以表示为:

Dp(fp)=d(p,l)p∈GR(11)

其中,l表示元素重心所拟合的直线。

平滑代价描述了集合所有元素之间的空间关系。当元素p,q之间的平滑代价值越小,空间 关系就越紧密,属于同一类格式塔心理学准则的可能性就越大。本发明主要考虑建筑轮廓之 间的空间位置关系,所以,

Vp,q(fp,fq)=d(p,q)-1(12)

这里,d(p,q)表示建筑元素p,q轮廓多边形之间的邻近距离,平滑代价的值也映射到0到1 之间。

标注代价表示对过于复杂分组的惩罚,一般表示为

Flabelcost=ΣlLhl·δl(f)---(13)

这里,L表示不同的格式塔心理学准则分组,hl表示标签label的标注代价,而δl(f)是一个 指示性函数:

δl(f)=1p:fp=l0otherwise---(14)

所有的标注代价映射到0到1之间。

对于邻近性,SConvexHull表示集合所有元素所构成凸包的面积,表示集合所有元素的面积之和。

对于相似性,hl表示为表示集合所有建筑元素 轮廓多边形面积与平均面积差值的方差。

对于规则性,hl表示为,var(d(p,l))表示集合内所有建筑轮廓多边形 重心到拟合直线的距离的方差。

步骤三:城市建筑轮廓的综合

对于规则分组,由步骤一获得满足潜在格式塔心理学准则性聚类分组GR1,由步骤二获 得多标签图割优化分类结果GR2;求出GR1、GR2两者交集GR=GR1∩GR2,在GR中,若存在两 个分组子集合Ga和Gb,则保留Gb,删除Ga,GR即为待合并的建筑轮廓的集合。

类似的,计算潜在格式塔相似性聚类分组GS和邻近性分组GP

在GR、GS、GP中,分别把每一子集中的轮廓合并成一个多边形。下面给出了任一子集 Gi轮廓的合并过程:

(1)计算Gi中所有轮廓的凸包多边形CH

(2)CH中,在相邻且属于不同原始轮廓的两个顶点间,插入离这两顶点连线距离最近的 原始轮廓的凸顶点。

(3)CH中,在相邻且属于同一原始轮廓的两个顶点间,插入原始轮廓顶点。

(4)重复上述步骤,直到完成所有多边形顶点的遍历,形成最后的合并结果。

给出不同的格式塔聚类阈值tp、ts、to,能够控制GR、GS、GP中每一子集的轮廓数量, 从而生成不同细节层次的建筑综合结果。

在二维轮廓图基础上加入对应建筑高度值,生成建筑的三维模型,合并后建筑的高度Z 为:

Z=Σi=1nSiS×hi---(15)

这里,Si为建筑i的面积,为合并组建筑的平均面积,hi为建筑i的高度。

步骤四:大规模城市建筑可视化

为了便于存储和渲染综合后的建筑多细节层次模型,本发明构建一种树结构SceneTree管 理和索引不同的细节层次模型。SceneTree的结点(Scene node)管理每个城市街区建筑轮廓 的细节层次模型。在图2中,Block ID记录了Block的编号,LOD Level存储了相应Block内轮廓 综合后的细节层次,LOD(i)polygons是第i层次的分组多边形,包含了两部分的内容:顶点数 组(VertexArray)和顶点索引(Vertex Index),顶点数组存储了建筑三维空间点的坐标;顶点索引 存储了对应顶点的索引编号,ptr_father是指向父节点的指针,ptr_child指向子节点。

SceneTree遍历时,先对索引树进行遍历,根据当前视点和视线方向确定Block ID,由Block ID检索出落在视景体内的街区。根据视点到视景体内每个街区的距离,计算对应街区的轮廓 综合后的LOD level,确定GR、GS、GP中相应待合并的子集,视线动态的渲染与可视化。

随着视点的变化,显示的数据不断的进行更新。利用并行计算技术,将数据绘制和数据 装载分离。包含两个主要的进程:一个负责绘制和SceneTree的遍历(简称为T1),一个负责I/O 管理和数据预取(简称为T2)。这两个进程异步并行运行。T1根据当前视点计算当前可见视域, 并遍历SceneTree,选择合适的结点,向T2发出命令。当T2接收到命令后,将在当前视域范围 内,但不在内存中的结点装载入内存。当一块数据被装载入内存后,T1即可其进行绘制,而 不需等待其他数据的装载。这种方法,在不影响对城市建筑群进行空间表达的前提下,有效 地提高了渲染的速度。

3、优点及功效:本发明提出了基于格式塔心理学与多标签图割实现大规模城市建筑综合的方 法,以实现兼顾仿真建筑与城市意象的三维城市场景可视化。已有方法相比,本发明对大规 模复杂城市建筑模型具有较高的综合精度和较高的可视化效率,并能保证标志性建筑在综合 前后结构形态不会发生改变。能够自动可视化大规模城市地图,满足城市导航的需要。

四、附图说明

图1大规模城市建筑可视化的流程

图2SceneTree结点的数据结构

图3(a)-(d)本发明的结果图

图3(e)-(h)是Zhang等人方法(2012)的综合结果图

五、具体实施方式

本发明涉及基于格式塔心理学准则和多标签图割最小化的城市建筑可视化的方法,该方 法的具体步骤如下:

步骤一:建筑潜在格式塔聚类的提取

对城市建筑分类时,依据格式塔心理学邻近性、相似性和规则性等准则提取城市建筑空 间分布模式,以实现满足格式塔心理学规则的建筑聚类。对街区内的建筑集合 G={B0,B1,...,Bn}依据上述格式塔心理学准则进行聚类分析,把属于同一类型的建筑划分为一 类。

邻近性:用建筑轮廓多边形之间的邻近距离衡量建筑间的邻近性程度,把距离小于阈值 tp的建筑归为一类GP

GP={(Bi,Bj)|d(Bi,Bj)<tp}(1)

式(1)中,d(Bi,Bj)为建筑Bi,Bj之间的欧氏距离,tp为距离阈值。

相似性:结合建筑轮廓的两个特点提取相似性特征:第一,两栋建筑轮廓多边形面积大 小比值应在阈值和之内,得到聚类GS1;如果它们最小外包矩形长和宽有较大差异,比 较两栋建筑轮廓多边形的最小外包矩形,得到聚类GS2。最后得到相似性的一类GS

GS1={(Bi,Bj)|ts2<SBi/SBj<ts1}---(2)

式(2)中,SBi,SBj分别为建筑Bi,Bj平面轮廓图的面积。

GS2={(Bi,Bj)|R(Hi,Hj)+R(Wi,Wj)2>ts}---(3)

GS=GS1∩GS2(4)

其中,W和H表示建筑轮廓图的最小外包矩形的宽度和高度,并且

R(a,b)=a/bifa<bb/aotherwise---(5)

对于规则性特征,侧重于建筑线状方向一致性。建筑主朝向一致,主朝向的夹角应小于 阈值to,得到聚类GR1;另外,计算聚类中初始建筑B0的中心到当前建筑Bi中心的方向di, 与建筑Bi的中心到聚类中上一幢建筑Bi-1重心的方向di-1,对于di和di-1间夹角θdir,只要 θdirdir,则认为Bi满足线状方向一致性,得到聚类GR2

GR1={(Bi,Bj)|θ(OBi,OBj)<to}---(6)

GR=GR1∩GR2(7)

以B0为起始建筑,初始时G={B0},将满足上述所定义的格式塔心理学邻近性、相似性 和规则性等准则的建筑B1,...,Bn加入G中,将集合G中元素个数不小于3的保留下来。

步骤二:基于多标签图割最小化的建筑格式塔聚类的优化

使用上述格式塔心理学准则之后,每个街区内的建筑被划分成不同的类别,有的建筑同 时满足不同的准则,这就产生了冲突。基于多标签图割对格式塔心理学准则的上述特征集成, 通过下面的能量函数最小化解决上述冲突。

E(f)=ΣpPD(p,f)+Σp,qNVp,q+ΣlLhl·δl(f)---(8)

在上式(8)中,δl(f)=1p:fp=l0otherwise,p表示某一栋建筑,fp表示一种格式塔心理学准 则,Dp(fp)表示数据代价,Vp,q(fp,fq)表示平滑代价,用来表示标注代价。

数据代价定义了一个标签赋给一个集合元素的代价值。一个元素被赋予的标签代价值越 大,属于这类标签所定义的元素可能性就越小。为了便于衡量不同格式塔心理学准则所描述 的代价值,把每个元素所对应的代价值统一映射到0-1之间。

格式塔心理学邻近性准则的数据代价:一个元素p到所属集合Gp元素的最小距离。通过 邻近距离反映建筑物在二维平面上的空间关系,因此,赋给一元素邻近性的代价值可表示为:

Dp(fp)=mind(p,q)p,q∈Gp(9)

格式塔心理学相似性准则的数据代价:获得相似性分组Gs中元素的平均面积计算 每个元素p或q的面积与平均面积差值的绝对值,数据代价可以表示为:

Dp(fp)=|SP-Sarea|pGs---(10)

其中,表示集合元素的平均面积。

格式塔心理学规则性准则的数据代价:计算规则性分组GR中所有建筑元素轮廓多边形的 重心,通过最小二乘法把这些坐标线性拟合为一条直线,求得重心坐标到拟合直线的距离, 作为对应建筑规则性格式塔分组的代价值,具体可以表示为:

Dp(fp)=d(p,l)p∈GR(11)

其中,l表示元素重心所拟合的直线。

平滑代价描述了集合所有元素之间的空间关系。当元素p,q之间的平滑代价值越小,空间 关系就越紧密,属于同一类格式塔心理学准则的可能性就越大。本发明主要考虑建筑轮廓之 间的空间位置关系,所以,

Vp,q(fp,fq)=d(p,q)-1(12)

这里,d(p,q)表示建筑元素p,q轮廓多边形之间的邻近距离,平滑代价的值也映射到0到1 之间。

标注代价表示对过于复杂分组的惩罚,一般表示为

Flabelcost=ΣlLhl·δl(f)---(13)

这里,L表示不同的格式塔心理学准则分组,hl表示标签label的标注代价,而δl(f)是一个 指示性函数:

δl(f)=1p:fp=l0otherwise---(14)

所有的标注代价映射到0到1之间。

对于邻近性,SConvexHull表示集合所有元素所构成凸包的面积,表示集合所有元素的面积之和。

对于相似性,hl表示为表示集合所有建筑元素 轮廓多边形面积与平均面积差值的方差。

对于规则性,hl表示为,var(d(p,l))表示集合内所有建筑轮廓多边形 重心到拟合直线的距离的方差。

步骤三:城市建筑轮廓的综合

对于规则分组,由步骤一获得满足潜在格式塔心理学准则性聚类分组GR1,由步骤二获 得多标签图割优化分类结果GR2;求出GR1、GR2两者交集GR=GR1∩GR2,在GR中,若存在两 个分组子集合Ga和Gb,则保留Gb,删除Ga,GR即为待合并的建筑轮廓的集合。

类似的,计算潜在格式塔相似性聚类分组GS和邻近性分组GP

在GR、GS、GP中,分别把每一子集中的轮廓合并成一个多边形。下面给出了任一子集 Gi轮廓的合并过程:

(1)计算Gi中所有轮廓的凸包多边形CH

(2)CH中,在相邻且属于不同原始轮廓的两个顶点间,插入离这两顶点连线距离最近的 原始轮廓的凸顶点。

(3)CH中,在相邻且属于同一原始轮廓的两个顶点间,插入原始轮廓顶点。

(4)重复上述步骤,直到完成所有多边形顶点的遍历,形成最后的合并结果。

给出不同的格式塔聚类阈值tp、ts、to,能够控制GR、GS、GP中每一子集的轮廓数量, 从而生成不同细节层次的建筑综合结果。

在二维轮廓图基础上加入对应建筑高度值,生成建筑的三维模型,合并后建筑的高度Z 为:

Z=Σi=1nSiS×hi---(15)

这里,Si为建筑i的面积,为合并组建筑的平均面积,hi为建筑i的高度。

步骤四:大规模城市建筑可视化

为了便于存储和渲染综合后的建筑多细节层次模型,本发明构建一种树结构SceneTree管 理和索引不同的细节层次模型。SceneTree的结点(Scene node)管理每个城市街区建筑轮廓 的细节层次模型。在图2中,Block ID记录了Block的编号,LOD Level存储了相应Block内轮廓 综合后的细节层次,LOD(i)polygons是第i层次的分组多边形,包含了两部分的内容:顶点数 组(Vertex Array)和顶点索引(Vertex Index),顶点数组存储了建筑三维空间点的坐标;顶点索引 存储了对应顶点的索引编号,ptr_father是指向父节点的指针,ptr_child指向子节点。

SceneTree遍历时,先对索引树进行遍历,根据当前视点和视线方向确定Block ID,由Block ID检索出落在视景体内的街区。根据视点到视景体内每个街区的距离,计算对应街区的轮廓 综合后的LOD level,确定GR、GS、GP中相应待合并的子集,视线动态的渲染与可视化。

随着视点的变化,显示的数据不断的进行更新。利用并行计算技术,将数据绘制和数据 装载分离。包含两个主要的进程:一个负责绘制和SceneTree的遍历(简称为T1),一个负责I/O 管理和数据预取(简称为T2)。这两个进程异步并行运行。T1根据当前视点计算当前可见视域, 并遍历SceneTree,选择合适的结点,向T2发出命令。当T2接收到命令后,将在当前视域范围 内,但不在内存中的结点装载入内存。当一块数据被装载入内存后,T1即可其进行绘制,而 不需等待其他数据的装载。这种方法,在不影响对城市建筑群进行空间表达的前提下,有效 地提高了渲染的速度。

实施例1:

为了验证本发明的效率和效果,采用VC++2010和OpenGL编程实现了城市三维建 筑可视化系统。实验所用计算机配置如下:Windows7专业版操作系统,3.2GHz Intel(R) Core(TM)i5-3470CPU,4GB的内存,Nvidia Geforce GT620显卡。实验所用的数据为北京地 区城市建筑群16,056栋建筑,243,687个建筑轮廓多边形顶点。

用本发明方法和Zhang等人(Zhang,M.,Zhang,L.,Takis Mathiopoulos,P.,Xie,W.,Ding,Y. and Wang,H.,2012.A geometry and texture coupled flexible generalization of urban building models.ISPRS Journal of Photogrammetry and Remote Sensing,70,1-14)的方法分别对相同区域 城市建筑轮廓图进行聚类和合并。Zhang等人首先采用图像处理方法提取道路,然后通过缓 冲区聚类简化和综合城市建筑轮廓。该方法需要将矢量轮廓转换成栅格图像,在城市路网或 者城市建筑分布比较复杂时,Zhang等人的方法难以精确的提取道路信息,造成了建筑轮廓 聚类精度不高。本发明的方法考虑了城市街区建筑间邻近性、相似性和规则性等格式塔特征, 并基于多标签图割优化了聚类结果,较好反映了城市建筑空间分布特征,使聚类和合并后的 结果更符合人类的视觉感知。从图3可以看出,与Zhang等的方法相比,本发明的轮廓聚类 能够很好的顾及城市意象,合并后的结果保持了轮廓间原有的空间邻接和关联关系。

在时间效率方面,对于图3所示的城市建筑群,Zhang对轮廓进行聚类合并所用的时间为 1180秒,而本发明处理同样的数据花费时间仅为126秒,时效提高了近90%。

本发明采用SceneTree管理城市综合后的多细节层次模型。根据视点和视线方向计算视景 体,通过视景体裁剪确定落在视域范围内SceneTree的结点。由视点到每个结点中心的距离, 确定相应结点选择的细节层次模型以实现该街区内轮廓的可视化。离视点近的区域,轮廓保 持了较高的细节,远离视点的区域细节层次较粗。由于轮廓在聚类和合并过程中采用了并行 算法和SceneTree的快速检索机制,所以随着视点的变化,能够动态生成大规模场景,并保 持了较高的帧速率,很好的维持了格式塔心理学原则和城市意象,符合人的视觉习惯。帧速 率一般保持在60帧/秒左右。

去获取专利,查看全文>

相似文献

  • 专利
  • 中文文献
  • 外文文献
获取专利

客服邮箱:kefu@zhangqiaokeyan.com

京公网安备:11010802029741号 ICP备案号:京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有
  • 客服微信

  • 服务号