healpy

    1熱度

    1回答

    當healpy工作,我能夠用繪製在Mollview一個Healpix地圖 import healpy map = 'filename.fits' healpy.visufunc.mollview(map) 或教程 >>> import numpy as np >>> import healpy as hp >>> NSIDE = 32 >>> m = np.arange(hp.nsi

    1熱度

    1回答

    我已經嘗試調試我的代碼,並且我已經意識到,當我嘗試將我的AltAz座標保存到.csv文件時,它最終會崩潰,因爲它不是一個numpy數組,它是一個SkyCoord對象。有人可能會提出一個簡單的方法來將大型赤道座標錶轉換爲AltAz,或者我可以如何將我的代碼保存到文件中。 # Get time now time = astropy.time.Time.now() time.delta_ut1_ut

    0熱度

    1回答

    我想從適合文件訪問極化並使用Healpy繪製它。我看到Healpix有一個名爲POLARIZATION的內置選項,但我沒有看到Healpy的相同選項。 如何繪製這些數據? 謝謝! 編輯: @AndreaZonca 這是我目前使用的繪製溫度代碼: t = "%s Field - WMAP at W Band (90 GHz) Foreground Reduced iquma"%(name) map

    2熱度

    1回答

    我剛剛開始與Healpy一起工作,並且注意到如果我使用地圖來獲取alm,然後使用這些alm來生成新地圖,則我不會獲得我開始使用的地圖。以下是我正在看的: import numpy as np import healpy as hp nside = 2 # healpix nside parameter m = np.arange(hp.nside2npix(nside)) # create

    0熱度

    1回答

    我是HEALPix的新手,也是Python的新手。我嘗試使用healpy將HEALPix索引轉換爲RA,Dec。 ,我感到我必須使用pix2ang,但找不出如何輸出θ,Φ轉換爲RA,十二月......我想這: import healpy as hp import numpy as np theta, phi = hp.pix2ang(256, 632668 ,nest=True) ra= p

    1熱度

    1回答

    我想用自定義顏色替換「hp.mollview」給出的顏色條。我特別感興趣的是: 旋轉90度的彩條(即由一個垂直替換水平) 使用兩個標籤(左和彩條的右邊) 設置自定義蜱 指示通過設置「cmap.set_over」來設置範圍(通過「max」參數)。 很少的代碼量: import numpy as np import healpy as hp m = np.arange(hp.nside2npix(

    0熱度

    1回答

    我有一系列具有兩個不同索引i和j的地圖。讓這個索引像map_series [i] [j]。 編輯1/21:的最小工作的例子是像 map_series=np.array([np.array([np.arange(12) + 0.1*(i+1) + 0.01*(j+1) for j in range(3)]) for i in range(5)]) 我想相同的掩模應用到每一個;如果map_series

    4熱度

    3回答

    我有一個HEALPix地圖,我已經閱讀過使用healpy,但它在銀河座標系中,我需要它在天體/赤道座標系中。有誰知道轉換地圖的簡單方法嗎? 我已經嘗試使用healpy.Rotator從(l,b)轉換爲(phi,theta),然後使用healpy.ang2pix重新排列像素,但該圖仍然看起來很奇怪。 如果有類似於Rotator的功能,您可以這樣稱呼它:map = AnotherRotator(map

    1熱度

    1回答

    我試圖在healpix地圖上使用healpy來產生一個波束。對於初學者,我希望能夠在mollweide投影中產生2D高斯,但我真的不知道從哪裏開始。 我可以定義一個二維高斯: import numpy as np def gaussian_2D(x,y,mu_x=0.,mu_y=0.,sig_x=1.,sig_y=1.): return np.exp(-0.5*(((x-mu_x)/s

    3熱度

    2回答

    我想對healpy.mollview 使用別人的colormaps我succeded與此代碼 from healpy import mollview from pylab import arange, show, cm m = arange(768) mollview(m, cmap=cm.bwr) show() ,但我得到一個意外的藍色背景,有沒有辦法,我可以將其設置爲白色