Changeset 20594

Timestamp:

Nov 8, 2008, 2:29:21 PM (18 years ago)

Author:

eugene

Message:

fix the PSF image mosaics

File:

• trunk/psphot/src/psphotVisual.c (modified) (3 diffs)

trunk/psphot/src/psphotVisual.c

@@ -691 +691 @@
     int DY = 64;
 
-    psImage *psfMosaic = psImageAlloc (3*DX, 3*DY, PS_TYPE_F32);
+    psImage *psfMosaic = psImageAlloc (5*DX, 5*DY, PS_TYPE_F32);
     psImageInit (psfMosaic, 0.0);
     
-    psImage *funMosaic = psImageAlloc (3*DX, 3*DY, PS_TYPE_F32);
+    psImage *funMosaic = psImageAlloc (5*DX, 5*DY, PS_TYPE_F32);
     psImageInit (funMosaic, 0.0);
 
-    psImage *resMosaic = psImageAlloc (3*DX, 3*DY, PS_TYPE_F32);
+    psImage *resMosaic = psImageAlloc (5*DX, 5*DY, PS_TYPE_F32);
     psImageInit (resMosaic, 0.0);
 
@@ -703 +703 @@
 
     // generate a fake model at each of the 3x3 image grid positions
-    for (int x = -1; x <= +1; x ++) {
-        for (int y = -1; y <= +1; y ++) {
+    for (int x = -2; x <= +2; x ++) {
+        for (int y = -2; y <= +2; y ++) {
             // use the center of the center pixel of the image
-            float xc = (int)((0.5 + 0.45*x)*readout->image->numCols) + readout->image->col0 + 0.5;
-            float yc = (int)((0.5 + 0.45*y)*readout->image->numRows) + readout->image->row0 + 0.5;
+            float xc = (int)((0.5 + 0.225*x)*readout->image->numCols) + readout->image->col0 + 0.5;
+            float yc = (int)((0.5 + 0.225*y)*readout->image->numRows) + readout->image->row0 + 0.5;
 
             // assign the x and y coords to the image center
@@ -723 +723 @@
             // no need to mask the source here
             // XXX should we measure this for the analytical model only or the full model?
-            pmModelAddWithOffset (psfMosaic, NULL, model, PM_MODEL_OP_FULL | PM_MODEL_OP_CENTER, 0, -x*DX, -y*DY);
-            pmModelAddWithOffset (funMosaic, NULL, model, PM_MODEL_OP_FUNC | PM_MODEL_OP_CENTER, 0, -x*DX, -y*DY);
-            pmModelAddWithOffset (resMosaic, NULL, model, PM_MODEL_OP_RES0 | PM_MODEL_OP_RES1 | PM_MODEL_OP_CENTER, 0, -x*DX, -y*DY);
-
+            pmModelAddWithOffset (psfMosaic, NULL, model, PM_MODEL_OP_FULL | PM_MODEL_OP_CENTER, 0, x*DX, y*DY);
+            pmModelAddWithOffset (funMosaic, NULL, model, PM_MODEL_OP_FUNC | PM_MODEL_OP_CENTER, 0, x*DX, y*DY);
+            pmModelAddWithOffset (resMosaic, NULL, model, PM_MODEL_OP_RES0 | PM_MODEL_OP_RES1 | PM_MODEL_OP_CENTER, 0, x*DX, y*DY);
             psFree (model);
         }