daloroat / first_breaks_picking Goto Github PK

I would like to ask how to use your code. I just don't kwon how to start it. Could you teach me?

Hi @DaloroAT

I was using the tool today and noticed this issue with manual time picks, which goes back to original position when amplitudes are adjusted.

Regards

It's a good tool to pick the first arrival time of seismic data. Thank you for your contribution. The input file for PyGIMLI need the first-arrival-time file, which contains three columns data: shot position, geophone position, travel-time. How can I obtain this formatted file from first_breaks_picking tools? If I can get this file, the workflow for tomography and full-waveform inversion can continue.
Yours Sincerely,
Li Jian
China University of Geosciences

I tried to run this on 2 different land datasets.

Fisrt one has pretty good signal/noise ratio and the algorithm works pretty good:

The second dataset is of poor quality: geology includes permafrost with high energy attenuation, surface noise is high:

Do you have ideas what we can do to improve the result of the second dataset?

Feature request described here

#21 (comment)

@DaloroAT when we do VSP time picking we prefer to have the traces plotted horizontally than vertically for accuracy in time picking purpose. Is it possible you can deploy it Visual Settings.

Described here

#21 (comment)

@DaloroAT , Google Colab is getting crashed when I am running the function,

export_image(task, image_filename,
             time_window=(0, 60),
             traces_window=(79.5, 90.5),
             show_processing_region=False,
             headers_total_pixels=80,
             height=500,
             width=700,
             hide_traces_axis=True
             )

Below is the Google Colab log file:

Timestamp,Level,Message
"Jun 6, 2024, 12:37:59 AM",WARNING,WARNING:root:kernel 7ef480cf-4b19-47a6-8eb2-038a24f6e581 restarted

"Jun 6, 2024, 12:37:59 AM",INFO,"KernelRestarter: restarting kernel (1/5), keep random ports"

"Jun 6, 2024, 12:37:58 AM",WARNING,"Available platform plugins are: eglfs, linuxfb, minimal, minimalegl, offscreen, vnc, wayland-egl, wayland, wayland-xcomposite-egl, wayland-xcomposite-glx, webgl, xcb."

"Jun 6, 2024, 12:37:58 AM",WARNING,This application failed to start because no Qt platform plugin could be initialized. Reinstalling the application may fix this problem.

