mask7odd.c

一个有关eigenface的matlab 例程

float mask7odd[529] = {
 0.000000,  0.002207,  0.005995,  0.010154,  0.011825,  0.007498, -0.004610, -0.022054, -0.037764, -0.043073, -0.032962, -0.009984,  0.016379,  0.034948,  0.039155,  0.030096,  0.014642,  0.000615, -0.007354, -0.008965, -0.006699, -0.003448, -0.000997,
-0.002207,  0.000000,  0.005904,  0.015266,  0.024614,  0.027290,  0.016474, -0.009640, -0.043911, -0.071582, -0.077728, -0.056630, -0.016332,  0.025499,  0.051799,  0.055250,  0.040430,  0.018726,  0.000750, -0.008523, -0.009892, -0.007036, -0.003448,
-0.005995, -0.005904,  0.000000,  0.014311,  0.035231,  0.054076,  0.057080,  0.032804, -0.018269, -0.079236, -0.122972, -0.127123, -0.088175, -0.024213,  0.035977,  0.069582,  0.070657,  0.049222,  0.021705,  0.000828, -0.008952, -0.009892, -0.006699,
-0.010154, -0.015266, -0.014311,  0.000000,  0.031440,  0.073685,  0.107672,  0.108198,  0.059200, -0.031380, -0.129584, -0.191461, -0.188426, -0.124425, -0.032530,  0.046006,  0.084713,  0.081893,  0.054312,  0.022800,  0.000828, -0.008523, -0.008965,
-0.011825, -0.024614, -0.035231, -0.031440,  0.000000,  0.062600,  0.139672,  0.194300,  0.185879,  0.096824, -0.048852, -0.192066, -0.270163, -0.253121, -0.159124, -0.039608,  0.053320,  0.093470,  0.086023,  0.054312,  0.021705,  0.000750, -0.007354,
-0.007498, -0.027290, -0.054076, -0.073685, -0.062600,  0.000000,  0.112963,  0.239945,  0.317769,  0.289407,  0.143520, -0.068926, -0.258004, -0.345497, -0.308168, -0.184431, -0.043705,  0.056008,  0.093470,  0.081893,  0.049222,  0.018726,  0.000615,
 0.004610, -0.016474, -0.057080, -0.107672, -0.139672, -0.112963,  0.000000,  0.184745,  0.373581,  0.471002,  0.408376,  0.192802, -0.088140, -0.314106, -0.400438, -0.340029, -0.193732, -0.043705,  0.053320,  0.084713,  0.070657,  0.040429,  0.014642,
 0.022054,  0.009640, -0.032804, -0.108198, -0.194300, -0.239945, -0.184745,  0.000000,  0.273830,  0.527145,  0.632711,  0.522254,  0.234733, -0.102151, -0.346578, -0.420628, -0.340029, -0.184431, -0.039608,  0.046006,  0.069582,  0.055250,  0.030096,
 0.037764,  0.043911,  0.018269, -0.059200, -0.185879, -0.317769, -0.373581, -0.273830, -0.000000,  0.367842,  0.674137,  0.770302,  0.605306,  0.259004, -0.107301, -0.346578, -0.400438, -0.308168, -0.159124, -0.032530,  0.035977,  0.051799,  0.039155,
 0.043073,  0.071582,  0.079236,  0.031380, -0.096824, -0.289407, -0.471002, -0.527145, -0.367842,  0.000000,  0.447832,  0.781337,  0.849940,  0.635827,  0.259004, -0.102151, -0.314106, -0.345497, -0.253121, -0.124425, -0.024213,  0.025499,  0.034948,
 0.032962,  0.077728,  0.122972,  0.129584,  0.048852, -0.143520, -0.408376, -0.632711, -0.674137, -0.447832, -0.000000,  0.494131,  0.820734,  0.849940,  0.605306,  0.234733, -0.088140, -0.258004, -0.270163, -0.188426, -0.088175, -0.016332,  0.016379,
 0.009984,  0.056630,  0.127123,  0.191461,  0.192066,  0.068926, -0.192802, -0.522254, -0.770302, -0.781337, -0.494131, -0.000000,  0.494131,  0.781338,  0.770302,  0.522254,  0.192802, -0.068926, -0.192066, -0.191461, -0.127123, -0.056630, -0.009984,
-0.016379,  0.016332,  0.088175,  0.188426,  0.270163,  0.258004,  0.088140, -0.234734, -0.605306, -0.849940, -0.820734, -0.494131,  0.000000,  0.447832,  0.674137,  0.632711,  0.408376,  0.143520, -0.048852, -0.129584, -0.122972, -0.077728, -0.032962,
-0.034948, -0.025499,  0.024213,  0.124425,  0.253121,  0.345497,  0.314106,  0.102151, -0.259004, -0.635827, -0.849940, -0.781338, -0.447832, -0.000000,  0.367842,  0.527145,  0.471002,  0.289407,  0.096824, -0.031380, -0.079236, -0.071582, -0.043073,
-0.039155, -0.051799, -0.035977,  0.032530,  0.159124,  0.308168,  0.400438,  0.346578,  0.107301, -0.259004, -0.605306, -0.770302, -0.674137, -0.367842,  0.000000,  0.273830,  0.373581,  0.317769,  0.185879,  0.059200, -0.018269, -0.043911, -0.037764,
-0.030096, -0.055250, -0.069582, -0.046006,  0.039608,  0.184431,  0.340029,  0.420628,  0.346578,  0.102151, -0.234734, -0.522254, -0.632711, -0.527145, -0.273830, -0.000000,  0.184745,  0.239945,  0.194300,  0.108198,  0.032804, -0.009640, -0.022054,
-0.014642, -0.040429, -0.070657, -0.084713, -0.053320,  0.043705,  0.193732,  0.340029,  0.400438,  0.314106,  0.088140, -0.192802, -0.408376, -0.471002, -0.373581, -0.184745, -0.000000,  0.112963,  0.139672,  0.107672,  0.057080,  0.016474, -0.004610,
-0.000615, -0.018726, -0.049222, -0.081893, -0.093470, -0.056008,  0.043705,  0.184431,  0.308168,  0.345497,  0.258004,  0.068926, -0.143520, -0.289407, -0.317769, -0.239945, -0.112963, -0.000000,  0.062600,  0.073685,  0.054076,  0.027290,  0.007498,
 0.007354, -0.000750, -0.021705, -0.054312, -0.086023, -0.093470, -0.053320,  0.039608,  0.159124,  0.253121,  0.270163,  0.192066,  0.048852, -0.096824, -0.185879, -0.194300, -0.139672, -0.062600, -0.000000,  0.031440,  0.035231,  0.024614,  0.011825,
 0.008965,  0.008523, -0.000828, -0.022800, -0.054312, -0.081893, -0.084713, -0.046006,  0.032530,  0.124425,  0.188426,  0.191461,  0.129584,  0.031380, -0.059200, -0.108198, -0.107672, -0.073685, -0.031440, -0.000000,  0.014311,  0.015266,  0.010154,
 0.006699,  0.009892,  0.008952, -0.000828, -0.021705, -0.049222, -0.070657, -0.069582, -0.035977,  0.024213,  0.088175,  0.127123,  0.122972,  0.079236,  0.018269, -0.032804, -0.057080, -0.054076, -0.035231, -0.014311, -0.000000,  0.005904,  0.005995,
 0.003448,  0.007036,  0.009892,  0.008523, -0.000750, -0.018726, -0.040430, -0.055250, -0.051799, -0.025499,  0.016332,  0.056630,  0.077728,  0.071582,  0.043911,  0.009640, -0.016474, -0.027290, -0.024614, -0.015266, -0.005904, -0.000000,  0.002207,
 0.000997,  0.003448,  0.006699,  0.008965,  0.007354, -0.000615, -0.014642, -0.030096, -0.039155, -0.034948, -0.016379,  0.009984,  0.032962,  0.043073,  0.037764,  0.022054,  0.004610, -0.007498, -0.011825, -0.010154, -0.005995, -0.002207, -0.000000
};


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