python 读取地震道头数据_python地震数据可视化详解

心已赠人 2023-01-08 08:06 143阅读 0赞

本文实例为大家分享了python地震数据可视化的具体代码，供大家参考，具体内容如下

参考源码：seisplot

准备工作：

在windows10下安装python3.7，下载参考源码到本地。

1. demo绘图测试

demo绘图指令

cmd> python seisplot.py --demo

问题1)缺少依赖包

File "D:/Desktop/python可视化/seisplot/seisplot.py", line 16, in

import yaml

ModuleNotFoundError: No module named 'yaml'

>pip3 install yaml

Could not find a version that satisfies the requirement yaml (from versions: )

No matching distribution found for yaml

查看environment.yml:

channels:

\- defaults

\- conda-forge

dependencies:

\- matplotlib

\- numpy

\- obspy

\- pillow

\- pyyaml

安装相关依赖包，如 obspy pillow pyyaml

cmd> pip3 install obspy

出现proxy error， 添加选项 --proxy server:port

安装ok

问题2)配置文件的编码异常

D:\\Desktop\\python可视化\\seisplot>python seisplot.py --demo

\[91m

Welcome to

┌─┐┌─┐┬┌─┐┌─┐┬ ┌─┐┌┬┐

└─┐├┤ │└─┐├─┘│ │ │ │

└─┘└─┘┴└─┘┴ ┴─┘└─┘ ┴

Good luck\[0m

Traceback (most recent call last):

File "seisplot.py", line 421, in

cfg = yaml.load(f)

File "D:\\ProgramFiles\\Python36\\lib\\site-packages\\yaml\\\_\_init\_\_.py", line 70, in load

loader = Loader(stream)

File "D:\\ProgramFiles\\Python36\\lib\\site-packages\\yaml\\loader.py", line 34, in \_\_init\_\_

Reader.\_\_init\_\_(self, stream)

File "D:\\ProgramFiles\\Python36\\lib\\site-packages\\yaml\\reader.py", line 85, in \_\_init\_\_

self.determine\_encoding()

File "D:\\ProgramFiles\\Python36\\lib\\site-packages\\yaml\\reader.py", line 124, in determine\_encoding

self.update\_raw()

File "D:\\ProgramFiles\\Python36\\lib\\site-packages\\yaml\\reader.py", line 178, in update\_raw

data = self.stream.read(size)

UnicodeDecodeError: 'gbk' codec can't decode byte 0x99 in position 2003: illegal multibyte sequence

分析源码中参数传递情况：

def load(stream, Loader=Loader):

"""

Parse the first YAML document in a stream

and produce the corresponding Python object.

"""

print(stream) \#\# debug by huanying03

loader = Loader(stream)

try:

return loader.get\_single\_data()

finally:

loader.dispose()

打印结果如下：

编码不对？

65001  ：UTF-8代码页

936 ：默认的GBK

437 ：是美国英语

仔细检查发现config.yml中有异常字符，更正config.yml中错误字符：ok

问题3)python调试打印语句

print ("value=%d"%value)

demo数据绘图结果

如下：

![Image 1][]

2. 本地数据绘图测试

应用指令

cmd > python seisplot.py data/my.sgy

问题1)数据道的采样值异常

filename data2/t10.sgy

seismic.py:104:dt=4000

seismic.py:104:ns=512

seismic.py:144:header=b'C 1 CLIENT'

n\_traces 500

n\_samples 512

dt 0.004

t\_start 0

t\_end 2.044

max\_val inf

min\_val -inf

clip\_val 180462501764003194804057887685476352

Read data in 0.4 s

\++++++++++++++++++++++++++++++++++++++++++++++++++++++++

PLOTTING

plot width 17.00 in

plot height 11.00 in

Traceback (most recent call last):

File "seisplot.py", line 442, in

main(t, cfg)

File "seisplot.py", line 214, in main

cfg)

File "D:\\Desktop\\python可视化\\seisplot\\plotter.py", line 123, in plot\_histogram

y, x, \_ = ax.hist(np.ravel(data), bins=int(100.0 / (clip\_val / largest)),

OverflowError: cannot convert float infinity to integer

解决办法：将数据做归一化处理

问题2)数据文件的卷头异常

filename data2/marmousi\_vel.segy.hdrs.sgy.cdp.sgy

seismic.py:104:dt=10000

seismic.py:104:ns=350

seismic.py:144:header=b'\\xc3@\\xf1@\\xc3\\xd3\\xc9\\xc5\\xd5\\xe3'

Traceback (most recent call last):

File "seisplot.py", line 442, in

main(t, cfg)

File "seisplot.py", line 40, in main

s = Seismic.from\_segy(target, params=\{'ndim': cfg\['ndim'\]\})

File "D:\\Desktop\\python可视化\\seisplot\\seismic.py", line 164, in from\_segy

return cls.from\_obspy(stream, params=params)

File "D:\\Desktop\\python可视化\\seisplot\\seismic.py", line 147, in from\_obspy

x =np.array(list(stream.textual\_file\_header.decode()))

UnicodeDecodeError: 'utf-8' codec can't decode byte 0xc3 in position 0: invalid continuation byte

解决办法：重写3200字节的segy卷头

marmousi相关数据绘图

速度模型绘图结果如下：

![Image 1][]

marmousi模型正演的炮集数据绘图结果如下：

![Image 1][]

测试备忘

cygwin执行指令，可以看到反馈信息为彩色，分析信息很方便。

相比之下，cmd终端执行指令后的反馈信息可读性较差。

本地数据出现异常后，通过cwp中的segyread, sushw, sugain, segywrite等程序进行修改，满足绘图程序要求。

1)3200字节的卷头必须规范

2)道头字中的cdpt取值必须规范(xline)

3)道采样值最好归一化处理。

以上就是本文的全部内容，希望对大家的学习有所帮助，也希望大家多多支持我们。

[Image 1]:

python 读取地震道头数据_python地震数据可视化详解

发表评论取消回复

还没有评论，来说两句吧...

相关阅读

相关分析获取地震台网地震数据

相关 Python的数据可视化

相关数据之美-Python数据可视化

相关【Python】Python库之数据可视化

相关 Python数据可视化图实现过程详解

相关 python 读取地震道头数据_python地震数据可视化详解

相关 Python数据分析可视化--Titanic

相关 python数据可视化

相关【python可视化系列】python数据可视化利器--pyecharts

相关 python实现数据可视化