16.6.2 MS-DOS-Stapelverarbeitungsdatei POLMAP.BAT

cls echo off echo =================================================== echo Die CCD-Daten werden aus Sicherheitsgründen ins Arbeitsverzeichnis echo kopiert und dargestellt. echo ===================================================

   copy   c:$\backslash$pcvista$\backslash$daten$\backslash$lkh198$_{-}$1.fts   c:$\backslash$pcvista$\backslash$lkh198$_{-}$1.fts 


copy   c:$\backslash$pcvista$\backslash$daten$\backslash$lkh198$_{-}$2.fts 		
      c:$\backslash$pcvista$\backslash$lkh198$_{-}$2.fts 


copy  c:$\backslash$pcvista$\backslash$daten$\backslash$lkh198$_{-}$3.fts 		
      c:$\backslash$pcvista$\backslash$lkh198$_{-}$3.fts 


copy  c:$\backslash$pcvista$\backslash$daten$\backslash$lkh198$_{-}$4.fts 		
      c:$\backslash$pcvista$\backslash$lkh198$_{-}$4.fts 


copy  c:$\backslash$pcvista$\backslash$daten$\backslash$darksig.fts 		
     c:$\backslash$pcvista$\backslash$darksig.fts 


copy  c:$\backslash$pcvista$\backslash$daten$\backslash$flatsig.fts 		
      c:$\backslash$pcvista$\backslash$flatsig.fts

rem ...Ändern Sie Pfade und Dateinamen bei Bedarf! rem Die Darstellung erfolgt mittels der Stapelverarbeitungsdatei view.bat rem oben links - lkh198$_{-}$1 (0$^\circ$), rem oben rechts - lkh198$_{-}$2 (45$^\circ$), rem mitte links - lkh198$_{-}$3 (90$^\circ$), rem mitte rechts - lkh198$_{-}$4 (135$^\circ$), rem unten links - darksig (enthält das beim Auslesen der CCD und das rem thermisch bedingte Dunkelsignal), rem unten rechts - flatsig.

   call  view

echo =================================================== echo Reduktion der durch nicht funktionsfähige Pixel und durch den echo Einfluß der kosmischen Strahlung hervorgerufene Fehler echo =================================================== rem Zunächst werden die fehlerhaften Ränder der CCD-Aufnahmen rem abgeschnitten.

    call  view   


box  1  cr=255  cc=164  nr=460  nc=280  show

rem ....... Wählen Sie bei Bedarf eine selbstgewählte Boxgröße!

    window  lkh198$_{-}$1 box=1 


window  lkh198$_{-}$2 		 box=1 


window  lkh198$_{-}$3 		 box=1 


window  lkh198$_{-}$4 		 box=1 


window  darksig 		  box=1 


window  flatsig 		 box=1 


cls 		 


call  view 		

rem Aufgrund der geringen Spanne in den Pixelwerten der Files dark- und rem flatsig ist diese Korrektur bei diesen Datenfiles leicht möglich. rem Dazu wird jeweils der Befehl sky (Max.- und Min.-Bestimmung) aus- rem geführt. Ein Vergleich mit den Nachbarwerten durch den Befehl rem cursor ermöglicht durch lineare Interpolation die Bestimmung der rem Werte, auf die die Intensitäten der betreffenden Pixel mit Hilfe rem des Befehls clip zurückgesetzt werden.

   sky  darksig

rem Der Maximumwert (375) unterscheidet sich nicht gravierend von den rem Werten der Nachbarpixel. Der Minimumwert könnte durch einen Fehler an rem der CCD hervorgerufen sein. Der Wert wird auf den durch Interpolation rem der vier Nachbarpixel ermittelten Wert von 339 angehoben.

   clip  darksig min=298 vmin=339 


sky  flatsig 		 		

