CN103778657B - 一种基于空间剖分的声线束追踪方法 - Google Patents
一种基于空间剖分的声线束追踪方法 Download PDFInfo
- Publication number
- CN103778657B CN103778657B CN201410071831.8A CN201410071831A CN103778657B CN 103778657 B CN103778657 B CN 103778657B CN 201410071831 A CN201410071831 A CN 201410071831A CN 103778657 B CN103778657 B CN 103778657B
- Authority
- CN
- China
- Prior art keywords
- bundle
- rays
- sound
- tetrahedron
- tetrahedral
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000005192 partition Methods 0.000 title abstract description 8
- 230000005540 biological transmission Effects 0.000 claims abstract description 21
- 238000005516 engineering process Methods 0.000 claims abstract description 8
- 230000008569 process Effects 0.000 claims description 13
- 238000013316 zoning Methods 0.000 claims description 4
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000012966 insertion method Methods 0.000 claims description 3
- 230000000644 propagated effect Effects 0.000 claims description 3
- 230000011514 reflex Effects 0.000 claims description 3
- 230000006854 communication Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 5
- 230000006872 improvement Effects 0.000 abstract description 2
- 230000004888 barrier function Effects 0.000 description 6
- 238000003780 insertion Methods 0.000 description 3
- 230000037431 insertion Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000000638 solvent extraction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Landscapes
- Image Generation (AREA)
Abstract
本发明公开一种基于空间剖分的声线束追踪方法,是对现有声线束追踪法技术上的改进。具体过程为:采用四面体剖分对声音传播的全空间进行四面体剖分形成多个四面体空间,多个四面体空间构成四面体网格,为剖分后的四面体空间编号,四面体剖分技术符合声线束传播的空间几何规律;对四面体网格进行约束处理,恢复障碍物边界条件,并存储四面体网格各个要素;在剖分生成四面体网格后,对四面体空间采用递归方式的构建声线束树形结构;采用声线束追踪法在声线束树形结构中搜索所有节点并生成声音传播路径;对生成的声音传播路径计算声音衰减。本发明应用于复杂障碍物情况下三维声场的计算,可有效的解决高密度障碍物大区域声音计算时间效率过低的问题。
Description
技术领域
本发明涉及室外声传播的技术领域,更具体地,涉及一种基于空间剖分的声线束追踪方法。
背景技术
声线束追踪法(Beam Tracing Method (BTM))是一种解决声波遇到障碍物反射和衍射问题的重要技术方法。算法追踪的是声线束(一簇声线),而并非是单条声线,通常用的比较多的是三棱锥形状的声线束。声源在某一时刻朝着某一方向发出声线束,当遇到障碍物时,声线束进行反射,声线沿新方向继续前进,直到声线能量低于设定的阈值时追踪完毕。声线束追踪法能达到很高的计算精度,声线束追踪法并不需要进行可见性判断,声线束本身就代表声源点能够到达的区域;高阶的反射声线束由低阶的声线束与障碍物可见平面部分细分而来。
当障碍物较多时,现有的声线束跟踪方法需要模拟声波向四周发散的过程,其计算量较大。
发明内容
本发明目的在于克服现有技术的缺点与不足, 提供一种计算效率高、处理方法简单的基于空间剖分的声线束追踪方法,用于众多障碍物遮挡后的接收点声音计算。
为解决上述技术问题,本发明采用的技术方案是:
一种基于空间剖分的声线束追踪方法,包含以下步骤:
(a)采用四面体剖分技术,对声音传播的全空间进行四面体剖分形成多个四面体空间,上述多个四面体空间构成四面体网格,为剖分后的四面体空间编号,上述四面体剖分技术符合声线束传播的空间几何规律;
(b)对所述四面体网格进行约束处理,恢复障碍物边界条件,并存储四面体网格各个要素;
(c)在剖分生成四面体网格后,对四面体空间采用递归方式的构建声线束树形结构;
(d)采用声线束追踪法在声线束树形结构中搜索所有节点并生成声音传播路径;对生成的声音传播路径计算声音衰减;
所述步骤(c)中采用递归方式构建声线束树形结构的具体方式为:采用逐点***的方法,对于任意的点,找到距离该约束最近的边,连接点和边构成一个平面,然后找到距离这个平面最近的点,形成一个初始四面体形,然后以此为基础向外扩展,直至所有点都被加入到构成的四面体网格中。
本发明首先将计算区域进行子空间划分,然后采用声线束跟踪法等几何声学办法,大大降低现有声线束追踪技术的计算量。发明提出将整个空间剖分成有序四面体网格的空间剖分技术,把需要全局搜索的声线束追踪转换为局部搜索,从而起到加速追踪声线束追踪的作用。在生成四面体网格后,采用递归的方式构建声线束树形结构,使得对每一对声源点和接收点,可以通过声线束追踪法在生成的声线束树形结构中搜索所有节点并生成声音传播路径。采用递归方式构建声线束树形结构可以有效的提高声线束追踪的速度;通过该技术过程得到声音传播路径后,应用声音传播衰减模型对接收点声音进行计算。
本发明是采用四面体剖分、声线束追踪和路径生成三个技术手段解决声音传播计算问题。
更进一步的,所述步骤(a)中,通过采用逐点***方法,将全空间剖分成多个四面体空间。
更进一步的,所述存储四面体网格的各个要素包括点、线和三角形、四面体空间。
更进一步的,所述声线束树形结构由节点元素构成,节点要素包括当前四面体编号、当前面编号、当前边编号、反射次数、声源点坐标和衍射点位置。
更进一步的,所述步骤(d)中,声线束追踪法在所述树形结构中进行传播和跟踪,记录声线束传播约束条件的编号并根据传播过程进行更新以生成声音传播路径。
更进一步的,所述步骤(d)中,声音传播路径生成过程中,对声线束追踪中经过的约束边或面以及影响到的四面体进行记录,用于寻找目标节点并生成一条从声源点到接收点的三维有效路径。
更进一步的,所述声线束追踪法的追踪过程在满足以下四个条件之一时停止:
1)声线束的长度达到自定义的长度阈值;
2)声线束的反射次数超过自定义的反射阈值;
3)声线束的衍射次数超过自定义的衍射阈值;
4)声线束到达计算区域的边界阈值。
更进一步的,所述声音传播路径的声线包括四类:直射声线、反射声线、衍射声线和包含反射和衍射的复合声线。
与现有技术相比,本发明有益效果是:在空间剖分技术下进行声线束追踪,使得声音在大区域建筑物群之间声音衰减计算时,其时间效率将会大大提高,其根本原因在于本发明提出的声线束复杂空间下的路径提取上。其具体表现为:
(1)在追踪时,现有的声线束追踪方法需要遍历全部空间,而应用空间剖分技术后,声线束仅需对其所经过的四面体进行求交点以及反射和衍射等处理,由全局搜索变为局部搜索,大大减少了声线束所需要求交点的平面个数;
(2)声线束树形结构体技术的应用,树形结构建立前,声线束所在空间是任意的,需要对全闭合空间内所有的四面体空间进行标记并判断声线束是否经过此四面体空间,而空间剖分后每个四面体空间均具有独立的编号,并可以通过树状结构和其他节点联系,树形结构节点存储的数据包括:当前四面体编号、当前边编号、当前面编号、反射次数、声源点坐标和衍射点位置,应用树形结构体技术可以引导路径是按照正确的路径生成的;
(3)路径生成技术包括了四面体空间标记的过程,形成一对声源点和接收点之间的传播路径后,四面体网格就标记了经过的节点信息,下一对声源点和接收点路径中重复的部分直接可以获取,省去了部分路径重合带来的重复追踪带来的时间和浪费。
附图说明
图1是声音具体传播路径算法示意图。
图2空间剖分过程示意图。
图3声线束追踪过程和结构树示意图。
具体实施方式
下面结合附图对本发明做进一步的描述,但本发明的实施方式并不限于此。
如图1所示,本发明提供一种基于空间剖分的声线束追踪下声衰减计算方法。通过以下步骤实现:
步骤1:如附图2(a)所示,***边界4个点,形成***边界。***四个边界点O、P、Q和R包围所有的散点集,形成边界封闭曲线OPQR。根据边界及内部点生成四面体。依次对所有的边按顺序从外到里循环上述过程,就可对该区域生成一个四面体网。
步骤2:如附图2(b)所示,为空间剖分后形成的四面体编号,编号顺序从外到里。
步骤3:如附图2(c)所示,依次***障碍物边,与障碍物相交的边所在的四面体形即为影响区域,在影响区域内重复步骤2,直到生成最优空间剖分结构,并将所有的点、线和三角形、四面体分别存入数据库。
步骤4:如附图2所示,采用逐点***的方法,对于任意的点,找到距离该约束最近的边,连接点和边构成一个平面,然后找到距离这个平面最近的点,形成一个初始四面体形,然后以此为基础向外扩展,直至所有点都被加入到构成的四面体网格中。
步骤5:如附图3(a)所示,枚举声线束树形结构体中所有的节点,找到符合追踪条件的目标节点直至追踪结束,生产声音传播路径。在路径生成的过程中,需要对声线束追踪中经过的约束边或面以及影响到的四面体进行记录。
步骤6:如附图3(b)所示,找到所有目标节点,生成一条从声源点到接收点的三维有效路径。
步骤7:对声音传播过程在生成的路径上进行计算。所有的有效路径中的声线可以分成4类:直射声线、反射声线、衍射声线和包含反射和衍射的复合声线,所有声线总的声能贡献即为接收点声能。
本发明结合空间剖分对声线束过程进行有目的性追踪,当声线束在某一四面体空间里进行追踪时,只需要与组成该四面体空间的面以及相邻的四面体空间进行路径判断以及数据交换,该四面体空间实时记录追踪信息。追踪过程在满足以下四个条件之一停止:1)声线束的长度达到用户定义的阈值,简称为长度阈值;2)声线束的反射次数超过用户自定义的阈值,简称为反射阈值;3)声线束的衍射次数超过用户自定义的阈值,简称为衍射阈值;4)声线束到达计算区域的边界,简称为边界阈值。
以上所述的本发明的实施方式,并不构成对本发明保护范围的限定。任何在本发明的精神原则之内所作出的修改、等同替换和改进等,均应包含在本发明的权利要求保护范围之内。
Claims (7)
1.一种基于空间剖分的声线束追踪方法,其特征在于,包含以下步骤:
(a)采用四面体剖分技术,对声音传播的全空间进行四面体剖分形成多个四面体空间,上述多个四面体空间构成四面体网格,为剖分后的四面体空间编号,上述四面体剖分技术符合声线束传播的空间几何规律;
(b)对所述四面体网格进行约束处理,恢复障碍物边界条件,并存储四面体网格各个要素;
(c)在剖分生成四面体网格后,对四面体空间采用递归方式的构建声线束树形结构;
(d)采用声线束追踪法在声线束树形结构中搜索所有节点并生成声音传播路径;对生成的声音传播路径计算声音衰减;
所述步骤(c)中采用递归方式构建声线束树形结构的具体方式为:采用逐点***的方法,对于任意的点,找到距离该点最近的边,连接点和边构成一个平面,然后找到距离这个平面最近的点,形成一个初始四面体形,然后以此为基础向外扩展,直至所有点都被加入到构成的四面体网格中。
2.根据权利要求1所述的基于空间剖分的声线束追踪方法,其特征在于,存储的所述四面体网格的各个要素包括点、线和三角形、四面体空间。
3.根据权利要求1所述的基于空间剖分的声线束追踪方法,其特征在于,所述声线束树形结构由节点要素构成,节点要素包括当前四面体编号、当前面编号、当前边编号、反射次数、声源点坐标和衍射点位置。
4.根据权利要求1所述的基于空间剖分的声线束追踪方法,其特征在于,所述步骤(d)中,声线束追踪法在所述声线树形结构中进行传播和跟踪,记录声线束传播所影响的节点要素,并根据传播过程进行更新以生成声音传播路径。
5.根据权利要求1所述的基于空间剖分的声线束追踪方法,其特征在于,所述步骤(d)中,声音传播路径生成过程中,对声线束追踪中经过的约束边或面以及影响到的四面体进行记录,用于寻找目标节点并生成一条从声源点到接收点的三维有效路径。
6.根据权利要求1所述的基于空间剖分的声线束追踪方法,其特征在于,所述声线束追踪法的追踪过程在满足以下四个条件之一时停止:
1)声线束的长度达到自定义的长度阈值;
2)声线束的反射次数超过自定义的反射阈值;
3)声线束的衍射次数超过自定义的衍射阈值;
4)声线束到达计算区域的边界阈值。
7.根据权利要求1所述的基于空间剖分的声线束追踪方法,其特征在于,所述声音传播路径的声线包括四类:直射声线、反射声线、衍射声线和包含反射和衍射的复合声线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410071831.8A CN103778657B (zh) | 2014-02-28 | 2014-02-28 | 一种基于空间剖分的声线束追踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410071831.8A CN103778657B (zh) | 2014-02-28 | 2014-02-28 | 一种基于空间剖分的声线束追踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103778657A CN103778657A (zh) | 2014-05-07 |
CN103778657B true CN103778657B (zh) | 2017-01-25 |
Family
ID=50570852
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410071831.8A Expired - Fee Related CN103778657B (zh) | 2014-02-28 | 2014-02-28 | 一种基于空间剖分的声线束追踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103778657B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105451151B (zh) * | 2014-08-29 | 2018-09-21 | 华为技术有限公司 | 一种处理声音信号的方法及装置 |
CN106096321B (zh) * | 2016-07-20 | 2019-02-01 | 中山大学 | 结合声线跟踪法与声束跟踪法的室内外声音传播模拟方法 |
CN109597032B (zh) * | 2019-01-31 | 2020-02-14 | 中国科学院深海科学与工程研究所 | 一种水下声学定位通信方法 |
CN114003855B (zh) * | 2021-10-20 | 2024-06-11 | 国网重庆市电力公司电力科学研究院 | 基于声线跟踪理论的声波障碍物绕射模拟方法及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6751322B1 (en) * | 1997-10-03 | 2004-06-15 | Lucent Technologies Inc. | Acoustic modeling system and method using pre-computed data structures for beam tracing and path generation |
CN101140612A (zh) * | 2006-09-05 | 2008-03-12 | 刘礼白 | 声波识别(swid)技术 |
CN102589676A (zh) * | 2011-12-21 | 2012-07-18 | 中山大学 | 一种应用于声线追踪的室内空间剖分方法 |
-
2014
- 2014-02-28 CN CN201410071831.8A patent/CN103778657B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6751322B1 (en) * | 1997-10-03 | 2004-06-15 | Lucent Technologies Inc. | Acoustic modeling system and method using pre-computed data structures for beam tracing and path generation |
CN101140612A (zh) * | 2006-09-05 | 2008-03-12 | 刘礼白 | 声波识别(swid)技术 |
CN102589676A (zh) * | 2011-12-21 | 2012-07-18 | 中山大学 | 一种应用于声线追踪的室内空间剖分方法 |
Non-Patent Citations (3)
Title |
---|
"The adaptive beam-tracing algorithm";I.A.Drumm,et al.;《The Journal of the Acoustical Society of America》;20000331;第107卷(第3期);全文 * |
"三维声线追踪的正三棱锥前向伸展算法";姜薇等;《声学学报》;20050930;第30卷(第5期);全文 * |
"基于参数模型两种剖分方式的射线追踪比较";成谷等;《中山大学学报(自然科学版)》;20080930;第47卷(第5期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103778657A (zh) | 2014-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102819865B (zh) | 一种大地电磁三维地质结构模型的建模方法 | |
KR101697238B1 (ko) | 영상 처리 장치 및 방법 | |
CN103778657B (zh) | 一种基于空间剖分的声线束追踪方法 | |
Kallmann et al. | Navigation meshes and real-time dynamic planning for virtual worlds | |
CN105844709A (zh) | 复杂河道地形流域洪水演进虚拟仿真的淹没线追踪方法 | |
EA201070615A8 (ru) | Формирование геологической модели | |
CN102057368B (zh) | 利用最大连续场在三维体积模型中分布性质 | |
CN102609982B (zh) | 空间地质数据非结构化模式的拓扑发现方法 | |
CN106600697A (zh) | 一种面向3d场景的导航网格地图表示方法 | |
CN105093319B (zh) | 基于三维地震数据的地面微地震静校正方法 | |
CN106845074A (zh) | 建立六边形格网模型的方法、洪水推演模拟方法及其*** | |
CN103105624A (zh) | 基于数据库技术的纵横波时差定位方法 | |
CN102446122A (zh) | 一种基于包围盒树的碰撞检测方法 | |
KR20130133778A (ko) | 트레이싱 정보에 기초한 충돌 회피 방식을 사용하여 캐릭터들을 애니메이팅하기 위한 방법 | |
Bosurgi et al. | A PSO highway alignment optimization algorithm considering environmental constraints. | |
CN104360396B (zh) | 一种海上井间tti介质三种初至波走时层析成像方法 | |
CN106096321B (zh) | 结合声线跟踪法与声束跟踪法的室内外声音传播模拟方法 | |
KR100903777B1 (ko) | 3차원 광선 추적 시스템에서 적응형 다중 계층 케이디-트리구조 알고리즘을 이용한 광선 추적 방법 및 장치 | |
US20130033500A1 (en) | Dynamic collision avoidance for crowd simulation over structured paths that intersect at waypoints | |
Chentanez et al. | Fast grid-free surface tracking | |
CN104376050A (zh) | 一种在三维场景动态调度并显示海量矢量数据的方法 | |
CN103698810A (zh) | 混合网最小走时射线追踪层析成像方法 | |
CN110349262A (zh) | 一种基于bim的隧道超前地质预报三维建模方法及*** | |
CN106373192B (zh) | 一种非拓扑一致性三维网格块体追踪算法 | |
CN103886129B (zh) | 将测井数据离散到储层网格模型的方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information |
Inventor after: Wang Haibo Inventor after: Cai Ming Inventor after: Yu Zhi Inventor before: Cai Ming Inventor before: Wang Haibo Inventor before: Yu Zhi |
|
COR | Change of bibliographic data | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170125 |
|
CF01 | Termination of patent right due to non-payment of annual fee |