在第5课《怎样修改保存二维和三维地震数据?》中,我们给大家讲解了怎样读取三维SEGY地震数据并保存为二进制类型数据,比如dat或npy格式的二进制数据。那如果我们对三维数据进行处理后,还希望保存回SEGY格式地震数据,又怎样来操作呢?
我们还是以1996年新西兰塔拉纳基盆地叠前数据Kerry3D.segy为例来进行讲解。这个数据概貌是这样的:
三维SEGY数据保存的方法主要代码是这样的:
其实这个代码和二维数据保存差别不大,有几个关键点讲一下:
一是主要的要素和二维程序类似,都是要指定的,比如sorting,format,samples,tracecount等等。如果变动不大,就可以直接用原来的数据取值。
二是要指定ilines和xline,这和二维程序不同。这里不太建议使用src.ilines这种方法,因为数据结构可能不规范,无法读出数值。
三是在open的参数里面,最好用strict=false,这样保证了数据结构不太规范也能读出来数值。
保存好数据之后,我们查看一个剖面,是可以正常显示的:
以上就是今天的更新内容。如果想了解更多的数据保存方法,请查看原课程《怎样修改保存二维和三维地震数据?》学习。如果对程序感兴趣,欢迎查看原课程下载,再见。
扩展阅读:
网友评论