mask7even.c

一个有关eigenface的matlab 例程

float mask7even[529] = {
 0.002508,  0.003579,  0.003012, -0.000887, -0.008800, -0.018858, -0.026197, -0.024803, -0.011326,  0.011514,  0.034892,  0.048171,  0.045118,  0.027685,  0.004616, -0.013994, -0.022207, -0.020285, -0.012775, -0.004849,  0.000381,  0.002407,  0.002302,
 0.003579,  0.007048,  0.009572,  0.007669, -0.002149, -0.020307, -0.041431, -0.054792, -0.049388, -0.021473,  0.020775,  0.059939,  0.078782,  0.070249,  0.041038,  0.006516, -0.018798, -0.028398, -0.024696, -0.014807, -0.005351,  0.000400,  0.002407,
 0.003012,  0.009572,  0.017946,  0.023204,  0.017699, -0.004721, -0.042471, -0.082493, -0.103860, -0.089126, -0.036893,  0.033972,  0.093320,  0.116771,  0.099128,  0.055131,  0.008335, -0.022884, -0.032914, -0.027250, -0.015553, -0.005351,  0.000381,
-0.000887,  0.007669,  0.023204,  0.041415,  0.050979,  0.037018, -0.009398, -0.080504, -0.148860, -0.178425, -0.145765, -0.057447,  0.050349,  0.131677,  0.156862,  0.126772,  0.067122,  0.009662, -0.025250, -0.034574, -0.027250, -0.014807, -0.004849,
-0.008800, -0.002149,  0.017699,  0.050979,  0.086619,  0.101505,  0.070171, -0.016957, -0.138296, -0.243453, -0.277801, -0.216061, -0.081068,  0.067630,  0.168393,  0.190974,  0.146933,  0.074063,  0.010150, -0.025250, -0.032914, -0.024696, -0.012775,
-0.018858, -0.020307, -0.004721,  0.037018,  0.101505,  0.164190,  0.183171,  0.120550, -0.027728, -0.215316, -0.360847, -0.391998, -0.290246, -0.103679,  0.082332,  0.195169,  0.210718,  0.154342,  0.074063,  0.009662, -0.022884, -0.028398, -0.020285,
-0.026197, -0.041431, -0.042471, -0.009398,  0.070171,  0.183171,  0.282065,  0.299568,  0.187692, -0.041094, -0.303820, -0.484735, -0.501307, -0.353368, -0.120170,  0.090842,  0.205009,  0.210718,  0.146933,  0.067122,  0.008335, -0.018798, -0.022207,
-0.024803, -0.054792, -0.082493, -0.080504, -0.016957,  0.120550,  0.299568,  0.439161,  0.444025,  0.264849, -0.055198, -0.388535, -0.590144, -0.581026, -0.389903, -0.126231,  0.090842,  0.195169,  0.190974,  0.126772,  0.055131,  0.006516, -0.013994,
-0.011326, -0.049388, -0.103860, -0.148860, -0.138296, -0.027728,  0.187693,  0.444025,  0.619684,  0.596473,  0.338705, -0.067198, -0.450318, -0.651156, -0.610324, -0.389903, -0.120170,  0.082332,  0.168393,  0.156862,  0.099128,  0.041038,  0.004616,
 0.011514, -0.021473, -0.089126, -0.178425, -0.243453, -0.215316, -0.041094,  0.264849,  0.596473,  0.792481,  0.726183,  0.392568, -0.074143, -0.473022, -0.651156, -0.581027, -0.353368, -0.103679,  0.067630,  0.131677,  0.116771,  0.070249,  0.027685,
 0.034892,  0.020775, -0.036893, -0.145766, -0.277801, -0.360847, -0.303820, -0.055198,  0.338705,  0.726183,  0.918503,  0.801261,  0.412362, -0.074143, -0.450317, -0.590145, -0.501307, -0.290246, -0.081068,  0.050349,  0.093320,  0.078782,  0.045118,
 0.048171,  0.059939,  0.033972, -0.057447, -0.216061, -0.391998, -0.484735, -0.388535, -0.067198,  0.392568,  0.801261,  0.964815,  0.801261,  0.392568, -0.067198, -0.388535, -0.484735, -0.391998, -0.216061, -0.057447,  0.033972,  0.059939,  0.048171,
 0.045118,  0.078782,  0.093320,  0.050349, -0.081068, -0.290246, -0.501307, -0.590145, -0.450318, -0.074143,  0.412362,  0.801261,  0.918502,  0.726183,  0.338705, -0.055198, -0.303820, -0.360847, -0.277801, -0.145765, -0.036893,  0.020775,  0.034892,
 0.027685,  0.070249,  0.116771,  0.131677,  0.067630, -0.103679, -0.353368, -0.581027, -0.651156, -0.473022, -0.074143,  0.392568,  0.726183,  0.792481,  0.596473,  0.264849, -0.041094, -0.215316, -0.243453, -0.178425, -0.089126, -0.021473,  0.011514,
 0.004616,  0.041038,  0.099128,  0.156862,  0.168393,  0.082332, -0.120170, -0.389903, -0.610323, -0.651156, -0.450318, -0.067198,  0.338705,  0.596473,  0.619684,  0.444025,  0.187692, -0.027728, -0.138296, -0.148860, -0.103860, -0.049388, -0.011326,
-0.013994,  0.006516,  0.055131,  0.126772,  0.190974,  0.195169,  0.090842, -0.126231, -0.389903, -0.581027, -0.590145, -0.388535, -0.055198,  0.264849,  0.444025,  0.439161,  0.299568,  0.120550, -0.016957, -0.080504, -0.082493, -0.054792, -0.024803,
-0.022207, -0.018798,  0.008335,  0.067122,  0.146933,  0.210718,  0.205009,  0.090842, -0.120170, -0.353368, -0.501307, -0.484735, -0.303820, -0.041094,  0.187693,  0.299568,  0.282065,  0.183171,  0.070171, -0.009398, -0.042471, -0.041431, -0.026197,
-0.020285, -0.028398, -0.022884,  0.009662,  0.074063,  0.154342,  0.210718,  0.195169,  0.082332, -0.103679, -0.290246, -0.391998, -0.360847, -0.215316, -0.027728,  0.120550,  0.183171,  0.164190,  0.101505,  0.037018, -0.004721, -0.020307, -0.018858,
-0.012775, -0.024696, -0.032914, -0.025250,  0.010150,  0.074063,  0.146933,  0.190974,  0.168393,  0.067630, -0.081067, -0.216061, -0.277801, -0.243453, -0.138296, -0.016957,  0.070171,  0.101505,  0.086619,  0.050979,  0.017699, -0.002149, -0.008800,
-0.004849, -0.014807, -0.027250, -0.034574, -0.025250,  0.009662,  0.067122,  0.126772,  0.156862,  0.131677,  0.050349, -0.057447, -0.145766, -0.178425, -0.148860, -0.080504, -0.009398,  0.037018,  0.050979,  0.041415,  0.023204,  0.007669, -0.000887,
 0.000381, -0.005351, -0.015553, -0.027250, -0.032914, -0.022884,  0.008335,  0.055131,  0.099128,  0.116771,  0.093320,  0.033972, -0.036893, -0.089126, -0.103860, -0.082493, -0.042471, -0.004721,  0.017699,  0.023204,  0.017946,  0.009572,  0.003012,
 0.002407,  0.000400, -0.005351, -0.014807, -0.024696, -0.028398, -0.018798,  0.006516,  0.041039,  0.070249,  0.078782,  0.059939,  0.020775, -0.021473, -0.049388, -0.054792, -0.041431, -0.020307, -0.002149,  0.007669,  0.009572,  0.007048,  0.003579,
 0.002302,  0.002407,  0.000381, -0.004849, -0.012775, -0.020285, -0.022207, -0.013994,  0.004616,  0.027685,  0.045118,  0.048171,  0.034892,  0.011514, -0.011326, -0.024803, -0.026197, -0.018858, -0.008800, -0.000887,  0.003012,  0.003579,  0.002508
};


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

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