Theli-Forum

von **mischa** » Sonntag 23. Oktober 2011, 23:02

Kann man aber das drizzeln irgendwie deaktivieren? Bzw kann man irgendwo einstellen, dass die (End-)Bilder in der gleichen Auflösung wie die Rohbilder gespeichert werden?

nein, abstellen in dem sinn kannst du das resampling nicht. aber du kannst die pixelskala angeben, welche das summenbild haben soll. und diese wird automatisch so eingestellt, dass sie der 'nativen' pixelskala deiner kamera entspricht. wenn das nicht der fall ist, dann hast du dein kamera-config verkorkst, oder versehentlich eine falsche pixelskala under coadd-config angegeben, oder aus irgendeinem grund ein anderes instrument mit deutlich unterschiedlichem abbildungsmasstab ausgewaehlt.

mischa

von **fehlfarbe** » Sonntag 23. Oktober 2011, 23:22

Hallo Mischa,

alles klar, die Pixel Scale steht bei coadd auf 1", obwohl sie ~3,3" beträgt.
Aber gut, dass jetzt dank eurer Hilfe schon mal tolle Bilder heraus kommen! Dauert nur um einiges länger, als das Ganze mal schnell 10 Minuten durch DSS zu schicken

Die Farben sehen aber nach dem kalibrieren ziemlich gut aus.

Grüße,
Marcus

von **fehlfarbe** » Montag 24. Oktober 2011, 14:21

So, das RGB stacken funktioniert wirklich einwandfrei, jetzt auch in der richtigen Größe

Beim H-Alpha Bild spinnt er aber immer noch. Kann es sein, dass der nicht damit klar kommt, dass der G und B Kanal kaum Informationen enthalten?

Mir sind beim RGB Bild außerdem noch mehrere "Löcher" im Bild aufgefallen. Haben da die Flats nicht gepasst oder woher kommen die Dinger?

Grüße,
Marcus

von **moos** » Montag 24. Oktober 2011, 15:39

Hi Markus,
die Löcher stammen vom Global weight.
Schau mal bitte mit dem Skycat nach, ob der globalweight_1.mc.fits dort die werte Null hat.
Dazu gibt es gerade ein Thema http://theli-forum.info/viewtopic.php?f=12&t=56

Demnach solltest du im der Konfiguration für die Global-weights die Parameter fürs defect detection alle leer lassen.

von **fehlfarbe** » Montag 24. Oktober 2011, 15:48

Danke Carsten,

ich werd mal die Standardwerte alle heraus nehmen und schauen, ob die Löcher weg sind.

von **spiegelei** » Montag 24. Oktober 2011, 16:06

Hallo,

mischa hat geschrieben:Da stimmt was mit dem Format oder dem numerischen wert deines "EXPOSURE" keywords nicht, sonst wuerde theli nicht "out of range" schreiben.

EZCAP schreibt die Zeiten in ms raus (warum auch immer). Da läuft dann der Zähler über. Vor kurzem gabs da mal einen Tip, an welcher Stelle in den Sourcen man das hochsetzen kann.

Schau mal im Vrzeichnis /<woauchimmer>/theli/theli/eclipsetools/tools, da gibt's eine Datei preflist_split.h, dort musst Du nach

{"EXPTIME", P_FLOAT, &prefs.exptime, 0, 0, 0.0, 999999.0,
{""}, 1, 1, NULL, 0},

suchen. Statt der 999999 steht da ein viel kleinerer Wert drin, den musst Du hochsetzen und dann alles neu kompilieren.

Gruesze,
Karsten

von **mischa** » Montag 24. Oktober 2011, 17:16

Schau mal im Vrzeichnis /<woauchimmer>/theli/theli/eclipsetools/tools, da gibt's eine Datei preflist_split.h

Das empfehle ich NICHT. das fuehrt zwar auch zu einem (demselben) ergebnis, kann aber unter umstaenden an anderer
stelle probleme machen, wenn man speziellere funktionen in theli benutzen will (die von sekunden ausgehen, nicht
millisekunden).
der richtige weg ist, das header-keyword im splitting-script zu korrigieren, und das wurde ja auch gemacht.

lg

mischa

von **Knickohr** » Montag 24. Oktober 2011, 19:02

Hallo Leute,

bitte nicht hier im Beitrag verzetteln !

Ursprünglich ging es mal um "ALCCD6 Fehler bei Preparation"

Damit das Forum übersichtlich bleibt, sollte nicht alles in einen einzigen Beitrag abgehandelt werden.

Danke,
Thomas

von **fehlfarbe** » Montag 24. Oktober 2011, 20:57

Hi Karsten,

ja das wäre auch eine Idee, aber wie Mischa ja schon sagte, gibt es da bestimmt Probleme mit anderen Funktionen bzw. Kameras. Mit der Ergänzung im Bashscript funktioniert das ja ganz gut, muss dann halt nur mal für jedes Script gemacht werden, aber das geht ja fix.

Die Löcher habe ich nun durch die gelöschten Standardwerte bei "defect detection" weg bekommen, super

Jetzt muss ich Theli nur noch irgendwie überreden den Rotkanal der Ha Bilder zu stacken.
Die "Weights" im Coadd Verzeichnis sehen so komisch aus, während die resamp.fits usw. alle ganz normal aussehen. Zum Vergleich mal die 2 Screens.

Grüße,
Marcus

von **schubi** » Donnerstag 15. Dezember 2011, 19:02

Hallo Marcus,
funktioniert die Bearbeitung der Rohbilder bei Dir ?
Ich musste auf einen neueren Treiber für die ALCCD6cPro umsteigen und dieser macht etwas anderes als die Vorgänger.
Wenn es damit getan ist am Header etwas zu ändern dann o.k.

Ich hänge derzeit mit meinem Theli fest.

Welche Chipgeometrie hast Du bei Deinem Theli eingestellt ?

vg
Christian

von **fehlfarbe** » Donnerstag 15. Dezember 2011, 21:13

Hallo Christian,

durch die Änderung des Scripts funktioniert die Bearbeitung nun einwandfrei. Es lag ja nur daran, dass die EZCap Software die Belichtungszeit in Millisekunden in den FITS Header schreibt.

Meine Einstellungen:

Overscan xmin/max: 0
cut on x/y: 1
x/y size: 3040 2016
CRPIX1/2: 1516 1004

Beim Splitting Script natürlich dann RGB Matrix, FITS und das erste Bayer Pattern.

Der Fehler bei den komischen Ha Bildern wie im letzten Screenshot lag glaube an der COSMICS Einstellung beim weighting. Bin mir aber nicht mehr sicher. Zumindest geht jetzt alles.

Grüße,
Marcus

von **schubi** » Donnerstag 15. Dezember 2011, 23:30

Hallo Marcus,
ich habe die Einstellungen für die Kamera von Dir übernommen. Dennoch funktioniert es nicht.
Hast Du noch manuelle Änderungen an dem Splitting-Skript durchgeführt oder nur das Splitting-Skript aus der GUI heraus neu erstellt ?
Hier das Ini-File:

Code: Alles auswählen: # We need to export these export INSTRUMENT export NCHIPS # Load directories for programs . progs.ini # Geographical latitude of the observatory OBSLAT=50 # Pixel scale in x and y [degrees] (CDELT1/2 keywords) PIXSCX=-0.000575 PIXSCY=0.000575 # Pixel scale in arcsec PIXSCALE=2.07 # Detector gain GAIN=1.0 # Overscan regions. The values are the x-range. # The min is given by OVSCANX1, the max by OVSCANX2. # We need separate regions for each chip. OVSCANX1=([1]=0) OVSCANX2=([1]=0) # Sections for the cutting. Select those pixels that receive light. # Give min values for x and y, and the extent in x and y. CUTX=([1]=1) CUTY=([1]=1) SIZEX=([1]=3040) SIZEY=([1]=2016) # Reference pixels for each chip. REFPIXX=([1]=1516) REFPIXY=([1]=1004)

und hier ein Auszug aus dem Log-File

