我正在尝试使用 python 绘制银河坐标图。假设我有这些数据:
数据 = [(0.261, -7.123, 13.03, '未识别'), (-0.326, 77, 13.03, '星系')]
其中每个元组的形式为(ra、dec、flux、type)。
我被要求使用astropy + matplotlib,所以:
c = SkyCoord(ra = ra*u.degree, dec = dec*u.degree)
galactic = c.galactic
这是我的问题出现的地方,我正在使用以下代码:
from mpl_toolkits.basemap import Basemap
import numpy as np
import matplotlib.pyplot as plt
# lon_0 is central longitude of projection.
# resolution = 'c' means use crude resolution coastlines.
m = Basemap(projection='hammer',lon_0=0,resolution='c')
m.drawcoastlines()
m.fillcontinents(color='coral',lake_color='aqua')
# draw parallels and meridians.
m.drawparallels(np.arange(-90.,120.,30.))
m.drawmeridians(np.arange(0.,420.,60.))
m.drawmapboundary(fill_color='aqua')
plt.title("Hammer Projection")
plt.show()
但是我无法在银河坐标中绘制数据,我也不知道为什么。我还需要每个点根据类型具有不同的颜色,并且根据通量的值具有不同的大小。我需要实现这样的目标(我对python有点陌生,我从未使用过astropy,我还没有找到好的例子):
希望有人能帮忙。