Modul 1 Auswahl der Sterne von M4 mit Python: mp_auswahl.py zeigt in pyplot nichts an

append delete 12 Jan ’22 · 12:51 Gab

Ich denke ich bin allen Anweisungen gefolgt, allerdings zeigt Pyplot beim Ausführen keine Daten an? Es ist nur ein leeres Fenster zu sehen.

Trotz dessen, dass Pyplot einen pmra von "-5" bis "-20", und pmdec von "-7" bis "-25" anzeigen sollte, zeigt es einen ganz anderen Bereich an. ([-0.04;0.04])
Auch beim Navigieren zum eigentlichen Bereich, ist nichts zu sehen.

Hier ist der Code:

import numpy as np
import matplotlib.pyplot as plt

katalog = np.genfromtxt("katalog_um_M4.tsv", delimiter="\t", names=True)

pmra_min = -5
katalog = np.extract(katalog["ra"] > pmra_min , katalog)
pmra_max = -20
katalog = np.extract(katalog["ra"] < pmra_max , katalog)
pmdec_min = -7
katalog = np.extract(katalog["dec"] > pmdec_min , katalog)
pmdec_max = -25
katalog = np.extract(katalog["dec"] < pmdec_max , katalog)

plt.plot(katalog["pmra"], katalog["pmdec"],
marker=".", markersize=1, linestyle="None")

plt.title("Gaia-Eigenbewegungen im Himmelsauschnitt um M4")
plt.xlabel("Eigenbewegung in RA [mas/year]")
plt.ylabel("Eigenbewegung in Dec [mas/year]")

plt.show()

Antworten

12 Jan ’22 · 13:29 #1. Gab

This post was deleted by a moderator

12 Jan ’22 · 14:55 #2. AndromedaMalte

This post was deleted by its owner

append delete 20 Jan ’22 · 13:46 #3. AndromedaMalte

Zusammenfassung des Threads: hier gab es zwei Probleme
Problem 1 ist, dass du katalog["ra"] und nicht katalog["pmra"] in dem np.extract() verwendest. Wir wollen die Auswahl ja in Eigenbewegung (also pm...) und nicht in Himmelskoordinaten treffen. Das gleiche gilt natürlich für katalog["dec"].
Problem 2 ist, dass pmra_min grösser ist als pmra_max (das gleiche bei pmdec). Die Bedingung, die du in np.extract() schreibst, ist eine Bedingung um den Stern in der Auswahl zu behalten, und nicht um ihn rauszuschmeissen. Daher muss katalog["pmra"] *größer* als pmra_min sein, und *kleiner* als pmra_max. Und das geht nur wenn pmra_min < pmra_max.

Nützlicher tipp: nach der Auswahl mal ein mal ein print(katalog) setzen (oder sogar ein print(len(katalog)), um zu sehen ob da noch etwas übrigbleibt. Dann kann man schon mal sehen, ob das Problem an der Auswahl oder am Diagram liegt.

Antworten