Skip to content

Instantly share code, notes, and snippets.

View braunfuss's full-sized avatar

Andreas Steinberg braunfuss

View GitHub Profile
@braunfuss
braunfuss / mtqt_source_plot.py
Last active January 4, 2024 15:08
mtqt_source_plot
from pyrocko.gf import meta
from pyrocko.gf.seismosizer import Source
from pyrocko import gf
from pyrocko.guts import Float
from pyrocko import moment_tensor as mtm
import numpy as num
pi = num.pi
pi4 = pi / 4.
km = 1000.
@braunfuss
braunfuss / seigerdown.py
Created July 14, 2021 14:46
Download data from the seismological resarch network Rhineland lower palatinate
#!/usr/bin/env pyrocko-python
import sys
import logging
import tempfile
import math
import os.path as op
import shutil
try:
from urllib.error import HTTPError
@braunfuss
braunfuss / gist:2b3fdb351b43abe40c6275ef9d11b308
Last active June 5, 2021 16:47
plot_street_orientation_city.py
import matplotlib.pyplot as plt
import numpy as np
import osmnx as ox
'''
Slimmed down example plot of city street network orientations for a single city,
instead of multiple as in the osmnx notebook.
Needs the awesome package osmnx:
https://github.com/gboeing/osmnx
@braunfuss
braunfuss / cake2locsat_shm.py
Created June 4, 2021 11:00
Convert pyrocko cake format travel time tables to locsat or Seismic Handler readable tables
import numpy as num
from pyrocko import cake
from silvertine.util.ref_mods import *
km = 1000.
def convert2shm_ttt(mod):
# define cake phase to be used
phases = cake.PhaseDef.classic('Pg')
@braunfuss
braunfuss / locsat2pyrocko.py
Created June 4, 2021 10:59
Convert locasat output to pyrocko Phasemarkers.
from pyrocko import model
import numpy as np
from pyrocko.gui.pile_viewer import PhaseMarker, EventMarker
from pyrocko import cake, util
from pyrocko.gui import marker as pym
from pyrocko.util import stt
def monthToNum(shortMonth):
@braunfuss
braunfuss / download_strassburg_event.py
Created May 26, 2021 16:50
Download data and visualize for the Strassburg event on 26.05.2021 using pyrocko
from pyrocko.client import fdsn
from pyrocko import util, io, trace, model
from pyrocko.io import quakeml
fdsn.g_timeout = 60.
# Start and end time
tmin = util.stt('2021-05-26 13:40:00')
tmax = util.stt('2021-05-26 14:22:00')
@braunfuss
braunfuss / snappy_proc.py
Last active July 15, 2022 20:13
Sentinel-1 processing with snappy (snaps python api) for one swath (https://step.esa.int/main/toolboxes/snap/)
import os
import snappy
from snappy import GPF
from snappy import ProductIO
from snappy import HashMap
from snappy import jpy
import subprocess
from time import *
# documentation: http://step.esa.int/docs/v2.0/apidoc/engine/overview-summary.html