Seite 1 von 2

DATE-OBS keyword wird verändert

BeitragVerfasst: Sonntag 22. Februar 2015, 11:45
von Michael
Hallo!

Ich möchte die proper-motion bei einem Kometenfoto berücksichtigen. Doch egal was ich dort eingebe: es wird ein normales Coadd auf die Sterne und nicht auf den Kometen gemacht.
Nachdem ich die mittlerweile exzellent gewordenen Theli-Doku ansehen habe, dachte ich es liegt an den DATE-OBS keywords im FITS-Header.
Und tatsächlich, so ist es!

Im Aufnahme.FITS steht zB:
DATE-OBS= '2015-02-16T20:26:28' / YYYY-MM-DDThh:mm:ss observation start, UT

und im *OFC.fits steht nur noch:
DATE-OBS= 2015

Hat vielleicht Jemand einen Tipp wie ich das DATE-OBS Keyword elegant per Terminal-Befehl von meinen Aufnahme.FITS in die 3 R/G/B_1OFC.fits kopieren kann?

btw: ich verwende ubuntheli mit der THELI GUI 2.9.5

Danke,
Michael

Re: DATE-OBS keyword wird verändert

BeitragVerfasst: Sonntag 22. Februar 2015, 15:23
von mischa
Hi Michael,

lad mal dein process_split_*.sh Skript hier hoch, und das logfile der "preparation".

mischa

Re: DATE-OBS keyword wird verändert

BeitragVerfasst: Sonntag 22. Februar 2015, 16:25
von Michael
Aber gerne, die Dateien findest du im Anahng.

Ich habe jetzt mühsam mittels "replacekey" den (meiner Meinung nach korrekten String) in jedes der *1OFC.fits geschrieben, doch leider bringt das nichts da ja zur Coadd die 1OFC.sub.fits herangezogen werden.
Also fix nochmal skysub laufen lassen und jetzt verschwindet der DATE-OBS Eintrag vom 1OFC.fits zum 1OFC.sub.fits :mrgreen:
Darum habe ich auch gleich das logfile vom skysub angehängt.

Header von einem *1OFC.fits:

