Changeset 36855 for trunk/ppStack/src/ppStackPrepare.c

Timestamp:

Jun 12, 2014, 11:52:06 PM (12 years ago)

Author:

mhuber

Message:

adding BSCALEOFFSET option hack to apply -0.5*bscale offset to each input warp, default is -0.0 offset

File:

• trunk/ppStack/src/ppStackPrepare.c (modified) (2 diffs)

trunk/ppStack/src/ppStackPrepare.c

@@ -148 +148 @@
     }
 
+    //MEH -- bscale offset hack for warps 
+    options->bscaleOffset = psVectorAlloc(options->num, PS_TYPE_F32);
+    psVectorInit(options->bscaleOffset, 0.0);
+    bool setbscaleoffset = psMetadataLookupBool(NULL,recipe, "BSCALEOFFSET");
+
     psArray *psfs = psArrayAlloc(num); // PSFs for PSF envelope
     int numCols = 0, numRows = 0;   // Size of image
@@ -162 +167 @@
         }
         options->sumExposure += options->exposures->data.F32[i];
+
+        //MEH -- hdu get redefn again below if conv (lots of this happening..) 
+        if (setbscaleoffset) {
+            pmHDU *hdu = pmHDUFromCell(cell);
+            assert(hdu && hdu->header);
+            float bscale = psMetadataLookupF32(NULL, hdu->header, "BSCALE");
+            options->bscaleOffset->data.F32[i] = 0.5*bscale;
+        }
+        psLogMsg("ppStack", PS_LOG_INFO, "bscaleOffset: %d %f", i,options->bscaleOffset->data.F32[i]);
 
         // Get list of PSFs, to determine target PSF