将excel数据转为地震数据

拿地震数据做网络训练时,需要生成相关的label.本文将excel表格数据转化为sac格式label。

import os
import pandas as pd
import numpy as np
from obspy import Trace
from tqdm import tqdm	
import matplotlib.pyplot as plt
import obspy
if __name__ == '__main__':
for i in tqdm(range(10000),desc = 'processing'):
    file = pd.read_excel('/Users/xichengpeng/work/data/xlsx/Z1.xlsx')
    win = np.zeros([len(file),2])
    win[:,0] = np.array(file[0])
    win[:,1] = np.array(file[1])
    pre = np.zeros([4992,520])
         
    for j in range(0,8):
        # i = 519
        pre[int(win[i*8+j,0]):int(win[i*8+j,1]),1] = 1
        time = np.arange(4992) * 0.0006
        sacfile = Trace()
        sacfile.data = pre[:,1]
        sacfile.write('/Users/xichengpeng/work/data/label/Z1_label/' + str(i) + '.sac')
        # plt.plot(time,pre,'r')
  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值