Проблемы с использованием Healpix - PullRequest
Новая
       29

Проблемы с использованием Healpix

0 голосов
/ 01 мая 2019

Я начинаю использовать healpix и пытаюсь построить источники из каталога планок. У меня есть файл FITS с информацией, такой как галактические координаты, и я хочу создать карту в проекции Моллвейд.

Сначала я использовал функцию «healpy.pixelfunc.ang2pix» для преобразования координат в пиксели. Один из параметров, который я до сих пор не понимаю: NSIDE, поэтому я просто поставил случайное число, кратное 2. Во-вторых, я использовал «healpy.visufunc.mollview» с параметром map, который представляет собой вектор пикселей, который я получил из предыдущего функция. Но я получаю это сообщение об ошибке "плохое количество пикселей". Я не знаю, что я делаю неправильно .. 

    import healpy as hp
    from astropy.io import fits
    import numpy as np
    from astropy.table import Table
    import copy as copy

    catalogue=fits.open("/home/luana/Documents/Mestrado/Planck_data/COM_PCCS_070_R2.04.fits")
    COORD=Table(catalogue[1].data)
    nside=np.power(2,4)
    GLON=copy.copy(COORD['GLON'])
    GLAT=copy.copy(COORD['GLAT'])
    hp_pixel_numbers=hp.pixelfunc.ang2pix(nside,GLON,GLAT,nest=True,lonlat=True)     
    hp.visufunc.mollview(map=hp_pixel_numbers, fig=None, rot=None, nest=True,cbar=True, cmap=None, notext=False, norm=None, hold=False, margins=None, sub=None, nlocs=2, return_projected_map=False)
...