SIMPLE = T / file does conform to FITS standard
BITPIX = 16 / number of bits per data pixel
NAXIS = 2 / number of data axes
NAXIS1 = 3036 / length of data axis 1
NAXIS2 = 2012 / length of data axis 2
EXTEND = T / FITS dataset may contain extensions
COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy
COMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H
AIRMASS = 1. / (average) Airmass during observation
BAYERPAT = 'RGGB '
CD1_1 = -0.00222222 / WCS Coordinate scale matrix
CD1_2 = 0. / WCS Coordinate scale matrix
CD2_1 = 0. / WCS Coordinate scale matrix
CD2_2 = 0.00222222 / WCS Coordinate scale matrix
COEFF = 1. / photometric extiction coefficient
CRPIX1 = 1520. / WCS Coordinate reference pixel
CRPIX2 = 1008. / WCS Coordinate reference pixel
CRVAL1 = 28.36250000
CRVAL2 = 50.78722222
CTYPE1 = 'RA---TAN' / WCS coordinate type
CTYPE2 = 'DEC--TAN' / WCS coordinate type
DATAMD5 = '0a93847c325b28861aa32c216f655b31' / MD5 checksum
DATE-OBS = '2015-02-16T19:26:44'
DUMMY11 = 0 / Dummy for adding new FITS cards later
DUMMY12 = 0 / Dummy for adding new FITS cards later
DUMMY13 = 0 / Dummy for adding new FITS cards later
DUMMY14 = 0 / Dummy for adding new FITS cards later
DUMMY15 = 0 / Dummy for adding new FITS cards later
DUMMY16 = 0 / Dummy for adding new FITS cards later
DUMMY17 = 0 / Dummy for adding new FITS cards later
DUMMY18 = 0 / Dummy for adding new FITS cards later
DUMMY19 = 0 / Dummy for adding new FITS cards later
DUMMY20 = 0 / Dummy for adding new FITS cards later
DUMMY21 = 0 / Dummy for adding new FITS cards later
DUMMY22 = 0 / Dummy for adding new FITS cards later
DUMMY23 = 0 / Dummy for adding new FITS cards later
DUMMY24 = 0 / Dummy for adding new FITS cards later
DUMMY25 = 0 / Dummy for adding new FITS cards later
DUMMY26 = 0 / Dummy for adding new FITS cards later
DUMMY27 = 0 / Dummy for adding new FITS cards later
DUMMY28 = 0 / Dummy for adding new FITS cards later
DUMMY29 = 0 / Dummy for adding new FITS cards later
DUMMY3 = 0 / Dummy for adding new FITS cards later
DUMMY30 = 0 / Dummy for adding new FITS cards later
DUMMY31 = 0 / Dummy for adding new FITS cards later
DUMMY32 = 0 / Dummy for adding new FITS cards later
DUMMY33 = 0 / Dummy for adding new FITS cards later
DUMMY34 = 0 / Dummy for adding new FITS cards later
DUMMY35 = 0 / Dummy for adding new FITS cards later
DUMMY36 = 0 / Dummy for adding new FITS cards later
DUMMY37 = 0 / Dummy for adding new FITS cards later
DUMMY38 = 0 / Dummy for adding new FITS cards later
DUMMY39 = 0 / Dummy for adding new FITS cards later
DUMMY4 = 0 / Dummy for adding new FITS cards later
DUMMY40 = 0 / Dummy for adding new FITS cards later
DUMMY41 = 0 / Dummy for adding new FITS cards later
DUMMY42 = 0 / Dummy for adding new FITS cards later
DUMMY43 = 0 / Dummy for adding new FITS cards later
DUMMY44 = 0 / Dummy for adding new FITS cards later
DUMMY45 = 0 / Dummy for adding new FITS cards later
DUMMY46 = 0 / Dummy for adding new FITS cards later
DUMMY47 = 0 / Dummy for adding new FITS cards later
DUMMY48 = 0 / Dummy for adding new FITS cards later
DUMMY49 = 0 / Dummy for adding new FITS cards later
DUMMY5 = 0 / Dummy for adding new FITS cards later
DUMMY50 = 0 / Dummy for adding new FITS cards later
DUMMY6 = 0 / Dummy for adding new FITS cards later
DUMMY7 = 0 / Dummy for adding new FITS cards later
DUMMY8 = 0 / Dummy for adding new FITS cards later
DUMMY9 = 0 / Dummy for adding new FITS cards later
EQUINOX = 2000. / Equinox of coordinates
EXPTIME = 180. / effective exposure time
FILTER = 'Blue'
GABODSID = 821 / Obs. date in days since 31/12/1998
IMAGEID = 1 / Chip Number
OBJECT = 'UNKNOWN ' / observed target
RADECSYS = 'FK5 ' / Coordinate system for equinox (FK4/FK5/GAPPT)
ZP = -1. / photometric zeropoint
INPUTFMT = 'FITS ' / Format of file from which image was read
SWOWNER = 'michi ' / Licensed owner of software
BSCALE = 1.0000000000000000 /physical = BZERO + BSCALE*array_value


und von dazugehörigen *1OFC.sub.fits

