Power spectral density multi-station

How to build multi-station power spectral density in Obspy with mseed files and StationXml files.
Look like this photo


I was able to build a PSD for one channel of one station

You could use this routine to extract a certain percentile and then make a custom plot with multiple stations from the data. Or you could find a psd at a certain time of interest by looking at ppsd.psd_values and finding the right point in time by comparing the index in the list with ppsd.times_processed.