График спектра мощности на основе набора данных с использованием healpy - PullRequest
Новая
       74

График спектра мощности на основе набора данных с использованием healpy

0 голосов
/ 19 сентября 2019

Я только начинаю использовать healpy в python.Вкратце, у меня есть данные о положении let 

x = np.array([7.47098722,  7.47805867,  7.47098779....12.57804461, 12.5809734])

y = np.array([58.32516929, 58.33223959, 58.33516828, ..., 63.41415155, 63.41707922])

, где phi (x, y) - потенциал вещества в этом месте.Теперь я хочу рассчитать спектр мощности для этого.Но я не понимаю, как установить мои данные для healpy, чтобы сделать карту, используя mollview и вычисляя спектр мощности, используя anafast.

1 Ответ

0 голосов
/ 23 сентября 2019

Вам просто нужно преобразовать тэты и phis в индексы пикселей, используя одну из healpy функций пикселей, healpy.pixelfunc.ang2pix.Вы можете прочитать больше об этом из https://healpy.readthedocs.io/en/latest/generated/healpy.pixelfunc.ang2pix.html

После того, как у вас есть индексы, у вас также есть потенциальные значения для этих индексов, другими словами, у вас есть соответствующая карта HEALPix (используйте mollview, чтобы получить mollview проекция этой карты).Наконец, используйте anafast, https://healpy.readthedocs.io/en/latest/healpy_spht.htm от Healpy, и вы получите спектр мощности.

Для получения дополнительной информации о healpy, пожалуйста, обратитесь к https://buildmedia.readthedocs.org/media/pdf/healpy/1.8.6/healpy.pdf

https://healpix.sourceforge.io/ также, посетите https://healpix.jpl.nasa.gov

...