Changeset 1531 in Sophya for trunk/ArchTOIPipe

Timestamp:

Jun 15, 2001, 10:12:53 AM (24 years ago)

Author:

ansari

Message:

Correction filtre gaussien - Reza 15/6/2001

Location:

trunk/ArchTOIPipe

Files:

• ProcWSophya/simtoipr.cc (modified) (1 diff)
• TestPipes/simtst.cc (modified) (2 diffs)

trunk/ArchTOIPipe/ProcWSophya/simtoipr.cc

@@ -378 +378 @@
     case GaussFilter :
       for(kk=-(wsz/2); kk<=(wsz/2); kk++) 
-        coef[kk+(wsz/2)] = a*exp(-(double)(kk*kk)/(s*s));
+        coef[kk+(wsz/2)] = a*exp(-(double)(kk*kk)/(2*s*s));
       break;
     case DiffFilter :

trunk/ArchTOIPipe/TestPipes/simtst.cc

@@ -295 +295 @@
 
     cout << " Creating a GaussianFilter SimpleFilter Object " << endl;
-    double G_sigma = 2.0;
+    double G_sigma = gfilt_sigma;
     double G_a = 1./(G_sigma*sqrt(M_PI*2.));
-    SimpleFilter filt(16, SimpleFilter::GaussFilter, G_a, G_sigma);
+    SimpleFilter filt(gfilt_wsz, SimpleFilter::GaussFilter, G_a, G_sigma);
 
     filt.addInput("in", toideglioff);
@@ -316 +316 @@
         if ( (kfj < 0) || (kfj >= fcoef.Size() ) )  continue;
         xfreq = (cfreq-kfj)/sigf_killfreq;
-        fcoef(kfj) -= exp(-xfreq*xfreq);
+        fcoef(kfj) -= exp(-xfreq*xfreq/2.0);
       }
     }