Changeset 490 in Sophya for trunk/SophyaLib/NTools/imageop.cc

Timestamp:

Oct 21, 1999, 5:25:53 PM (26 years ago)

Author:

ansari

Message:

Merge avec PEIDA++ (~V 3.8) et nettoyage pour nouveau PPersist Reza+cmv 21/10/99

File:

• trunk/SophyaLib/NTools/imageop.cc (modified) (1 diff)

trunk/SophyaLib/NTools/imageop.cc

@@ -528 +528 @@
 
 }
+
+/* Nouvelle-Fonction */
+//++
+ImageR4* NoiseFiltImage(ImageR4& img, ImageR4& filtre, DynCCD& dynccd)
+//      Calcule une image de bruit, a partir de "img" et de "dynccd"
+//      prenant en compte le filtre de convolution "filtre".
+//      L'image de bruit renvoyee contient les correlations de bruit entre
+//      pixels, du au filtrage (convolution) de l'image "img" par le 
+//      filtre "filtre".
+//--
+{
+  ImageR4* ImgBrt = NoiseImage(&img,&dynccd);
+
+  int i,j;  
+  for (i=0; i<ImgBrt->XSize(); i++)
+    for (j=0; j<ImgBrt->YSize(); j++)  (*ImgBrt)(i,j) =(*ImgBrt)(i,j) *(*ImgBrt)(i,j);
+
+  ImageR4 Filter2(filtre.XSize(), filtre.YSize()) ;
+  for (i=0; i<filtre.XSize(); i++)
+    for (j=0; j<filtre.YSize(); j++) (Filter2)(i,j) = filtre(i,j) * filtre(i,j);  
+ 
+  ImageR4* out2;
+  out2 = new ImageR4(ImgBrt->XSize(), ImgBrt->YSize());
+  FilterImage(ImgBrt, out2, &Filter2);   
+  for (i=0; i<out2->XSize(); i++)
+    for (j=0; j<out2->YSize(); j++)  (*out2)(i,j) =sqrt((*out2)(i,j));
+  delete ImgBrt;
+  return(out2);
+}
+
 
 //////////////////////////////////////////////////////////////////////////////