rem Der Maximumwert (12257) unterscheidet sich erheblich von den Werten rem der Nachbarpixel. Die Interpolation mit den vier Nachbarpixeln er- rem gibt einen Wert von 8809. Auffällig ist die Zeile, die auch den rem Minimumwert beinhaltet. Ursache dafür liegt offensichtlich in der rem CCD. Die Werte werden auf 9200 gesetzt, was dem ungefähren Mittel- rem wert der Nachbarzeilen entspricht.

   clip  flatsig max=12256  vmax=8809  min=5000  vmin=9200

rem Die "`Pixelkorrektur"' in den 4 Bilddatenfiles ist wegen des erhöh- rem ten Aufwandes im vorliegenden Versuch nicht unbedingt durchzuführen.
echo =================================================== echo Reduktion des Dunkelsignals, das durch das Auslesen und die ther- echo mische Ladungsträgergeneration entstanden ist echo =================================================== rem ....... Führen Sie die Dunkelsignalreduktion durch! ................................... ................................... ................................... ................................... ...................................
echo =================================================== echo Flat-Field-Reduktion der vier Bilddatenfiles zum Ausgleich der echo unterschiedlichen Pixelempfindlichkeiten: File * mn(flatsig)/flatsig echo =================================================== rem ....... Bilden Sie den Mittelwert des Datenfiles flatsig und führen rem ....... Sie im Anschluß daran die flat-field-Reduktion durch! ................................... ................................... ................................... ................................... ...................................
echo =================================================== echo Kalibration und Zentrierung echo =================================================== rem Zum Ausgleich unterschiedlicher Beobachtungsbedingungen und zum Aus- rem schluß der interstellaren Polarisation werden die vier Bild-Files rem zueinander kalibriert. Zur Ermittelung der Wichtungsfaktoren werden rem die Intensitätswerte von Sternen in den vier verschiedenen Aufnahmen rem verglichen (stars-Befehl, Sterne dürfen nicht durch den Nebel be- rem einflußt werden und müssen in etwa in Nebelentfernung liegen). rem Die Anwendung der Prozedur stars benötigt den durch rem die Prozedur sky bereitgestellten Parameter sky. rem Der stars-Befehl liefert gleichzeitig die exakten Ortskoordi- rem naten dieser Sterne, die für die nachfolgende Zentrierung erforder- rem lich sind.

   box  2  cr=60   cc=230  nr=100  nc=80  show 


sky  lkh198$_{-}$1 		 box=2 


let  -e  a=sky 		 


sky  lkh198$_{-}$2 		  box=2 


let  -e  b=sky 		 


sky  lkh198$_{-}$3 		 box=2 


let  -e  c=sky 		 


sky  lkh198$_{-}$4 		 box=2 


let  -e  d=sky 		 


stars  lkh198$_{-}$1 		 box=2  sky=a  minsig=6 


stars  lkh198$_{-}$2 		 box=2  sky=b  minsig=6 


stars  lkh198$_{-}$3 		 box=2  sky=c  minsig=6 


stars  lkh198$_{-}$4 		 box=2  sky=d  minsig=6 

rem Aus den 4 Intensitätswerten des hellsten Sternes ergeben sich fol- rem gende Wichtungsfaktoren: w(45$^\circ$) = 8883/9234 = 0,96, rem w(90$^\circ$) = 8883/8715 = 1,02, w(135$^\circ$) = 8883/9780 = 0,91. rem Da die Multiplikation mit diesen Faktoren im Integer-Zahlenbereich rem sehr ungünstig ist, erfolgt die Kalibrierung erst nach Umschreibung rem ins ASCII-Datenformat. rem ....... Weisen Sie den Variablen w45, w90 und w135 die ermittelten rem ....... Kalibrationsfaktoren zu. ................................... ................................... ................................... rem Zentrierung: die nachfolgende Verschiebungen werden durchgeführt: rem lkh198$_{-}$1 : nr= 1 nc=-1 rem lkh198$_{-}$2 : nr= 1 nc= 0 rem lkh198$_{-}$3 : nr= 0 nc= 0 rem lkh198$_{-}$4 : nr= 2 nc=-1 rem ....... Zentrieren Sie die vier Bild-Files zueinander! rem ....... Beachten Sie dabei den Verlauf der Spalten und Zeilen. ................................... ................................... ...................................
echo =================================================== echo Reduktion des Himmelshintergrundsignals echo ================================================== rem Hierzu wird ein Feld (box 3) genutzt, daß relativ stern- und rem staubfrei ist.

   tv  lkh198$_{-}$1   


