我有一个星系列表可以绘制到healpix地图上(我使用healpy来做)每个星系都有一个通量集,我需要让它们以这样一种方式绘图:通量对于每个星系在地图上是守恒的。 这是我的代码: import numpy as np
import matplotlib.pyplot as plt
import healpy as hp
pi = np.pi
nside = 8
xsize = 100
我有一个银河坐标的地图,我需要将它保存在另一个文件的赤道坐标中。我知道我可以使用: import healpy as hp
map=hp.read_map('file.fits')
map_rot=hp.mollview(map, coord=['G','C'], return_projected_map=True)
,应返回存储在map_rot二维数组numpy。但是当我读到map_ro