Astronomie Logo
Diskussionsforen Spektrum Hintergrund abziehen wie | Spektroskopie
Optionen
Thema bewerten
Pauldd Offline
Mitglied

Registriert: 13/05/2012
Beiträge: 148
Ort: Dresden
#1272625 - 03/03/2017 09:20 Spektrum Hintergrund abziehen wie


Hallo

Ich versuchs einfach mal auf Gut-Glück ob mir hier jemand weiterhelfen kann.

Ich habe ein Tageslichtspektrum aufgenommen und wollte es gerne analysieren.



Man erkennt sehr gut die Frauhnhofer Linien und auch die Natrium Doppellinie. Da ich ausschließlich mit linuxeigenen Mitteln arbeite habe ich mir überlegt wie ich das Spektrum in Zahlen verwandeln kann um damit zu arbeiten. Dazu habe ich das Spektrum mit "convert" auf ein minimales Maß ge-croppt.



Dann habe ich das Spektrum mit "convert spektrum.png spektrum.txt" in Zahlen zerlegt. Nun habe ich für jedes Pixel ein RGB Wert.

Dann habe ich ein Script geschrieben das jeweils zu jedem X-Pixel einen Durchschnitt aller Y-Pixel berechnet, um sozusagen den Mittelwert in der Vertikalen zu bekommen.
Damit bekomme ich eine X Reihe mit gemittelten Farbwerten.

Diese X-Reihe Plotte ich dann mit Gnuplot. Der Einfachheit halber habe ich mal nur den roten Kanal geplottet.



Man erkennt einwandfrei die zwei Natrium Linien (D) und rechts daneben Ha(C) und O2(B).

Meine eigentliche Frage ist nun wie kann ich aus dem Plot folgenden Plot machen der mir nur die Linien anzeigt, also der Rest rausgerechnet wird? Ich habe das mal schwarz nachgezeichnet wie ich das meine:



Könnte man das nicht so machen indem man einen gleitenden Durchschnitt zu der ganzen Reihe berechnet und dann jedem Wert zum Durchschnitt vergleicht und ab einer bestimmten Abweichung vom Durchschnitt den Wert extrahiert? Oder habt ihr eine Andere Idee?
_________________________
--Mensch ohne Windows--

[zum Seitenanfang]
mkoenig66 Offline
Mitglied

Registriert: 12/04/2003
Beiträge: 1539
Ort: Rimbach Odw.
#1276135 - 25/03/2017 10:07 Re: Spektrum Hintergrund abziehen wie [Re: Pauldd]


Hallo Paul,

du bist auf dem richtigen Weg, allerdings ist das Abziehen eines "running mean" nicht sinnvoll, denn damit läufst du Gefahr zu viele Details/Linien mit zu entfernen, wobei dies von der Breite deines Mittelungsfensters ist.

Am besten wäre es, wenn du ein Kontinuum anfittest und dieses dann abziehst, dann bleiben "die Linien übrig". Die Modellierung des Kontinuums gelingt - je nachdem, welche Tools du zur Verfügung hast - ganz gut mit einem Polynom n-ten Grades. Hier gibts es in den Mathematik-Bibliotheken meist entsprechende Fit-Funktionen.

Schöne Grüsse,
Michael
_________________________
www.astro-images.de

[zum Seitenanfang]
Pauldd Offline
Mitglied

Registriert: 13/05/2012
Beiträge: 148
Ort: Dresden
#1316868 - 28/12/2017 13:01 Re: Spektrum Hintergrund abziehen wie [Re: Pauldd]


So, nach längerem rumtrödeln (9 Monate biggrin ) in Libreoffice Calc und dessem Polynom bin ich zu Python gewechselt.

Ich habe jetzt die Rohdaten in eine Liste gelesen, dann das Polynom 10ten Grades erstellt und eine zweite Liste mit dessen y Werten gefüllt. Dieser Liste habe ich dann die Rohdaten abgezogen und herausgekommen ist das:



Ich habe keine Ahnung ob das so gut aussieht. Irgendwie konnte das Ergebnis noch ein bisschen deutlicher nach oben skaliert werden.

Dann werde ich mal versuchen das ganze Spektrum auszugeben und die X-Achse richtig mit Wellenlänge anstatt nur der x-pixel position.


Bearbeitet von Pauldd (28/12/2017 13:07)
_________________________
--Mensch ohne Windows--

[zum Seitenanfang]
TorstenHansen Offline
Mitglied

Registriert: 24/03/2007
Beiträge: 1371
Ort: Unterallgäu
#1317163 - 30/12/2017 15:39 Re: Spektrum Hintergrund abziehen wie [Re: Pauldd]


Hallo Paul,

wenn Du Deine "raw data" durch Dein "10th deg. poly." dividierst, bleiben das (Pseudo-)Kontinuum und die Linien als Absorptionslinien (wie im Rohspektrum) übrig.

Grüße
Torsten

[zum Seitenanfang]
Pauldd Offline
Mitglied

Registriert: 13/05/2012
Beiträge: 148
Ort: Dresden
#1318021 - 06/01/2018 10:53 Re: Spektrum Hintergrund abziehen wie [Re: Pauldd]


Danke Torsten, für den Tipp.
Ich hätte noch eine Frage. Hier mal ein RGB Plot der Rohdaten.



Wie man sieht hat sich im blauem Bereich (Richtung UV) noch ein Rest Rot/Grünanteil und im rotem Bereich (Richtung IR) ein Rest Blau/Grün versteckt. Sehe ich das richtig das ich diese Rester einfach rauschneiden kann so dass nur noch die drei großen RGB Keulen übrig bleiben?
_________________________
--Mensch ohne Windows--

[zum Seitenanfang]
Pauldd Offline
Mitglied

Registriert: 13/05/2012
Beiträge: 148
Ort: Dresden
#1318410 - 09/01/2018 08:44 Re: Spektrum Hintergrund abziehen wie [Re: Pauldd]


Und gleich noch eine Frage. Kann ich eigentlich das ganze RGB Spektrum nicht in Mein Monochombild umwandeln? Oder gegen dabei Informationen verloren? Eigentlich müssten doch die Absorptionslinien dabei erhalten bleiben.
_________________________
--Mensch ohne Windows--

[zum Seitenanfang]


Moderator:  Bernd_Hoffmann 
Aktuelle Fotos
Hier die Milchstraße mit 35mm
Neues von La Palma
Messier 64 - NGC 4826
Cirrus-Nebel / veil nebula
Aufnahmeplanung La Palma
Astrourlaub La Palma
Sonne am 14.6.18
Wer ist Online
63 Mitglieder (Auriga_HH, astra1400, Alex1981, Andi8307, backyardastro, 17 unsichtbar), 511 Gäste und 5 Suchmaschinen online.
Status: Admin, GlobalModerator, Moderator
Forum Statistik
27710 Mitglieder
68 Foren
254907 Themen
1328047 Beiträge

Besucher Rekord: 1886 @ 17/10/2009 22:24