Code: Alles auswählen: #!/bin/bash -xv # *************************************** # ***** INSTRUMENT splitting script ***** # ***** created by THELI ***** # *************************************** # splits QHY8Pro@TMB130 images # $1: main directory # $2: science directory INSTDIR=instruments_professional + INSTDIR=instruments_professional if [ ! -f ${INSTDIR}/${INSTRUMENT}.ini ]; then INSTDIR=instruments_commercial fi + '[' '!' -f instruments_professional/QHY8Pro@TMB130.ini ']' + INSTDIR=instruments_commercial if [ ! -f ${INSTDIR}/${INSTRUMENT}.ini ]; then INSTDIR=~/.theli/instruments_user/ fi + '[' '!' -f instruments_commercial/QHY8Pro@TMB130.ini ']' + INSTDIR=/home/schubi/.theli/instruments_user/ . ${INSTDIR}/${INSTRUMENT:?}.ini + . /home/schubi/.theli/instruments_user//QHY8Pro@TMB130.ini # ************************************************************** # Instrument config file created by THELI # ************************************************************** # Name of the instrument (must be unique!) INSTRUMENT=QHY8Pro@TMB130 ++ INSTRUMENT=QHY8Pro@TMB130 # Number of chips NCHIPS=1 ++ NCHIPS=1 # We need to export these export INSTRUMENT ++ export INSTRUMENT export NCHIPS ++ export NCHIPS # Load directories for programs . progs.ini ++ . progs.ini # # the following language setting ensures that # we do not have to deal with country specific # settings influencing programs like gawk (giving # a ',' instead of '.' for the decimal point in numbers # within a german environment). export LANG=C +++ export LANG=C +++ LANG=C # # how many processors do we have # used directories PIPESOFT=/home/schubi/THELI/theli/ +++ PIPESOFT=/home/schubi/THELI/theli/ BIN=${PIPESOFT}/bin/Linux_64/ +++ BIN=/home/schubi/THELI/theli//bin/Linux_64/ CONF=${PIPESOFT}/gui/reduction; export CONF +++ CONF=/home/schubi/THELI/theli//gui/reduction +++ export CONF STANDARDSTARSCAT= +++ STANDARDSTARSCAT= PHOTSTANDARDSTARSCAT= +++ PHOTSTANDARDSTARSCAT= SCRIPTS=${PIPESOFT}/scripts/reduction +++ SCRIPTS=/home/schubi/THELI/theli//scripts/reduction UTILSCRIPTS=${PIPESOFT}/scripts/Linux_64/ +++ UTILSCRIPTS=/home/schubi/THELI/theli//scripts/Linux_64/ SMMACROS=${PIPESOFT}/scripts/sm_macros +++ SMMACROS=/home/schubi/THELI/theli//scripts/sm_macros export TEMPDIR=~/.theli/tmp/ +++ export TEMPDIR=/home/schubi/.theli/tmp/ +++ TEMPDIR=/home/schubi/.theli/tmp/ # for compatibility I need the following DATACONF=${CONF}; export DATACONF +++ DATACONF=/home/schubi/THELI/theli//gui/reduction +++ export DATACONF # # define here the location of script interpreters: P_PERL=/usr/bin/perl +++ P_PERL=/usr/bin/perl P_PYTHON=/usr/bin/python +++ P_PYTHON=/usr/bin/python # # here are the used programs listed # # LDAC programs # P_ALBUM=${BIN}/album +++ P_ALBUM=/home/schubi/THELI/theli//bin/Linux_64//album P_ANALYSELDAC=${BIN}/analyseldac +++ P_ANALYSELDAC=/home/schubi/THELI/theli//bin/Linux_64//analyseldac P_ANISOTROPY=${BIN}/anisotropy +++ P_ANISOTROPY=/home/schubi/THELI/theli//bin/Linux_64//anisotropy P_ASCTOLDAC=${BIN}/asctoldac +++ P_ASCTOLDAC=/home/schubi/THELI/theli//bin/Linux_64//asctoldac P_ASSOCIATE=${BIN}/associate +++ P_ASSOCIATE=/home/schubi/THELI/theli//bin/Linux_64//associate P_DECIMALTOHMS=${BIN}/decimaltohms +++ P_DECIMALTOHMS=/home/schubi/THELI/theli//bin/Linux_64//decimaltohms P_DECIMALTODMS=${BIN}/decimaltodms +++ P_DECIMALTODMS=/home/schubi/THELI/theli//bin/Linux_64//decimaltodms P_DMSTODECIMAL=${BIN}/dmstodecimal +++ P_DMSTODECIMAL=/home/schubi/THELI/theli//bin/Linux_64//dmstodecimal P_HMSTODECIMAL=${BIN}/hmstodecimal +++ P_HMSTODECIMAL=/home/schubi/THELI/theli//bin/Linux_64//hmstodecimal P_IC=${BIN}/ic +++ P_IC=/home/schubi/THELI/theli//bin/Linux_64//ic P_IMCOMBFLAT_IMCAT=${BIN}/imcombflat +++ P_IMCOMBFLAT_IMCAT=/home/schubi/THELI/theli//bin/Linux_64//imcombflat P_IMSTATS=${BIN}/imstats +++ P_IMSTATS=/home/schubi/THELI/theli//bin/Linux_64//imstats P_LDACADDMASK=${BIN}/ldacaddmask +++ P_LDACADDMASK=/home/schubi/THELI/theli//bin/Linux_64//ldacaddmask P_LDACCONV=${BIN}/ldacconv +++ P_LDACCONV=/home/schubi/THELI/theli//bin/Linux_64//ldacconv P_LDACPASTE=${BIN}/ldacpaste +++ P_LDACPASTE=/home/schubi/THELI/theli//bin/Linux_64//ldacpaste P_LDACCALC=${BIN}/ldaccalc +++ P_LDACCALC=/home/schubi/THELI/theli//bin/Linux_64//ldaccalc P_LDACADDKEY=${BIN}/ldacaddkey +++ P_LDACADDKEY=/home/schubi/THELI/theli//bin/Linux_64//ldacaddkey P_LDACRENKEY=${BIN}/ldacrenkey +++ P_LDACRENKEY=/home/schubi/THELI/theli//bin/Linux_64//ldacrenkey P_LDACADDTAB=${BIN}/ldacaddtab +++ P_LDACADDTAB=/home/schubi/THELI/theli//bin/Linux_64//ldacaddtab P_LDACRENTAB=${BIN}/ldacrentab +++ P_LDACRENTAB=/home/schubi/THELI/theli//bin/Linux_64//ldacrentab P_LDACFILTER=${BIN}/ldacfilter +++ P_LDACFILTER=/home/schubi/THELI/theli//bin/Linux_64//ldacfilter P_LDACTESTEXIST=${BIN}/ldactestexist +++ P_LDACTESTEXIST=/home/schubi/THELI/theli//bin/Linux_64//ldactestexist P_LDACTOASC=${BIN}/ldactoasc +++ P_LDACTOASC=/home/schubi/THELI/theli//bin/Linux_64//ldactoasc P_LDACTOSKYCAT=${BIN}/ldactoskycat +++ P_LDACTOSKYCAT=/home/schubi/THELI/theli//bin/Linux_64//ldactoskycat P_LDACJOINKEY=${BIN}/ldacjoinkey +++ P_LDACJOINKEY=/home/schubi/THELI/theli//bin/Linux_64//ldacjoinkey P_MAKEJOIN=${BIN}/make_join +++ P_MAKEJOIN=/home/schubi/THELI/theli//bin/Linux_64//make_join P_MAKESSC=${BIN}/make_ssc +++ P_MAKESSC=/home/schubi/THELI/theli//bin/Linux_64//make_ssc P_MAKESUBIMAGE=${BIN}/makesubimage +++ P_MAKESUBIMAGE=/home/schubi/THELI/theli//bin/Linux_64//makesubimage P_PREANISOTROPY=${BIN}/preanisotropy +++ P_PREANISOTROPY=/home/schubi/THELI/theli//bin/Linux_64//preanisotropy P_SELECTOVERLAPS=${BIN}/select_overlaps +++ P_SELECTOVERLAPS=/home/schubi/THELI/theli//bin/Linux_64//select_overlaps P_STATS=${BIN}/stats +++ P_STATS=/home/schubi/THELI/theli//bin/Linux_64//stats P_TRANSFORMIMAGE=${BIN}/transformimage +++ P_TRANSFORMIMAGE=/home/schubi/THELI/theli//bin/Linux_64//transformimage # # Terapix programs # P_SEX=${BIN}/sex_theli +++ P_SEX=/home/schubi/THELI/theli//bin/Linux_64//sex_theli P_WW=${BIN}/ww_theli +++ P_WW=/home/schubi/THELI/theli//bin/Linux_64//ww_theli P_SWARP=${BIN}/swarp_theli +++ P_SWARP=/home/schubi/THELI/theli//bin/Linux_64//swarp_theli # # programs from the eclipse package # P_DFITS=${BIN}/dfits +++ P_DFITS=/home/schubi/THELI/theli//bin/Linux_64//dfits P_FITSORT=${BIN}/fitsort +++ P_FITSORT=/home/schubi/THELI/theli//bin/Linux_64//fitsort P_HIERARCH=${BIN}/hierarch28 +++ P_HIERARCH=/home/schubi/THELI/theli//bin/Linux_64//hierarch28 P_REPLACEKEY=${BIN}/replacekey +++ P_REPLACEKEY=/home/schubi/THELI/theli//bin/Linux_64//replacekey P_IOFITS=${BIN}/iofits +++ P_IOFITS=/home/schubi/THELI/theli//bin/Linux_64//iofits # # own programs developed with # eclipse and qfits # P_FITSSPLIT_ECL=${BIN}/mefsplit +++ P_FITSSPLIT_ECL=/home/schubi/THELI/theli//bin/Linux_64//mefsplit P_FIXCR_ECL=${BIN}/fixcr +++ P_FIXCR_ECL=/home/schubi/THELI/theli//bin/Linux_64//fixcr P_IMRED_ECL=${BIN}/preprocess +++ P_IMRED_ECL=/home/schubi/THELI/theli//bin/Linux_64//preprocess # # For THELI modified versions of the ASTROMETRIX and # PHOTOMETRIX scripts: # S_ASTROMETRIX=/home/schubi/THELI/theli-1.5.1/scripts/astrometrix/astrom_theli +++ S_ASTROMETRIX=/home/schubi/THELI/theli-1.5.1/scripts/astrometrix/astrom_theli S_PHOTOMETRIX=/home/schubi/THELI/theli-1.5.1/scripts/astrometrix/photom_theli +++ S_PHOTOMETRIX=/home/schubi/THELI/theli-1.5.1/scripts/astrometrix/photom_theli # # other utility programs and scripts # P_READLINK=${BIN}/readlink +++ P_READLINK=/home/schubi/THELI/theli//bin/Linux_64//readlink P_AIRMASS=${BIN}/airmass +++ P_AIRMASS=/home/schubi/THELI/theli//bin/Linux_64//airmass P_NIGHTID=${BIN}/nightid +++ P_NIGHTID=/home/schubi/THELI/theli//bin/Linux_64//nightid P_CALDATE=${BIN}/caldate +++ P_CALDATE=/home/schubi/THELI/theli//bin/Linux_64//caldate S_SCAMPCAT=${UTILSCRIPTS}/scampcat.py +++ S_SCAMPCAT=/home/schubi/THELI/theli//scripts/Linux_64//scampcat.py S_PHOTOABS=${UTILSCRIPTS}/photo_abs.py +++ S_PHOTOABS=/home/schubi/THELI/theli//scripts/Linux_64//photo_abs.py S_APLASTROMSCAMP=${UTILSCRIPTS}/aplastrom_scamp.py +++ S_APLASTROMSCAMP=/home/schubi/THELI/theli//scripts/Linux_64//aplastrom_scamp.py # # Programs and scripts not included in the pipeline distribution # P_FIND=/usr/bin/find +++ P_FIND=/usr/bin/find P_GAWK=/usr/bin/gawk +++ P_GAWK=/usr/bin/gawk P_SM= +++ P_SM= P_SORT=/usr/bin/sort +++ P_SORT=/usr/bin/sort # # this is needed for the GUI # P_PIPELIBS=/home/schubi/THELI/PipeLibs//pipe.csh +++ P_PIPELIBS=/home/schubi/THELI/PipeLibs//pipe.csh P_GETPMMOBJECTS=${BIN}/get_pmm_objects +++ P_GETPMMOBJECTS=/home/schubi/THELI/theli//bin/Linux_64//get_pmm_objects P_LDACDESC=${BIN}/ldacdesc +++ P_LDACDESC=/home/schubi/THELI/theli//bin/Linux_64//ldacdesc P_LDACDELTAB=${BIN}/ldacdeltab +++ P_LDACDELTAB=/home/schubi/THELI/theli//bin/Linux_64//ldacdeltab P_LDACDELKEY=${BIN}/ldacdelkey +++ P_LDACDELKEY=/home/schubi/THELI/theli//bin/Linux_64//ldacdelkey P_LDACPUTXY=${BIN}/ldacputxy +++ P_LDACPUTXY=/home/schubi/THELI/theli//bin/Linux_64//ldacputxy P_FITSFLIP=${BIN}/fitsflip +++ P_FITSFLIP=/home/schubi/THELI/theli//bin/Linux_64//fitsflip P_DEBLOOM=${BIN}/debloom +++ P_DEBLOOM=/home/schubi/THELI/theli//bin/Linux_64//debloom P_EQUALISEBAYERFLAT=${BIN}/equalisebayerflat +++ P_EQUALISEBAYERFLAT=/home/schubi/THELI/theli//bin/Linux_64//equalisebayerflat P_FITSCOLLAPSE=${BIN}/fitscollapse +++ P_FITSCOLLAPSE=/home/schubi/THELI/theli//bin/Linux_64//fitscollapse P_FITSWRAP=${BIN}/fitswrap +++ P_FITSWRAP=/home/schubi/THELI/theli//bin/Linux_64//fitswrap P_FITSCUBESLICE=${BIN}/fitscubeslice +++ P_FITSCUBESLICE=/home/schubi/THELI/theli//bin/Linux_64//fitscubeslice P_FITSADDKEY=${BIN}/fitsaddkey +++ P_FITSADDKEY=/home/schubi/THELI/theli//bin/Linux_64//fitsaddkey P_FITSCUT=${BIN}/fitscut +++ P_FITSCUT=/home/schubi/THELI/theli//bin/Linux_64//fitscut P_FITSBIN=${BIN}/fitsbin +++ P_FITSBIN=/home/schubi/THELI/theli//bin/Linux_64//fitsbin P_FITSBLOOMDETECT=${BIN}/fitsbloomdetect +++ P_FITSBLOOMDETECT=/home/schubi/THELI/theli//bin/Linux_64//fitsbloomdetect P_FITSFILL=${BIN}/fitsfill +++ P_FITSFILL=/home/schubi/THELI/theli//bin/Linux_64//fitsfill P_FITSGRAD=${BIN}/fitsgrad +++ P_FITSGRAD=/home/schubi/THELI/theli//bin/Linux_64//fitsgrad P_FITSDEMOSAICBAYER=${BIN}/fitsdemosaicbayer +++ P_FITSDEMOSAICBAYER=/home/schubi/THELI/theli//bin/Linux_64//fitsdemosaicbayer P_FITS3MAX=${BIN}/fits3max +++ P_FITS3MAX=/home/schubi/THELI/theli//bin/Linux_64//fits3max P_FITSGAUSS=${BIN}/fitsgauss +++ P_FITSGAUSS=/home/schubi/THELI/theli//bin/Linux_64//fitsgauss P_FITS2TIFF=${BIN}/fits2tiff +++ P_FITS2TIFF=/home/schubi/THELI/theli//bin/Linux_64//fits2tiff P_TIFF2FITS=${BIN}/tiff2fits +++ P_TIFF2FITS=/home/schubi/THELI/theli//bin/Linux_64//tiff2fits P_CORRLIRIS=${BIN}/corrliris +++ P_CORRLIRIS=/home/schubi/THELI/theli//bin/Linux_64//corrliris P_CROSSTALK=${BIN}/crosstalk +++ P_CROSSTALK=/home/schubi/THELI/theli//bin/Linux_64//crosstalk P_FINDGSC=${BIN}/findgsc +++ P_FINDGSC=/home/schubi/THELI/theli//bin/Linux_64//findgsc P_FINDUSNOB1=${BIN}/findusnob1 +++ P_FINDUSNOB1=/home/schubi/THELI/theli//bin/Linux_64//findusnob1 P_FIND2MASS=${BIN}/find2mass +++ P_FIND2MASS=/home/schubi/THELI/theli//bin/Linux_64//find2mass P_FINDSDSS7=${BIN}/findsdss7 +++ P_FINDSDSS7=/home/schubi/THELI/theli//bin/Linux_64//findsdss7 P_FINDUCAC3=${BIN}/finducac3 +++ P_FINDUCAC3=/home/schubi/THELI/theli//bin/Linux_64//finducac3 P_FINDNOMAD1=${BIN}/findnomad1 +++ P_FINDNOMAD1=/home/schubi/THELI/theli//bin/Linux_64//findnomad1 P_VIZQUERY=${BIN}/vizquery +++ P_VIZQUERY=/home/schubi/THELI/theli//bin/Linux_64//vizquery P_ACLIENT=${BIN}/aclient +++ P_ACLIENT=/home/schubi/THELI/theli//bin/Linux_64//aclient P_SCAMP=${BIN}/scamp +++ P_SCAMP=/home/schubi/THELI/theli//bin/Linux_64//scamp P_RENAME=${BIN}/rename +++ P_RENAME=/home/schubi/THELI/theli//bin/Linux_64//rename P_MISSFITS=${BIN}/missfits +++ P_MISSFITS=/home/schubi/THELI/theli//bin/Linux_64//missfits P_SWARPFILTER=${BIN}/swarpfilter +++ P_SWARPFILTER=/home/schubi/THELI/theli//bin/Linux_64//swarpfilter P_FITSBLOCKEDIT=${BIN}/fitsblockedit +++ P_FITSBLOCKEDIT=/home/schubi/THELI/theli//bin/Linux_64//fitsblockedit P_FITSMASK=${BIN}/fitsmask +++ P_FITSMASK=/home/schubi/THELI/theli//bin/Linux_64//fitsmask P_FITSMEDIAN=${BIN}/fitsmedian +++ P_FITSMEDIAN=/home/schubi/THELI/theli//bin/Linux_64//fitsmedian P_FITSNANMASK=${BIN}/fitsnanmask +++ P_FITSNANMASK=/home/schubi/THELI/theli//bin/Linux_64//fitsnanmask P_FITSSPLITBAYER=${BIN}/fitssplitbayer +++ P_FITSSPLITBAYER=/home/schubi/THELI/theli//bin/Linux_64//fitssplitbayer P_FITSGUESSGEOM=${BIN}/fitsguessgeom +++ P_FITSGUESSGEOM=/home/schubi/THELI/theli//bin/Linux_64//fitsguessgeom P_SPLITVISIRFITSCUBE=${BIN}/split_visirfitscube +++ P_SPLITVISIRFITSCUBE=/home/schubi/THELI/theli//bin/Linux_64//split_visirfitscube P_SPLITSUPRIMECAM=${BIN}/splitsuprimecam +++ P_SPLITSUPRIMECAM=/home/schubi/THELI/theli//bin/Linux_64//splitsuprimecam P_PASTEMOSAICII=${BIN}/paste_mosaicII +++ P_PASTEMOSAICII=/home/schubi/THELI/theli//bin/Linux_64//paste_mosaicII P_PASTESOI=${BIN}/paste_soi +++ P_PASTESOI=/home/schubi/THELI/theli//bin/Linux_64//paste_soi P_GETSHIFT=${BIN}/get_shift +++ P_GETSHIFT=/home/schubi/THELI/theli//bin/Linux_64//get_shift P_FITSNONLINEARITY=${BIN}/fitsnonlinearity +++ P_FITSNONLINEARITY=/home/schubi/THELI/theli//bin/Linux_64//fitsnonlinearity P_FITSPIC=${BIN}/fitspic +++ P_FITSPIC=/home/schubi/THELI/theli//bin/Linux_64//fitspic P_DCRAW=${BIN}/dcraw +++ P_DCRAW=/home/schubi/THELI/theli//bin/Linux_64//dcraw P_DECSEXCONV=${BIN}/decsexconv +++ P_DECSEXCONV=/home/schubi/THELI/theli//bin/Linux_64//decsexconv P_STIFF=${BIN}/stiff +++ P_STIFF=/home/schubi/THELI/theli//bin/Linux_64//stiff P_GIFSICLE=${BIN}/gifsicle +++ P_GIFSICLE=/home/schubi/THELI/theli//bin/Linux_64//gifsicle P_XY2SKY=${BIN}/xy2sky +++ P_XY2SKY=/home/schubi/THELI/theli//bin/Linux_64//xy2sky P_SKY2XY=${BIN}/sky2xy +++ P_SKY2XY=/home/schubi/THELI/theli//bin/Linux_64//sky2xy P_GETKEY=${BIN}/getkey +++ P_GETKEY=/home/schubi/THELI/theli//bin/Linux_64//getkey P_SKY2XY_HELPER=${BIN}/sky2xy_helper +++ P_SKY2XY_HELPER=/home/schubi/THELI/theli//bin/Linux_64//sky2xy_helper P_GET_ROTIMSIZE=${BIN}/get_rotimsize +++ P_GET_ROTIMSIZE=/home/schubi/THELI/theli//bin/Linux_64//get_rotimsize P_CLEANSTRING=${BIN}/cleanstring +++ P_CLEANSTRING=/home/schubi/THELI/theli//bin/Linux_64//cleanstring P_NUMCPU=${BIN}/numcpu +++ P_NUMCPU=/home/schubi/THELI/theli//bin/Linux_64//numcpu P_GETSTATS=${BIN}/getstats +++ P_GETSTATS=/home/schubi/THELI/theli//bin/Linux_64//getstats P_CDELT2CD=${BIN}/cdelt2cd +++ P_CDELT2CD=/home/schubi/THELI/theli//bin/Linux_64//cdelt2cd P_GETCDMATRIX=${BIN}/get_cdmatrix +++ P_GETCDMATRIX=/home/schubi/THELI/theli//bin/Linux_64//get_cdmatrix P_GETPIXSCALE=${BIN}/get_pixscale +++ P_GETPIXSCALE=/home/schubi/THELI/theli//bin/Linux_64//get_pixscale P_GETPOSANGLE=${BIN}/get_posangle +++ P_GETPOSANGLE=/home/schubi/THELI/theli//bin/Linux_64//get_posangle P_ROTATECDMATRIX=${BIN}/rotate_cdmatrix +++ P_ROTATECDMATRIX=/home/schubi/THELI/theli//bin/Linux_64//rotate_cdmatrix P_MJD=${BIN}/mjd +++ P_MJD=/home/schubi/THELI/theli//bin/Linux_64//mjd P_LST=${BIN}/lst +++ P_LST=/home/schubi/THELI/theli//bin/Linux_64//lst P_IMALYZER=${BIN}/imalyzer +++ P_IMALYZER=/home/schubi/THELI/theli//bin/Linux_64//imalyzer P_FITSSMOOTHEDGE=${BIN}/fitssmoothedge +++ P_FITSSMOOTHEDGE=/home/schubi/THELI/theli//bin/Linux_64//fitssmoothedge P_FITSEXTRACTEXT=${BIN}/fitsextractext +++ P_FITSEXTRACTEXT=/home/schubi/THELI/theli//bin/Linux_64//fitsextractext P_ERRTEST=${BIN}/errtest +++ P_ERRTEST=/home/schubi/THELI/theli//bin/Linux_64//errtest # # load the parameter files # . ~/.theli/param_set1.ini +++ . /home/schubi/.theli/param_set1.ini NPARA=7 ++++ NPARA=7 GUIVERSION=2.3.10 ++++ GUIVERSION=2.3.10 NFRAMES=14 ++++ NFRAMES=14 V_DO_BIAS=N ++++ V_DO_BIAS=N V_DO_FLAT=N ++++ V_DO_FLAT=N V_NONLINCORR=N ++++ V_NONLINCORR=N V_DO_SUPERFLAT=N ++++ V_DO_SUPERFLAT=N V_ILLUMFRINGEBACKSIZE=256 ++++ V_ILLUMFRINGEBACKSIZE=256 V_SMOOTHFRINGEBACKSIZE= ++++ V_SMOOTHFRINGEBACKSIZE= V_SF_SMOOTH=SMOOTHED ++++ V_SF_SMOOTH=SMOOTHED V_FRINGESCALE=Y ++++ V_FRINGESCALE=Y V_ADJUSTGAINS=N ++++ V_ADJUSTGAINS=N V_RADIUS=40 ++++ V_RADIUS=40 V_AP_REFRA=nora ++++ V_AP_REFRA=nora V_AP_REFDEC=nodec ++++ V_AP_REFDEC=nodec . ~/.theli/param_set2.ini +++ . /home/schubi/.theli/param_set2.ini V_PRE_RENAME_CHECKED=0 ++++ V_PRE_RENAME_CHECKED=0 V_RENAME_FITSKEY= ++++ V_RENAME_FITSKEY= V_SORT_FITSKEY= ++++ V_SORT_FITSKEY= V_SORT_BIAS= ++++ V_SORT_BIAS= V_SORT_DARK= ++++ V_SORT_DARK= V_SORT_DOMEFLAT= ++++ V_SORT_DOMEFLAT= V_SORT_SKYFLAT= ++++ V_SORT_SKYFLAT= V_SORT_STD= ++++ V_SORT_STD= V_PRE_XTALK_NOR_CHECKED=0 ++++ V_PRE_XTALK_NOR_CHECKED=0 V_PRE_XTALK_ROW_CHECKED=0 ++++ V_PRE_XTALK_ROW_CHECKED=0 V_PRE_XTALK_MULTI_CHECKED=0 ++++ V_PRE_XTALK_MULTI_CHECKED=0 V_PRE_XTALK_NOR_AMPLITUDE= ++++ V_PRE_XTALK_NOR_AMPLITUDE= V_PRE_XTALK_ROW_AMPLITUDE= ++++ V_PRE_XTALK_ROW_AMPLITUDE= V_PRE_XTALK_NOR_BUTTONID=0 ++++ V_PRE_XTALK_NOR_BUTTONID=0 V_PRE_XTALK_ROW_BUTTONID=0 ++++ V_PRE_XTALK_ROW_BUTTONID=0 V_PRE_XTALK_MULTI_BUTTONID=0 ++++ V_PRE_XTALK_MULTI_BUTTONID=0 V_PRE_XTALK_MULTI_NSECTION= ++++ V_PRE_XTALK_MULTI_NSECTION= V_CAL_OVSCANNLOW=0 ++++ V_CAL_OVSCANNLOW=0 V_CAL_OVSCANNHIGH=1 ++++ V_CAL_OVSCANNHIGH=1 V_CAL_BIASNLOW=0 ++++ V_CAL_BIASNLOW=0 V_CAL_BIASNHIGH=1 ++++ V_CAL_BIASNHIGH=1 V_CAL_DARKNLOW=0 ++++ V_CAL_DARKNLOW=0 V_CAL_DARKNHIGH=1 ++++ V_CAL_DARKNHIGH=1 V_CAL_FLATNLOW=0 ++++ V_CAL_FLATNLOW=0 V_CAL_FLATNHIGH=1 ++++ V_CAL_FLATNHIGH=1 V_CAL_SCIENCENLOW=0 ++++ V_CAL_SCIENCENLOW=0 V_CAL_SCIENCENHIGH=1 ++++ V_CAL_SCIENCENHIGH=1 V_CAL_DETECTTHRESH=1.3 ++++ V_CAL_DETECTTHRESH=1.3 V_CAL_DETECTMINAREA=5 ++++ V_CAL_DETECTMINAREA=5 V_CAL_SEXBACKSIZE=256 ++++ V_CAL_SEXBACKSIZE=256 V_CAL_WINDOWSIZE=0 ++++ V_CAL_WINDOWSIZE=0 V_CAL_WINDOWMETHOD=0 ++++ V_CAL_WINDOWMETHOD=0 V_CAL_COMBINEMETHOD=MEDIAN ++++ V_CAL_COMBINEMETHOD=MEDIAN V_CAL_EXCLUDECENTERIMAGE=TRUE ++++ V_CAL_EXCLUDECENTERIMAGE=TRUE V_CAL_SEXFILTER=N ++++ V_CAL_SEXFILTER=N V_CAL_USEDARK=FALSE ++++ V_CAL_USEDARK=FALSE V_COLLDETECTTHRESH=1.0 ++++ V_COLLDETECTTHRESH=1.0 V_COLLDETECTMINAREA=5 ++++ V_COLLDETECTMINAREA=5 V_COLLREJECTTHRESH=1.5 ++++ V_COLLREJECTTHRESH=1.5 V_COLLDIRECTION=x ++++ V_COLLDIRECTION=x V_COLLXMIN= ++++ V_COLLXMIN= V_COLLXMAX= ++++ V_COLLXMAX= V_COLLYMIN= ++++ V_COLLYMIN= V_COLLYMAX= ++++ V_COLLYMAX= V_COLLMASKACTION=0 ++++ V_COLLMASKACTION=0 V_COLLAUTOTHRESHOLD=0 ++++ V_COLLAUTOTHRESHOLD=0 V_WEIGHTBINMIN=-100 ++++ V_WEIGHTBINMIN=-100 V_WEIGHTBINMAX=500 ++++ V_WEIGHTBINMAX=500 V_WEIGHTBINSIZE=4 ++++ V_WEIGHTBINSIZE=4 V_WEIGHTLOWTHRESHOLD= ++++ V_WEIGHTLOWTHRESHOLD= V_WEIGHTHIGHTHRESHOLD= ++++ V_WEIGHTHIGHTHRESHOLD= V_GLOBWFLATMIN=0.5 ++++ V_GLOBWFLATMIN=0.5 V_GLOBWFLATMAX=1.6 ++++ V_GLOBWFLATMAX=1.6 V_GLOBWDARKMIN= ++++ V_GLOBWDARKMIN= V_GLOBWDARKMAX= ++++ V_GLOBWDARKMAX= V_GLOBWDARKDIR= ++++ V_GLOBWDARKDIR= V_COSMICSTHRESHOLD=0.1 ++++ V_COSMICSTHRESHOLD=0.1 V_COSMICDT=6 ++++ V_COSMICDT=6 V_COSMICDMIN=1 ++++ V_COSMICDMIN=1 V_DEFECT_KERNELSIZE= ++++ V_DEFECT_KERNELSIZE= V_DEFECT_ROWTOL= ++++ V_DEFECT_ROWTOL= V_DEFECT_COLTOL= ++++ V_DEFECT_COLTOL= V_DEFECT_CLUSTOL= ++++ V_DEFECT_CLUSTOL= V_DEFECT_KERNELSIZE_SF= ++++ V_DEFECT_KERNELSIZE_SF= V_DEFECT_ROWTOL_SF= ++++ V_DEFECT_ROWTOL_SF= V_DEFECT_COLTOL_SF= ++++ V_DEFECT_COLTOL_SF= V_DEFECT_CLUSTOL_SF= ++++ V_DEFECT_CLUSTOL_SF= V_BLOOMLOWLIMIT=20000 ++++ V_BLOOMLOWLIMIT=20000 V_BLOOMMINAREA=100 ++++ V_BLOOMMINAREA=100 V_BLOOMWIDTH=5000 ++++ V_BLOOMWIDTH=5000 V_WEIGHT_BINOUTLIER=FALSE ++++ V_WEIGHT_BINOUTLIER=FALSE V_GLOBW_UNIFORMWEIGHT=FALSE ++++ V_GLOBW_UNIFORMWEIGHT=FALSE V_MASKBLOOMSPIKE=0 ++++ V_MASKBLOOMSPIKE=0 V_AP_DETECTTHRESH=10.0 ++++ V_AP_DETECTTHRESH=10.0 V_AP_DETECTMINAREA=10 ++++ V_AP_DETECTMINAREA=10 V_DEBLENDMINCONT=0.0005 ++++ V_DEBLENDMINCONT=0.0005 V_AP_LOWNUM=1 ++++ V_AP_LOWNUM=1 V_SEXCATMINFWHM=1.5 ++++ V_SEXCATMINFWHM=1.5 V_SEXCATFLAG=0 ++++ V_SEXCATFLAG=0 V_SEXBACKLEVEL= ++++ V_SEXBACKLEVEL= V_AP_FILTER=N ++++ V_AP_FILTER=N V_SCAMP_MAXPOSANGLE=2.0 ++++ V_SCAMP_MAXPOSANGLE=2.0 V_SCAMP_MAXOFFSET=2.0 ++++ V_SCAMP_MAXOFFSET=2.0 V_SCAMP_MAXSCALE=1.05 ++++ V_SCAMP_MAXSCALE=1.05 V_SCAMP_SNLOW=5 ++++ V_SCAMP_SNLOW=5 V_SCAMP_SNHIGH=20 ++++ V_SCAMP_SNHIGH=20 V_SCAMP_FGROUPRADIUS= ++++ V_SCAMP_FGROUPRADIUS= V_SCAMP_CROSSIDRADIUS= ++++ V_SCAMP_CROSSIDRADIUS= V_SCAMP_ASTREFWEIGHT= ++++ V_SCAMP_ASTREFWEIGHT= V_SCAMP_ASTRINSTRUKEY= ++++ V_SCAMP_ASTRINSTRUKEY= V_SCAMP_PHOTINSTRUKEY= ++++ V_SCAMP_PHOTINSTRUKEY= V_SCAMP_POLYORDER=3 ++++ V_SCAMP_POLYORDER=3 V_SCAMP_POLYORDER2= ++++ V_SCAMP_POLYORDER2= V_SCAMP_DISTORTGROUPS= ++++ V_SCAMP_DISTORTGROUPS= V_SCAMP_DISTORTKEYS= ++++ V_SCAMP_DISTORTKEYS= V_SCAMP_STABILITY=INSTRUMENT ++++ V_SCAMP_STABILITY=INSTRUMENT V_SCAMP_MOSAICTYPE=UNCHANGED ++++ V_SCAMP_MOSAICTYPE=UNCHANGED V_SCAMP_FOCALPLANE=DEFAULT ++++ V_SCAMP_FOCALPLANE=DEFAULT V_SCAMP_MATCHFLIPPED=Y ++++ V_SCAMP_MATCHFLIPPED=Y V_ASTRX_KANGLE=0.0 ++++ V_ASTRX_KANGLE=0.0 V_ASTRX_PSORDOLOFF=2 ++++ V_ASTRX_PSORDOLOFF=2 V_ASTRX_PSORDOLON=3 ++++ V_ASTRX_PSORDOLON=3 V_ABSPHOT_FILTER=B ++++ V_ABSPHOT_FILTER=B V_ABSPHOT_COLORFIXED=0.1 ++++ V_ABSPHOT_COLORFIXED=0.1 V_ABSPHOT_EXTINCTION=-0.1 ++++ V_ABSPHOT_EXTINCTION=-0.1 V_ABSPHOT_ZPMIN=24 ++++ V_ABSPHOT_ZPMIN=24 V_ABSPHOT_ZPMAX=25 ++++ V_ABSPHOT_ZPMAX=25 V_ABSPHOT_STDMINMAG=0.0 ++++ V_ABSPHOT_STDMINMAG=0.0 V_ABSPHOT_MAXPHOTERR=0.05 ++++ V_ABSPHOT_MAXPHOTERR=0.05 V_ABSPHOT_STDCAT_DIRECT=SDSS ++++ V_ABSPHOT_STDCAT_DIRECT=SDSS V_ABSPHOT_STDCAT_INDIRECT=LANDOLT ++++ V_ABSPHOT_STDCAT_INDIRECT=LANDOLT V_ABSPHOT_STDFILTER=U ++++ V_ABSPHOT_STDFILTER=U V_ABSPHOT_STDCOLOR=UmB ++++ V_ABSPHOT_STDCOLOR=UmB V_ABSPHOT_STDFILTERM2=u ++++ V_ABSPHOT_STDFILTERM2=u V_ABSPHOT_FITTINGMETHODM2=CHIP ++++ V_ABSPHOT_FITTINGMETHODM2=CHIP V_COADD_REFRA= ++++ V_COADD_REFRA= V_COADD_REFDEC= ++++ V_COADD_REFDEC= V_COADD_IDENT= ++++ V_COADD_IDENT= V_COADD_SEEING= ++++ V_COADD_SEEING= V_COADD_RZP= ++++ V_COADD_RZP= V_COADD_PIXSCALE= ++++ V_COADD_PIXSCALE= V_COADD_SIZEX= ++++ V_COADD_SIZEX= V_COADD_SIZEY= ++++ V_COADD_SIZEY= V_COADD_SKYPOSANGLE= ++++ V_COADD_SKYPOSANGLE= V_COADD_CHIPS="" ++++ V_COADD_CHIPS= V_COADD_FILTERTHRESHOLD= ++++ V_COADD_FILTERTHRESHOLD= V_COADD_FILTERCLUSTERSIZE= ++++ V_COADD_FILTERCLUSTERSIZE= V_COADD_FILTERBORDERWIDTH= ++++ V_COADD_FILTERBORDERWIDTH= V_COADD_SMOOTHEDGE= ++++ V_COADD_SMOOTHEDGE= V_COADD_PROPMOTIONRA= ++++ V_COADD_PROPMOTIONRA= V_COADD_PROPMOTIONDEC= ++++ V_COADD_PROPMOTIONDEC= V_COADD_KERNEL=LANCZOS3 ++++ V_COADD_KERNEL=LANCZOS3 V_COADD_FILTER=All ++++ V_COADD_FILTER=All V_COADD_PROJECTION=TAN ++++ V_COADD_PROJECTION=TAN V_COADD_CELESTIALTYPE=EQUATORIAL ++++ V_COADD_CELESTIALTYPE=EQUATORIAL V_COADD_COMBINETYPE=WEIGHTED ++++ V_COADD_COMBINETYPE=WEIGHTED . ~/.theli/param_set3.ini +++ . /home/schubi/.theli/param_set3.ini export PYTHONPATH=/home/schubi/THELI/gui-2.4.1/packages//python/ +++ export PYTHONPATH=/home/schubi/THELI/gui-2.4.1/packages//python/ +++ PYTHONPATH=/home/schubi/THELI/gui-2.4.1/packages//python/ # Geographical latitude of the observatory OBSLAT=50 ++ OBSLAT=50 # Pixel scale in x and y [degrees] (CDELT1/2 keywords) PIXSCX=-0.000575 ++ PIXSCX=-0.000575 PIXSCY=0.000575 ++ PIXSCY=0.000575 # Pixel scale in arcsec PIXSCALE=2.07 ++ PIXSCALE=2.07 # Detector gain GAIN=1.0 ++ GAIN=1.0 # Overscan regions. The values are the x-range. # The min is given by OVSCANX1, the max by OVSCANX2. # We need separate regions for each chip. OVSCANX1=([1]=0) ++ OVSCANX1=([1]=0) OVSCANX2=([1]=0) ++ OVSCANX2=([1]=0) # Sections for the cutting. Select those pixels that receive light. # Give min values for x and y, and the extent in x and y. CUTX=([1]=1) ++ CUTX=([1]=1) CUTY=([1]=1) ++ CUTY=([1]=1) SIZEX=([1]=3040) ++ SIZEX=([1]=3040) SIZEY=([1]=2016) ++ SIZEY=([1]=2016) # Reference pixels for each chip. REFPIXX=([1]=1516) ++ REFPIXX=([1]=1516) REFPIXY=([1]=1004) ++ REFPIXY=([1]=1004) . bash.include + . bash.include function get_wcs { local k k=1 cd11="" cd12="" cd21="" cd22="" cpx1="" cpx2="" while [ "${k}" -le "${NCHIPS}" ] do if [ "${NCHIPS}" = "1" ]; then ext=0 else ext=${k} fi # test if the CD matrix is present testcd=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CD1_1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CDELT1 | ${P_GAWK} '(NR==2) {print $2}'` testcdelt2=`echo ${testcdelt} | ${P_GAWK} '{if ($1>=1) print "bad"; else print "good"}'` if [ "${testcd}"_A != "_A" ] && [ "${testcd}" != "KEY_N/A" ]; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CD1_1 CD1_2 CD2_1 CD2_2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3, $4, $5}'` else if [ "${testcdelt}"_A != "_A" ] && [ "${testcdelt}" != "KEY_N/A" ] && [ "${testcdelt2}" = "good"] ; then cdmatrix=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CDELT1 CDELT2 | \ ${P_GAWK} '($1!="FILE") {print $2, "0.0", "0.0", $3}'` else cdmatrix="${PIXSCX} 0.0 0.0 ${PIXSCY}" fi fi # test if the reference pixel is present testcrpix=`${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CRPIX1 | ${P_GAWK} '(NR==2) {print $2}'` if [ "${testcrpix}"_A != "_A" ] && [ "${testcrpix}" != "KEY_N/A" ]; then refpixel=`${P_DFITS} -x ${ext} ${FILE} | \ ${P_FITSORT} CRPIX1 CRPIX2 | \ ${P_GAWK} '($1!="FILE") {print $2, $3}'` else rpix1=${REFPIXX[${k}]} rpix2=${REFPIXY[${k}]} if [ "${rpix1}_A" = "_A" ]; then rpix1=1 fi if [ "${rpix2}_A" = "_A" ]; then rpix2=1 fi refpixel="${rpix1} ${rpix2}" fi c11=`echo ${cdmatrix} | ${P_GAWK} '{print $1}'` c12=`echo ${cdmatrix} | ${P_GAWK} '{print $2}'` c21=`echo ${cdmatrix} | ${P_GAWK} '{print $3}'` c22=`echo ${cdmatrix} | ${P_GAWK} '{print $4}'` cr1=`echo ${refpixel} | ${P_GAWK} '{print $1}'` cr2=`echo ${refpixel} | ${P_GAWK} '{print $2}'` if [ "${k}" = "1" ]; then cd11="${c11}" cd12="${c12}" cd21="${c21}" cd22="${c22}" cpx1="${cr1}" cpx2="${cr2}" else cd11="${cd11},${c11}" cd12="${cd12},${c12}" cd21="${cd21},${c21}" cd22="${cd22},${c22}" cpx1="${cpx1},${cr1}" cpx2="${cpx2},${cr2}" fi k=$(( $k + 1 )) done } function get_statminsection { local minx local miny k=1 statminsection="" minx=${SIZEX[1]} miny=${SIZEY[1]} while [ "${k}" -le "${NCHIPS}" ] do if [ ${SIZEX[k]} -lt ${minx} ]; then minx=${SIZEX[k]} fi if [ ${SIZEY[k]} -lt ${miny} ]; then miny=${SIZEY[k]} fi k=$((${k}+1)) done minx=$((${minx}-1)) miny=$((${miny}-1)) statminsection="2,${minx},2,${miny}" } function p_ls { ls -l | grep ^- | ${P_GAWK} '{print $NF}' } function get_key { ${P_DFITS} $1 | ${P_FITSORT} "$2" | ${P_GAWK} '(NR==2) {if ($2!="KEY_N/A") print $2}' } # rename the images to the value of a FITS key (if requested) if [ "${V_PRE_RENAME_CHECKED}" = "1" ]; then cd /$1/$2 mkdir ORIGINALS p_ls > ${TEMPDIR}/fitskey_$$ cat ${TEMPDIR}/fitskey_$$ |\ { while read FILE do newname=`get_key ${FILE} ${V_RENAME_FITSKEY}` if [ "${newname}" != "" ]; then mv ${FILE} ${newname} else echo ERROR: Could not find ${V_RENAME_FITSKEY} in the FITS header mv ${FILE} ORIGINALS/ fi done } fi + '[' 0 = 1 ']' cd /$1/$2 + cd //home/schubi/THELI/REDU/Frames/20111003_M1//Bias FILES=`p_ls` p_ls ++ p_ls ++ ls -l ++ grep '^-' ++ /usr/bin/gawk '{print $NF}' + FILES='Bias_NewDiverOpt_001B.fit Bias_NewDiverOpt_002B.fit Bias_NewDiverOpt_003B.fit Bias_NewDiverOpt_004B.fit Bias_NewDiverOpt_005B.fit Bias_NewDiverOpt_006B.fit Bias_NewDiverOpt_007B.fit Bias_NewDiverOpt_008B.fit Bias_NewDiverOpt_009B.fit Bias_NewDiverOpt_010B.fit Bias_NewDiverOpt_011B.fit Bias_NewDiverOpt_012B.fit Bias_NewDiverOpt_013B.fit Bias_NewDiverOpt_014B.fit Bias_NewDiverOpt_015B.fit Bias_NewDiverOpt_016B.fit Bias_NewDiverOpt_017B.fit Bias_NewDiverOpt_018B.fit Bias_NewDiverOpt_019B.fit Bias_NewDiverOpt_020B.fit' mkdir ORIGINALS + mkdir ORIGINALS for FILE in ${FILES} do SUFFIX=`echo $FILE | ${P_GAWK} 'BEGIN {FS="."} {print "."$NF}'` BASE=`basename $FILE ${SUFFIX}` mv ${FILE} ${BASE}.fits FILE=${BASE}.fits RA=`get_key ${FILE} "RA"` if [ "${RA}" = "" ]; then RA=0.001 fi DEC=`get_key ${FILE} "DEC"` if [ "${DEC}" = "" ]; then DEC=0.0 fi EXPTIME=`get_key ${FILE} "EXPTIME"` if [ "${EXPTIME}" = "" ]; then EXPTIME=`get_key ${FILE} "EXPOSURE"` fi if [ "${EXPTIME}" = "" ]; then EXPTIME=1.0 fi OBJECT=`get_key ${FILE} "OBJECT"` if [ "${OBJECT}" = "" ]; then OBJECT=UNKNOWN fi EQUINOX=`get_key ${FILE} "EQUINOX" | sed 's/[^0-9.]//g'` if [ "${EQUINOX}" = "" ]; then EQUINOX=2000.0 fi MJD=`get_key ${FILE} "MJD-OBS"` if [ "${MJD}" = "" ]; then MJD=52000.0 fi FILTER=`get_key ${FILE} "FILTER" | ${P_CLEANSTRING}` if [ "${FILTER}" = "" ]; then FILTER=RGB else FILTER=RGB_${FILTER} fi LST=`get_key ${FILE} "LST"` if [ "${LST}" = "" ]; then LST=0.0 fi AIRMASS=`get_key ${FILE} "AIRMASS"` if [ "${AIRMASS}" = "" ]; then AIRMASS=1.0 fi GABODSID=`${P_NIGHTID} -t 16:00:00 -d 31/12/1998 -m ${MJD} |\ ${P_GAWK} ' ($1 ~ /Days/) {print $6}' | ${P_GAWK} 'BEGIN{ FS="."} {print $1}'` # Preserve the CD matrix if present. # If this does not work, replace the 'get_wcs' function below # with the explicit code in 'bash.include', and then modify it # according to your needs. get_wcs ${P_FITSSPLIT_ECL} \ -CRPIX1 "${cpx1}" \ -CRPIX2 "${cpx2}" \ -CD11 "${cd11}" \ -CD12 "${cd12}" \ -CD21 "${cd21}" \ -CD22 "${cd22}" \ -CRVAL1 ${RA} -CRVAL2 ${DEC} -EXPTIME ${EXPTIME} -AIRMASS ${AIRMASS} \ -GABODSID ${GABODSID} -FILTER ${FILTER} -EQUINOX ${EQUINOX} -OBJECT ${OBJECT} \ -M11 "1" \ -M12 "0" \ -M21 "0" \ -M22 "1" \ ${FILE} k=1 while [ "${k}" -le "${NCHIPS}" ] do tmp=`get_key ${FILE} "DATE-OBS"` ${P_REPLACEKEY} ${BASE}_${k}.fits "DATE-OBS= ${tmp}" DUMMY1 k=$(( $k + 1 )) done k=1 while [ "${k}" -le "${NCHIPS}" ] do ${P_REPLACEKEY} ${BASE}_${k}.fits "BAYERPAT= 'RGGB'" DUMMY10 k=$(( $k + 1 )) done mv ${FILE} ORIGINALS/${BASE}${SUFFIX} done + for FILE in '${FILES}' echo $FILE | ${P_GAWK} 'BEGIN {FS="."} {print "."$NF}' ++ echo Bias_NewDiverOpt_001B.fit ++ /usr/bin/gawk 'BEGIN {FS="."} {print "."$NF}' + SUFFIX=.fit basename $FILE ${SUFFIX} ++ basename Bias_NewDiverOpt_001B.fit .fit + BASE=Bias_NewDiverOpt_001B + mv Bias_NewDiverOpt_001B.fit Bias_NewDiverOpt_001B.fits + FILE=Bias_NewDiverOpt_001B.fits get_key ${FILE} "RA" ++ get_key Bias_NewDiverOpt_001B.fits RA ++ /home/schubi/THELI/theli//bin/Linux_64//dfits Bias_NewDiverOpt_001B.fits ++ /home/schubi/THELI/theli//bin/Linux_64//fitsort RA ++ /usr/bin/gawk '(NR==2) {if ($2!="KEY_N/A") print $2}' + RA= + '[' '' = '' ']' + RA=0.001 get_key ${FILE} "DEC" ++ get_key Bias_NewDiverOpt_001B.fits DEC ++ /home/schubi/THELI/theli//bin/Linux_64//dfits Bias_NewDiverOpt_001B.fits ++ /home/schubi/THELI/theli//bin/Linux_64//fitsort DEC ++ /usr/bin/gawk '(NR==2) {if ($2!="KEY_N/A") print $2}' + DEC= + '[' '' = '' ']' + DEC=0.0 get_key ${FILE} "EXPTIME" ++ get_key Bias_NewDiverOpt_001B.fits EXPTIME ++ /home/schubi/THELI/theli//bin/Linux_64//dfits Bias_NewDiverOpt_001B.fits ++ /home/schubi/THELI/theli//bin/Linux_64//fitsort EXPTIME ++ /usr/bin/gawk '(NR==2) {if ($2!="KEY_N/A") print $2}' + EXPTIME=1.000000000000E-003 + '[' 1.000000000000E-003 = '' ']' + '[' 1.000000000000E-003 = '' ']' get_key ${FILE} "OBJECT" ++ get_key Bias_NewDiverOpt_001B.fits OBJECT ++ /home/schubi/THELI/theli//bin/Linux_64//dfits Bias_NewDiverOpt_001B.fits ++ /home/schubi/THELI/theli//bin/Linux_64//fitsort OBJECT ++ /usr/bin/gawk '(NR==2) {if ($2!="KEY_N/A") print $2}' + OBJECT=Bias_ab + '[' Bias_ab = '' ']' get_key ${FILE} "EQUINOX" | sed 's/[^0-9.]//g' ++ get_key Bias_NewDiverOpt_001B.fits EQUINOX ++ sed 's/[^0-9.]//g' ++ /home/schubi/THELI/theli//bin/Linux_64//dfits Bias_NewDiverOpt_001B.fits ++ /home/schubi/THELI/theli//bin/Linux_64//fitsort EQUINOX ++ /usr/bin/gawk '(NR==2) {if ($2!="KEY_N/A") print $2}' + EQUINOX= + '[' '' = '' ']' + EQUINOX=2000.0 get_key ${FILE} "MJD-OBS" ++ get_key Bias_NewDiverOpt_001B.fits MJD-OBS ++ /home/schubi/THELI/theli//bin/Linux_64//dfits Bias_NewDiverOpt_001B.fits ++ /home/schubi/THELI/theli//bin/Linux_64//fitsort MJD-OBS ++ /usr/bin/gawk '(NR==2) {if ($2!="KEY_N/A") print $2}' + MJD= + '[' '' = '' ']' + MJD=52000.0 get_key ${FILE} "FILTER" | ${P_CLEANSTRING} ++ get_key Bias_NewDiverOpt_001B.fits FILTER ++ /home/schubi/THELI/theli//bin/Linux_64//cleanstring ++ /home/schubi/THELI/theli//bin/Linux_64//dfits Bias_NewDiverOpt_001B.fits ++ /home/schubi/THELI/theli//bin/Linux_64//fitsort FILTER ++ /usr/bin/gawk '(NR==2) {if ($2!="KEY_N/A") print $2}' + FILTER= + '[' '' = '' ']' + FILTER=RGB get_key ${FILE} "LST" ++ get_key Bias_NewDiverOpt_001B.fits LST ++ /home/schubi/THELI/theli//bin/Linux_64//dfits Bias_NewDiverOpt_001B.fits ++ /home/schubi/THELI/theli//bin/Linux_64//fitsort LST ++ /usr/bin/gawk '(NR==2) {if ($2!="KEY_N/A") print $2}' + LST= + '[' '' = '' ']' + LST=0.0 get_key ${FILE} "AIRMASS" ++ get_key Bias_NewDiverOpt_001B.fits AIRMASS ++ /home/schubi/THELI/theli//bin/Linux_64//dfits Bias_NewDiverOpt_001B.fits ++ /home/schubi/THELI/theli//bin/Linux_64//fitsort AIRMASS ++ /usr/bin/gawk '(NR==2) {if ($2!="KEY_N/A") print $2}' + AIRMASS= + '[' '' = '' ']' + AIRMASS=1.0 ${P_NIGHTID} -t 16:00:00 -d 31/12/1998 -m ${MJD} | ${P_GAWK} ' ($1 ~ /Days/) {print $6}' | ${P_GAWK} 'BEGIN{ FS="."} {print $1}' ++ /home/schubi/THELI/theli//bin/Linux_64//nightid -t 16:00:00 -d 31/12/1998 -m 52000.0 ++ /usr/bin/gawk ' ($1 ~ /Days/) {print $6}' ++ /usr/bin/gawk 'BEGIN{ FS="."} {print $1}' + GABODSID=821 + get_wcs + local k + k=1 + cd11= + cd12= + cd21= + cd22= + cpx1= + cpx2= + '[' 1 -le 1 ']' + '[' 1 = 1 ']' + ext=0 ${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CD1_1 | ${P_GAWK} '(NR==2) {print $2}' ++ /home/schubi/THELI/theli//bin/Linux_64//dfits -x 0 Bias_NewDiverOpt_001B.fits ++ /home/schubi/THELI/theli//bin/Linux_64//fitsort CD1_1 ++ /usr/bin/gawk '(NR==2) {print $2}' + testcd=KEY_N/A ${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CDELT1 | ${P_GAWK} '(NR==2) {print $2}' ++ /home/schubi/THELI/theli//bin/Linux_64//dfits -x 0 Bias_NewDiverOpt_001B.fits ++ /home/schubi/THELI/theli//bin/Linux_64//fitsort CDELT1 ++ /usr/bin/gawk '(NR==2) {print $2}' + testcdelt=KEY_N/A echo ${testcdelt} | ${P_GAWK} '{if ($1>=1) print "bad"; else print "good"}' ++ echo KEY_N/A ++ /usr/bin/gawk '{if ($1>=1) print "bad"; else print "good"}' + testcdelt2=bad + '[' KEY_N/A_A '!=' _A ']' + '[' KEY_N/A '!=' KEY_N/A ']' + '[' KEY_N/A_A '!=' _A ']' + '[' KEY_N/A '!=' KEY_N/A ']' + cdmatrix='-0.000575 0.0 0.0 0.000575' ${P_DFITS} -x ${ext} ${FILE} | ${P_FITSORT} CRPIX1 | ${P_GAWK} '(NR==2) {print $2}' ++ /home/schubi/THELI/theli//bin/Linux_64//dfits -x 0 Bias_NewDiverOpt_001B.fits ++ /home/schubi/THELI/theli//bin/Linux_64//fitsort CRPIX1 ++ /usr/bin/gawk '(NR==2) {print $2}' + testcrpix=KEY_N/A + '[' KEY_N/A_A '!=' _A ']' + '[' KEY_N/A '!=' KEY_N/A ']' + rpix1=1516 + rpix2=1004 + '[' 1516_A = _A ']' + '[' 1004_A = _A ']' + refpixel='1516 1004' echo ${cdmatrix} | ${P_GAWK} '{print $1}' ++ echo -0.000575 0.0 0.0 0.000575 ++ /usr/bin/gawk '{print $1}' + c11=-0.000575 echo ${cdmatrix} | ${P_GAWK} '{print $2}' ++ echo -0.000575 0.0 0.0 0.000575 ++ /usr/bin/gawk '{print $2}' + c12=0.0 echo ${cdmatrix} | ${P_GAWK} '{print $3}' ++ echo -0.000575 0.0 0.0 0.000575 ++ /usr/bin/gawk '{print $3}' + c21=0.0 echo ${cdmatrix} | ${P_GAWK} '{print $4}' ++ echo -0.000575 0.0 0.0 0.000575 ++ /usr/bin/gawk '{print $4}' + c22=0.000575 echo ${refpixel} | ${P_GAWK} '{print $1}' ++ echo 1516 1004 ++ /usr/bin/gawk '{print $1}' + cr1=1516 echo ${refpixel} | ${P_GAWK} '{print $2}' ++ echo 1516 1004 ++ /usr/bin/gawk '{print $2}' + cr2=1004 + '[' 1 = 1 ']' + cd11=-0.000575 + cd12=0.0 + cd21=0.0 + cd22=0.000575 + cpx1=1516 + cpx2=1004 + k=2 + '[' 2 -le 1 ']' + /home/schubi/THELI/theli//bin/Linux_64//mefsplit -CRPIX1 1516 -CRPIX2 1004 -CD11 -0.000575 -CD12 0.0 -CD21 0.0 -CD22 0.000575 -CRVAL1 0.001 -CRVAL2 0.0 -EXPTIME 1.000000000000E-003 -AIRMASS 1.0 -GABODSID 821 -FILTER RGB -EQUINOX 2000.0 -OBJECT Bias_ab -M11 1 -M12 0 -M21 0 -M22 1 Bias_NewDiverOpt_001B.fits > WARNING: default.split not found, using internal defaults reading main header [1merror: reading main header from [Bias_NewDiverOpt_001B.fits] [0m+ k=1 + '[' 1 -le 1 ']' get_key ${FILE} "DATE-OBS" ++ get_key Bias_NewDiverOpt_001B.fits DATE-OBS ++ /home/schubi/THELI/theli//bin/Linux_64//dfits Bias_NewDiverOpt_001B.fits ++ /home/schubi/THELI/theli//bin/Linux_64//fitsort DATE-OBS ++ /usr/bin/gawk '(NR==2) {if ($2!="KEY_N/A") print $2}' + tmp=2011-10-02T12:10:40 + /home/schubi/THELI/theli//bin/Linux_64//replacekey Bias_NewDiverOpt_001B_1.fits 'DATE-OBS= 2011-10-02T12:10:40' DUMMY1 cannot open Bias_NewDiverOpt_001B_1.fits: aborting error getting FITS header size for Bias_NewDiverOpt_001B_1.fits error during placeholder replacement: aborting + k=2 + '[' 2 -le 1 ']' + k=1 + '[' 1 -le 1 ']' + /home/schubi/THELI/theli//bin/Linux_64//replacekey Bias_NewDiverOpt_001B_1.fits 'BAYERPAT= '\''RGGB'\''' DUMMY10 cannot open Bias_NewDiverOpt_001B_1.fits: aborting error getting FITS header size for Bias_NewDiverOpt_001B_1.fits error during placeholder replacement: aborting + k=2 + '[' 2 -le 1 ']' + mv Bias_NewDiverOpt_001B.fits ORIGINALS/Bias_NewDiverOpt_001B.fit + for FILE in '${FILES}' echo $FILE | ${P_GAWK} 'BEGIN {FS="."} {print "."$NF}' ++ echo Bias_NewDiverOpt_002B.fit ++ /usr/bin/gawk 'BEGIN {FS="."} {print "."$NF}' + SUFFIX=.fit basename $FILE ${SUFFIX}

