Changeset 3426

Timestamp:

Mar 15, 2005, 1:18:29 PM (21 years ago)

Author:

Paul Price

Message:

Put image combination, subtraction in main part of document instead of appendix (whoops!).

File:

• trunk/doc/modules/ModulesSDRS.tex (modified) (3 diffs)

trunk/doc/modules/ModulesSDRS.tex

@@ -1 @@
-%%% $Id: ModulesSDRS.tex,v 1.33 2005-03-15 23:13:16 price Exp $
+%%% $Id: ModulesSDRS.tex,v 1.34 2005-03-15 23:18:29 price Exp $
 \documentclass[panstarrs]{panstarrs}
 
@@ -1535 +1535 @@
 \end{verbatim}
 
-\appendix
-
-\section{Pseudo-C PSPhot}
-\label{psphot}
-
-\begin{verbatim}
-# include <pslib.h>
-# include <psmodule.h>
-
-main () {
-
-  psMetadata *header;
-  psImage *image;
-
-  fd = psFitsOpen (argv[1]);
-  md = psFitsReadHeader (fd);
-  image = psFitsReadImage (fd, md);
-
-  stats = psImageStats (NULL, image);
-
-  RDNOISE  = psMetadataLookup (md, "RDNOISE");
-  GAIN     = psMetadataLookup (md, "GAIN");
-  INNER    = psMetadataLookup (config, "INNER_RADIUS");
-  OUTER    = psMetadataLookup (config, "OUTER_RADIUS");
-  SATURATE = psMetadataLookup (config, "SATURATE");
-  NSIGMA   = psMetadataLookup (config, "PSF_PEAK_THRESHOLD");
-  RADIUS   = psMetadataLookup (config, "PSF_MOMENTS_RADIUS");
-  XBORDER  = psMetadataLookup (config, "XBORDER");
-  YBORDER  = psMetadataLookup (config, "YBORDER");
-
-  keep = psRegionAlloc (XBORDER, image->nCol - YBORDER, 
-                        YBORDER, image->nRow - YBORDER);
-  
-  Sky = stats->median;
-  Sig = sqrt(Sky/GAIN + SQ(RDNOISE));
-
-  kernel = psKernelParts ();
-  smooth = psImageConvolve (NULL, image, kernel, PS_PARTS);
-
-  peaks = pmFindImagePeaks (smooth, NSIGMA*Sig + Sky);
-
-  peaks = pmCullImagePeaks (peaks, SATURATE, keep);
-  
-  sources = pmSourceLocalSky (image, peaks, INNER, OUTER);
-
-  sources = pmSourceMoments (image, sources, RADIUS);
-  
-  sources = pmSourceRoughClassify (sources, SATURATE, MIN_SN_LIM, keep);
-
-  stars = pmSourceSelectBrightStars (sources);
-
-  stars = pmSourceFitModel (
-}
-
-
-
-psArray *pmFindImagePeaks (psImage *image, float threshold) {
-
-  psVector *row;
-
-  row = psVectorAlloc (image[0].Ncol);
-
-  /* find peaks in each row */
-  for (i = 0; i < image[0].Nrow; i++) {
-    rowpeaks = pmFindVectorPeaks (row, threshold);
-    peaks.x = rowpeaks;
-    peaks.y = i;
-    peaks.z = image (i, x);
-  }
-
-  /* drop non-local peaks (peaks with neighbors) */
-  for (n = 0; n < peaks.n; n++) {
-    if (!local_peak) {
-      drop_peak;
-    }
-  }
-  return (peaks);
-}
-\end{verbatim}
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -2184 +2105 @@
     psFree(kernels);
     psFree(kernelCoeffs);
-\ebd{verbatim}
+\end{verbatim}
+
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+\appendix
+
+\section{Pseudo-C PSPhot}
+\label{psphot}
+
+\begin{verbatim}
+# include <pslib.h>
+# include <psmodule.h>
+
+main () {
+
+  psMetadata *header;
+  psImage *image;
+
+  fd = psFitsOpen (argv[1]);
+  md = psFitsReadHeader (fd);
+  image = psFitsReadImage (fd, md);
+
+  stats = psImageStats (NULL, image);
+
+  RDNOISE  = psMetadataLookup (md, "RDNOISE");
+  GAIN     = psMetadataLookup (md, "GAIN");
+  INNER    = psMetadataLookup (config, "INNER_RADIUS");
+  OUTER    = psMetadataLookup (config, "OUTER_RADIUS");
+  SATURATE = psMetadataLookup (config, "SATURATE");
+  NSIGMA   = psMetadataLookup (config, "PSF_PEAK_THRESHOLD");
+  RADIUS   = psMetadataLookup (config, "PSF_MOMENTS_RADIUS");
+  XBORDER  = psMetadataLookup (config, "XBORDER");
+  YBORDER  = psMetadataLookup (config, "YBORDER");
+
+  keep = psRegionAlloc (XBORDER, image->nCol - YBORDER, 
+                        YBORDER, image->nRow - YBORDER);
+  
+  Sky = stats->median;
+  Sig = sqrt(Sky/GAIN + SQ(RDNOISE));
+
+  kernel = psKernelParts ();
+  smooth = psImageConvolve (NULL, image, kernel, PS_PARTS);
+
+  peaks = pmFindImagePeaks (smooth, NSIGMA*Sig + Sky);
+
+  peaks = pmCullImagePeaks (peaks, SATURATE, keep);
+  
+  sources = pmSourceLocalSky (image, peaks, INNER, OUTER);
+
+  sources = pmSourceMoments (image, sources, RADIUS);
+  
+  sources = pmSourceRoughClassify (sources, SATURATE, MIN_SN_LIM, keep);
+
+  stars = pmSourceSelectBrightStars (sources);
+
+  stars = pmSourceFitModel (
+}
+
+
+
+psArray *pmFindImagePeaks (psImage *image, float threshold) {
+
+  psVector *row;
+
+  row = psVectorAlloc (image[0].Ncol);
+
+  /* find peaks in each row */
+  for (i = 0; i < image[0].Nrow; i++) {
+    rowpeaks = pmFindVectorPeaks (row, threshold);
+    peaks.x = rowpeaks;
+    peaks.y = i;
+    peaks.z = image (i, x);
+  }
+
+  /* drop non-local peaks (peaks with neighbors) */
+  for (n = 0; n < peaks.n; n++) {
+    if (!local_peak) {
+      drop_peak;
+    }
+  }
+  return (peaks);
+}
+\end{verbatim}
 
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%