Hi everyone,
I am looking for obspy function to plot the amplitude spectra of the vertical waveform.
Here is my script
import numpy as np
from obspy import read
from obspy.clients.fdsn import Client
from obspy.geodetics.base import gps2dist_azimuth
from obspy import UTCDateTime
from datetime import datetime
from datetime import timedelta
import pandas as pd
client = Client("IRIS")
T_P=events["P_Time"]
ev_lat=events["lat"]
ev_long=events["long"]
ev_dep=events["depth"]
station_latitude=45.90
station_longitude=-130
PGV_Z=[]
for a in range (110, 111,1):
time=T_P[a]
time=UTCDateTime(time)
depth=ev_dep[a]
starttime=time-timedelta(minutes=10)
endtime= time+timedelta(minutes=30)
net = "OO"
sta = "AXEC2"
loc = "*"
chan = "HHZ"
st = client.get_waveforms(net, sta, loc, chan, starttime, endtime, attach_response= True);
st_rem=st.copy()
st_tr=st_rem.remove_response(output = 'VEL',plot=False);
st_filt=st_tr.filter('bandpass', freqmin =0.001, freqmax=0.1)
st_f=st_filt.detrend('demean').detrend('linear').taper(0.05)
st_f.plot(equal_scale=False, automerge=False);
data=st_f[0].data
pgv=max(data)
PGV_Z.append(float(pgv))
print(PGV_Z)
Here is what I required.