Changeset 36307 for branches/eam_branches/ipp-20130904/psLib/src/imageops/psImageConvolve2dCache.c

Timestamp:

Nov 23, 2013, 12:09:42 PM (13 years ago)

Author:

eugene

Message:

add gauss to 2d cache smooth; some test code in pmPCM_Minimize

File:

• branches/eam_branches/ipp-20130904/psLib/src/imageops/psImageConvolve2dCache.c (modified) (1 diff)

branches/eam_branches/ipp-20130904/psLib/src/imageops/psImageConvolve2dCache.c

@@ -125 +125 @@
 }
 
+// generate a 2D smoothing kernel for supplied sigma & kappa (PS1_V1 profile).  
+bool psImageSmooth2dCacheKernel_Gauss (psImageSmooth2dCacheData *smdata, float sigma) {
+
+    // check for NULL structure elements?
+    int Ns = (int)(smdata->Nsigma * sigma);
+    Ns = PS_MAX (3, PS_MIN (Ns, 10));
+    smdata->Ns = Ns;
+
+    int Ns2 = Ns * Ns;
+
+    // we are going to use a hard-wired set of radial points
+    smdata->radflux = psAlloc(sizeof(float)*NRAD_MAX);
+
+    float sum = 0.0;
+    for (int i = 0; i < NRAD_MAX; i++) {
+        if (radii2[i] > Ns2) {
+            smdata->radflux[i] = 0.0;
+            continue;
+        }
+        float z = 0.5 * radii2[i] / PS_SQR(sigma);
+        smdata->radflux[i] = exp(-z);
+        sum += radiiN[i] * smdata->radflux[i];
+    }
+    for (int i = 0; i < NRAD_MAX; i++) {
+        smdata->radflux[i] = smdata->radflux[i] / sum;
+    }
+
+    return true;
+}
+
 // we can use the same DATA structure on multiple images of the same size
 bool psImageSmooth2dCache_F32(psImage *image, psImageSmooth2dCacheData *smdata)