DATE-OBS keyword wird verändert

[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

[mischa]

Hi Michael,

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

mischa

[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
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

[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

[Michael]

Hi Mischa,

ja bitte das wäre spitze.

Danke,
Michael

[mischa]

Hi Michael,

Probier's mal damit?

mischa

[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?

[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?

[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

[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!


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

[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

[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
Die korrekten Werte habe ich, aber auch damit "bewegt" sich nichts.

[Michael]

Habe jetzt nochmal bei Null angefangen:
DATE-OBS ist nach der Preparation noch ok.
DATE-OBS ist nach der Calibration auf 2015 beschnitten.

[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$


[mischa]

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