fsmooth.c

image processing including fourier,wavelet,segmentation etc.

/*----------------------Megawave2 Module----------------------*/
/*mwcommand
  name = {fsmooth};
  version = {"1.1"};
  author = {"Jacques Froment, Yann Guyonvarc'h"};
  function = {"Smoothing with a binomial filter"};
  usage = {
            'S':[S=2]->S  "standard deviation (default 2)",
            'W':[W=1]->W  "Weight of the considered pixel (default 1)",
            in->in        "Input Fimage",
            out<-out      "Output Fimage"
          };
*/
/*----------------------------------------------------------------------
 v1.1: bug nrow/ncol fixed (L.Moisan)
----------------------------------------------------------------------*/

#include <stdio.h>
#include <math.h>
#include "mw.h"

#define OUT(i,j) (ptrOUT[(j)*ncOUT + (i)])

void fsmooth(S,W,in,out)

int *S,*W;
Fimage in,out;

{

 register int i,j,ncOUT;
 register float *ptrOUT;
 int nrOUT,n,N;
 
 N=((2+(*W))*(*S)*(*S));

 out=mw_change_fimage(out,in->nrow,in->ncol);
 if (!out) mwerror(FATAL,1,"not enough memory\n"); 

 mw_copy_fimage(in,out);
 ptrOUT = out->gray;
 ncOUT = out->ncol;
 nrOUT = out->nrow;

 for (n=1;n<=N;n++)
   { 

     /*horizontal*/

     for (j=0;j<(nrOUT);j++)
       for (i=1;i<(ncOUT)-1;i++)
	 OUT(i,j) = (OUT((i-1),j)+OUT(i+1,j)+(*W)*OUT(i,j))/(2+(*W));

     for (j=0;j<(nrOUT);j++)
       OUT(ncOUT-1,j) = (OUT(((ncOUT)-2),j)+OUT((ncOUT)-1,j)+(*W)*OUT(((ncOUT)-1),j))/(2+(*W));

     for (j=0;j<(nrOUT);j++)
       OUT(0,j)=(OUT(0,j)+OUT(1,j)+(*W)*OUT(0,j))/(2+(*W));
	

     /*Vertical*/

     for (i=0;i<(ncOUT);i++)
       for (j=1;j<(nrOUT)-1;j++)
	 OUT(i,j) = (OUT(i,j-1)+OUT(i,j+1)+(*W)*OUT(i,j))/(2+(*W));
     
  
     for (i=0;i<(ncOUT);i++)
        OUT(i,nrOUT-1)=(OUT(i,((nrOUT)-2))+OUT(i,(nrOUT)-1)+(*W)*OUT(i,(nrOUT)-1))/(2+(*W));
     
     for (i=0;i<(ncOUT);i++)
       OUT(i,0) = (OUT(i,0)+OUT(i,1)+(*W)*OUT(i,0))/(2+(*W));
 
   }
}