11,846
社区成员




我用segyio读取地震数据的道头字,inline、xline、offset正常,azimuth为零0。用其他商业软件查看,azimuth正常。请问如何解决?部分代码如下:
……
with segyio.open(segy_path, "r", strict=False) as segy:
trace_count = segy.tracecount
sample_count = segy.samples.size # 总样本数据长度
#trace_length = 240 + sample_count*4 # 道头240B + 样本数据
sample_interval = segy.samples[1] - segy.samples[0]
for i in range(trace_count):
header = segy.header[i]
# 读取动态道头字段
inline = header[line_field] #正常
xline = header[cdp_field] #正常
offset = header[offset_field] #正常
azimuth = header[azimuth_field] #为0,异常
……