python地震數(shù)據(jù)可視化詳解

更新時(shí)間：2019年06月18日 14:12:58 作者：幻影03

這篇文章主要為大家詳細(xì)介紹了python地震數(shù)據(jù)可視化的相關(guān)資料，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下

本文實(shí)例為大家分享了python地震數(shù)據(jù)可視化的具體代碼，供大家參考，具體內(nèi)容如下

參考源碼：seisplot

準(zhǔn)備工作：

在windows10下安裝python3.7，下載參考源碼到本地。

1. demo繪圖測試

demo繪圖指令

cmd> python seisplot.py --demo

問題1）缺少依賴包

 File "D:/Desktop/python可視化/seisplot/seisplot.py", line 16, in <module>
 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

安裝相關(guān)依賴包，如 obspy pillow pyyaml

cmd> pip3 install obspy

出現(xiàn)proxy error，添加選項(xiàng) --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 <module>
 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

分析源碼中參數(shù)傳遞情況：

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()

打印結(jié)果如下：

<_io.TextIOWrapper name='config.yml' mode='r' encoding='cp936'>

編碼不對？

65001 ：UTF-8代碼頁

936 ：默認(rèn)的GBK

437 ：是美國英語

仔細(xì)檢查發(fā)現(xiàn)config.yml中有異常字符，更正config.yml中錯誤字符：ok

問題3）python調(diào)試打印語句

print ("value=%d"%value)

demo數(shù)據(jù)繪圖結(jié)果

如下：

2. 本地?cái)?shù)據(jù)繪圖測試

應(yīng)用指令

cmd > python seisplot.py data/my.sgy

問題1）數(shù)據(jù)道的采樣值異常

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 <module>
 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

解決辦法：將數(shù)據(jù)做歸一化處理

問題2）數(shù)據(jù)文件的卷頭異常

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 <module>
 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字節(jié)的segy卷頭

marmousi相關(guān)數(shù)據(jù)繪圖

速度模型繪圖結(jié)果如下：