mask10odd.c

一个有关eigenface的matlab 例程

-0.005241, -0.012582, -0.021299, -0.029568, -0.034619, -0.033129, -0.022002,  0.000560,  0.033912,  0.074149,  0.113973,  0.143571,  0.152536,  0.132551,  0.080200, -0.000916, -0.099569, -0.198303, -0.276769, -0.316367, -0.305006, -0.240670, -0.132657, -0.000000,  0.132657,  0.240670,  0.305006,  0.316367,  0.276769,  0.198303,  0.099569,  0.000916, -0.080199, -0.132551, -0.152536, -0.143571, -0.113973, -0.074149, -0.033912, -0.000560,  0.022002,  0.033129,  0.034619,  0.029568,  0.021299,  0.012582,  0.005241,
-0.004436, -0.010650, -0.018029, -0.025029, -0.029305, -0.028044, -0.018624,  0.000474,  0.028706,  0.062766,  0.096477,  0.121531,  0.129121,  0.112203,  0.067888, -0.000775, -0.084285, -0.167862, -0.234283, -0.267802, -0.258185, -0.203725, -0.112293, -0.000000,  0.112293,  0.203725,  0.258185,  0.267802,  0.234283,  0.167862,  0.084285,  0.000775, -0.067888, -0.112203, -0.129121, -0.121531, -0.096477, -0.062766, -0.028706, -0.000474,  0.018624,  0.028044,  0.029305,  0.025029,  0.018029,  0.010650,  0.004436,
-0.003709, -0.008905, -0.015074, -0.020927, -0.024502, -0.023448, -0.015572,  0.000396,  0.024001,  0.052479,  0.080665,  0.101613,  0.107959,  0.093814,  0.056762, -0.000648, -0.070471, -0.140350, -0.195885, -0.223911, -0.215870, -0.170336, -0.093889, -0.000000,  0.093889,  0.170335,  0.215870,  0.223911,  0.195885,  0.140350,  0.070471,  0.000648, -0.056762, -0.093814, -0.107958, -0.101613, -0.080665, -0.052479, -0.024001, -0.000396,  0.015572,  0.023448,  0.024502,  0.020927,  0.015074,  0.008905,  0.003709,
-0.003063, -0.007354, -0.012449, -0.017282, -0.020235, -0.019364, -0.012860,  0.000327,  0.019822,  0.043340,  0.066617,  0.083917,  0.089157,  0.077476,  0.046877, -0.000535, -0.058198, -0.115908, -0.161771, -0.184916, -0.178276, -0.140671, -0.077538, -0.000000,  0.077538,  0.140671,  0.178276,  0.184916,  0.161771,  0.115908,  0.058198,  0.000535, -0.046877, -0.077476, -0.089157, -0.083917, -0.066617, -0.043340, -0.019822, -0.000327,  0.012860,  0.019364,  0.020235,  0.017282,  0.012449,  0.007354,  0.003063,
-0.002499, -0.005999, -0.010155, -0.014097, -0.016506, -0.015796, -0.010490,  0.000267,  0.016169,  0.035353,  0.054340,  0.068452,  0.072727,  0.063198,  0.038238, -0.000437, -0.047473, -0.094548, -0.131959, -0.150839, -0.145422, -0.114748, -0.063249, -0.000000,  0.063249,  0.114748,  0.145422,  0.150839,  0.131959,  0.094548,  0.047473,  0.000437, -0.038238, -0.063198, -0.072727, -0.068452, -0.054340, -0.035353, -0.016169, -0.000267,  0.010490,  0.015796,  0.016506,  0.014097,  0.010155,  0.005999,  0.002499,
-0.002013, -0.004833, -0.008182, -0.011358, -0.013299, -0.012727, -0.008452,  0.000215,  0.013027,  0.028484,  0.043782,  0.055153,  0.058597,  0.050919,  0.030809, -0.000352, -0.038249, -0.076178, -0.106321, -0.121532, -0.117168, -0.092453, -0.050960, -0.000000,  0.050960,  0.092453,  0.117168,  0.121532,  0.106321,  0.076178,  0.038250,  0.000352, -0.030809, -0.050919, -0.058597, -0.055153, -0.043782, -0.028484, -0.013027, -0.000215,  0.008452,  0.012727,  0.013299,  0.011358,  0.008182,  0.004833,  0.002013,
-0.001602, -0.003846, -0.006511, -0.009039, -0.010583, -0.010128, -0.006726,  0.000171,  0.010367,  0.022668,  0.034843,  0.043892,  0.046633,  0.040523,  0.024518, -0.000280, -0.030440, -0.060624, -0.084612, -0.096718, -0.093245, -0.073576, -0.040555, -0.000000,  0.040555,  0.073576,  0.093245,  0.096718,  0.084612,  0.060624,  0.030440,  0.000280, -0.024518, -0.040523, -0.046633, -0.043892, -0.034843, -0.022668, -0.010367, -0.000171,  0.006726,  0.010128,  0.010583,  0.009039,  0.006511,  0.003846,  0.001602,
-0.001259, -0.003024, -0.005118, -0.007105, -0.008319, -0.007961, -0.005287,  0.000134,  0.008149,  0.017819,  0.027389,  0.034501,  0.036656,  0.031853,  0.019273, -0.000220, -0.023927, -0.047654, -0.066510, -0.076026, -0.073296, -0.057835, -0.031879, -0.000000,  0.031879,  0.057835,  0.073296,  0.076026,  0.066510,  0.047654,  0.023927,  0.000220, -0.019273, -0.031853, -0.036656, -0.034501, -0.027389, -0.017819, -0.008149, -0.000134,  0.005287,  0.007961,  0.008319,  0.007105,  0.005118,  0.003024,  0.001259,
-0.000978, -0.002348, -0.003974, -0.005517, -0.006459, -0.006181, -0.004105,  0.000104,  0.006327,  0.013835,  0.021265,  0.026787,  0.028460,  0.024731,  0.014964, -0.000171, -0.018578, -0.036999, -0.051639, -0.059028, -0.056908, -0.044904, -0.024751, -0.000000,  0.024751,  0.044904,  0.056908,  0.059028,  0.051639,  0.036999,  0.018578,  0.000171, -0.014964, -0.024731, -0.028460, -0.026787, -0.021265, -0.013835, -0.006327, -0.000104,  0.004105,  0.006181,  0.006459,  0.005517,  0.003974,  0.002348,  0.000978,
-0.000750, -0.001800, -0.003048, -0.004231, -0.004953, -0.004740, -0.003148,  0.000080,  0.004852,  0.010610,  0.016308,  0.020543,  0.021826,  0.018966,  0.011475, -0.000131, -0.014247, -0.028374, -0.039602, -0.045268, -0.043642, -0.034436, -0.018981, -0.000000,  0.018981,  0.034436,  0.043642,  0.045268,  0.039602,  0.028374,  0.014247,  0.000131, -0.011475, -0.018966, -0.021826, -0.020543, -0.016308, -0.010610, -0.004852, -0.000080,  0.003148,  0.004740,  0.004953,  0.004231,  0.003048,  0.001800,  0.000750,
-0.000568, -0.001364, -0.002308, -0.003205, -0.003752, -0.003591, -0.002385,  0.000061,  0.003676,  0.008037,  0.012353,  0.015561,  0.016533,  0.014366,  0.008692, -0.000099, -0.010792, -0.021493, -0.029997, -0.034289, -0.033058, -0.026085, -0.014378, -0.000000,  0.014378,  0.026085,  0.033058,  0.034289,  0.029997,  0.021493,  0.010792,  0.000099, -0.008692, -0.014366, -0.016533, -0.015561, -0.012353, -0.008037, -0.003676, -0.000061,  0.002385,  0.003591,  0.003752,  0.003205,  0.002308,  0.001364,  0.000568
};


float apply_mask10odd( 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 += 
mask10odd[h+k* 47] * image[imagexoffset+h+imagew*(imageyoffset+k)]; 
} 
} 
return returnvalue; 
}