Forum der Vereinigung der Sternfreunde

Hallo zusammen,
der Befehl compute/barycorr hat mir heute bei der gleichen Beobachtung zwei leicht verschiedene Korrekturwerte für Total barycentric RV correction geliefert. Lasse ich Datum und Uhrzeit aus der bdf-Datei lesen, so ist der Wert um 0,4 km/s verschieden, als wenn ich Datum und Uhrzeit von Hand in den Befehl einsetze.

Anzeige bei Datum und Uhrzeit aus der Datei (Stationskoordinaten sind in rvcorr.prg enthalten):
Midas 005> (==>)(==>) rvcorr 18,13,45.8 -21,03,32
current dir. now: /home/jottschi/spektro/analyse
no CDELTi nor CDi_j matrix found... STEPi set to 1.0
FITS file: tmp00001.fts converted to: barycorr.bdf
Barycentric correction time: 5.817411185035E-03 day
Heliocentric correction time: 5.816522465156E-03 day

Total barycentric RV correction: -3.89611 km/s
Total heliocentric RV correction: -3.906156 km/s
(incl. diurnal RV correction of -0.0776 km/s)

Descriptor O_TIME of image barycorr.bdf used for date and UT.
current dir. now: /home/jottschi/spektro

Anzeige bei händischer Eingabe:
Midas 003> compute/barycorr 2011,07,03 23,16,53.1 18,13,45.8 -21,03,32 +07,59,33 +47,07,23
Barycentric correction time: 5.804280576029E-03 day
Heliocentric correction time: 5.803359576328E-03 day

Total barycentric RV correction: -4.324439 km/s
Total heliocentric RV correction: -4.33448 km/s
(incl. diurnal RV correction of -0.027816 km/s)

Ich habe mir dann den O_TIME descriptor angesehen. Bei der Wandlung von fts nach bdf wird aus 2011-07-03T23-16-53.096 der Wert 2011.5013698630

Gebe ich diesen Wert bei der händischen Eingabe ein, so erhalte ich wieder den Korrekturwert, als würde compute/barycorr den Wert selbst aus dem header der bdf-Datei lesen.

Läuft da bei der Wandlung was schief?

Ratlos,
Jörg!

_________________
Wer hat's gemacht? .... Die Schweizer ...

Hallo Jörg,

On 12/15/2011 08:50 PM, Jörg Schirmer wrote: Da scheint mit den Zeiten im Header was nicht zu stimmen.
0.5013698630*365 = 183. Das sind also nur die vollen Tage. Die Stunden
müssen noch woanders im o_time header stehen und sind vermutlich bei der
Umwandlung verloren gegangen. Das erste Ergebnis passt auch zu UT=0h.

Midas 019> compute/barycorr 2011,07,03 0,0,0 18,13,45.8 -21,03,32
+07,59,33 +47,07,23
Barycentric correction time: 5.817411185035E-03 day
Heliocentric correction time: 5.816522465156E-03 day

Total barycentric RV correction: -3.89611 km/s
Total heliocentric RV correction: -3.906156 km/s
(incl. diurnal RV correction of -0.0776 km/s)

Es sieht also so aus, als ob Midas nur das Datum, aber nicht die Uhrzeit
korrekt aus dem Header liest. Wie sieht denn der FITS-Header genau aus?
Herzliche Grüße,
Otmar

Sali Otmar,

und vielen Dank. Tja, mein altes Leiden ist immer noch aktiv. Wieder einmal nicht genau genug hingeschaut, um die Ursache selbst finden zu können.
Gut, bleibt aber noch die Tatsache, dass die Beobachtungszeit nicht vollständig übernommen wird. Deshalb anhängend der Fits-Header als Textausgabe.

Schöne Grüsse,
Jörg!