box  3  cr=60 		 cc=250  nr=100  nc=40 show 


sky  lkh198$_{-}$1 		 box=3 


sub  lkh198$_{-}$1 		 const=sky 


sky  lkh198$_{-}$2 		 box=3 


sub  lkh198$_{-}$2 		 const=sky 


sky  lkh198$_{-}$3 		 box=3 


sub  lkh198$_{-}$3 		 const=sky 


sky  lkh198$_{-}$4 		 box=3 


sub  lkh198$_{-}$4 		 const=sky

rem Aufgrund des Signalrauschens kann es passieren, daß nach Abzug des rem Hintergrundsignals einige Pixel negative Werte aufweisen. Dies wird rem durch ein generelles Anheben aller negativen Werte auf Null rem korrigiert.

     clip  lkh198$_{-}$1 min=0 


clip  lkh198$_{-}$2 		 min=0 


clip  lkh198$_{-}$3 		 min=0 


clip  lkh198$_{-}$4 		 min=0

echo =================================================== echo Ausschneiden der als Polarisationskarten auszuschneidenden Gebiete echo =================================================== rem Die Zeilen- und Spaltenanzahl der dazu nötigen quadratischen Boxen rem muß durch drei teilbar sein. rem ....... Schneiden Sie die entsprechenden Gebiete heraus! rem ....... Beachten Sie dabei, daß die Files des Gesamtbildes dabei rem ....... nicht überschrieben werden. rem ....... Veranschaulichen Sie sich zunächst die auszuschneidenden rem ....... Gebiete im Gesamtbild mittels der Befehle box (Option show) rem ....... und tv. ................................... ................................... ...................................

   copy  c:$\backslash$pcvista$\backslash$lkh198$_{-}$1.fts
               c:$\backslash$pcvista$\backslash$v376$_{-}$1.fts 


copy  c:$\backslash$pcvista$\backslash$lkh198$_{-}$2.fts
 		  c:$\backslash$pcvista$\backslash$v376$_{-}$2.fts 


copy  c:$\backslash$pcvista$\backslash$lkh198$_{-}$3.fts
 		  c:$\backslash$pcvista$\backslash$v376$_{-}$3.fts 


copy  c:$\backslash$pcvista$\backslash$lkh198$_{-}$4.fts
 		  c:$\backslash$pcvista$\backslash$v376$_{-}$4.fts

................................... ................................... ................................... ................................... ................................... ................................... ................................... ................................... ................................... ................................... ................................... ................................... ................................... ................................... ................................... ...................................
echo =================================================== echo Umschreiben der CCD-Daten in das ASCII-Format, Ausführung der echo Kalibration (w45, w90, w135) und Erstellung der Polarisationskarten. echo ===================================================

  readfits  lkh198$_{-}$1  lkh198$_{-}$2  lkh198$_{-}$3                  lkh198$_{-}$4   


readfits  v376$_{-}$1  v376$_{-}$2  v376$_{-}$3                  v376$_{-}$4 


kalibrat  lkh198$_{-}$2  lkh198$_{-}$3  lkh198$_{-}$4                  0.96  1.02  0.91 


kalibrat  v376$_{-}$2  v376$_{-}$3  v376$_{-}$4  0.96                   1.02  0.91 


linpol  lkh198$_{-}$1  lkh198$_{-}$2  lkh198$_{-}$3                   lkh198$_{-}$4  100  0 


linpol  v376$_{-}$1  v376$_{-}$2  v376$_{-}$3                  v376$_{-}$4   100  0