2015-10-13 81 views
6

私は銀河座標のプロットをパイソンを使って作ろうとしています。のは、私はこのデータを持っているとしましょう:Astrophy matplotlibとプロット銀河座標

データ= [(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を使用したことがない、私は良い例が見つかっていない)、このような何かを達成する必要があります:

enter image description here

誰かが助けることができる希望。

+0

を「私は銀河の座標でデータをプロットすることができず、なぜわからないのですか?何を試しましたか?エラーメッセージは何ですか、またはポイントが表示されませんでしたか? – keflavich

答えて

6

http://www.astropy.org/astropy-tutorials/plot-catalog.htmlの下部にある例を見てください。私はプロット銀河座標ときに実行する一般的な問題は、あなたが-180から+180にプロットしたいということですが、デフォルトではあなたがwrap_atでこれを変更することができ、0から360までの座標、例えば与えることです:

plot(galactic.l.wrap_at(180*u.deg), galactic.b.wrap_at(180*u.deg)) 
関連する問題