SIMPLE = T / file does conform to FITS standard
BITPIX = 16 / number of bits per data pixel
NAXIS = 2 / number of data axes
NAXIS1 = 3036 / length of data axis 1
NAXIS2 = 2012 / length of data axis 2
EXTEND = T / FITS dataset may contain extensions
COMMENT FITS (Flexible Image Transport System) format is defined in 'Astronomy
COMMENT and Astrophysics', volume 376, page 359; bibcode: 2001A&A...376..359H
AIRMASS = 1. / (average) Airmass during observation
BAYERPAT = 'RGGB '
CD1_1 = -0.00222222 / WCS Coordinate scale matrix
CD1_2 = 0. / WCS Coordinate scale matrix
CD2_1 = 0. / WCS Coordinate scale matrix
CD2_2 = 0.00222222 / WCS Coordinate scale matrix
COEFF = 1. / photometric extiction coefficient
CRPIX1 = 1520. / WCS Coordinate reference pixel
CRPIX2 = 1008. / WCS Coordinate reference pixel
CRVAL1 = 28.3625
CRVAL2 = 50.78722222
CTYPE1 = 'RA---TAN' / WCS coordinate type
CTYPE2 = 'DEC--TAN' / WCS coordinate type
DATAMD5 = '0a93847c325b28861aa32c216f655b31' / MD5 checksum
DUMMY11 = 0 / Dummy for adding new FITS cards later
DUMMY12 = 0 / Dummy for adding new FITS cards later
DUMMY13 = 0 / Dummy for adding new FITS cards later
DUMMY14 = 0 / Dummy for adding new FITS cards later
DUMMY15 = 0 / Dummy for adding new FITS cards later
DUMMY16 = 0 / Dummy for adding new FITS cards later
DUMMY17 = 0 / Dummy for adding new FITS cards later
DUMMY18 = 0 / Dummy for adding new FITS cards later
DUMMY19 = 0 / Dummy for adding new FITS cards later
DUMMY20 = 0 / Dummy for adding new FITS cards later
DUMMY21 = 0 / Dummy for adding new FITS cards later
DUMMY22 = 0 / Dummy for adding new FITS cards later
DUMMY23 = 0 / Dummy for adding new FITS cards later
DUMMY24 = 0 / Dummy for adding new FITS cards later
DUMMY25 = 0 / Dummy for adding new FITS cards later
DUMMY26 = 0 / Dummy for adding new FITS cards later
DUMMY27 = 0 / Dummy for adding new FITS cards later
DUMMY28 = 0 / Dummy for adding new FITS cards later
DUMMY29 = 0 / Dummy for adding new FITS cards later
DUMMY3 = 0 / Dummy for adding new FITS cards later
DUMMY30 = 0 / Dummy for adding new FITS cards later
DUMMY31 = 0 / Dummy for adding new FITS cards later
DUMMY32 = 0 / Dummy for adding new FITS cards later
DUMMY33 = 0 / Dummy for adding new FITS cards later
DUMMY34 = 0 / Dummy for adding new FITS cards later
DUMMY35 = 0 / Dummy for adding new FITS cards later
DUMMY36 = 0 / Dummy for adding new FITS cards later
DUMMY37 = 0 / Dummy for adding new FITS cards later
DUMMY38 = 0 / Dummy for adding new FITS cards later
DUMMY39 = 0 / Dummy for adding new FITS cards later
DUMMY4 = 0 / Dummy for adding new FITS cards later
DUMMY40 = 0 / Dummy for adding new FITS cards later
DUMMY41 = 0 / Dummy for adding new FITS cards later
DUMMY42 = 0 / Dummy for adding new FITS cards later
DUMMY43 = 0 / Dummy for adding new FITS cards later
DUMMY44 = 0 / Dummy for adding new FITS cards later
DUMMY45 = 0 / Dummy for adding new FITS cards later
DUMMY46 = 0 / Dummy for adding new FITS cards later
DUMMY47 = 0 / Dummy for adding new FITS cards later
DUMMY48 = 0 / Dummy for adding new FITS cards later
DUMMY49 = 0 / Dummy for adding new FITS cards later
DUMMY5 = 0 / Dummy for adding new FITS cards later
DUMMY50 = 0 / Dummy for adding new FITS cards later
DUMMY6 = 0 / Dummy for adding new FITS cards later
DUMMY7 = 0 / Dummy for adding new FITS cards later
DUMMY8 = 0 / Dummy for adding new FITS cards later
DUMMY9 = 0 / Dummy for adding new FITS cards later
EQUINOX = 2000. / Equinox of coordinates
EXPTIME = 180. / effective exposure time
FILTER = 'Blue '
GABODSID = 821 / Obs. date in days since 31/12/1998
IMAGEID = 1 / Chip Number
OBJECT = 'UNKNOWN ' / observed target
RADECSYS = 'FK5 ' / Coordinate system for equinox (FK4/FK5/GAPPT)
ZP = -1. / photometric zeropoint
HISTORY /home/theli/THELI/theli//bin/Linux//subsky -i //home/theli/Lovejoy//Ligh
HISTORY t/Lovejoy_001_B_1OFC.fits -s /home/theli/Lovejoy//Light/Lovejoy_001_B_1O
HISTORY FC_seg.fits -o /home/theli/Lovejoy//Light/Lovejoy_001_B_1OFC.sub.fits -n
HISTORY 5000 -c 1 -f 256 -t -1e8 1e8 -w //home/theli/Lovejoy//WEIGHTS/Lovejoy_0
HISTORY 01_B_1OFC.weight.fits -d
INPUTFMT = 'FITS ' / Format of file from which image was read
SWOWNER = 'michi ' / Licensed owner of software
BSCALE = 1.0000000000000000 /physical = BZERO + BSCALE*array_value

