mask8even.c

一个有关eigenface的matlab 例程

 0.005025,  0.006634,  0.008651,  0.011142,  0.014175,  0.017812,  0.022107,  0.027101,  0.032816,  0.039249,  0.046367,  0.054103,  0.062356,  0.070986,  0.079819,  0.088649,  0.097249,  0.105374,  0.112777,  0.119219,  0.124483,  0.128385,  0.130785,  0.131594,  0.130785,  0.128385,  0.124483,  0.119219,  0.112777,  0.105374,  0.097249,  0.088649,  0.079819,  0.070986,  0.062356,  0.054103,  0.046367,  0.039249,  0.032816,  0.027101,  0.022107,  0.017812,  0.014175,  0.011142,  0.008651,  0.006634,  0.005025,
 0.007928,  0.010466,  0.013647,  0.017578,  0.022362,  0.028099,  0.034875,  0.042754,  0.051769,  0.061917,  0.073146,  0.085350,  0.098369,  0.111983,  0.125917,  0.139848,  0.153415,  0.166232,  0.177911,  0.188074,  0.196378,  0.202533,  0.206318,  0.207596,  0.206318,  0.202533,  0.196378,  0.188074,  0.177911,  0.166232,  0.153415,  0.139848,  0.125917,  0.111983,  0.098369,  0.085350,  0.073146,  0.061917,  0.051769,  0.042754,  0.034875,  0.028099,  0.022362,  0.017578,  0.013647,  0.010466,  0.007928,
 0.008694,  0.011477,  0.014966,  0.019275,  0.024522,  0.030813,  0.038243,  0.046883,  0.056770,  0.067898,  0.080211,  0.093594,  0.107871,  0.122800,  0.138080,  0.153356,  0.168233,  0.182288,  0.195095,  0.206240,  0.215346,  0.222095,  0.226246,  0.227647,  0.226246,  0.222095,  0.215346,  0.206240,  0.195095,  0.182288,  0.168233,  0.153356,  0.138080,  0.122800,  0.107871,  0.093594,  0.080211,  0.067898,  0.056770,  0.046883,  0.038243,  0.030813,  0.024522,  0.019275,  0.014966,  0.011477,  0.008694,
 0.007795,  0.010290,  0.013418,  0.017282,  0.021986,  0.027627,  0.034289,  0.042035,  0.050900,  0.060877,  0.071917,  0.083917,  0.096717,  0.110102,  0.123802,  0.137499,  0.150837,  0.163440,  0.174922,  0.184914,  0.193079,  0.199131,  0.202852,  0.204108,  0.202852,  0.199131,  0.193079,  0.184914,  0.174922,  0.163440,  0.150837,  0.137499,  0.123802,  0.110102,  0.096717,  0.083917,  0.071917,  0.060877,  0.050900,  0.042035,  0.034289,  0.027627,  0.021986,  0.017282,  0.013418,  0.010290,  0.007795,
 0.005894,  0.007782,  0.010147,  0.013069,  0.016626,  0.020891,  0.025929,  0.031787,  0.038490,  0.046035,  0.054384,  0.063458,  0.073138,  0.083260,  0.093619,  0.103977,  0.114064,  0.123593,  0.132276,  0.139832,  0.146007,  0.150583,  0.153397,  0.154347,  0.153397,  0.150583,  0.146007,  0.139832,  0.132276,  0.123593,  0.114064,  0.103977,  0.093619,  0.083260,  0.073138,  0.063458,  0.054384,  0.046035,  0.038490,  0.031787,  0.025929,  0.020891,  0.016626,  0.013069,  0.010147,  0.007782,  0.005894,
 0.003654,  0.004824,  0.006290,  0.008102,  0.010307,  0.012951,  0.016075,  0.019706,  0.023862,  0.028539,  0.033714,  0.039340,  0.045340,  0.051615,  0.058038,  0.064459,  0.070712,  0.076620,  0.082003,  0.086687,  0.090515,  0.093352,  0.095096,  0.095685,  0.095096,  0.093352,  0.090515,  0.086687,  0.082003,  0.076620,  0.070712,  0.064459,  0.058038,  0.051615,  0.045340,  0.039340,  0.033714,  0.028539,  0.023862,  0.019706,  0.016075,  0.012951,  0.010307,  0.008102,  0.006290,  0.004824,  0.003654,
 0.001595,  0.002106,  0.002746,  0.003536,  0.004499,  0.005653,  0.007016,  0.008601,  0.010415,  0.012457,  0.014716,  0.017171,  0.019790,  0.022529,  0.025332,  0.028135,  0.030864,  0.033443,  0.035792,  0.037837,  0.039508,  0.040746,  0.041508,  0.041765,  0.041508,  0.040746,  0.039508,  0.037837,  0.035792,  0.033443,  0.030864,  0.028135,  0.025332,  0.022529,  0.019790,  0.017171,  0.014716,  0.012457,  0.010415,  0.008601,  0.007016,  0.005653,  0.004499,  0.003536,  0.002746,  0.002106,  0.001595,
 0.000031,  0.000041,  0.000054,  0.000069,  0.000088,  0.000110,  0.000137,  0.000168,  0.000203,  0.000243,  0.000287,  0.000335,  0.000387,  0.000440,  0.000495,  0.000550,  0.000603,  0.000653,  0.000699,  0.000739,  0.000772,  0.000796,  0.000811,  0.000816,  0.000811,  0.000796,  0.000772,  0.000739,  0.000699,  0.000653,  0.000603,  0.000550,  0.000495,  0.000440,  0.000387,  0.000335,  0.000287,  0.000243,  0.000203,  0.000168,  0.000137,  0.000110,  0.000088,  0.000069,  0.000054,  0.000041,  0.000031,
