Noise Power Spectrum

Hello,

I am looking for a way to plot the noise power spectrum for a Catphan504 (though I assume the process is similar for other CatPhans). I know that the results_data() method provides values like the maximum and the mean, but I am specifically interested in generating the graph.

I attempted the method described in the documentation, but it didn’t work as expected. Is there an alternative way to access and plot this graph?

Regards,
Safir

from matplotlib import pyplot as plt

from pylinac import CatPhan504
from pylinac.core.nps import plot_nps1d

ct = CatPhan504.from_demo_images()
ct.analyze()
plot_nps1d(ct.ctp486.power_spectrum_1d)
plt.show()

Thank you you’re the best!