"Jun 6, 2024, 12:37:58 AM",WARNING,"qt.qpa.plugin: Could not load the Qt platform plugin ""xcb"" in """" even though it was found."

"Jun 6, 2024, 12:36:34 AM",WARNING,"�[0;93m2024-06-05 19:06:34.678717360 [W:onnxruntime:Default, onnxruntime_pybind_state.cc:870 CreateExecutionProviderInstance] Failed to create CUDAExecutionProvider. Please reference https://onnxruntime.ai/docs/execution-providers/CUDA-ExecutionProvider.html#requirementsto ensure all dependencies are met.�[m"

"Jun 6, 2024, 12:36:34 AM",WARNING,�[m

"Jun 6, 2024, 12:36:34 AM",WARNING,"�[1;31m2024-06-05 19:06:34.678683952 [E:onnxruntime:Default, provider_bridge_ort.cc:1744 TryGetProviderInfo_CUDA] /onnxruntime_src/onnxruntime/core/session/provider_bridge_ort.cc:1426 onnxruntime::Provider& onnxruntime::ProviderLibrary::Get() [ONNXRuntimeError] : 1 : FAIL : Failed to load library libonnxruntime_providers_cuda.so with error: libcublasLt.so.11: cannot open shared object file: No such file or directory"

"Jun 6, 2024, 12:36:33 AM",WARNING,�[m

"Jun 6, 2024, 12:36:33 AM",WARNING,"�[1;31m2024-06-05 19:06:33.857753107 [E:onnxruntime:Default, provider_bridge_ort.cc:1744 TryGetProviderInfo_CUDA] /onnxruntime_src/onnxruntime/core/session/provider_bridge_ort.cc:1426 onnxruntime::Provider& onnxruntime::ProviderLibrary::Get() [ONNXRuntimeError] : 1 : FAIL : Failed to load library libonnxruntime_providers_cuda.so with error: libcublasLt.so.11: cannot open shared object file: No such file or directory"

For saving the data, I get this error.

np.save(args.data_dir / f'model_{k}.npy', model_and_picking)
  File "<__array_function__ internals>", line 200, in save
  File "env/lib/python3.9/site-packages/numpy/lib/npyio.py", line 521, in save
    arr = np.asanyarray(arr)
ValueError: setting an array element with a sequence. The requested array has an inhomogeneous shape after 1 dimensions. The detected shape was (2,) + inhomogeneous part.

To avoid this problem, the output of model_generator.get_random_model() should be converted to ndarray with dtype="object"

Hi,

Are there any reasons why it doesn't work with python > 3.10?
For example when trying to run pip install first-breaks-picking-gpu I get error:

ERROR: Ignored the following versions that require a different python version: 0.3.100 Requires-Python <=3.10,>=3.8; 0.4.0a0 Requires-Python <=3.10,>=3.8; 0.4.0a2 Requires-Python <=3.10,>=3.8; 0.5.1 Requires-Python <=3.10,>=3.8; 0.6.0 Requires-Python <=3.10,>=3.8
ERROR: Could not find a version that satisfies the requirement first-breaks-picking-gpu (from versions: none)
ERROR: No matching distribution found for first-breaks-picking-gpu

[notice] A new release of pip available: 22.2.2 -> 23.3.1
[notice] To update, run: python.exe -m pip install --upgrade pip

Hi @HHTC0.

I added the ability to refine the first breaks manually and export the pick to SGY. Could you check it?

Described here

#21 (comment)

I'm a beginner in deep learning.My data is in IBM format.I got a error when I used 'sgy reader.py' to read my data.Could you teach me how to use your program in my real data?I want to get the time of firs break.Thanks!

Feature request described here

#21 (comment)

Right now They are stored in the last 4 bytes (236 byte if counting from 0) to be decoded as an unsigned integer. Values are in microseconds.

Hello @kerim371

I added GPU support to pick first breaks as we discussed earlier. Could you check it pls? Installation guide is here.

Hi @DaloroAT,

In VSP view, first break pick looks much smoother when initially picked (red line). However, only can see it unsmooth when loaded back (blue line).

If you spare sometime to fix this over weekend, that would be awesome as I am in middle of project.

Regards

Hi @DaloroAT,

Hope that you are well. Hope you had good break from my requests :)

Is it possible to modify the spectrum functionality so that it give individual spectrum for each trace in loaded gather? Currently it gives average out spectrum for selected window.

Here is an example I could find https://www.researchgate.net/figure/Frequency-content-of-signal-Example-of-CRP-gather-of-14-position-in2-2-seismic-line-a_fig3_348750792

Regards

          @DaloroAT if you planning to develop it further. Having a BPF would be helpful for time picking.

Originally posted by @HHTC0 in #29 (comment)

Hi,

I would like to try this package but it is difficult to start.

Is there any examples how to use it? Especially I'm looking for the python commands...

Hi,

Thank you for sharing this with us.

I receive the error **** after I selected SGY file.

Also, it breaks my Anaconda installation after I install the library and I have reinstalled.

Some other error message I received when I installed the library are below. I installed these libraries afterwards and then reinstalled the main library. But still the same error message as above.

ERROR: pip's dependency resolver does not currently take into account all the packages that are installed. This behaviour is the source of the following dependency conflicts.
gensim 4.3.0 requires FuzzyTM>=0.4.0, which is not installed.
conda-repo-cli 1.0.41 requires requests_mock, which is not installed.
numba 0.56.4 requires numpy<1.24,>=1.18, but you have numpy 1.24.2 which is incompatible.
conda-repo-cli 1.0.41 requires clyent==1.2.1, but you have clyent 1.2.2 which is incompatible.
conda-repo-cli 1.0.41 requires nbformat==5.4.0, but you have nbformat 5.7.0 which is incompatible.
conda-repo-cli 1.0.41 requires requests==2.28.1, but you have requests 2.28.2 which is incompatible.

ERROR: pip's dependency resolver does not currently take into account all the packages that are installed. This behaviour is the source of the following dependency conflicts.
conda-repo-cli 1.0.41 requires requests==2.28.1, but you have requests 2.28.2 which is incompatible.

#9 (comment)

#9 (comment)

could you help me out , how to pick first break? I have read all of the codes. but only functions are defined. Didn't understand what are the functions are to be called ?