ldgyb的博客

python segy格式地震数据读写包segyio学习笔记（二）最近大致搞明白了segyio读取叠后和叠前segy数据的方法，以及内部存储结构，以两段代码为例：叠后数据读取。这是一个从给定时窗内提取每道的均方根振幅，然后进行拟合求取tar因子的简单代码：import segyioimport numpy as npimport matplotlib.pyplot as pltfile = ‘/home/geou6/0922_j301_down_test.segy’data =np.loadt

今天安装segyio包后读取了一块叠后数据，运行程序后报错unable to find sorting，使用github提供的测试数据可以正常读取，一度怀疑是程序包本身的问题，上网查资料有网友反映也遇到过这个问题，segyio 1.60以及以前的版本可以正常运行，于是去pypi翻出了1.6版本，发现只能在python3.6上运行，于是去清华大学的镜像网站下找到历史版本，在工作站上登录一个新的用户重新安装python3.6，然后pip安装上segyio1.6，折腾一下午，再运行，还是报错！strict还是不能

python中对序列FFT后结果不是共轭对称的问题我对雷克子波进行了fft后，发现结果并不是共轭对称的。这导致我把fft后结果取前一半区间，让后共轭对称到后一半，再反变换回去的结果就有点异常了。如果直接把fft后的全序列做ifft，结果就是对的。第二张图是fft后，把前一半区间取共轭对称，再和后一半区间相减，可以看到误差还是挺大的，而且越到后面越大。我在别的论坛看到有人用matlab也遇到过这个问题。请问这是什么原因？附上代码：import numpy as npfrom scipy imp

用普通账户在linunx环境下安装anaconda3（1）从官网下载anaconda的安装文件Aanconda******.sh，放到个人账号的路径/home/w022/下（w022是一个个人账户）。执行命令： bash Aanconda******.sh然后按照提示按回车和yes就可以了。（2）安装完成后，anaconda3装在了/home/w022/下。（3）在当前目录下vi .ba...

python和matlab一样，可以创建函数的句柄。比如：window = turtle.Screen()window.bgcolor(‘blue’)window.exitonclick()window就是一个句柄。句柄就是将直接调用改为间接调用，以避免每次调用的时候在库里搜寻，提高调用速度。可以这样理解吗？...

今天发了一个有意思的事，sum（）ndarray数组求和时，默认是按axis=0轴聚合的。比如，import numpy as np 后a = np.array([1, 2, 3])sum(a)结果自然是6a = np.array([[1,2,3],[4,5,6],[7,8,9]])此时a为[[1 2 3][4 5 6][7,8,9]]sum(a)结果成了array[12,15...

做了一道九宫格的题，用了itertools中的排列方法。

今天用函数的时候遇到了一个问题，后来总结了下：函数中若对变量进行了初始化，则该变量为局部变量；若没有初始化而直接对该变量操作，则该变量为全局变量。一般这种情况出现在列表和字典的使用中。如：a = 1b = 1list = []dict = {}def ch1():b = a + 1def ch2():list.append(10)def ch3():dict[‘x’] = 1...