CN103544675B - 地震波传播展示方法 - Google Patents
地震波传播展示方法 Download PDFInfo
- Publication number
- CN103544675B CN103544675B CN201310501172.2A CN201310501172A CN103544675B CN 103544675 B CN103544675 B CN 103544675B CN 201310501172 A CN201310501172 A CN 201310501172A CN 103544675 B CN103544675 B CN 103544675B
- Authority
- CN
- China
- Prior art keywords
- data
- gis
- earthquake intensity
- seismic
- earth
- 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.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及计算机信息技术领域，特别涉及地震波传播展示方法。地震波传播展示方法，包括：预先获取地震波数据，根据地震波数据获取地震波位移信息及烈度信息；在地理信息系统中，根据地震波位移信息对应修改地理信息系统的地表地形网格的顶点坐标，且根据烈度信息对应在地理信息系统的地表影像中渲染烈度标识，以在地理信息系统进行地震波传播展示。本发明提供的地震波传播展示方法能够直观可视的展示地震波在地表的传播状况。
Description
技术领域
本发明涉及计算机信息技术领域，具体而言，涉及地震波传播展示方法。
背景技术
相关技术中，为展示地震波在地表的传播，所采用的方法主要为，在三维空间数据场中利用地震波数据构造中间几何图元(如曲面、平面、切片等)，然后根据计算机图形学绘制展示面，从绘制的展示面中抽取地震波传播等值面，利用抽取的地震波传播等值面进行地震波在地表的传播展示。
但，相关技术中的地震波传播展示方法，主要以构建等值面的波形为主，并未考虑地形地貌等地理因素对地震波传播的影响，导致无法直观可视地展示地震波在地表的传播状况。
发明内容
本发明的目的在于提供地震波传播展示方法，以解决上述的问题。
在本发明的实施例中提供了地震波传播展示方法，包括：
在地震场模拟计算模型中，预先获取地震波数据，根据所述地震波数据获取地震波位移信息及烈度信息；
在地理信息系统中，根据所述地震波位移信息对应修改所述地理信息系统的地表地形网格的顶点坐标，且根据所述烈度信息对应在所述地理信息系统的地表影像中渲染烈度标识，以在所述地理信息系统进行地震波传播展示；
其中，所述根据所述地震波数据获取地震波位移信息及烈度信息，包括：
所述地震波数据为地震波速度三分量数据；
将所述地震波速度三分量数据与该数据采集的时间值进行运算，获取地震波位移信息；
根据所述地震波速度三分量数据获取对应的地震波峰值速率；
根据所述地震波峰值速率获取对应的烈度信息。
本发明上述实施例的地震波传播展示方法中，根据地震波数据获取地震波位移信息及烈度信息，并且在地理信息系统中，通过地震波位移信息对地表地形网格的顶点坐标进行修改，地理信息系统中的地表地形层模拟实际的地理地形地貌，通过对地表地形网格顶点坐标的修改能够展示地震波的传播状态，以及地震时地表的起伏状态，通过烈度信息在地理信息系统的地表影像中渲染烈度标识，能够直观展示地震时地表烈度，因此本发明实施例的地震波传播展示方法能够直观可视的展示地震波在地表的传播状况。
附图说明
图1示出了本发明实施例的地震波传播展示方法的流程图。
具体实施方式
下面通过具体的实施例子并结合附图对本发明做进一步的详细描述。
本发明实施例提供的地震波传播展示方法，如图1所示，主要处理过程包括：
步骤S11：预先获取地震波数据，根据地震波数据获取地震波位移信息及烈度信息；
步骤S12：在地理信息系统中，根据地震波位移信息对应修改地理信息系统的地表地形网格的顶点坐标，且根据烈度信息对应在地理信息系统的地表影像中渲染烈度标识，以在地理信息系统进行地震波传播展示。
获取地震波数据，包括：在地震场模拟计算模型中，设置地表等间隔采样点，并按预设时间间隔获取每个地表等间隔采样点的地震波数据。
优选地，在地震场模拟计算模型中，模拟产生地震发生到地震结束过程中所产生的、以震源为中心的地震波数据。
具体地，以预设的震源为中心在地震场模拟计算模型中设置地表等间隔采样点，并按一定时间间隔获取每个地表等间隔采样点的地震波数据。
优选地，时间间隔之间的、每个地表等间隔采样点的地震波数据可以利用线性插值，根据渲染时间点来获取。
本发明实施例中，根据地震波数据获取地震波位移信息及烈度信息，包括：地震波数据为地震波速度三分量数据；将地震波速度三分量数据与该数据采集的时间值进行运算，获取地震波位移信息；根据地震波速度三分量数据获取对应的地震波峰值速率；根据地震波峰值速率获取对应的烈度信息。
预先获取地震波数据，根据地震波数据获取地震波位移信息及烈度信息之后，在地理信息系统中，根据地震波位移信息对应修改地理信息系统的地表地形网格的顶点坐标之前，还包括：将地震波位移信息及烈度信息均转换为对应的纹理数据。
在地理信息系统中，根据地震波位移信息对应修改地理信息系统的地表地形网格的顶点坐标，且根据烈度信息对应在地理信息系统的地表影像中渲染烈度标识，包括：纹理数据为离散采样数据；通过对离散采样数据进行插值运算，得到任意时间值的纹理数据；通过多线程调度任意时间值的纹理数据，在地理信息系统中，实现随设定时间值的不同，对应修改地理信息系统的地表地形网格的顶点坐标，且对应在地理信息系统的地表影像中渲染烈度标识。
进一步地，将地震波位移信息及烈度信息均转换为对应的纹理数据，包括：地震波位移信息包括三个位移子分量；将三个位移子分量及烈度信息分别对应存入RGBA的四个色彩通道中，并对应输出地震波纹理位移数据及烈度纹理坐标。纹理数据包含四个通道RGBA，其中RGB存储位移信息，A存储烈度信息。
本发明实施例中，烈度通过颜色表现，烈度纹理数据保存了不同烈度所对应的颜色数据值，通过烈度纹理坐标即可计算出不同烈度所对应的颜色数据值。因此烈度信息分为两部分保存，分别为烈度的颜色纹理图片和采样点的烈度纹理坐标。
本发明实施例中，将获取的地震波数据转换为纹理数据，达到简化地震波数据的目的，转化后的数据量为原始数据量的三分之一。
进一步，利用显卡里面的图形处理器(GPU)接收产生的地震波纹理位移数据及烈度纹理坐标，并分别对应提取地震波纹理位移数据的三个位移分量及烈度纹理坐标对应的烈度颜色信息，从而实现将原来的一个地震波数据场优化为可满足实时计算的一帧一帧的数据。
根据地震波位移信息对应修改地理信息系统的地表地形网格的顶点坐标，包括：获取地理信息系统的地表地形网格的顶点经纬度坐标；将顶点经纬度坐标转换为顶点纹理坐标；地震波纹理位移数据中包括纹理位移数据，提取地震波纹理位移数据中的纹理位移数据；根据纹理位移数据，对应计算顶点纹理坐标的变化差值；根据变化差值对应修改地理信息系统的地表地形网格的顶点坐标。
优选地，利用显卡的GPU，将地理信息系统中的地表地形网格上顶点经纬度坐标从内存传入到GPU显存，通过显卡双线性插值，滤波操作等，将地形网格上顶点的经纬度坐标转换为顶点纹理坐标。利用GPU的顶点着色程序(VertexShader)读取地震波纹理中保存的纹理位移数据，根据纹理位移数据，对应计算顶点纹理维度坐标值的变化差值；根据变化差值对应修改地理信息系统的地表地形网格的顶点坐标。从而达到根据地震波的位移信息随时间变化修改地形网格上的顶点坐标，通过对离散的纹理数据做插值实现地理信息系统的地表在空间上和时间上的随地震波得动态的连续起伏变化。
根据烈度信息对应在地理信息系统的地表影像中渲染烈度标识，包括：通过图像硬件光栅化插值，获取地理信息系统的地表地形网格的顶点的像素所对应的片段的经纬度；将片段的经纬度转化为片段纹理坐标；根据片段纹理坐标获取烈度颜色数据；获取片段在地理信息系统的地表影像颜色数据；按预设规则对相对应的烈度颜色数据及地表影像颜色数据进行加权求和运算，得到烈度标识数据；根据烈度标识数据对应在地理信息系统的地表影像中渲染烈度标识。
由此可以实现，烈度的颜色和影像的颜色的融合，从而实现地理信息系统中地表影像与烈度信息动态融合。
本发明实施例中，在地理信息系统中，根据地震波位移信息对应修改地理信息系统的地表地形网格的顶点坐标，且根据烈度信息对应在地理信息系统的地表影像中渲染烈度标识，包括：纹理数据为离散采样数据；通过对离散采样数据进行插值运算，得到任意时间值的纹理数据；通过多线程调度任意时间值的纹理数据，在地理信息系统中，实现随设定时间值的不同，对应修改地理信息系统的地表地形网格的顶点坐标，且对应在地理信息系统的地表影像中渲染烈度标识。
将离散的采样时间的位移和烈度的信息表达的纹理数据，通过插值来获取任意时间的位移和烈度。采用多线程调度纹理数据，载入内存的纹理数据保存在纹理队列中，从而保证纹理数据的渲染流畅。为了防止多线程访问冲突，在访问线程间的共享数据时加锁。从而实现基于地理信息系统的地震波纹理位移数据调度。
从地震发生到地震消失这个时间段内，地震波数据量大，本发明实施例中将地震波数据转换为纹理数据，解决了地震波数据的优化组织问题，缓解内存、外存调度实现的困难。
本发明实施例的数据波传播展示方法，采用图形处理器(GPU)，加快了地震波数据的处理，实现了纹理数据调度可视化，提高了地震波模拟的效率。
本发明实施例中，将地震波与地理信息系统中的地形、地貌、地理要素等信息相融合进行展示，增强地震波传播展示效果。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (7)
1.地震波传播展示方法，其特征在于，包括：
在地震场模拟计算模型中，预先获取地震波数据，根据所述地震波数据获取地震波位移信息及烈度信息；
在地理信息系统中，根据所述地震波位移信息对应修改所述地理信息系统的地表地形网格的顶点坐标，且根据所述烈度信息对应在所述地理信息系统的地表影像中渲染烈度标识，以在所述地理信息系统进行地震波传播展示；
其中，所述根据所述地震波数据获取地震波位移信息及烈度信息，包括：
所述地震波数据为地震波速度三分量数据；
将所述地震波速度三分量数据与该数据采集的时间值进行运算，获取地震波位移信息；
根据所述地震波速度三分量数据获取对应的地震波峰值速率；
根据所述地震波峰值速率获取对应的烈度信息。
2.根据权利要求1所述的方法，其特征在于，所述获取地震波数据，包括：
在地震场模拟计算模型中，设置地表等间隔采样点，并按预设时间间隔获取每个所述地表等间隔采样点的地震波数据。
3.根据权利要求2所述的方法，其特征在于，所述预先获取地震波数据，根据所述地震波数据获取地震波位移信息及烈度信息之后，在地理信息系统中，根据所述地震波位移信息对应修改所述地理信息系统的地表地形网格的顶点坐标之前，还包括：
将所述地震波位移信息及所述烈度信息均转换为对应的纹理数据。
4.根据权利要求3所述的方法，其特征在于，所述在地理信息系统中，根据所述地震波位移信息对应修改所述地理信息系统的地表地形网格的顶点坐标，且根据所述烈度信息对应在所述地理信息系统的地表影像中渲染烈度标识，包括：
所述纹理数据为离散采样数据；
通过对所述离散采样数据进行插值运算，得到任意时间值的纹理数据；
通过多线程调度任意时间值的纹理数据，在地理信息系统中，实现随设定时间值的不同，对应修改所述地理信息系统的地表地形网格的顶点坐标，且对应在所述地理信息系统的地表影像中渲染烈度标识。
5.根据权利要求3所述的方法，其特征在于，所述将所述地震波位移信息及所述烈度信息均转换为对应的纹理数据，包括：
所述地震波位移信息包括三个位移子分量；
将三个所述位移子分量及所述烈度信息分别对应存入RGBA的四个色彩通道中，并对应输出地震波纹理位移数据及烈度纹理坐标。
6.根据权利要求5所述的方法，其特征在于，所述根据所述地震波位移信息对应修改所述地理信息系统的地表地形网格的顶点坐标，包括：
获取所述地理信息系统的地表地形网格的顶点经纬度坐标；
将所述顶点经纬度坐标转换为顶点纹理坐标；
所述地震波纹理位移数据中包括纹理位移数据，提取所述地震波纹理位移数据中的纹理位移数据；
根据所述纹理位移数据，对应计算所述顶点纹理坐标的变化差值；
根据所述变化差值对应修改所述地理信息系统的地表地形网格的顶点坐标。
7.根据权利要求4所述的方法，其特征在于，所述根据所述烈度信息对应在所述地理信息系统的地表影像中渲染烈度标识，包括：
通过图像硬件光栅化插值，获取所述地理信息系统的地表地形网格的顶点的像素所对应的片段的经纬度；
将所述片段的经纬度转化为片段纹理坐标；
根据所述片段纹理坐标获取烈度颜色数据；
获取所述片段在所述地理信息系统的地表影像颜色数据；
按预设规则对相对应的所述烈度颜色数据及所述地表影像颜色数据进行加权求和运算，得到烈度标识数据；
根据所述烈度标识数据对应在所述地理信息系统的地表影像中渲染烈度标识。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310501172.2A CN103544675B (zh) | 2013-10-22 | 2013-10-22 | 地震波传播展示方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310501172.2A CN103544675B (zh) | 2013-10-22 | 2013-10-22 | 地震波传播展示方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103544675A CN103544675A (zh) | 2014-01-29 |
| CN103544675B true CN103544675B (zh) | 2016-06-15 |
Family
ID=49968100
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310501172.2A Active CN103544675B (zh) | 2013-10-22 | 2013-10-22 | 地震波传播展示方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103544675B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105096385B (zh) * | 2015-07-20 | 2018-06-12 | 北京工业大学 | 一种二维地震剖面三维显示方法 |
| CN107527383B (zh) * | 2017-08-30 | 2020-12-25 | 北京市地震局 | 一种地震影响场三维漫延可视化方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101545982A (zh) * | 2009-05-05 | 2009-09-30 | 江西省电力信息通讯有限公司 | 地震计算模型的动态拟合方法 |
| CN101980226A (zh) * | 2010-11-19 | 2011-02-23 | 中国地震灾害防御中心 | 城市震害预测虚拟仿真系统 |
| CN201780711U (zh) * | 2009-10-10 | 2011-03-30 | 嘉兴市科技馆 | 地震震级与烈度演示展板 |
| CN102183790A (zh) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | 基于时空双变网格的弹性波正演模拟技术 |
-
2013
- 2013-10-22 CN CN201310501172.2A patent/CN103544675B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101545982A (zh) * | 2009-05-05 | 2009-09-30 | 江西省电力信息通讯有限公司 | 地震计算模型的动态拟合方法 |
| CN201780711U (zh) * | 2009-10-10 | 2011-03-30 | 嘉兴市科技馆 | 地震震级与烈度演示展板 |
| CN101980226A (zh) * | 2010-11-19 | 2011-02-23 | 中国地震灾害防御中心 | 城市震害预测虚拟仿真系统 |
| CN102183790A (zh) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | 基于时空双变网格的弹性波正演模拟技术 |
Non-Patent Citations (2)
| Title |
|---|
| 肖金城.大规模地形场景三维实时漫游显示技术研究.《遥感信息》.2002, * |
| 褚春雷等.非规则三角网格有限差分法地震正演模拟.《中国海洋大学学报》.2005,第35卷(第1期), * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103544675A (zh) | 2014-01-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104200506B (zh) | 三维gis海量矢量数据渲染方法及装置 | |
| CN108470369B (zh) | 一种水面渲染方法及装置 | |
| CN109102560B (zh) | 三维模型渲染方法及装置 | |
| Marston et al. | Improving the representation of major landforms in analytical relief shading | |
| CN103606192A (zh) | 一种基于三维虚拟地球的风场可视化展示方法 | |
| CN113034656B (zh) | 游戏场景中光照信息的渲染方法、装置及设备 | |
| CN110706326B (zh) | 数据展示方法及装置 | |
| CN110796732A (zh) | 基于WebGL与C#的矿山三维展示方法 | |
| CN103530901A (zh) | 一种匹配真实三维环境的动态水域可视化方法及系统 | |
| CN112328929A (zh) | 基于矢量瓦片地图发布技术的人流密度动态可视化方法 | |
| Du et al. | Spatio-temporal visualization of air–sea CO2 flux and carbon budget using volume rendering | |
| CN114820990A (zh) | 一种基于数字孪生的流域防洪可视化方法及系统 | |
| CN104463937A (zh) | 基于不规则物体的动画生成方法及装置 | |
| CN103544675B (zh) | 地震波传播展示方法 | |
| CN105427378B (zh) | 基于数字多媒体技术的文化古城展示方法及系统 | |
| Semmo et al. | Real-time rendering of water surfaces with cartography-oriented design | |
| CN109064534A (zh) | 一种基于网页的全球海浪真实感绘制方法 | |
| CN113144613A (zh) | 基于模型的体积云生成的方法 | |
| CN106873931A (zh) | 栅格数据可视化装置、方法及计算机系统 | |
| Whelan et al. | Formulated silhouettes for sketching terrain | |
| Döllner | Geovisualization and real-time 3D computer graphics | |
| CN103413356A (zh) | 一种地形地貌三维动态模拟显示方法和系统 | |
| CN109658495A (zh) | 环境光遮蔽效果的渲染方法、装置及电子设备 | |
| CN111243099B (zh) | 一种处理图像的方法和装置以及在ar设备中显示图像的方法和装置 | |
| She et al. | Rendering 2d lines on 3d terrain model with optimization in visual quality and running performance |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |