CATALOG.FIT を healpix フィット画像に変換しようとしています。私は適合カタログを読み、ra と dec を theta と phi に変換することができます。しかし、ang2pix を使用すると、エラーがスローされます」
Error encountered at /home/faisal/healpy/healpy-healpy-aa5f605/hpbeta
/Healpix_cxx/healpix_base.cc, line 829
(function I T_Healpix_Base<I>::loc2pix(double, double, double, bool) const [with I =
long int])
起こしてはならない
「PlanckError」のインスタンスをスローした後に呼び出される終了
-------------------
import numpy as np
import string as s
from scipy.interpolate import interp1d
import matplotlib.pyplot as plt
import healpy as hp
from numpy import sin, cos, arccos, arcsin, arctan2
import astro_util as au
#start healpy---------------------------
nvss_ra=np.arange(1,100,1)
nvss_dec=np.arange(1,100,1)
print nvss_ra
theta,phi= au.euler(nvss_ra,nvss_dec,1)
print theta
theta1=theta/(57.3248)
phi1=phi/(57.3248)+((3.14159265359)/2)
nvssmap=np.zeros(hp.nside2npix(512))
#Working till this point
pix=hp.pixelfunc.ang2pix(512,theta1, phi1)
nvssmap[pix]=1
hp.fitsfunc.write_map("nvss_map.fits",nvssmap)
誰かが同様の問題に直面した場合は返信してください。
ありがとう
ファイサル