Ejemplos de parse_NSHA2018_catalogue en Python

Ejemplo n.º 1

#from hmtk.parsers.catalogue.csv_catalogue_parser import CsvCatalogueParser
#from hmtk.seismicity.utils import haversine
from catalogue.parsers import parse_NSHA2018_catalogue
from misc_tools import dictlist2array, timedelta2days_hours_minutes, toYearFraction, ymd2doy
from numpy import array, where, hstack, delete
import matplotlib.pyplot as plt
from datetime import datetime as dt

#hmtk_csv = '/nas/gemd/ehp/georisk_earthquake/modelling/sandpits/tallen/NSHA2018/catalogue/data/AUSTCAT_V0.12_hmtk_deblast.csv'
nsha_csv = 'data/NSHA18CAT.MW.V0.1.csv'

# parse HMTK csv
#parser = CsvCatalogueParser(hmtk_csv)
#nshacat = parser.read_file()

nshacat = parse_NSHA2018_catalogue(nsha_csv)
mx_orig = dictlist2array(nshacat, 'mx_orig')
mx_rev_ml = dictlist2array(nshacat, 'mx_rev_ml')
mw_pref = dictlist2array(nshacat, 'prefmag')
evdt = dictlist2array(nshacat, 'datetime')
ev_type = dictlist2array(nshacat, 'ev_type')
lat = dictlist2array(nshacat, 'lat')
lon = dictlist2array(nshacat, 'lon')

# get indexes to delete
delidx = where((ev_type == 'blast') | (ev_type == 'coal'))[0]
delidx = hstack((delidx, where(lat > -12)[0]))
delidx = hstack((delidx, where(lon < 135)[0]))  # just eastern events
#delidx = hstack((delidx, where(lon > 134)[0])) # just western events
datelim = dt(1900, 1, 1)
delidx = hstack((delidx, where(evdt < datelim)[0]))

Ejemplo n.º 2

Archivo: augment_nsha18_cat.py Proyecto: treviallen/manuscripts

from misc_tools import dictlist2array
from numpy import array, where, unique
from datetime import datetime
from os import path, mkdir, getcwd

################################################################################
# parse nsha catalogue
################################################################################

if getcwd(),startswith('/nas'):
    catfile = '/nas/active/ops/community_safety/ehp/georisk_earthquake/modelling/sandpits/tallen/NSHA2018/catalogue/data/NSHA18CAT.MW.V0.1.csv'
else:
    catfile = '/Users/trev/Documents/Geoscience_Australia/NSHA2018/catalogue/data/NSHA18CAT.MW.V0.1.csv'

# parse catalogue
cat = parse_NSHA2018_catalogue(catfile)

# get data arrays
nsha18_mw = dictlist2array(cat, 'prefmag')
mw_src = dictlist2array(cat, 'mw_src')
ml_region = dictlist2array(cat, 'ml_region')
nsah18_auth = dictlist2array(cat, 'auth')
nsah18_lats = dictlist2array(cat, 'lat')
nsah18_lons = dictlist2array(cat, 'lon')
nsah18_deps = dictlist2array(cat, 'dep')
nsah18_dt = dictlist2array(cat, 'datetime')

################################################################################
# parse updated GA catalogue
################################################################################