Wie zeichne ich galaktische Koordinaten mit Matplotlib und Astropie in Python?

Question

view on stackexchange narkive permalink

Ich möchte eine Projektion von Millisekundenpulsaren auf der galaktischen Ebene machen, ähnlich wie diese von Sala et al. 2004:

Ich habe verschiedene Methoden ausprobiert und bin nirgendwo hingekommen. Dies ist mein aktueller Code, zusammen mit dem, was er produziert:

  importiere Mathematik als mimport numpy als npimport csvimport matplotlib.pyplot als pltimport pandas als pdimport astropy.koordinaten als coordimport astropy.units als ufrom astropy .io importieren asciifrom astropy.coordinates importieren SkyCoorddata = pd.read_csv ('galacticwperiod.csv') xarr = np.array (data.iloc [:, 0]) yarr = np.array (data.iloc [:, 1]) eq = SkyCoord (xarr [:], yarr [:], unit = u.deg) gal = eq.galactic # print (xarr) #xarr = np.array (df.iloc [:, 0]) # yarr = np .array (df.iloc [:, 1]) # zarr = np.array (df.iloc [:, 2]) # ra = coord.Angle (xarr [:], unit = u.hour) # ra.degree #ra = ra.wrap_at (180 * u.degree) #dec = coord.Angle (yarr [:], unit = u.deg) #print (ra) plt.figure (figsize = (6,5)) fig = plt.figure () ax = fig.add_subplot (111, projection = "aitoff") plt.plot (gal.l.wrap_at (180 * u.deg), gal.b.wrap_at (180 * u.deg), linestyle = 'Keine') ax.scatter (gal.l, gal.b, linestyle = 'Keine') # ax.set_facecolor ('xkcd: Schlachtschiff grau') # fig.patch.set_facecolor ('weiß') # ax.tick_params (Achse = 'beide', w hich = 'major', labelize = 10) # ax.grid (color = 'b', linestyle = 'solid') fig.show () # plt.savefig ('millisecondcoloraitoff.png', dpi = 600)  

Hier einige Zeilen der Eingabedatei 'galacticwperiod.csv':

  Gl, Gb111.383, -54.849305.898, -44.888305.913, - 44.877  

Dies ist das Bild, das es erzeugt:

Ich bin fast sicher, dass dies falsch ist, da sie nicht entlang der verteilt sind galaktische Ebene, die sie sein sollten. Die Daten, die ich verwende, stammen aus dem ATNF Pulsar Catalog.

Dies sind die Websites, die ich mir bereits als Referenz angesehen habe:

https://stackoverflow.com/questions/33105898/astropy-matplotlib-and-plot- Galaktische Koordinaten https://astropy4scipy2014.readthedocs.io/_static/notebooks/06_Celestial_Coordinates_solutions.html

Jede Hilfe dabei wäre sehr dankbar.

Answer 1

view on stackexchange narkive permalink

Dieser Code liest Koordinaten als äquatorial (ra, dec) und wandelt sie in galaktische (l, b) um:

  eq = SkyCoord (xarr [:], yarr [:], Einheit = u.deg) gal = eq.galactic  

Der Inhalt von 'galacticwperiod.csv' ist bereits in galaktischen Koordinaten angegeben und sollte nicht transformiert werden. So etwas kann zu besseren Ergebnissen führen:

  gal = SkyCoord (xarr [:], yarr [:], frame = 'galactic', unit = u.deg)  

Das andere Problem ist das Die geografischen Projektionen von pyplot scheinen Winkel im Bogenmaß zu erwarten. Dieser Code:

  plt.subplot (111, projection = 'aitoff') plt.grid (True) plt.scatter (gal.l.wrap_at ('180d'). Radian, gal.b.radian)  

erzeugt dieses Diagramm:

Answer 2

view on stackexchange narkive permalink

Wenn ich mir die Handlung ansehe, stelle ich fest, dass sich an den "Polen" eine Konzentration von Punkten befindet. Die Konzentration sollte entlang des Äquators und insbesondere in Richtung des Zentrums der Galaxie liegen.

Ich denke, Sie haben Ihre Achsen umgeschaltet, sodass Punkte entlang des galaktischen Äquators entlang des vertikalen "Nullmeridians" gezeichnet werden.

Answer 3

view on stackexchange narkive permalink

  aus astropy.io importieren passt aus astropy importieren wcs importieren matplotlib.pyplot als pltx = fit.open ('f160w_noopt_drz.fits') ['SCI'] y = wcs.WCS (x.header) plt. Unterzeichnungen (figsize = (8, 8), subplot_kw = {'Projektion': y}) plt.imshow (x.data, cmap = 'viridis', origin = 'lower', vmin = 0, vmax = 10) plt. show ()