[PYTHON] Un moyen facile de tracer un diagramme HR à l'aide d'une requête
Diagramme HR et astroquery
Présentation d'un tracé simple utilisant la méthode de tracé du diagramme Hertzsprung-Russell (communément appelé diagramme HR) en utilisant la base de données astronomique spatiale appelée astroquery. Faire. Les gens qui ne sont pas bons dans ce domaine peuvent se familiariser avec le traçage pour le moment.
Pour ceux qui n'ont besoin que du code, veuillez consulter la page google Colab.
installez astroquery
python
pip install astroquery
Ensuite, entrez astroquery.
Obtenir des données à l'aide d'une requête
Les données proviennent d'une base de données appelée Vizir
- https://vizier.u-strasbg.fr/viz-bin/VizieR?-source=1239/hip_main
Utilisez les données du satellite Hipparcos dans.
python
from astroquery.vizier import Vizier
v = Vizier(catalog="I/239/hip_main",columns=['HIP',"Vmag","B-V","Plx"],row_limit=-1)
data = v.query_constraints()
vmag = data[0]["Vmag"]
bv = data[0]["B-V"]
plx = data[0]["Plx"]
Vous n'avez besoin que de trois.
--Vmag: luminosité de la bande V --BV: Rapport de luminosité entre la bande B et la bande V (notez qu'il s'agit d'un rapport car la note est log bien qu'elle soit écrite en B-V.) --plx: dioptrie annuelle. J'ai juste besoin d'informations sur la distance. Puisque la distance des étoiles proches est mesurée par triangulation, convertissez la parallaxe annuelle en distance.
Est.
La signification de l'option est de spécifier le catalogue avec catalog = "I / 239 / hip_main", de spécifier les colonnes avec les colonnes = ['HIP', "Vmag", "BV", "Plx"] et row_limit = - Obtenez toutes les données en supprimant la limite supérieure avec 1.
Créer un diagramme RH
Tracez une étoile avec une dilatation annuelle de 20 à 25 mili arcsec (c'est-à-dire une distance d'environ 40 à 50 pc).
def to_parsec(marcsec):
return 1./(1e-3*marcsec) # pc
pmin=20 # m arcsec
pmax=25 # m arcsec
dmax = to_parsec(pmin)
dmin = to_parsec(pmax)
dlabel = str("%3.1f" % dmin) + " pc to " + str("%3.1f" % dmax) + " pc"
print (dlabel)
vmag_cut = vmag[ (plx > pmin) & (plx < pmax)]
bv_cut = bv [ (plx > pmin) & (plx < pmax)]
from matplotlib import pyplot as plt
plt.title(dlabel)
plt.xlim(-0.1,2)
plt.ylim(-13,0)
plt.xlabel("B-V")
plt.ylabel("Vmag")
plt.scatter(bv_cut, -vmag_cut, s=1)
Puis
De cette manière, le diagramme HR est appliqué. Et si je ne coupe pas à distance? Veuillez essayer. B-V est
python
B-V \equiv log(Luminosité de la bande B)-log(Luminosité de la bande V) = log(Luminosité de la bande B/Luminosité de la bande V)
Parce que cela signifie, faisons attention au rapport de luminosité.
Page connexe
Recommended Posts