1

我正在使用 python lasio.py 库，我必须在文档中缺少某些内容。

是否有一种本地方法可以在输出文件的 ~ASCII (~A) 部分中添加标题列。

我需要输出类似

~A  DEPTH     DT    RHOB        NPHI   SFLU    SFLA      ILM      ILD
1670.000   123.450 2550.000    0.450  123.450  123.450  110.200  105.600
1669.875   123.450 2550.000    0.450  123.450  123.450  110.200  105.600
1669.750   123.450 2550.000    0.450  123.450  123.450  110.200  105.600

我实际上得到的地方

~ASCII -----------------------------------------------------
   1670     123.45       2550       0.45     123.45     123.45      110.2      105.6
 1669.9     123.45       2550       0.45     123.45     123.45      110.2      105.6
 1669.8     123.45       2550       0.45     123.45     123.45      110.2      105.6

写入我的输出文件。

我觉得我错过了一些非常明显的东西。

编辑：文件的输出非常简单：

las.write(outfilename + '.las', fmt='%.3f', version=2)

干杯

4

1 回答 1

1

对于遇到此问题的任何人，这里有一个解决方法。需要修改 lasio 包中的 writer.py 文件。使用曲线助记符创建一个列表，将其转换为 np 数组，使用 np.vstack()将所述数组堆叠在data_arr之上

data_arr = las.data    # original code

# adding curve names as headers - insert these 5 lines in writer.py
hdr_list = []
for hdr_item in las.curves:
    hdr_list.append(hdr_item.mnemonic)
hdrs = np.array(hdr_list).reshape(1, -1)
data_arr = np.vstack((hdrs, data_arr))

nrows, ncols = data_arr.shape   # original code
于 2021-04-04T08:44:54.753 回答