将一炮地震数据标记的label画在一张图上

朋哥的求职之路

于 2020-03-15 20:52:29 发布

阅读量185

点赞数 1

本文链接：https://blog.csdn.net/weixin_46518295/article/details/104885747

版权

这篇博客主要介绍了如何使用Python的库，如obspy和matplotlib，来读取地震数据并将其标记的label绘制在同一张图上，便于数据理解与分析。

摘要由CSDN通过智能技术生成

import os
import pandas as pd
import numpy as np
from obspy import read
import matplotlib.pyplot as plt
from tqdm import tqdm

if name == ‘main’:

for z in tqdm(range(37960),desc='Processing'):
    datapath = '/Users/xichengpeng/work/data/train_data/syn/R/'+str(z)+'.sac'
  
    win = np.array(pd.read_excel('/Users/xichengpeng/work/data/train_data/xlsx/R.xlsx'))[:,1:3]
    selectwin = np.zeros((300,10,2))
   
    for i in range(300):
        selectwin[i] = win[i*10:(i+1)*10]
   
    del win
    os.chdir(datapath)
    fig = plt.figure()
    ax = fig.add_subplot()
   
    for i in range(300):
        filename = str(i) + ".sac"
        sac = read(filename)
        data = sac[0].data
        data = data / np.max(np.abs(data))
        data = data + i * 1
        time = np.arange(len(dat