equinor / segyio

Fast Python library for SEGY files.

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

I want to know which data represent final_clusters? Thanks in anticipation

farooqad0 opened this issue · comments

I want to know which data represent final_clusters? Thanks in anticipation

spec = segyio.spec()
file = 'segy_attribute_classifications_run2.segy'
spec.format = 1
spec.xline = 193
spec.iline = 189
spec.offsets = [0]
spec.samples = np.unique(final_values[:,2])
spec.ilines = np.unique(final_values[:,0])
spec.xlines = np.unique(final_values[:,1])
with segyio.create(file, spec) as f:
    ix = len_ilines*len_xlines
    f.trace[:ix] = final_clusters[:ix]
    hdx = 0
    for ilno in spec.ilines:
        for xlno in spec.xlines:
            f.header[hdx] = {segyio.TraceField.INLINE_3D: int(ilno),
                             segyio.TraceField.CROSSLINE_3D: int(xlno),
                             segyio.TraceField.offset:0}
            hdx = hdx+1

c = segyio.open(file, 'r+')

Is this code from nitishkannan23's issue?

Closing due to inactivity