mask8odd.c

一个有关eigenface的matlab 例程

-0.012353, -0.016308, -0.021265, -0.027389, -0.034843, -0.043782, -0.054340, -0.066617, -0.080665, -0.096477, -0.113973, -0.132989, -0.153275, -0.174488, -0.196199, -0.217906, -0.239044, -0.259016, -0.277213, -0.293049, -0.305988, -0.315578, -0.321476, -0.323467, -0.321476, -0.315578, -0.305988, -0.293049, -0.277213, -0.259016, -0.239044, -0.217906, -0.196199, -0.174488, -0.153275, -0.132989, -0.113973, -0.096477, -0.080665, -0.066617, -0.054340, -0.043782, -0.034843, -0.027389, -0.021265, -0.016308, -0.012353,
-0.008037, -0.010610, -0.013835, -0.017819, -0.022668, -0.028484, -0.035353, -0.043340, -0.052479, -0.062766, -0.074149, -0.086521, -0.099718, -0.113519, -0.127644, -0.141766, -0.155518, -0.168512, -0.180350, -0.190652, -0.199071, -0.205310, -0.209147, -0.210442, -0.209147, -0.205310, -0.199071, -0.190652, -0.180350, -0.168512, -0.155518, -0.141766, -0.127644, -0.113519, -0.099718, -0.086521, -0.074149, -0.062766, -0.052479, -0.043340, -0.035353, -0.028484, -0.022668, -0.017819, -0.013835, -0.010610, -0.008037,
-0.003676, -0.004852, -0.006327, -0.008149, -0.010367, -0.013027, -0.016169, -0.019822, -0.024001, -0.028706, -0.033912, -0.039570, -0.045606, -0.051918, -0.058378, -0.064837, -0.071127, -0.077069, -0.082483, -0.087195, -0.091045, -0.093899, -0.095654, -0.096246, -0.095654, -0.093899, -0.091045, -0.087195, -0.082483, -0.077069, -0.071127, -0.064837, -0.058378, -0.051918, -0.045606, -0.039570, -0.033912, -0.028706, -0.024001, -0.019822, -0.016169, -0.013027, -0.010367, -0.008149, -0.006327, -0.004852, -0.003676,
-0.000061, -0.000080, -0.000104, -0.000134, -0.000171, -0.000215, -0.000267, -0.000327, -0.000396, -0.000474, -0.000560, -0.000653, -0.000753, -0.000857, -0.000963, -0.001070, -0.001174, -0.001272, -0.001361, -0.001439, -0.001502, -0.001549, -0.001578, -0.001588, -0.001578, -0.001549, -0.001502, -0.001439, -0.001361, -0.001272, -0.001174, -0.001070, -0.000963, -0.000857, -0.000753, -0.000653, -0.000560, -0.000474, -0.000396, -0.000327, -0.000267, -0.000215, -0.000171, -0.000134, -0.000104, -0.000080, -0.000061,
 0.002385,  0.003148,  0.004105,  0.005287,  0.006726,  0.008452,  0.010490,  0.012860,  0.015572,  0.018624,  0.022002,  0.025673,  0.029589,  0.033684,  0.037875,  0.042065,  0.046146,  0.050001,  0.053514,  0.056571,  0.059069,  0.060920,  0.062059,  0.062443,  0.062059,  0.060920,  0.059069,  0.056571,  0.053514,  0.050001,  0.046146,  0.042065,  0.037875,  0.033684,  0.029589,  0.025673,  0.022002,  0.018624,  0.015572,  0.012860,  0.010490,  0.008452,  0.006726,  0.005287,  0.004105,  0.003148,  0.002385,
 0.003591,  0.004740,  0.006181,  0.007961,  0.010128,  0.012727,  0.015796,  0.019364,  0.023448,  0.028044,  0.033129,  0.038657,  0.044554,  0.050720,  0.057031,  0.063340,  0.069485,  0.075290,  0.080580,  0.085183,  0.088944,  0.091732,  0.093446,  0.094025,  0.093446,  0.091732,  0.088944,  0.085183,  0.080580,  0.075290,  0.069485,  0.063340,  0.057031,  0.050720,  0.044554,  0.038657,  0.033129,  0.028044,  0.023448,  0.019364,  0.015796,  0.012727,  0.010128,  0.007961,  0.006181,  0.004740,  0.003591,
 0.003752,  0.004953,  0.006459,  0.008319,  0.010583,  0.013299,  0.016506,  0.020235,  0.024502,  0.029305,  0.034619,  0.040395,  0.046557,  0.053000,  0.059595,  0.066188,  0.072609,  0.078675,  0.084202,  0.089012,  0.092943,  0.095856,  0.097647,  0.098252,  0.097647,  0.095856,  0.092943,  0.089012,  0.084202,  0.078675,  0.072609,  0.066188,  0.059595,  0.053000,  0.046557,  0.040395,  0.034619,  0.029305,  0.024502,  0.020235,  0.016506,  0.013299,  0.010583,  0.008319,  0.006459,  0.004953,  0.003752,
 0.003205,  0.004231,  0.005517,  0.007105,  0.009039,  0.011358,  0.014097,  0.017282,  0.020927,  0.025029,  0.029568,  0.034501,  0.039764,  0.045267,  0.050899,  0.056531,  0.062014,  0.067196,  0.071916,  0.076025,  0.079381,  0.081869,  0.083400,  0.083916,  0.083400,  0.081869,  0.079381,  0.076025,  0.071916,  0.067196,  0.062014,  0.056531,  0.050899,  0.045267,  0.039764,  0.034501,  0.029568,  0.025029,  0.020927,  0.017282,  0.014097,  0.011358,  0.009039,  0.007105,  0.005517,  0.004231,  0.003205,
 0.002308,  0.003048,  0.003974,  0.005118,  0.006511,  0.008182,  0.010155,  0.012449,  0.015074,  0.018029,  0.021299,  0.024853,  0.028644,  0.032608,  0.036665,  0.040721,  0.044672,  0.048404,  0.051805,  0.054764,  0.057182,  0.058974,  0.060076,  0.060448,  0.060076,  0.058974,  0.057182,  0.054764,  0.051805,  0.048404,  0.044672,  0.040721,  0.036665,  0.032608,  0.028644,  0.024853,  0.021299,  0.018029,  0.015074,  0.012449,  0.010155,  0.008182,  0.006511,  0.005118,  0.003974,  0.003048,  0.002308,
 0.001364,  0.001800,  0.002348,  0.003024,  0.003846,  0.004833,  0.005999,  0.007354,  0.008905,  0.010650,  0.012582,  0.014681,  0.016921,  0.019262,  0.021659,  0.024055,  0.026389,  0.028594,  0.030602,  0.032351,  0.033779,  0.034838,  0.035489,  0.035708,  0.035489,  0.034838,  0.033779,  0.032351,  0.030602,  0.028594,  0.026389,  0.024055,  0.021659,  0.019262,  0.016921,  0.014681,  0.012582,  0.010650,  0.008905,  0.007354,  0.005999,  0.004833,  0.003846,  0.003024,  0.002348,  0.001800,  0.001364,
 0.000568,  0.000750,  0.000978,  0.001259,  0.001602,  0.002013,  0.002499,  0.003063,  0.003709,  0.004436,  0.005241,  0.006115,  0.007048,  0.008024,  0.009022,  0.010020,  0.010992,  0.011911,  0.012747,  0.013476,  0.014071,  0.014512,  0.014783,  0.014874,  0.014783,  0.014512,  0.014071,  0.013476,  0.012747,  0.011911,  0.010992,  0.010020,  0.009022,  0.008024,  0.007048,  0.006115,  0.005241,  0.004436,  0.003709,  0.003063,  0.002499,  0.002013,  0.001602,  0.001259,  0.000978,  0.000750,  0.000568
};


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