Re: DATE-OBS keyword wird verändert

BeitragVerfasst: Sonntag 22. Februar 2015, 17:53
von mischa
Hi Michael,

das DATE-OBS keyword ist irrelevant. Ausschlaggebend ist nur MJD (modified julian date), welches beim splitting aus dem DATE-OBS berechnet wird.Vermutlich ist letzteres in deinem splitting script nicht aktiviert. Ich kann dir morgen eine aktualisierte version schicken.

mischa

Re: DATE-OBS keyword wird verändert

BeitragVerfasst: Sonntag 22. Februar 2015, 18:24
von Michael
Hi Mischa,

ja bitte das wäre spitze.

Danke,
Michael

Re: DATE-OBS keyword wird verändert

BeitragVerfasst: Montag 23. Februar 2015, 15:23
von mischa
Hi Michael,

Probier's mal damit?

mischa

Re: DATE-OBS keyword wird verändert

BeitragVerfasst: Montag 23. Februar 2015, 17:00
von Michael
Bin mir nicht zu 100% sicher ob das jetzt stimmt.
Im Header finde ich jetzt nach der Preparation folgende Header-Einträge:

DATE-OBS= 2015-02-16T20:00:28
MJD-OBS = 2015-02-16T20:00:28

Stimmt das Format von MJD-OBS? Das ist doch kein modified julian date?

Re: DATE-OBS keyword wird verändert

BeitragVerfasst: Montag 23. Februar 2015, 23:47
von Michael
nach der Calibration habe ich in den 1OFC.fits files folgendes im Fits-Header stehen:

DATE-OBS= 2015
MJD-OBS = 2015


Das verunsichert mich jetzt eigentlich noch mehr. Glaubst du es macht Sinn so die Daten weiter zu verarbeiten?

Re: DATE-OBS keyword wird verändert

BeitragVerfasst: Dienstag 24. Februar 2015, 13:08
von mischa
hi michael,

sorry, war ein schluderer meinerseits. Ersetze bitte in Zeile 144

Code: Alles auswählen
${tmp}


durch

Code: Alles auswählen
${MJD}


mischa

Re: DATE-OBS keyword wird verändert

BeitragVerfasst: Mittwoch 25. Februar 2015, 14:26
von Michael
Ausgezeichnet, jetzt funktioniert das mit dem MJD-OBS Header-Eintrag.
Auch in den Skysubs ist der richtige Wert eingetragen und er wird auch bei jeder Aufnahme größer.
Vielen Dank! :D


Alleine das Coadd ist noch immer auf Sterne gestackt und nicht auf den Kometenkern.
Ich habe schon viele unterschiedliche proper-motion Werte ausprobiert, es wird aber immer auf die Sterne gestackt.

Hast du vielleicht eine Erklärung dafür?
Im Anhang das Coadd-Log.


Danke,
Michael

Re: DATE-OBS keyword wird verändert

BeitragVerfasst: Mittwoch 25. Februar 2015, 19:09
von mischa
Hi Michael,

ich seh grad dass in dem beetreffenden script sowohl DATE-OBS als auch MJD verwendet wird. Kannst du mal in deinen headern nachschauen, ob das DATE_OBS ebenfalls richtig gesetzt wurde? mir scheint, da steht immer noch nur 2015 drin.
Ausserdem: Du hast fuer die Eigenbewegung in RA und DEC jeweils -2.0000 eingetragen, das sieht mir etwas unwahrscheinlich aus?

