Frequency plot and analysis

Dear all,
I want to produce a frequency plot and analysis the frequency content of my seismogram as below but I am not idea how to do that.

Look here: obspy.imaging - Plotting routines for ObsPy — ObsPy 1.4.0.post0+118.gebdcda16a9.obspy.master documentation

Or you could do it via PPSD and it’s plot spectrogram routine (bit of a different approach): obspy.signal.spectral_estimation.PPSD.plot_spectrogram — ObsPy 1.4.0.post0+118.gebdcda16a9.obspy.master documentation

Example plots with links to code in gallery for both:
https://docs.obspy.org/master/gallery.html#gallery