Bei einem Bias Bild bleibt wahrscheinlich nichts anderes übrig als die Belichtungszeit in 10-3 s anzugeben.
Kann ich am Header etwas ändern ?

Theli kann den Fitsheaader lesen. Skycat und Fitswork können das Bild anzeigen.

Momentan bin ich ratlos.

viele Grüsse
Christian

von **fehlfarbe** » Freitag 16. Dezember 2011, 23:07

Hallo Christian,

deine Chipgeometrie etc scheint ja zu stimmen. Bei mir gab es ja das Problem, dass die Belichtungszeit im FITS Header in ms angegeben war.

Durch anpassen der EXPTIME Zeile(n) im Splitting Script wird die Zeit durch 1000 geteilt und funktioniert dann auch wieder.

Code: Alles auswählen: EXPTIME=`get_key ${FILE} "EXPOSURE"` if [ "${EXPTIME}" = "" ]; then EXPTIME=`get_key ${FILE} "EXPTIME"` fi EXPTIME=`echo ${EXPTIME} | ${P_GAWK} '{print $1/1000.}'`

Also schau mal nach, ob es vllt am Header liegt.

Grüße,
Marcus

von **schubi** » Samstag 17. Dezember 2011, 00:33

Hallo Marcus,
ich habe Aufnahmen mit dem alten Treiber( ASCOM) und dem neuen Treiber (StarSenseSci) gemacht.
Die Dateien wurden danach mit Fitswork geöffnet und die Header verglichen.
Theli kann die mit dem neuen Treiber erstellten Dateien öffnen und den FitsHeader auch anzeigen.