mischa

Re: DATE-OBS keyword wird verändert

BeitragVerfasst: Mittwoch 25. Februar 2015, 19:13
von Michael
In DATE-OBS steht überall 2015 drinnen, da wird das Datum anscheinend bei irgendeinem Schritt beschnitten.
In den Ausgangs.fits ist das Datum korrekt eingetragen.

Die 2 in RA/DEC waren nur einer der verzweifelten und den Tränen nahen Versuche die Sterne in irgendein Eck zu bewegen :lol:
Die korrekten Werte habe ich, aber auch damit "bewegt" sich nichts.

Re: DATE-OBS keyword wird verändert

BeitragVerfasst: Donnerstag 26. Februar 2015, 11:22
von Michael
Habe jetzt nochmal bei Null angefangen:
DATE-OBS ist nach der Preparation noch ok.
DATE-OBS ist nach der Calibration auf 2015 beschnitten.

Re: DATE-OBS keyword wird verändert

BeitragVerfasst: Donnerstag 26. Februar 2015, 13:27
von Michael
Mein Linux-Know-How ist sehr bescheiden, aber ich denke ich habe den Fehler eingegrenzt:

Sobald ich die Bilder "debayere"
Code: Alles auswählen
${P_FITSDEMOSAICBAYER} -i ${FILE} -p ${BAYERPATTERN} -q 2
steht in DATE-OBS nur noch "2015" drinnen.

Lasse ich mir im Terminal den DATE-OBS Wert eines Fits files mit getkey anzeigen so erhalte ich nur 2015 obwohl der Wert definitiv ein komplettes Datum+Uhrzeit ist.
Vielleicht verwendet ja fitsdemosaicbayer intern getkey?

Code: Alles auswählen
theli@draco:~/Lovejoy/Light$ dfits Lovejoy_001_1.fits | fitsort DATE-OBS
FILE                 DATE-OBS              
Lovejoy_001_1.fits   2015-02-16T19:26:44   
theli@draco:~/Lovejoy/Light$ getkey Lovejoy_001_1.fits DATE-OBS
2015
theli@draco:~/Lovejoy/Light$

Re: DATE-OBS keyword wird verändert

BeitragVerfasst: Donnerstag 26. Februar 2015, 13:29
von mischa
Michael hat geschrieben:Habe jetzt nochmal bei Null angefangen:
DATE-OBS ist nach der Preparation noch ok.
DATE-OBS ist nach der Calibration auf 2015 beschnitten.


Du meinst nach dem du "Calibrate data" gemacht hast, und bevor du was anderes tust? Kannst du mir mal dein ~/.theli/scripts/theli_cal_comm.sh.log schicken?

Die Verwendung von DATE-OBS bei der Koaddition ist ueberfluessig. Ersetze mal in
Code: Alles auswählen
gui/scripts/prepare_coadd_swarp.sh
die Zeilen 220-223 wie folgt:

ALT:
Code: Alles auswählen
timezero=`${P_DFITS} ${firstimage} | ${P_FITSORT} -d DATE-OBS | ${P_GAWK} '{print $2}'`
timenow=`${P_DFITS} ${file} | ${P_FITSORT} -d DATE-OBS | ${P_GAWK} '{print $2}'`
timezero=`${P_MJD} -f ${timezero} | ${P_GAWK} '{print $7}'`
timenow=`${P_MJD} -f ${timenow} | ${P_GAWK} '{print $7}'`


NEU:
Code: Alles auswählen
timezero=`${P_DFITS} ${firstimage} | ${P_FITSORT} -d MJD-OBS | ${P_GAWK} '{print $2}'`
timenow=`${P_DFITS} ${file} | ${P_FITSORT} -d MJD-OBS | ${P_GAWK} '{print $2}'`
# timezero=`${P_MJD} -f ${timezero} | ${P_GAWK} '{print $7}'`
# timenow=`${P_MJD} -f ${timenow} | ${P_GAWK} '{print $7}'`


Damit sollte dann der Eigenbewegungsvektor korrekt ins koaddierte Bild uebertragen werden.

mischa