わかりました-これを行うには2つの方法があります。1つは、距離に応じて目立たない数のビンを用意することです(d <10pc、10pc <d <20pc、d> 20pcなど)。これは比較的簡単です。必要なのはいくつかのループだけです。これは3の例です。
raclose = []
ramid = []
rafar = []
decdlose = []
decmid = []
decfar = []
for ii in range(len(dist)):
if dist[ii] < 10.:
raclose.append(ra[ii])
decclose.append(dec[ii])
elif dist[ii] > 20.:
rafar.append(ra[ii])
decfar.append(dec[ii])
else:
ramid.append(ra[ii])
decmid.append(dec[ii])
plt.clf
ax1 = scatter(raclose, decclose, marker='o', s=20, color="darkgreen", alpha=0.6)
ax2 = scatter(ramid, decmid, marker='o', s=20, color="goldenrod", alpha=0.6)
ax3 = scatter(rafar, decfar, marker='o', s=20, color="firebrick", alpha=0.6)
line1 = Line2D(range(10), range(10), marker='o', color="darkgreen")
line2 = Line2D(range(10), range(10), marker='o',color="goldenrod")
line3 = Line2D(range(10), range(10), marker='o',color="firebrick")
plt.legend((line1,line2,line3),('d < 10pc','20pc > d > 10pc', 'd > 20pc'),numpoints=1, loc=3)
show()
または、x軸にRA、y軸にDecを指定し、プロットに距離を入力するように、等高線図を作成することもできます。RAとDecはどちらも、それぞれの座標を持つ1D配列です。次に、距離を使用して2D配列を作成します。距離の中央値/平均値を決定し、2D配列をその値で割って正規化します。最後に、次のような等高線図(contourfまたはimshowを使用)を使用してプロットします。
import matplotlib.pyplot as plt
from matplotlib import cm
ax = pylab.contourf(RA,Dec,dists, levels=[1, 5, 10, 15], cmap=plt.cm.spectral)
cbar=pylab.colorbar()