# Power spectral density multi-station 

**URL:** https://discourse.obspy.org/t/power-spectral-density-multi-station/2185
**Category:** Uncategorized
**Created:** [March 14, 2025, 7:50am UTC](https://discourse.obspy.org/t/power-spectral-density-multi-station/2185 "2025-03-14T07:50:45Z")
**Posts on this page:** 2
**Page:** 1

<div class="post-metadata">

### Author: ![Ulug\_bek\_Abdujabbaro](https://discourse.obspy.org/user_avatar/discourse.obspy.org/ulug_bek_abdujabbaro/32/948_2.png) [@Ulug\_bek\_Abdujabbaro](https://discourse.obspy.org/u/Ulug_bek_Abdujabbaro)
#### Post date: [March 14, 2025, 7:50am UTC](https://discourse.obspy.org/t/power-spectral-density-multi-station/2185/1 "2025-03-14T07:50:45Z")

</div>

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

 ![20250305_185524](https://discourse.obspy.org/uploads/default/original/1X/44a637834b382fd71893dc348f42c673673dff01.jpeg)  
I was able to build a PSD for one channel of one station  
 ![IMG_20250311_144617_497](https://discourse.obspy.org/uploads/default/original/1X/0975e598107c81f15299b1c5e43279e67b0796bd.png)

---

<div class="post-metadata">

### Author: ![megies](https://discourse.obspy.org/user_avatar/discourse.obspy.org/megies/32/9_2.png) [@megies](https://discourse.obspy.org/u/megies)
#### Post date: [March 14, 2025, 8:08am UTC](https://discourse.obspy.org/t/power-spectral-density-multi-station/2185/2 "2025-03-14T08:08:57Z")

</div>

You could use [this routine](https://docs.obspy.org/master/packages/autogen/obspy.signal.spectral_estimation.PPSD.get_percentile.html) 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`.