home *** CD-ROM | disk | FTP | other *** search
/ ftp.disi.unige.it / 2015-02-11.ftp.disi.unige.it.tar / ftp.disi.unige.it / pub / .person / BarlaA / sw / OLD / Simo / SVM_mono / kernel.c < prev    next >
C/C++ Source or Header  |  2002-06-25  |  4KB  |  176 lines
  1. #include <math.h>
  2. #include <stdio.h>
  3. #include <string.h>
  4. #include "nrutil.h"
  5. #include "patt_rec.h"
  6.  
  7. extern int    dim;
  8. extern int    degree;
  9. extern double    sigma, normalization;
  10. extern int      dim_r,dim_c, delta_g, delta_r, delta_c;
  11.  
  12. double linear_kernel(double *v1, double *v2)
  13. {
  14.     register  i;
  15.     double aux = 0.;
  16.  
  17.     for (i = dim; i > 0; i--)
  18.         aux += *v1++ * *v2++;
  19. /*    fprintf(stderr, "%lf\n", aux/normalization);*/
  20.     return aux/normalization;
  21. }
  22.  
  23. double negative_kernel(double *v1, double *v2)
  24. {
  25.     register  i;
  26.     double aux = 0., tmp;
  27.  
  28.     
  29.     for (i = dim; i > 0; i--) {
  30.         tmp = *v1++ - *v2++;
  31.             aux -= tmp * tmp;
  32.       }
  33.     
  34.     /*
  35.     for (i=dim; i>0; i--) {
  36.       tmp = - (*v1)*(*v1) - (*v2)*(*v2) + 2*(*v1++)(*v2++);
  37.       aux += tmp;
  38.     }
  39.     */
  40.  
  41.  
  42.     return (aux/normalization);
  43. }
  44.  
  45. double poly_kernel(double *v1, double *v2)
  46. {
  47.     register  i;
  48.     double aux = 0.;
  49.     
  50.     for (i = dim; i > 0; i--)
  51.         aux += *v1++ * *v2++;
  52. /*    fprintf(stderr, "%lf\n", pow((aux/normalization + 1.),(double)degree)-1.);*/
  53.     return(pow((aux/normalization + 1.),(double)degree)-1.);
  54. }
  55.  
  56. double rbf_kernel(double *v1, double *v2)
  57. {
  58.     register  i;
  59.     double aux, res = 0.;
  60.     
  61.     for (i = dim; i > 0; i--)
  62.     {
  63.         aux = *v1++ - *v2++;
  64.         res += aux * aux;
  65.     }
  66. /*    fprintf(stderr, "%lf\n", exp(-res/(2*sigma*sigma)));*/
  67.     return exp(-res / (2*sigma*sigma));
  68. }
  69.  
  70. double mv_kernel(double *mod,double *imm)
  71. {
  72.   register i;
  73.   double normal_imm = 0.; 
  74.   double normal_mod = 0.; 
  75.   double mod2[dim],imm2[dim];    /*vettori ausiliarii*/
  76.   double inter = 0.;
  77.   
  78.   for(i=0; i<dim; i++){ 
  79.     //fprintf(stdout,"uffa %lf %lf \n",imm[i],mod[i]);
  80.     normal_imm += imm[i];
  81.     normal_mod += mod[i];
  82.   }
  83.   
  84.   //fprintf(stdout,"%lf %lf\n",normal_imm,normal_mod);
  85.   for(i=0; i<dim; i++){
  86.     imm2[i] = imm[i]/normal_imm;
  87.     mod2[i] = mod[i]/normal_mod;
  88.     inter +=  ilmiomin(imm2[i],mod2[i]);
  89.   }
  90.   
  91.   
  92.   return(inter); 
  93.  
  94.   
  95. } 
  96.  
  97. double haus_kernel(double *v1, double *v2)
  98. {
  99.   register int big_g;
  100.   int i,j,k,ii,jj,kk,h;
  101.   int max_val, counter, tmp_counter;
  102.   int st_k;
  103.   int check=0;
  104.  
  105.   int verbose, occlusion, max_noise;
  106.  
  107.   int type;
  108.   int w_r, w_c;
  109.   double *data_i, *model_data;
  110.   double *data_s, *base, *startp;
  111.   byte pop[ 256 ], *grey;
  112.   double **mat_1, **mat_2;
  113.   int flag=0;
  114.  
  115.   double  *image_s;
  116.   long *addr,*next;
  117.   int **surface, **row_surface, **col_surface, **and,*buffer;
  118.  
  119.   long tot=0;
  120.   
  121.   model_data= v1; //model first
  122.   data_i = v2;    //image second
  123.  
  124.   w_r= dim_r;
  125.   w_c = dim_c;
  126.   counter=0;
  127.  
  128.   /*PROVA DI CAMBIAMENTO CODICE*/
  129.  
  130.   for (i=delta_r; i<(dim_r-delta_r); i++)
  131.     for (j=delta_c; j<(dim_c-delta_c); j++){
  132.       flag=1;
  133.       for (k=-delta_r; k<= delta_r; k++)
  134.     for(h=-delta_c; h<= delta_c; h++)
  135.         if ((fabs(model_data[i*dim_c+j]-data_i[(i+k)*dim_c+j+h]))<=delta_g){
  136.         counter++;
  137.         flag=0; k=delta_r+1; h=delta_c+1;
  138.         }
  139.     }
  140.   return((double)counter/normalization);
  141. }
  142.  
  143.  
  144.  
  145. static double (*kernel_v[])() = {linear_kernel, negative_kernel,poly_kernel, rbf_kernel, mv_kernel, haus_kernel};
  146.  
  147. double kernel(int k, double *v1, double *v2)
  148. {
  149.     switch (k)
  150.     {
  151.         case LINEAR_KERNEL:
  152.             return linear_kernel(v1, v2);
  153.         case NEGATIVE_KERNEL:
  154.             return negative_kernel(v1, v2);
  155.         case POLY_KERNEL:
  156.             return poly_kernel(v1, v2);
  157.         case RBF_KERNEL:
  158.             return rbf_kernel(v1, v2);
  159.         case MV_KERNEL:
  160.             return mv_kernel(v1, v2);
  161.         case HAUS_KERNEL:
  162.           {
  163.             return (haus_kernel(v1, v2)+haus_kernel(v2, v1))/2;
  164.           }
  165.                 case HAUS_KERNEL_SIMM:
  166.                   {
  167.                     return (haus_kernel(v1, v2)+haus_kernel(v2, v1))/2;
  168.                    }
  169.                 fprintf(stderr, "Unknown kernel\n");
  170.             exit(-1);
  171.     }
  172.     return (*kernel_v[k])(v1, v2);
  173. }
  174.  
  175.  
  176.