SIMPLE = T/CCDSOFT-SOFTWARE BISQUE 3
BITPIX = 16
NAXIS = 2
NAXIS1 = 1530
NAXIS2 = 1020
BSCALE = +1.000000000000e+000
BZERO = +3.276800000000e+004
BIAS = 100
FOCALLEN= +0.000000000000e+000
APTAREA = +0.000000000000e+000
APTDIA = +0.000000000000e+000
DATE-OBS= '2011-07-03T23:16:53.096'
TIME-OBS= '23:16:53.096 '
SWCREATE= 'CCDSoft Version 5.00.186'
SET-TEMP= -1.000000000000e+001
COLORCCD= 0
DISPCOLR= 1
IMAGETYP= 'Light Frame '
CCDSFPT = 1
XORGSUBF= 0
YORGSUBF= 0
CCDSUBFL= 0
CCDSUBFT= 0
XBINNING= 1
CCDXBIN = 1
YBINNING= 1
CCDYBIN = 1
EXPSTATE= 293
CCD-TEMP= -9.787817267812e+000
TEMPERAT= -9.787817267812e+000
INSTRUME= 'SBIG ST-8 3 CCD Camera'
EGAIN = +2.470000000000e+000
E-GAIN = +2.470000000000e+000
XPIXSZ = +9.000000000000e+000
YPIXSZ = +9.000000000000e+000
SBIGIMG = 6
USER_2 = 'SBIG ST-8 3 CCD Camera'
DATAMAX = 65535
SBSTDVER= 'SBFITSEXT Version 1.0'
EXPTIME = +6.000000000000e+002
EXPOSURE= +6.000000000000e+002
CBLACK = 47
CWHITE = 698
HISTORY Darkframe subtracted.
END

_________________
Wer hat's gemacht? .... Die Schweizer ...

Hallo Jörg,

2011/12/16 Jörg Schirmer <fg-spek-comp@vdsastro.de>: Midas hätte gerne einen Header TM-START mit der Startzeit in Sekunden.
Ich hänge mal ein kleines Perl-Skript an, dass diesen Header einfügt,
modifiziert aus meinem Baukasten

Das Skript z.B. als fixheader5.pl speichern, ausführbar machen (chmod
+x) und dann aufrufen mit

fixheader5.pl file*fits

file*fits ist die Liste der Files, oder der Name eines Files. Bitte
vorher die Files sichern! Das Skript überschreibt die Originale! Damit
sich die Länge des Headers nicht ändert, überschreibe ich die Zeile
FOCALLEN, in der offensichtlich nichts vernünftiges steht. Das kannst
Du auch ersetzen.

Herzliche Grüße,
Otmar

Am Samstag, den 17.12.2011, 10:02 +0100 schrieb Otmar Stahl: Das ist auch in SMS eingebaut.

günter

_________________
Im längsten Frieden spricht der Mensch nicht so viel Unsinn und Unwahrheit als im kürzesten Kriege. (Jean Paul)

Hallo Otmar,

vielen Dank für die rasche Hilfe. Als ich deinen Beitrag las, fiel mir genau das ein, was Günter in seinem Beitrag mitteilt und was auch in seiner Anleitung zu sms steht. Und siehe da, als ich compute/barycorr auf eine von sms bearbeitete Datei anwendete, bekam ich genau den Wert, als hätte ich Datum und Zeit von Hand eingegeben.
Gut, dann will ich mal meine RV-Werte korrigieren lassen!

Schöne Grüsse,
Jörg!

_________________
Wer hat's gemacht? .... Die Schweizer ...

Hallo an die Experten,

wie kann ich denn dann das Ergebnis aus compute/bary auf das Spektrum anwenden? Ist es ok, wenn man dies nach dem rebinning macht rebin/line macht?

Danke und Gruß,
Daniel

_________________
Daniel P. Sablowski

https://www.kunstmann.de/buch/axel_hack ... 42000/t-2/
https://dpsablowski.wordpress.com
https://github.com/DPSablowski
https://www.aip.de/mitglieder/dsablowski/

Hallo Daniel,
rebin/lin ist ok, wenn der Wellenlängenbereich klein ist. Ansonsten musst Du nicht-linear rebinnen. Ich habe dazu mal das Skript unten benutzt.

Herzliche Grüße,
Otmar

!
! rebin file to heliocentric velocity, keep number of pixels
!
set/format i5 g15.8
define/local cvel/d/1/1 2.99792458d5
!
define/local fact/d/1/1 0.0
define/local n/i/1/1 0

defi/para p1 ? ? "input file"
defi/para p2 ? ? "output file"
defi/para p3 ? ? "helvel"
!
!
fact = 1.0/(1.0+{p3}/cvel)           ! factor (1+v/c)
!
write/out rebininng file {p1} with velocity {p3} to new file {p2}
! 
rebi/ii {p2} {p1} {p1} lin 0,{fact} spg             ! rebin

Hi Otmar,

danke für die schnelle Antwort!

Hat geholfen. Aber ich denke ich baue das dann direkt in mein Tool ein, MIDAS ist da doch bissl umständlich . . .

LG,
Dani

_________________
Daniel P. Sablowski

https://www.kunstmann.de/buch/axel_hack ... 42000/t-2/
https://dpsablowski.wordpress.com
https://github.com/DPSablowski
https://www.aip.de/mitglieder/dsablowski/