mask9even.c

一个有关eigenface的matlab 例程

-0.012442, -0.014063, -0.013855, -0.010708, -0.003502,  0.008653,  0.026121,  0.048430,  0.074029,  0.100176,  0.123044,  0.138088,  0.140665,  0.126850,  0.094310,  0.043068, -0.024024, -0.101265, -0.180616, -0.252706, -0.308197, -0.339288, -0.341047, -0.312322, -0.256031, -0.178755, -0.089703,  0.000739,  0.082612,  0.147827,  0.191193,  0.210887,  0.208325,  0.187516,  0.154102,  0.114277,  0.073799,  0.037265,  0.007701, -0.013497, -0.026366, -0.031993, -0.032069, -0.028477, -0.022957, -0.016911, -0.011304,
-0.009017, -0.008994, -0.007038, -0.002330,  0.005829,  0.017815,  0.033441,  0.051752,  0.070900,  0.088167,  0.100176,  0.103313,  0.094324,  0.070998,  0.032825, -0.018538, -0.079111, -0.142854, -0.202354, -0.249854, -0.278476, -0.283397, -0.262751, -0.218070, -0.154143, -0.078313,  0.000653,  0.073925,  0.133925,  0.175365,  0.195832,  0.195855,  0.178481,  0.148500,  0.111490,  0.072894,  0.037265,  0.007797, -0.013834, -0.027361, -0.033612, -0.034111, -0.030666, -0.025030, -0.018666, -0.012632, -0.007563,
-0.005696, -0.004513, -0.001513,  0.003831,  0.011854,  0.022528,  0.035296,  0.048956,  0.061635,  0.070900,  0.074029,  0.068427,  0.052145,  0.024408, -0.013956, -0.060296, -0.110232, -0.158084, -0.197616, -0.222990, -0.229749, -0.215657, -0.181208, -0.129677, -0.066702,  0.000563,  0.064539,  0.118372,  0.156925,  0.177416,  0.179641,  0.165739,  0.139611,  0.106118,  0.070244,  0.036356,  0.007701, -0.013834, -0.027701, -0.034452, -0.035398, -0.032218, -0.026623, -0.020101, -0.013772, -0.008348, -0.004163,
-0.002823, -0.000958,  0.002457,  0.007695,  0.014806,  0.023486,  0.032980,  0.042037,  0.048956,  0.051752,  0.048430,  0.037365,  0.017707, -0.010250, -0.044835, -0.082985, -0.120487, -0.152488, -0.174205, -0.181714, -0.172688, -0.146905, -0.106435, -0.055426,  0.000474,  0.054970,  0.102074,  0.136999,  0.156813,  0.160751,  0.150153,  0.128053,  0.098542,  0.066039,  0.034604,  0.007421, -0.013497, -0.027361, -0.034452, -0.035838, -0.033023, -0.027628, -0.021118, -0.014649, -0.008990, -0.004539, -0.001391,
-0.000592,  0.001537,  0.004874,  0.009493,  0.015246,  0.021675,  0.027971,  0.032980,  0.035296,  0.033441,  0.026121,  0.012532, -0.007345, -0.032526, -0.060949, -0.089592, -0.114796, -0.132774, -0.140218, -0.134907, -0.116191, -0.085228, -0.044934,  0.000389,  0.045678,  0.085873,  0.116687,  0.135222,  0.140340,  0.132716,  0.114588,  0.089275,  0.060572,  0.032134,  0.006977, -0.012847, -0.026366, -0.033612, -0.035398, -0.033023, -0.027971, -0.021646, -0.015202, -0.009445, -0.004828, -0.001498,  0.000622,
 0.000938,  0.003011,  0.005939,  0.009656,  0.013898,  0.018158,  0.021675,  0.023486,  0.022528,  0.017815,  0.008653, -0.005134, -0.023020, -0.043673, -0.064995, -0.084313, -0.098728, -0.105559, -0.102822, -0.089657, -0.066582, -0.035539,  0.000312,  0.037031,  0.070482,  0.096962,  0.113760,  0.119532,  0.114443,  0.100038,  0.078908,  0.054203,  0.029112,  0.006399, -0.011930, -0.024788, -0.031993, -0.034111, -0.032218, -0.027628, -0.021646, -0.015391, -0.009681, -0.005010, -0.001574,  0.000661,  0.001878,
 0.001815,  0.003624,  0.005966,  0.008694,  0.011500,  0.013898,  0.015246,  0.014806,  0.011854,  0.005829, -0.003502, -0.015896, -0.030531, -0.046000, -0.060415, -0.071622, -0.077529, -0.076457, -0.067495, -0.050746, -0.027423,  0.000244,  0.029289,  0.056439,  0.078607,  0.093370,  0.099327,  0.096279,  0.085206,  0.068043,  0.047320,  0.025731,  0.005727, -0.010808, -0.022736, -0.029709, -0.032069, -0.030666, -0.026623, -0.021118, -0.015202, -0.009681, -0.005072, -0.001613,  0.000686,  0.001973,  0.002476,
 0.002158,  0.003597,  0.005306,  0.007106,  0.008694,  0.009656,  0.009493,  0.007695,  0.003831, -0.002330, -0.010708, -0.020823, -0.031763, -0.042234, -0.050691, -0.055553, -0.055466, -0.049573, -0.037734, -0.020645,  0.000186,  0.022601,  0.044091,  0.062172,  0.074766,  0.080524,  0.079022,  0.070803,  0.057243,  0.040304,  0.022188,  0.004999, -0.009553, -0.020345, -0.026915, -0.029414, -0.028477, -0.025030, -0.020101, -0.014649, -0.009445, -0.005010, -0.001613,  0.000695,  0.002022,  0.002569,  0.002569,
 0.002115,  0.003159,  0.004283,  0.005306,  0.005966,  0.005939,  0.004874,  0.002457, -0.001513, -0.007038, -0.013855, -0.021398, -0.028805, -0.035002, -0.038836, -0.039256, -0.035521, -0.027374, -0.015163,  0.000138,  0.017014,  0.033605,  0.047975,  0.058409,  0.063688,  0.063277,  0.057399,  0.046983,  0.033491,  0.018666,  0.004258, -0.008237, -0.017762, -0.023789, -0.026321, -0.025799, -0.022957, -0.018666, -0.013772, -0.008990, -0.004828, -0.001574,  0.000686,  0.002022,  0.002601,  0.002633,  0.002326,
 0.001835,  0.002519,  0.003159,  0.003597,  0.003624,  0.003011,  0.001537, -0.000958, -0.004513, -0.008994, -0.014063, -0.019167, -0.023580, -0.026487, -0.027106, -0.024832, -0.019374, -0.010865,  0.000100,  0.012496,  0.024988,  0.036116,  0.044518,  0.049144,  0.049433,  0.045399,  0.037622,  0.027151,  0.015321,  0.003538, -0.006930, -0.015128, -0.020514, -0.022979, -0.022803, -0.020543, -0.016911, -0.012632, -0.008348, -0.004539, -0.001498,  0.000661,  0.001973,  0.002569,  0.002633,  0.002355,  0.001901,
 0.001445,  0.001835,  0.002115,  0.002158,  0.001815,  0.000938, -0.000592, -0.002823, -0.005696, -0.009017, -0.012442, -0.015497, -0.017624, -0.018260, -0.016936, -0.013378, -0.007595,  0.000071,  0.008954,  0.018128,  0.026526,  0.033103,  0.036997,  0.037677,  0.035031,  0.029391,  0.021474,  0.012268,  0.002868, -0.005688, -0.012571, -0.017258, -0.019572, -0.019663, -0.017935, -0.014947, -0.011304, -0.007563, -0.004163, -0.001391,  0.000622,  0.001878,  0.002476,  0.002569,  0.002326,  0.001901,  0.001419
};


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