Changeset 112

Timestamp:

Oct 8, 2003, 12:00:41 AM (23 years ago)

Author:

eugene

Message:

added code to update DATASEC keyword
added / improved keyword entries & comments

File:

• trunk/Ohana/src/perl/src/fr.smooth (modified) (8 diffs)

trunk/Ohana/src/perl/src/fr.smooth

@@ -24 +24 @@
 $CCDKEYWORD = `gconfig CCDNUM-KEYWORD`; chop $CCDKEYWORD;
 $dbmode = `gconfig DETREND-DB-MODE`; chop $dbmode;
+
 @frstats = `cat $frstats`;
 for ($i = 0; $i < @frstats; $i++) {
-    if ($frstats[$i] =~ m|^\s*\#|) { splice @frstats, $i, 1; $i--; }
+    if ($frstats[$i] =~ m|^\s*\#|) { 
+        ($tmp, $entry, $value) = split (" ", $frstats[$i]);
+        if ($entry eq "REF_CHIP") { $RefChip = $ccds[$value]; }
+        splice @frstats, $i, 1; $i--; 
+    }
 }
 ($C0, $C1, $Fs, $Fr) = split (" ", $frstats[$ccd]);
+
+# find fringe-point file
+$frpts  = `detsearch -quiet -image $infile $ccds[$ccd] split -type frpts`; chop $frpts;
+if ($?) { &escape ("can't get fringe points for image"); }
+($name, $Fs, $Fr) = split (" ", `getfringe $infile $ccds[$ccd] $frpts $binning`);
 
 $mask = `detsearch -ccd $ccds[$ccd] -type mask -time $start -select -quiet`; chop $mask;
@@ -35 +45 @@
 # unbin and insert header keywords
 open (MANA, "|mana --norc");
+
+# convert BIASSEC string to x1, x2, y1, y2 vars
+print MANA "macro parsesec\n";
+print MANA "  \$f = \$1\n";
+print MANA "  getchr \$f : n\n";
+print MANA "  substr \$f 1 {\$n-1} x1\n";
+print MANA "  strlen \$f N\n";
+print MANA "  substr \$f \$n {\$N-\$n} f\n";
+  
+print MANA "  getchr \$f , n\n";
+print MANA "  substr \$f 1 {\$n-1} x2\n";
+print MANA "  strlen \$f N\n";
+print MANA "  substr \$f \$n {\$N-\$n} f\n";
+  
+print MANA "  getchr \$f : n\n";
+print MANA "  substr \$f 1 {\$n-1} y1\n";
+print MANA "  strlen \$f N\n";
+print MANA "  substr \$f \$n {\$N-\$n} f\n";
+  
+print MANA "  getchr \$f ] n\n";
+print MANA "  substr \$f 1 {\$n-1} y2\n";
+print MANA "  strlen \$f N\n";
+print MANA "  substr \$f \$n {\$N-\$n} f\n";
+print MANA "end\n";
+
 print MANA "list kernel\n";
 print MANA "  0.025\n";
@@ -47 +82 @@
 print MANA "end\n";
 print MANA "macro go\n";
-print MANA " rd a $infile\n";
-print MANA " rebin a A -$binning\n";
 print MANA " \$CCDKEYWORD = $CCDKEYWORD\n";
+# read in mask, rebin
 if ($dbmode eq "MEF") {
     print MANA " rd mask $mask -n $ccds[$ccd]\n";
@@ -55 +89 @@
     print MANA " rd mask $mask\n";
 }
+print MANA " keyword mask NAXIS1 nx2\n";
+print MANA " keyword mask NAXIS2 ny2\n";
+print MANA " rebin mask binmask $binning -norm\n";
+print MANA " clip binmask 0.9 0.0 1.1 0.0\n";
+
+print MANA " rd a $infile\n";
+print MANA " set a = a * binmask\n";
+# print MANA " set a = a * binmask + $Fs*not(binmask)\n";
+print MANA " rebin a A -$binning\n";
 print MANA " keyword A NAXIS1 nx1\n";
 print MANA " keyword A NAXIS2 ny1\n";
-print MANA " keyword mask NAXIS1 nx2\n";
-print MANA " keyword mask NAXIS2 ny2\n";
 print MANA " delete a\n";
-print MANA " extract A a 0 0 \$nx1 \$ny1 0 0 \$nx2 \$ny2\n";
-print MANA " kern a kernel \n";
-print MANA " set b = a*mask + $Fs*not(mask)\n";
+
+print MANA " extract A b 0 0 \$nx1 \$ny1 0 0 \$nx2 \$ny2\n";
+print MANA " keyword b DATASEC line\n";
+print MANA " parsesec \$line\n";
+print MANA " sprint newline \"[%d:%d,%d:%d]\" {\$x1*2} {\$x2*2} {\$y1*2} {\$y2*2}\n";
+print MANA " keyword b DATASEC -w \$newline\n";
+
+@time  = localtime;
+$now   = sprintf "%04d-%02d-%02dT%02d:%02d:%02d", $time[5]+1900, $time[4]+1, $time[3], $time[2], $time[1], $time[0];
+$elsys = `gconfig ELIXIR_VERSION`; chop ($elsys);
+
+@word = split ("/", $frpts);
+$frpts = $word[-1];
 
 # add keywords for fringe amplitude, fit, etc
@@ -69 +120 @@
 print MANA " keyword b FRNG_SKY -wf $Fs\n";
 print MANA " keyword b FRNG_RNG -wf $Fr\n";
-print MANA " keyword b FRNG_REF -wd 11\n";
-print MANA " keyword b TVSTART  -w $start\n";
-print MANA " keyword b TVSTOP   -w $stop\n";
+print MANA " keyword b FRNG_REF -w  $RefChip\n";
+print MANA " keyword b FRNG_PTS -w  \"$frpts\"\n";
+print MANA " keyword b TVSTART  -w  $start\n";
+print MANA " keyword b TVSTOP   -w  $stop\n";
+print MANA " keyword b EL_SYS   -w  \"$elsys\"\n";
+print MANA " keyword b EL_SMTH  -w  \"$version\"\n";
+print MANA " keyword b REL_DATE -w  \"$now\"\n";
 
 # add comments for the new keywords
@@ -79 +134 @@
 print MANA " keyword b FRNG_RNG -wc \"Fringe Frame fringe strength\"\n";
 print MANA " keyword b FRNG_REF -wc \"Fringe Cross Correlation reference chip\"\n";
+print MANA " keyword b FRNG_PTS -wc \"Fringe Point File\"\n";
 print MANA " keyword b TVSTART  -wc \"data validity start time\"\n";
 print MANA " keyword b TVSTOP   -wc \"data validity stop time\"\n";
 print MANA " keyword b OBSTYPE  -w  \"FRINGE\"\n";
 print MANA " keyword b EXPTYPE  -w  \"FRINGE\"\n";
+print MANA " keyword b EL_SYS   -wc \"Elixir System Version\"\n";
+print MANA " keyword b EL_SMTH  -wc \"Elixir Smooth Version\"\n";
+print MANA " keyword b REL_DATE -wc \"UTC Release Date\"\n";
 
 print MANA " keyword b COMMENT  -ws \"Fringe (NN) = C0 + C1 * Fringe (REF)\"\n";
@@ -91 +150 @@
 print MANA " keyword b COMMENT  -ws \"Elixir smoothing: 3x3 kernel (3): 0.025, 0.100, 0.025\"\n";
 
-print MANA " wd b $outfile\n";
+print MANA " stat b - - - -\n";
+print MANA " \$range = \$MAX - \$MIN + 4\n";
+print MANA " if (\$range < 100)\n";
+print MANA "  \$range = 100\n";
+print MANA " end\n";
+print MANA " sprintf bscale %9.4f {\$range / 2^16}\n";
+print MANA " \$bzero = \$bscale * int ((\$MIN + \$range / 2) / \$bscale)\n";
+
+print MANA " wd b $outfile -bzero \$bzero -bscale \$bscale \n";
 print MANA " exit 0\n";
 
@@ -126 +193 @@
 }
 
+#print MANA " kern a kernel \n";
+#print MANA " set b = a*mask + $Fs*not(mask)\n";