CN102589676A - 一种应用于声线追踪的室内空间剖分方法 - Google Patents

一种应用于声线追踪的室内空间剖分方法 Download PDF

Info

Publication number
CN102589676A
CN102589676A CN2011104354980A CN201110435498A CN102589676A CN 102589676 A CN102589676 A CN 102589676A CN 2011104354980 A CN2011104354980 A CN 2011104354980A CN 201110435498 A CN201110435498 A CN 201110435498A CN 102589676 A CN102589676 A CN 102589676A
Authority
CN
China
Prior art keywords
room
interior space
metope
door
sound ray
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.)
Granted
Application number
CN2011104354980A
Other languages
English (en)
Other versions
CN102589676B (zh
Inventor
蔡铭
陈志斌
李峰
余志�
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sun Yat Sen University
National Sun Yat Sen University
Original Assignee
National Sun Yat Sen University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National Sun Yat Sen University filed Critical National Sun Yat Sen University
Priority to CN201110435498.0A priority Critical patent/CN102589676B/zh
Publication of CN102589676A publication Critical patent/CN102589676A/zh
Application granted granted Critical
Publication of CN102589676B publication Critical patent/CN102589676B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明提出一种应用于声线追踪的室内空间剖分方法,该方法引入“房间”的概念,首先将室内空间的墙面要素存储到相应的数据结构中;然后根据空间的布局,将其剖分成若干个“房间”,并将墙面要素存储至对应的“房间”中;接着设置门窗、各“房间”之间的连通关系;最后将声场接收点放入在相应的“房间”里。本发明应用于声线跟踪在三维室内空间的声场计算,有效提高算法的效率。

Description

一种应用于声线追踪的室内空间剖分方法
技术领域
本发明涉及空间几何处理领域,特别的涉及一种应用于声线追踪的室内空间剖分方法。
背景技术
声线追踪方法作为一种几何声学模拟方法,原理是借鉴几何光学理论,假设声音沿直线传播,并忽略其波动特性,通过计算声音传播中能量的变化及反射到达的区域进行声场模拟。
该方法从声源向各方向发射“声粒子”,追踪它们的传播路径。“声粒子”因反射不断地失去能量,并按入射角等于反射角确定新的传播方向。可应用于三维空间中室内外声源对室内声场分布的计算。
在跟踪某一根声线时,确定声线与墙面的唯一碰撞点是关键,方法通常是先求出该声线与所有墙面的交点,然后找出离声线出发点距离最近的交点,再对该交点进行有效性判断(即判断该交点是否在实际墙面范围内)。如果是有效点,则该交点就是唯一碰撞点,如果不是,则找出离出发点距离第二近的交点,再进行有效性判断,如此循环,直至找出唯一碰撞点。常用的判断方法有夹角和检验法、叉积检验法和交点计数判别法3种,在三维空间内这3种方法都需要记录墙面的顶点集以及将三维空间中的交点和墙面做适当的投影,转化为判断二维空间中的点与多边形的位置关系,整个过程较为繁琐,计算量较大,效率不高。如果声线跟踪在一个封闭的凸多面体空间内进行,那么离声线出发点距离最近的交点必然落在实际墙面范围内,这时可以省略有效性判断的步骤。 
发明内容
本发明解决的技术问题是提供一种计算效率高、处理方法简单的应用于声线追踪的室内空间剖分方法。
为解决上述技术问题,本发明采用的技术方案是:提供一种应用于声线追踪的室内空间剖分方法,包含以下步骤:
(a)将室内空间的墙面要素存储成相应的数据结构;
(b)根据室内空间布局, 将室内空间剖分成1个以上的“房间”;
(c)存储各“房间”的墙面要素;
(d)存储各“房间”的门窗要素,以及各“房间”之间的连通关系;
(e)在“房间”中加入声场接收点。
  所述“房间”是指将室内空间(尤其针对凹多面体室内空间)按原有布局划分成的多个区域,其中每个区域为凸多面体空间。
进一步地,所述墙面要素由平面数学方程式中的参数表示。
进一步地,所述步骤(b)中,通过添加虚拟门,将室内空间剖分成1个以上的“房间”。
进一步地,门窗要素包括门窗所在“房间”编号以及所连通“房间”编号。每个墙面上的门窗属于且仅属于1个“房间”,连通且仅连通1个“房间”。
 
与现有技术相比,有益效果是:一是在跟踪时,剖分前的声线追踪需要遍历室内所有墙面进行求交点处理,而剖分后声线仅需对其所在的“房间”的墙面进行求交点处理,由全局搜索变为局部搜索,大大减少了声线所需要求交点的平面个数;二是室内空间进行“房间”剖分前,由于声线所在空间不一定是一个封闭的凸多面体空间,需要进行有效性判断,而剖分后每个“房间”均为凸多面体封闭空间,声线在“房间”内追踪,离声线出发点距离最近的交点必然落在实际墙面,这样就可以省去繁琐的有效性判断过程,大大提高了效率;三是将接收点按照所在“房间”进行归类,仅遍历搜索声线所在“房间”的所有接收点,对接收到声线的接收点进行计算,省去了剖分前对每根声线进行接收点的全局搜索判断过程。
附图说明
图1是除去顶面的室内空间示意图;
图2是进行以“房间”为单元的室内空间剖分后的示意图;
图3是若干个接收点在室内空间的位置示意图。
具体实施方式
为了表达方便,在下面的陈述中,图中所示的墙面用P加数字符号以及方向箭头表示,图中所示的“房间”用R加数字符号表示,图中所示的门用D加数字符号表示,图中所示的窗用W加数字符号表示,图中所示的虚拟门窗用V加数字符号表示。
在引入实施实例之前,先介绍门窗、墙面、“房间”的数据结构表:
表1:门窗要素的数据结构“Polygon”
数据类型 参数名称 备注
List<Point> Vertex 门窗的顶点链表
int now_room 所在“房间”的索引
int next_room 连通“房间”的索引
int flag 门窗的类别标志(门:1;窗:0;虚拟门:-1)
其中:Point为空间点的坐标。
表2:墙面要素的数据结构“Plane”
数据类型 参数名称 备注
double a 平面方程中的参数a
double b 平面方程中的参数b
double c 平面方程中的参数c
double d 平面方程中的参数d
List<Polygon> DW 该墙面上所包含的门窗链表
其中:a,b,c,d为数学上表示平面的方程:ax+by+cz+d=0的四个参数。
表3:“房间”要素的数据结构“Room”
数据类型 参数名称 备注
List<Plane> List_plane 组成该“房间”的墙面链表
List<Point> Recevers 该“房间”内的接收点
上述要素结构中,同种要素,数据结构相同,且要素之间为线性关系。
不同要素之间,关系如下所述:每个“房间”均由若干个墙面组成,包含“房间”里的接收点,而门窗则包含在所在的墙面上,对相应的“房间”进行连通。例如,“房间1”由墙面1、2、3、4、5、6组成,包含接收点p1、p2;而“门1”在墙面2上,其“now_room”属性为1,即其所在“房间”为“房间1”;其“next_room”属性为2,即其连通“房间”为“房间2”;当声线在“房间1”追踪时,仅需对“房间1”的6个组成墙面进行求交点处理,并对接收点p1、p2进行遍历计算”,若某次追踪求得的交点落在“墙面2”的“门1”上,则声线从“房间1”进入到“房间2”,再同理进行追踪。
下面结合附图所示的具体实例,来对剖分方法做具体说明。附图1是除去顶面的室内空间示意图,为所要进行剖分的室内空间。
整个剖分方法分为以下步骤:
步骤1:如附图2所示,将室内空间的所有墙面存入相应的数据结构数组中,如P0表示底面方向朝下的墙面,P1表示底面方向朝上的墙面,而未在图中显示的顶面,其方向朝下的墙面为P18,方向朝上的墙面为P19。其余各墙面的表示如附图2所示。
步骤2:根据室内空间布局, 将室内空间剖分成若干“房间”。如附图2所示,在厅的右侧添加多边形虚拟门V0,这样就将室内空间分成了4个“房间”:R1、R2、R3、R4,而室外空间则定义为R0。
步骤3:存储各“房间”的墙面要素;如附图2所示, “房间”R1由墙面P1、P3、P4、P14、P16、P18组成,“房间”R2由墙面P1、P4、P6、P16、P15、P18组成,“房间”R3由墙面P1、P3、P17、P6、P12、P18组成,“房间”R4由墙面P1、P10、P13、P6、P8、P18组成;而室外空间R0由墙面P0、P2、P5、P7、P9、P11、P13、P19组成,
步骤4:存储各“房间”的门窗要素,以及各“房间”之间的连通关系。如附图2所示,把窗户W0同时添加至墙面P2和P3上,不同的是,P2上的门窗W0所在的“房间”为R0,连通的“房间”为R1,即now_room=0,next_room=1;而P3上的门窗W0所在的“房间”为R1,连通的“房间”为R0,即now_room=1,next_room=0。由于W0是窗户,所以其flag属性为1。同理,其它门窗亦同样表示。
步骤5:在“房间”中加入声场接收点。如附图3所示,将接收点p0,p1,p2,p3,p4分别添加至“房间”R1,R2,R3,R4和室外空间R0的属性Recevers中。
本发明结合室内空间的布局进行以“房间”为单元的剖分后,当声线在某一“房间”追踪时,只需要与组成该“房间”的墙面求交点,离声线出发点距离最近的交点即为有效交点;判断该交点是否落在该“房间”的门窗上,若是,则进入对应的门窗的连通“房间”追踪;在每次追踪的过程中,仅遍历搜索该“房间”的所有接收点,对接收到声线的接收点进行计算;最后,当声线能量衰减至一定大小或声线从室内射到室外,即R0时,结束追踪。

Claims (4)

1.一种应用于声线追踪的室内空间剖分方法,其特征在于,包含以下步骤:
(a)将室内空间的墙面要素存储成相应的数据结构;
(b)根据室内空间布局,将室内空间剖分成1个以上的“房间”;
(c)存储各“房间”的墙面要素;
(d)存储各“房间”的门窗要素,以及各“房间”之间的连通关系;
(e)在“房间”中加入声场接收点。
2.根据权利要求1所述的室内空间剖分方法,其特征在于,所述墙面要素由平面数学方程式中的参数表示。
3.根据权利要求1所述的室内空间剖分方法,其特征在于:所述步骤(b)中,通过添加虚拟门,将室内空间剖分成1个以上的“房间”。
4.根据权利要求1所述的室内空间剖分方法,其特征在于:所述步骤(d)中,门窗要素包括门窗所在“房间”编号以及所连通“房间”编号。
CN201110435498.0A 2011-12-21 2011-12-21 一种基于室内空间剖分的声线追踪方法 Active CN102589676B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110435498.0A CN102589676B (zh) 2011-12-21 2011-12-21 一种基于室内空间剖分的声线追踪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110435498.0A CN102589676B (zh) 2011-12-21 2011-12-21 一种基于室内空间剖分的声线追踪方法

Publications (2)

Publication Number Publication Date
CN102589676A true CN102589676A (zh) 2012-07-18
CN102589676B CN102589676B (zh) 2015-08-12

Family

ID=46478642

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110435498.0A Active CN102589676B (zh) 2011-12-21 2011-12-21 一种基于室内空间剖分的声线追踪方法

Country Status (1)

Country Link
CN (1) CN102589676B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103778657A (zh) * 2014-02-28 2014-05-07 中山大学 一种基于空间剖分的声线束追踪方法
CN106096321A (zh) * 2016-07-20 2016-11-09 中山大学 结合声线跟踪法与声束跟踪法的室内外声音传播模拟方法
WO2023202551A1 (zh) * 2022-04-19 2023-10-26 北京字跳网络技术有限公司 声透射方法、装置和非易失性计算机可读存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60173423A (ja) * 1984-02-20 1985-09-06 Matsushita Electric Ind Co Ltd 音場解析法
JPH01248083A (ja) * 1988-03-29 1989-10-03 Toshiba Corp 波動源探査装置
JPH09166482A (ja) * 1995-12-14 1997-06-24 Nippon Hoso Kyokai <Nhk> 高精度音線追跡装置および高精度音線追跡方法
US5752167A (en) * 1995-01-23 1998-05-12 Advantest Corporation Radio propagation simulation method, wave field strength inference method and three-dimensional delay spread inference method
JP2004077352A (ja) * 2002-08-21 2004-03-11 Yamaha Corp 幾何音響シミュレーション装置、幾何音響シミュレーション方法及びプログラム
JP2006208273A (ja) * 2005-01-31 2006-08-10 Asahi Kasei Homes Kk 音場シミュレーションシステム

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60173423A (ja) * 1984-02-20 1985-09-06 Matsushita Electric Ind Co Ltd 音場解析法
JPH01248083A (ja) * 1988-03-29 1989-10-03 Toshiba Corp 波動源探査装置
US5752167A (en) * 1995-01-23 1998-05-12 Advantest Corporation Radio propagation simulation method, wave field strength inference method and three-dimensional delay spread inference method
JPH09166482A (ja) * 1995-12-14 1997-06-24 Nippon Hoso Kyokai <Nhk> 高精度音線追跡装置および高精度音線追跡方法
JP2004077352A (ja) * 2002-08-21 2004-03-11 Yamaha Corp 幾何音響シミュレーション装置、幾何音響シミュレーション方法及びプログラム
JP3900041B2 (ja) * 2002-08-21 2007-04-04 ヤマハ株式会社 幾何音響シミュレーション装置、幾何音響シミュレーション方法及びプログラム
JP2006208273A (ja) * 2005-01-31 2006-08-10 Asahi Kasei Homes Kk 音場シミュレーションシステム

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A. KROKSTAD 等: "Calculating the acoustical room response by the use of a ray tracing technique", 《J. SOUND VIB》 *
曾向阳 等: "计算房间脉冲响应的经典方法-声线跟踪法研究", 《电声技术》 *
王海涛 等: "复杂场景中的多重声反射跟踪算法研究", 《电声基础》 *
谭同德 等: "基于声线跟踪法的室外声场仿真", 《计算机工程与应用》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103778657A (zh) * 2014-02-28 2014-05-07 中山大学 一种基于空间剖分的声线束追踪方法
CN103778657B (zh) * 2014-02-28 2017-01-25 中山大学 一种基于空间剖分的声线束追踪方法
CN106096321A (zh) * 2016-07-20 2016-11-09 中山大学 结合声线跟踪法与声束跟踪法的室内外声音传播模拟方法
CN106096321B (zh) * 2016-07-20 2019-02-01 中山大学 结合声线跟踪法与声束跟踪法的室内外声音传播模拟方法
WO2023202551A1 (zh) * 2022-04-19 2023-10-26 北京字跳网络技术有限公司 声透射方法、装置和非易失性计算机可读存储介质

Also Published As

Publication number Publication date
CN102589676B (zh) 2015-08-12

Similar Documents

Publication Publication Date Title
Aartsen et al. Search for time-independent neutrino emission from astrophysical sources with 3 yr of IceCube data
CN103247041B (zh) 一种基于局部采样的多几何特征点云数据的分割方法
WO2008067483A1 (en) Ray tracing a three dimensional scene using a grid
US20180060463A1 (en) Hybrid Monte Carlo and Deterministic Particle Transport Method Based on Transition Area
CN102589676A (zh) 一种应用于声线追踪的室内空间剖分方法
CN105389850A (zh) 一种大规模三维场景的新型可见性生成方法
CN101533102B (zh) 二维复杂结构三角网射线追踪全局方法
CN106096321B (zh) 结合声线跟踪法与声束跟踪法的室内外声音传播模拟方法
Dong et al. An accelerated SBR for EM scattering from the electrically large complex objects
Samsing et al. Active galactic nuclei as factories for eccentric black hole mergers
CN104570080A (zh) 一种海量数据叠前逆时偏移多gpu卡协同快速计算方法
CN103698810A (zh) 混合网最小走时射线追踪层析成像方法
CN106250101A (zh) 基于MapReduce的叠前偏移并行处理方法和装置
CN102722653A (zh) 一种基于MapReduce的射线跟踪加速算法
CN109218961A (zh) 一种基于虚拟节点的多站协同干扰定位方法及系统
Cao et al. Design for additive manufacturing: Optimization of piping network in compact system with enhanced path-finding approach
CN104361215A (zh) 一种基于标记信息的蒙特卡罗聚变堆重复结构处理方法
CN103778657B (zh) 一种基于空间剖分的声线束追踪方法
de Berg et al. Kinetic convex hulls and Delaunay triangulations in the black-box model
CN103049614B (zh) 一种控制智慧楼宇中声波运动轨迹的方法、装置
CN102930102A (zh) 一种微波部件二次电子倍增仿真中粒子合并方法
CN106651979A (zh) 基于投影的空间曲面四边形网格剖分方法
CN103391139A (zh) 一种电波传播损耗的快速预测方法
Saiguran et al. Evaluation of Markov chains to describe movements on tiling
Abachi et al. The D0 upgrade

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant