segy地震数据的读取python_Python-segy格式地震资料segyio读写包说明（二）,pythonsegy,数据,学习,笔记...

python segy格式地震数据读写包segyio学习笔记(二)

最近大致搞明白了segyio读取叠后和叠前segy数据的方法，以及内部存储结构，以两段代码为例：

叠后数据读取。这是一个从给定时窗内提取每道的均方根振幅，然后进行拟合求取tar因子的简单代码：

import segyio

import numpy as np

import matplotlib.pyplot as plt

file = ‘/home/geou6/0922_j301_down_test.segy’

data =np.loadtxt(’/home/geou6/0922_j301_op_ampana.txt’

r = np.array([])

with segyio.open(file, iline=189, xline=21) as f:

for l in f.iline:

for t in l:

amp = np.sqrt(np.mean(t ** 2))

r = np.append(r, amp)

lnt = np.log(data[:, 2])

lnr = np.log®

f_1 = np.polyfit(lnt, lnr, 1)

p_1 = np.poly1d(f_1)

print(f_1)

print(p_1)

plt.plot(data[:, 2], r)

plt.show()

plt.plot(lnt, lnr)

plt.show()

读取叠后数据时，第一个循环l in f.iline，这时l是一条illine线，第二重循环t in l中，t才是一个地震道，是一个一维的numpy数组。注意f.iline和并不是ndarray。他们的存储结构我还没有搞太明白，感觉是个类似生成器或者类的东西？

读取叠前数据时：

file = ‘/home/geou6/0922_j301_down_gather_test1.segy’

with segyio.open(file, iline=189, xline=21) as f:

print(f.gather[1 ,1][0:2]) #这时打印的1线1点的第一道

for m in f.gather[:, :]: #m是个道集和非一条illine线

for n in m: #n是地震道

print(n)

for k in f.gather[1, 1]:

for j in k: # k 是地震道, j 就是样点了

print(j)

注意f.gather[1, 1][0 :2]不能写成f.gather[1,1,1]或f.gather[1,1][0],f.gather[:, :]也不能写成f.gather，这个我还要继续去搞明白

segy地震数据的读取python_Python-segy格式地震资料segyio读写包说明（二）,pythonsegy,数据,学习,笔记...相关推荐

matlab的seg是什么,地震数据文件读取——以SEG-Y格式为例
1.引言地震勘探是利用地下介质弹性和密度的差异,通过观测和分析大地对人工激发地震波的响应,推断地下岩层的性质和形态的地球物理勘探方法,是钻探前勘测石油与天然气资源的重要手段,在煤田和工程地质勘查. ...
PE文件数字签名信息读取存储及格式具体解释图之上(历史代码，贴出学习)
// 注意下图PE文件格式具体解释图中的 // IMAGE_NT_HEADERS------->OptionalHeader------>DataDirectory[IMAGE_DIREC ...
02_clickhouse安装，卸载，启动/关闭服务，交互式命令（数据库创建，数据导入，查询），批模式数据导入，MySQL接口操作ClickHouse，DBeaver可视化工具安装与使用（学习笔记）
1 ClickHouse安装安装文件清单 clickhouse-client-${version}.noarch.rpm clickhouse-common-static-dbg-${version ...
python怎么筛选excel数据_懂Excel也能轻松入门Python数据分析包pandas(二)：高级筛选(上)-excel筛选...
转发本文并私信我"python",即可获得Python资料以及各种心得(持续更新的) 系列文章: 前言经常听别人说 Python 在数据领域有多厉害,结果学了很长时间,连数据处理 ...
音视频封装格式转换器(支持avi格式转换)，基于FFmpeg4.1实现（音视频学习笔记二）
之前参照雷霄骅博士的最简单的基于FFMPEG的封装格式转换器(无编解码)的博客和FFmpeg官网的example,实现一个简单的封装格式转换器.但是后来我发现我想从mp4格式转换成avi格式的时候会报 ...
3-7 pandas数据的读取与存储
数据分析工具pandas 7. 数据的读取与存储 7.1 读操作 7.2 写操作 7.3 JSON格式 7.4 分块读取大文件 Pandas是一个强大的分析结构化数据的工具集,基于NumPy构建,提供 ...
matlab读取grib2数据,matlab读取grib2数据
2018 基于通用模型的 GRIB 格式数据读取技术王兵,李杰 (1. 南京航空航天大学民航学院,江苏南京 211106;2. 国家空管飞行流量技术重点实验室,江苏南京 211106) 摘要: ...
Python文件和数据格式化学习笔记
一.文件的使用 1. 文件的定义: 文件是数据的抽象和集合:文件是存储在辅助存储器上的数据序列:文件是数据存储的一种形式: 2. 文件的展现方式:文本文件和二进制文件(两种) 文本文件和二进制文件只是 ...
数据访问模式二：数据集和数据适配器（传统的数据访问模式）
上一篇文章介绍了使用DataSource控件访问数据库的过程,本节介绍利用数据适配集/数据适配器的访问数据库.这两种设计模式的差别,使得GridView的设计即要支持DataSource控件的数据绑定 ...

segy地震数据的读取python_Python-segy格式地震资料segyio读写包说明（二）,pythonsegy,数据,学习,笔记...

segy地震数据的读取python_Python-segy格式地震资料segyio读写包说明（二）,pythonsegy,数据,学习,笔记...相关推荐

最新文章

热门文章