mask9odd.c

一个有关eigenface的matlab 例程

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


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