-0.000929, -0.001226, -0.001599, -0.002059, -0.002620, -0.003292, -0.004085, -0.005008, -0.006065, -0.007253, -0.008569, -0.009998, -0.011523, -0.013118, -0.014751, -0.016382, -0.017972, -0.019473, -0.020841, -0.022032, -0.023005, -0.023726, -0.024169, -0.024319, -0.024169, -0.023726, -0.023005, -0.022032, -0.020841, -0.019473, -0.017972, -0.016382, -0.014751, -0.013118, -0.011523, -0.009998, -0.008569, -0.007253, -0.006065, -0.005008, -0.004085, -0.003292, -0.002620, -0.002059, -0.001599, -0.001226, -0.000929,
-0.001335, -0.001762, -0.002298, -0.002960, -0.003765, -0.004731, -0.005872, -0.007198, -0.008716, -0.010425, -0.012316, -0.014371, -0.016563, -0.018855, -0.021201, -0.023546, -0.025831, -0.027989, -0.029955, -0.031666, -0.033064, -0.034101, -0.034738, -0.034953, -0.034738, -0.034101, -0.033064, -0.031666, -0.029955, -0.027989, -0.025831, -0.023546, -0.021201, -0.018855, -0.016563, -0.014371, -0.012316, -0.010425, -0.008716, -0.007198, -0.005872, -0.004731, -0.003765, -0.002960, -0.002298, -0.001762, -0.001335,
-0.001329, -0.001755, -0.002288, -0.002947, -0.003749, -0.004711, -0.005847, -0.007168, -0.008679, -0.010381, -0.012263, -0.014310, -0.016492, -0.018775, -0.021111, -0.023446, -0.025721, -0.027870, -0.029828, -0.031532, -0.032924, -0.033956, -0.034591, -0.034805, -0.034591, -0.033956, -0.032924, -0.031532, -0.029828, -0.027870, -0.025721, -0.023446, -0.021111, -0.018775, -0.016492, -0.014310, -0.012263, -0.010381, -0.008679, -0.007168, -0.005847, -0.004711, -0.003749, -0.002947, -0.002288, -0.001755, -0.001329
};


float apply_mask8even( float* image, int imagew, int imageh, int x, int y ) { 
float returnvalue; 
int mx1 = x-( 47/2) >= 0 ? 0 : ( 47/2)-x; 
int mx2 = x+( 47/2) < imagew ? 47-1 : 47-1 - (x+( 47/2) - (imagew-1)); 
int my1 = y-( 47/2) >= 0 ? 0 : ( 47/2)-y; 
int my2 = y+( 47/2) < imageh ? 47-1 : 47-1 - (y+( 47/2) - (imageh-1)); 
int imagexoffset = x - 47/2; 
int imageyoffset = y - 47/2; 

int h , k; 
returnvalue = 0; 
for(k = my1; k <= my2; k++){ 
for(h = mx1; h <= mx2; h++){ 
returnvalue += 
mask8even[h+k* 47] * image[imagexoffset+h+imagew*(imageyoffset+k)]; 
} 
} 
return returnvalue; 
}