mask11odd.c

一个有关eigenface的matlab 例程

 0.007077,  0.004900,  0.000107, -0.007963, -0.019662, -0.034795, -0.052407, -0.070642, -0.086744, -0.097261, -0.098456, -0.086916, -0.060285, -0.017974,  0.038291,  0.104251,  0.173335,  0.237362,  0.287648,  0.316370,  0.317928,  0.290033,  0.234307,  0.156215,  0.064335, -0.030922, -0.118909, -0.190344, -0.238629, -0.260650, -0.256949, -0.231273, -0.189661, -0.139264, -0.087150, -0.039310, -0.000000,  0.028517,  0.045865,  0.053169,  0.052530,  0.046469,  0.037453,  0.027562,  0.018305,  0.010593,  0.004800,
 0.008381,  0.007908,  0.005408,  0.000117, -0.008575, -0.020913, -0.036556, -0.054384, -0.072407, -0.087822, -0.097261, -0.097248, -0.084796, -0.058093, -0.017108,  0.035999,  0.096809,  0.158986,  0.215040,  0.257401,  0.279630,  0.277559,  0.250099,  0.199567,  0.131421,  0.053460, -0.025380, -0.096399, -0.152418, -0.188738, -0.203626, -0.198271, -0.176269, -0.142780, -0.103554, -0.064008, -0.028517, -0.000000,  0.020183,  0.032063,  0.036713,  0.035827,  0.031304,  0.024921,  0.018114,  0.011883,  0.006792,
 0.008471,  0.009251,  0.008622,  0.005824,  0.000125, -0.009009, -0.021703, -0.037470, -0.055060, -0.072407, -0.086744, -0.094889, -0.093712, -0.080711, -0.054616, -0.015887,  0.033019,  0.087705,  0.142268,  0.190067,  0.224717,  0.241128,  0.236406,  0.210405,  0.165833,  0.107866,  0.043340, -0.020323, -0.076244, -0.119072, -0.145637, -0.155198, -0.149263, -0.131071, -0.104867, -0.075123, -0.045865, -0.020183, -0.000000,  0.013937,  0.021868,  0.024732,  0.023839,  0.020574,  0.016178,  0.011615,  0.007526,
 0.007670,  0.009236,  0.009963,  0.009171,  0.006119,  0.000129, -0.009234, -0.021972, -0.037470, -0.054384, -0.070642, -0.083590, -0.090318, -0.088103, -0.074949, -0.050094, -0.014393,  0.029547,  0.077519,  0.124203,  0.163897,  0.191399,  0.202857,  0.196444,  0.172693,  0.134440,  0.086374,  0.034279, -0.015876, -0.058833, -0.090753, -0.109638, -0.115402, -0.109627, -0.095085, -0.075142, -0.053169, -0.032063, -0.013937, -0.000000,  0.009389,  0.014551,  0.016255,  0.015475,  0.013192,  0.010246,  0.007266,
 0.006328,  0.008260,  0.009824,  0.010467,  0.009517,  0.006272,  0.000131, -0.009234, -0.021703, -0.036556, -0.052407, -0.067238, -0.078587, -0.083870, -0.080809, -0.067901, -0.044827, -0.012721,  0.025795,  0.066846,  0.105788,  0.137884,  0.159045,  0.166499,  0.159257,  0.138284,  0.106332,  0.067477,  0.026451, -0.012100, -0.044291, -0.067483, -0.080525, -0.083719, -0.078553, -0.067297, -0.052530, -0.036713, -0.021868, -0.009389, -0.000000,  0.006170,  0.009446,  0.010423,  0.009801,  0.008252,  0.006331,
 0.004767,  0.006731,  0.008678,  0.010194,  0.010728,  0.009635,  0.006272,  0.000129, -0.009009, -0.020913, -0.034795, -0.049270, -0.062438, -0.072081, -0.075983, -0.072312, -0.060015, -0.039135, -0.010970,  0.021970,  0.056236,  0.087906,  0.113171,  0.128937,  0.133324,  0.125960,  0.108031,  0.082050,  0.051429,  0.019913, -0.008998, -0.032530, -0.048955, -0.057700, -0.059252, -0.054915, -0.046469, -0.035827, -0.024732, -0.014551, -0.006170, -0.000000,  0.003956,  0.005982,  0.006520,  0.006056,  0.005037,
 0.003239,  0.005008,  0.006985,  0.008895,  0.010321,  0.010728,  0.009517,  0.006119,  0.000125, -0.008575, -0.019662, -0.032311, -0.045191, -0.056566, -0.064502, -0.067159, -0.063130, -0.051752, -0.033333, -0.009229,  0.018256,  0.046157,  0.071265,  0.090621,  0.101980,  0.104155,  0.097196,  0.082338,  0.061769,  0.038242,  0.014625, -0.006527, -0.023309, -0.034648, -0.040337, -0.040914, -0.037453, -0.031304, -0.023839, -0.016255, -0.009446, -0.003956, -0.000000,  0.002475,  0.003696,  0.003979,  0.003651,
 0.001916,  0.003362,  0.005133,  0.007072,  0.008895,  0.010194,  0.010467,  0.009171,  0.005824,  0.000117, -0.007963, -0.018034, -0.029272, -0.040439, -0.049997, -0.056311, -0.057912, -0.053770, -0.043538, -0.027698, -0.007575,  0.014800,  0.036960,  0.056365,  0.070795,  0.078691,  0.079384,  0.073171,  0.061225,  0.045367,  0.027743,  0.010480, -0.004620, -0.016295, -0.023925, -0.027511, -0.027562, -0.024921, -0.020574, -0.015475, -0.010423, -0.005982, -0.002475, -0.000000,  0.001511,  0.002228,  0.002369,
 0.000885,  0.001964,  0.003403,  0.005133,  0.006985,  0.008678,  0.009824,  0.009963,  0.008622,  0.005408,  0.000107, -0.007214, -0.016138, -0.025873, -0.035304, -0.043113, -0.047962, -0.048720, -0.044681, -0.035735, -0.022455, -0.006066,  0.011706,  0.028874,  0.043493,  0.053958,  0.059240,  0.059029,  0.053741,  0.044416,  0.032508,  0.019635,  0.007326, -0.003190, -0.011113, -0.016117, -0.018305, -0.018114, -0.016178, -0.013192, -0.009801, -0.006520, -0.003696, -0.001511, -0.000000,  0.000899,  0.001311,
 0.000164,  0.000896,  0.001964,  0.003362,  0.005008,  0.006731,  0.008260,  0.009236,  0.009251,  0.007908,  0.004900,  0.000096, -0.006376, -0.014088, -0.022311, -0.030070, -0.036270, -0.039855, -0.039988, -0.036223, -0.028615, -0.017760, -0.004739,  0.009033,  0.022007,  0.032743,  0.040122,  0.043510,  0.042822,  0.038508,  0.031436,  0.022725,  0.013558,  0.004997, -0.002149, -0.007395, -0.010593, -0.011883, -0.011615, -0.010246, -0.008252, -0.006056, -0.003979, -0.002228, -0.000899, -0.000000,  0.000522,
-0.000275,  0.000164,  0.000885,  0.001916,  0.003239,  0.004767,  0.006328,  0.007670,  0.008471,  0.008381,  0.007077,  0.004331,  0.000084, -0.005498, -0.012000, -0.018769, -0.024987, -0.029769, -0.032310, -0.032020, -0.028649, -0.022354, -0.013705, -0.003612,  0.006800,  0.016364,  0.024048,  0.029107,  0.031177,  0.030308,  0.026920,  0.021706,  0.015499,  0.009134,  0.003325, -0.001412, -0.004800, -0.006792, -0.007526, -0.007266, -0.006331, -0.005037, -0.003651, -0.002369, -0.001311, -0.000522, -0.000000
};


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