ich habe keinen wesentlichen Unterschied gefunden.
Beispielsweise sind Belichtungszeiten in beiden Fällen in Sekunden angegeben.
Die Treibernamen sind im FitsHeader eingetragen und sind unterschiedlich voneinander.
Der egain Wert ist statt 1.0000 jetzt 0.5...... Das kann doch nicht die Ursache sein?

Es kann sein, daß eine Ausleseroutine mit dem Header zurechtkommt. Das braucht aber nicht für das Splitting-Skript zu gelten.

Bisher hatte ich im Zusammenspiel von Nutzdaten und Theli keine solchen Probleme.
Mit den Treiber für die ALCCD6cPro gab es schon mal Probleme unter Windows.

viele Grüsse

von **moos** » Samstag 17. Dezember 2011, 10:55

Hi Christian,

EZCAP wohl in Millisekunden

Das geht schief.
Beim Bias wäre das in einem noch möglichen Wertebereich. Bei den Lights aber sicher nicht.

Es gab dazu eine Mail von Mischa auf der Liste:

hi carsten,

> EXPOSURE= +3.600000000000e+003
> Theli meckert: EXPTIME keyword out of range
> Ist das für theli wirklich zu viel??

du hast einfach eine willkuerlich gewaehlte interne parameter-obergrenze
von 3000s ueberschritten.
das kannst du in

