| #!/bin/csh | |
| # ++++++++++++++++++++++ | |
| # segy2ascii | |
| # Matt Hall, October 2009 | |
| # Converts SEGY file to a TIFF, one sample per pixel | |
| # First, export a SEGY file from Poststack, putting Trace in byte 13 | |
| # Note the name of the project and the name of the file | |
| # Run this file and follow the instructions | |
| # ++++++++++++++++++++++ | |
| # |
| """ | |
| Run as a script: | |
| python test_check.py /path/to/myfile.seg2 | |
| Adapted from https://github.com/agile-geoscience/jeepr | |
| """ | |
| from struct import unpack | |
| import sys |
| import re | |
| import urllib | |
| from io import BytesIO | |
| # If the fname you gave me looks like a URL... | |
| if re.match(r'https?://.+\..+/.+?', fname) is not None: | |
| try: | |
| data = urllib.request.urlopen(fname).read() | |
| except urllib.error.HTTPError as e: |
| #!/usr/bin/env python | |
| """ | |
| Reads special core analysis and core description from SPWLA formatted files. | |
| Usage: | |
| ./spwla.py infile.spwla outfile.csv | |
| Author: | |
| Matt Hall, ca. 2021 | |
| © agilescientific.com |
You'll need to install Python somehow before using this script.
I recommend installing Miniconda
Then do this in the terminal (or Anaconda prompt on Windows):
conda create -n segy2wav python=3.8 scipy
conda activate segy2wav
pip install click segyio
python segy2wav.py --help
| from scipy.signal import remez, minimum_phase | |
| import matplotlib.pyplot as plt | |
| # Based on this example from the scipy docs: | |
| # https://docs.scipy.org/doc/scipy/reference/generated/scipy.signal.minimum_phase.html | |
| freq = [0, 0.2, 0.3, 1.0] | |
| desired = [1, 0] | |
| h_linear = remez(151, freq, desired, Hz=2.0) | |
| h_min_hil = minimum_phase(h_linear, method='hilbert') |