Changeset 10181 for trunk/psModules/src/objects/models/pmModel_GAUSS.c

Timestamp:

Nov 24, 2006, 10:56:20 AM (19 years ago)

Author:

magnier

Message:

guess is off by sqrt(2)

File:

• trunk/psModules/src/objects/models/pmModel_GAUSS.c (modified) (3 diffs)

trunk/psModules/src/objects/models/pmModel_GAUSS.c

@@ -14 +14 @@
  * PM_PAR_XPOS 2  - X center of object
  * PM_PAR_YPOS 3  - Y center of object
- * PM_PAR_SXX 4   - X^2 term of elliptical contour (sqrt(2) / SigmaX)
- * PM_PAR_SYY 5   - Y^2 term of elliptical contour (sqrt(2) / SigmaY)
+ * PM_PAR_SXX 4   - X^2 term of elliptical contour (sqrt(2) * SigmaX)
+ * PM_PAR_SYY 5   - Y^2 term of elliptical contour (sqrt(2) * SigmaY)
  * PM_PAR_SXY 6   - X*Y term of elliptical contour
  *****************************************************************************/
@@ -95 +95 @@
 bool PM_MODEL_GUESS (pmModel *model, pmSource *source)
 {
-
     pmMoments *moments = source->moments;
     psF32     *PAR  = model->params->data.F32;
+
+    # if (0)
+
+        psEllipseMoments emoments;
+    emoments.x2 = moments->Sx;
+    emoments.y2 = moments->Sx;
+    emoments.xy = moments->Sxy;
+
+    psEllipseAxes axes = psEllipseMomentsToAxes (emoments);
+    psEllipseShape shape = psEllipseAxesToShape (axes);
+    # endif
 
     PAR[PM_PAR_SKY] = moments->Sky;
@@ -103 +113 @@
     PAR[PM_PAR_XPOS] = moments->x;
     PAR[PM_PAR_YPOS] = moments->y;
-    PAR[PM_PAR_SXX] = PS_MAX(0.5, moments->Sx);
-    PAR[PM_PAR_SYY] = PS_MAX(0.5, moments->Sy);
+    PAR[PM_PAR_SXX] = PS_MAX(0.5, M_SQRT2*moments->Sx);
+    PAR[PM_PAR_SYY] = PS_MAX(0.5, M_SQRT2*moments->Sy);
     PAR[PM_PAR_SXY] = 0.0;
+    # if (0)
+
+        PAR[PM_PAR_SXX] = PS_MAX(0.5, M_SQRT2*shape.sx);
+    PAR[PM_PAR_SYY] = PS_MAX(0.5, M_SQRT2*shape.sy);
+    PAR[PM_PAR_SXY] = shape.sxy;
+    # endif
+
     return(true);
 }