THELI/theli/eclipsetools/tools/preflist_split.h, zeile 73, hochsetzen.
danach den pipeline-kern neu kompilieren, und das gui ebenfalls.

hast du wirklich 3600s einzelbelichtungen?

mischa

Wenn Du willst, schreibe ich dir eine zeil für dein splitting script. Aber erst am Montag.
Dazu schickst Du mir am besten dein jetziges splitting script.
Brauche ich ja gar nicht, hat der Markus ja schon gemacht ;-)

Theli-Forum

ALCCD6 Fehler bei Preparation

Re: ALCCD6 Fehler bei Preparation

Re: ALCCD6 Fehler bei Preparation

Re: ALCCD6 Fehler bei Preparation

Re: ALCCD6 Fehler bei Preparation

Re: ALCCD6 Fehler bei Preparation

Re: ALCCD6 Fehler bei Preparation

Re: ALCCD6 Fehler bei Preparation

Re: ALCCD6 Fehler bei Preparation

Re: ALCCD6 Fehler bei Preparation

Re: ALCCD6 Fehler bei Preparation

Re: ALCCD6 Fehler bei Preparation

Re: ALCCD6 Fehler bei Preparation

Re: ALCCD6 Fehler bei Preparation

Re: ALCCD6 Fehler bei Preparation

Re: ALCCD6 Fehler bei Preparation

Wer ist online?