#include <stdlib.h>
#include <stdio.h>
#include <string.h>

#include "ctimer.h"
#include "imageop.h"

#include "pidemoup.h"

FitsImageR4 * PIUserProc_1(FitsImageR4 * myimg, int xp, int yp, 
                           char * arg1, char * arg2)
{
printf(" PIUserProc_1(%d %d %s %s)  Filtrage passe-bas 3*3 \n", xp, yp, arg1, arg2);
ImageR4 filt(3,3);
FitsImageR4 * out;

// Definition du filtre passe-bas
int i,j;
for(i=0; i<3; i++)
  for(j=0; j<3; j++)  filt(i,j) = 0.06;
filt(1,1) = 0.52;

filt.PrintPave();
out = new FitsImageR4(myimg->XSize(), myimg->YSize());
TIMEF ;
FilterImage((ImageR4 *) myimg, (ImageR4 *)out, &filt); 
return(out);
}
 
FitsImageR4 * PIUserProc_2(FitsImageR4 * myimg, int xp, int yp, 
                           char * arg1, char * arg2)
{
printf(" PIUserProc_2(%d %d %s %s)  Ombre 3*3 \n", xp, yp, arg1, arg2);
ImageR4 filt(3,3);
FitsImageR4 * out;

// Definition du filtre ombre
int i,j;
filt(0,0) = filt(1,1) = filt(2,2) = 0.;
filt(0,1) = filt(1,2) = -1.;
filt(1,0) = filt(2,1) = 1.;
filt(0,2) = -2.;
filt(2,0) = 2.;
filt.PrintPave();
out = new FitsImageR4(myimg->XSize(), myimg->YSize());
TIMEF ;
FilterImage((ImageR4 *) myimg, (ImageR4 *)out, &filt); 
return(out);
}
 
FitsImageR4 * PIUserProc_3(FitsImageR4 * /*myimg*/, int xp, int yp, 
                           char * arg1, char * arg2)
{
printf(" PIUserProc_3(%d %d %s %s) ne fait rien \n", xp, yp, arg1, arg2);
return(NULL);
}