Changeset 28676 for branches/eam_branches/ipp-20100621/psModules/src/objects/pmSourceFitModel.c

Timestamp:

Jul 15, 2010, 6:30:19 AM (16 years ago)

Author:

eugene

Message:

add pmSourceCopyData vs pmSourceCopy

File:

• branches/eam_branches/ipp-20100621/psModules/src/objects/pmSourceFitModel.c (modified) (2 diffs)

branches/eam_branches/ipp-20100621/psModules/src/objects/pmSourceFitModel.c

@@ -86 +86 @@
     psVector *yErr = psVectorAllocEmpty(nPix, PS_TYPE_F32);
 
+    // XXX for a test, skip the central pixel in the sersic fit
+    bool skipCenter = false && (model->type == pmModelClassGetType("PS_MODEL_SERSIC"));
+    float Xo = model->params->data.F32[PM_PAR_XPOS];
+    float Yo = model->params->data.F32[PM_PAR_YPOS];
+
     // fill in the coordinate and value entries
     nPix = 0;
@@ -109 +114 @@
             }
 
-            psVector *coord = psVectorAlloc(2, PS_TYPE_F32);
-
             // Convert i/j to image space:
             // 0.5 PIX: the coordinate values must be in pixel coords, not index            
-            coord->data.F32[0] = (psF32) (j + 0.5 + source->pixels->col0);
-            coord->data.F32[1] = (psF32) (i + 0.5 + source->pixels->row0);
+            float Xv = (psF32) (j + 0.5 + source->pixels->col0);
+            float Yv = (psF32) (i + 0.5 + source->pixels->row0);
+
+            // XXX possible skip of center pixel:
+            if (skipCenter) {
+                float r = hypot(Xv - Xo, Yv - Yo);
+                if (r < 0.75) {
+                    continue;
+                }
+            }
+
+            psVector *coord = psVectorAlloc(2, PS_TYPE_F32);
+            coord->data.F32[0] = Xv;
+            coord->data.F32[1] = Yv;
             x->data[nPix] = (psPtr *) coord;
             y->data.F32[nPix] = source->pixels->data.F32[i][j];