CN103986536B - 一种正向射线追踪等分隔射线束的生成方法 - Google Patents
一种正向射线追踪等分隔射线束的生成方法 Download PDFInfo
- Publication number
- CN103986536B CN103986536B CN201410234745.4A CN201410234745A CN103986536B CN 103986536 B CN103986536 B CN 103986536B CN 201410234745 A CN201410234745 A CN 201410234745A CN 103986536 B CN103986536 B CN 103986536B
- Authority
- CN
- China
- Prior art keywords
- triangle
- generation method
- source
- regular polygon
- center
- 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
Landscapes
- Image Generation (AREA)
Abstract
本发明涉及一种正向射线追踪等分隔射线束的生成方法,该方法包括:将源置于非正多面体中心;将所述非正多面体每个面的中心与所述每个面的顶点连接;对所述每个面上的三角形进行几何分形,直至所述源与每个三角形顶点形成的射线束符合预定条件。本发明利用非多面体来近似生成等分割射线束,能够更快速地实现正向射线追踪,并能够减少等分割射线束之间的角度误差。
Description
技术领域
本发明涉及正向射线跟踪,尤其涉及一种正向射线追踪等分隔射线束的生成方法。
背景技术
正向射线跟踪是一种广泛应用的预测无线电波传播特性的技术,它的原理是由源点向空间均匀发出大量的射线束,然后跟踪每条射线的路径并用接收球判定其对接收点的场强是否有贡献。
在正向射线跟踪的过程中,产生等间隔的射线束是整个过程中非常关键的关键。在现有技术中,一般通过正多面体来产生等间隔的射线束。现有技术中,仅存在五种正多面体,而正二十面体为面数最多,所以通过正二十面体来产生等间隔的射线束。
正二十面体如图1所示,每个面均为等边三角形,为了提高射线束密度,需要将每个面进一步细分为更多的小三角形。具体方法是将每条边N等分,然后分别过各个分点作其他边的平行线,从而每个面得到N2个小三角形。将源点置于正二十面体的中心,连接源点和各个三角形的顶点便得到10N2+2条射线束。图2所示为N=4时的正二十面体的一个面的细分情况。以AB边四等分为例,可得到M1,M2,M3三个分点,然后过这三个点作关于其它两个边的平行线,对其它两个边采用同样的方法,最后便得到图2的结果。
另一种方法较简便,用如下关于角度的循环公式来近似生成:
θ=iδi=1…int(180/δ)
φ=jδ/sinθj=1…int(360sinθ/δ)
其中δ为射线束夹角,θ为各射线束与z轴夹角,φ为射线束与x轴夹角,根据不同的(θ,φ)来确定每一条射线束的方向。
但这种方法同样不能得到真正的等分隔射线束,在求取小角度间隔的射线束时会产生较大的误差。
发明内容
为了解决上述的技术问题,提供了一种正向射线追踪等分隔射线束的生成方法,该方法包括:将源置于非正多面体中心;将所述非正多面体每个面的中心与所述每个面的顶点连接;对所述每个面上的三角形进行几何分形,直至所述源与每个三角形顶点形成的射线束符合预定条件。
优选地,所述非正多面体为C60分子结构多面体。
优选地,在进行几何分形时,通过连接三角形每条边的中点得到4个分形后的三角形。
优选地,所述预定条件为射线束角度间隔达到预定值。
本发明利用非多面体来近似生成等分割射线束。对于现有的角度循环公式法,其不能产生等分割,无论在射线较少或较多时均有很大误差,在靠近z轴时射线密集,其他部分稀疏,其缺点在发射机旋转时体现得尤为明显,发射机旋转会导致空间射线的整体稀疏和密集程度区域分布产生很大变化,加剧了射线追踪的误差;对于基于正二十面体的射线束产生方法,其空间对称性很好,但是当射线较多时,细剖分每个面三角形的产生射线的方法会使误差逐渐增大;本发明基于C60分子的射线束产生方法,将基础面增加至32个,在剖分每个面来加密射线时有更好的均匀性,尤其是在射线束要求很多时,与前两种已有方法相比更能体现出优势,而且在相同的分割率条件下,本发明相对于正二十面体的射线束产生方法提高为约9倍的细分效率。因此,对于射线密度要求较高的射线追踪算法来说,本发明相对于现有技术能够更快速地实现正向射线追踪,并能够减少等分割射线束之间的角度误差。
附图说明
图1为正二十面体空间结构示意图;
图2为为N=4时的正二十面体的一个面的细分情况;
图3为C60分子多面体的空间结构示意图;
图4为基于C60分子空间结构的近似等分隔射线束生成的流程图;
图5为谢尔宾斯基三角形分形方法;
图6和图7分别为正五边形面和正六边形面的细分方法示意图。
具体实施方式
下面结合附图,对本发明做进一步的详细描述。
为了得到分割效果更好的射线群束,本发明期望通过非正多面体来实现等分割射线束。优选地,本发明通过C60分子多面体来实现等分割射线束的近似生成。由于非正多面体的所有面并非完全相同,因此所生成的射线束之间并非完全等分,但可以满足正向射线跟踪的要求。
C60分子如图3所示,有60个顶点、90条棱和32个面,其中12个为正五边形,20个为正六边形,是一种高度对称的完美分子结构。本发明中,通过将源置于C60分子多面体的中心,并利用每个面的中心和顶点形成多个三角形,之后对每个三角形进行分形,从而得到符合要求的射线束。在三角形分割率N确定的情况下,传统的基于正二十面体方法产生射线的条数为10N2+2,而本发明的基于C60分子结构方法可以生成90N2+2条射线(N=2k,k≥1)。
作为示例,本发明提供的基于C60分子多面体的生成等分隔射线束的方法如图4中的流程图所示。为了提高细分效率,在将每个面细分为小三角形的过程中引入了改进的谢尔宾斯基三角形分形方法。谢尔宾斯基三角形是由波兰数学家谢尔宾斯基在1915年提出的,方法如图5所示。对于一个三角形,连接各个边的中点从而将其分为四个小三角形,然后挖去中间的一个小三角形,将剩下的三个小三角形按照上述方法继续分割并舍弃中间的,不断重复分割与舍弃的过程便得到谢尔宾斯基三角形。本发明中采用该分割方法,但并不舍弃分割后的中间的小三角形,即在每个过程中对所有的三角形进行分割。
作为示例,本发明提供的正向射线追踪等分割线束的生成方法的具体步骤如下:
步骤一、将辐射源放置在C60分子结构模型的中心,连接源点和模型的60个顶点,得到初始的60条射线束。
步骤二、C60分子结构模型共有32个面,连接源点和每个面的中心,得到32条补充射线束。
步骤三、对于12个正五边形面,将每个面上的中心和五个顶点相连接得到五个三角形,如图6所示;对于20个正六边形面,将每个面的中心和六个顶点连接得到六个正三角形,如图7所示。
步骤四、在每个三角形中,连接源和每条边的中点得到细分射线束。如图6所示,对于ΔAOE,A1,B1,C1分别为三条边的中点,与源点连接便得到1级细分射线束;对于公共边,如AO既为ΔAOE的边又为ΔAOB的边,则连接源点和AO边中点B1的射线束只取一次。设源点为O',即如果在ΔAOE中已经得到经过源点与三条边中点的射线O'A1,O'B1与O'C1,则在ΔAOB中则只连接源与AB,OB的中点得到射线O'M,O'N即可,而不需再连接源与AO的中点B1,从而在进行正向射线跟踪时避免对射线O'B1重复跟踪。
步骤五、采用谢尔宾斯基三角形分割方法,依次连接三角形各边的中点,从而将一个三角形细分为四个小三角形。如图3,以ΔAOE为例,分别连接三个边中点A1B1,B1C1,A1C1,则将ΔAOE划分为四个小三角形ΔAA1B1,ΔA1EC1,ΔA1B1C1和ΔC1OB1。
步骤六、重复执行步骤四、五k次,直到射线束角度间隔达到预定要求为止。每次重复执行步骤四、五则将每个三角形细分一次,最终得到符合角度间隔要求的k级细分射线束。作为示例,图6中示出了生成2级细分射线束的示例三角形A2B2C2,同样地,可以对三角形AA1B1以及三角形B1C1O进行分形,也得到2级细分射线束。
本发明利用非多面体来近似生成等分割射线束,能够更快速地实现正向射线追踪,并能够减少等分割射线束之间的角度误差。
本领域的技术人员在不脱离权利要求书确定的本发明的精神和范围的条件下,还可以对以上内容进行各种各样的修改。因此本发明的范围并不仅限于以上的说明,而是由权利要求书的范围来确定的。
Claims (4)
1.一种正向射线追踪等分隔射线束的生成方法,其特征在于,该方法包括:
将源置于非正多面体中心;
将所述非正多面体每个面的中心与所述每个面的顶点连接;
对所述每个面上的三角形进行几何分形,直至所述源与每个三角形顶点形成的射线束符合预定条件。
2.根据权利要求1所述的生成方法,其特征在于,所述非正多面体为C60分子结构多面体。
3.根据权利要求1所述的生成方法,其特征在于,在进行几何分形时,通过连接三角形每条边的中点得到4个分形后的三角形。
4.根据权利要求1-3任意一项所述的生成方法,其特征在于,所述预定条件为射线束角度间隔达到预定值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410234745.4A CN103986536B (zh) | 2014-05-29 | 2014-05-29 | 一种正向射线追踪等分隔射线束的生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410234745.4A CN103986536B (zh) | 2014-05-29 | 2014-05-29 | 一种正向射线追踪等分隔射线束的生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103986536A CN103986536A (zh) | 2014-08-13 |
CN103986536B true CN103986536B (zh) | 2016-04-13 |
Family
ID=51278367
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410234745.4A Expired - Fee Related CN103986536B (zh) | 2014-05-29 | 2014-05-29 | 一种正向射线追踪等分隔射线束的生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103986536B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107300645B (zh) * | 2017-06-06 | 2019-07-23 | 华中科技大学 | 一种快速的射线追踪方法及系统 |
CN107333278A (zh) * | 2017-06-30 | 2017-11-07 | 北京世纪德辰通信技术有限公司 | 一种公众移动通信信号覆盖可视化实现方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102036712A (zh) * | 2008-05-08 | 2011-04-27 | 约翰霍普金斯大学 | 卷积/迭加剂量计算方法的、使用图形处理单元加速的用于放射治疗的实时剂量计算 |
CN103067931A (zh) * | 2011-10-20 | 2013-04-24 | 中国移动通信集团广东有限公司 | 一种定位无线网络中干扰源的方法及装置 |
-
2014
- 2014-05-29 CN CN201410234745.4A patent/CN103986536B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102036712A (zh) * | 2008-05-08 | 2011-04-27 | 约翰霍普金斯大学 | 卷积/迭加剂量计算方法的、使用图形处理单元加速的用于放射治疗的实时剂量计算 |
CN103067931A (zh) * | 2011-10-20 | 2013-04-24 | 中国移动通信集团广东有限公司 | 一种定位无线网络中干扰源的方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103986536A (zh) | 2014-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107039885A (zh) | 应用于3d成像的激光阵列 | |
Aneesh et al. | Accordiohedra as positive geometries for generic scalar field theories | |
CN103986536B (zh) | 一种正向射线追踪等分隔射线束的生成方法 | |
US10174593B2 (en) | System and method for coarsening in reservoir simulation system | |
CN104063903A (zh) | 三维实体模型的四面体网格生成方法和装置 | |
CN104809293B (zh) | 一种声束在相贯线内传播的截面仿真方法和装置 | |
CN105223480B (zh) | 天线阵列时差法定位变电站局部放电源的定位误差仿真方法 | |
CN106251400B (zh) | 一种基于多四边形网格构建地形图的方法及装置 | |
CN102629390A (zh) | 海量机载LiDAR点云Delaunay三角网并行构建方法和装置 | |
CN101976468A (zh) | 一种多分辨率动态地形可视化方法及系统 | |
CN1885349A (zh) | 三维扫描的点云孔洞填补方法 | |
TW201827860A (zh) | 自由曲面照明系統的設計方法 | |
CN103425816B (zh) | 快速获取金属旋转对称体电磁散射特性的矩阵抽取方法 | |
TWI463435B (zh) | 點雲三角網格面構建系統及方法 | |
Feng et al. | An improved ray-tracing algorithm for SBR-based EM scattering computation of electrically large targets | |
CN111199109B (zh) | 基于计盒维数与聚类分析用于划分岩体均质区的耦合方法 | |
CN106023317B (zh) | 一种用于大数据测试的加权Voronoi图生成方法 | |
Meng et al. | An accelerated SBR method for RCS prediction of electrically large target | |
CN107506437A (zh) | 一种基于crushmap结构的OSD选取方法及装置 | |
CN104392027B (zh) | 一种涡轮叶片扰流柱的参数化造型方法 | |
CN107545596B (zh) | 一种点云模型最优切割平面的提取方法 | |
CN106651979A (zh) | 基于投影的空间曲面四边形网格剖分方法 | |
CN105427330A (zh) | 一种矢量多边形的任意多比例分割方法 | |
TW201827863A (zh) | 自由曲面照明系統 | |
CN106557602B (zh) | 一种基于gis的空间距离校核三维仿真方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160413 Termination date: 20170529 |