mask10even.c

一个有关eigenface的matlab 例程

-0.012263, -0.012316, -0.008569,  0.000287,  0.014716,  0.033714,  0.054384,  0.071917,  0.080211,  0.073146,  0.046367, -0.000838, -0.064141, -0.133847, -0.196043, -0.235245, -0.238085, -0.197148, -0.113854,  0.000616,  0.126860,  0.241139,  0.320704,  0.349206,  0.320704,  0.241139,  0.126860,  0.000616, -0.113853, -0.197148, -0.238085, -0.235245, -0.196043, -0.133847, -0.064141, -0.000838,  0.046367,  0.073146,  0.080211,  0.071917,  0.054384,  0.033714,  0.014716,  0.000287, -0.008569, -0.012316, -0.012263,
-0.010381, -0.010425, -0.007253,  0.000243,  0.012457,  0.028539,  0.046035,  0.060877,  0.067898,  0.061917,  0.039249, -0.000709, -0.054295, -0.113300, -0.165948, -0.199132, -0.201537, -0.166884, -0.096376,  0.000521,  0.107386,  0.204122,  0.271473,  0.295600,  0.271473,  0.204122,  0.107386,  0.000521, -0.096376, -0.166884, -0.201537, -0.199132, -0.165948, -0.113300, -0.054295, -0.000709,  0.039249,  0.061917,  0.067898,  0.060877,  0.046035,  0.028539,  0.012457,  0.000243, -0.007253, -0.010425, -0.010381,
-0.008679, -0.008716, -0.006065,  0.000203,  0.010415,  0.023862,  0.038490,  0.050900,  0.056770,  0.051769,  0.032816, -0.000593, -0.045396, -0.094731, -0.138750, -0.166496, -0.168506, -0.139533, -0.080581,  0.000436,  0.089786,  0.170667,  0.226980,  0.247153,  0.226980,  0.170667,  0.089786,  0.000436, -0.080580, -0.139533, -0.168506, -0.166496, -0.138750, -0.094731, -0.045396, -0.000593,  0.032816,  0.051769,  0.056770,  0.050900,  0.038490,  0.023862,  0.010415,  0.000203, -0.006065, -0.008716, -0.008679,
-0.007168, -0.007198, -0.005008,  0.000168,  0.008601,  0.019706,  0.031787,  0.042035,  0.046883,  0.042754,  0.027101, -0.000490, -0.037490, -0.078233, -0.114587, -0.137500, -0.139160, -0.115233, -0.066547,  0.000360,  0.074149,  0.140945,  0.187451,  0.204111,  0.187451,  0.140945,  0.074149,  0.000360, -0.066547, -0.115233, -0.139160, -0.137500, -0.114587, -0.078233, -0.037490, -0.000490,  0.027101,  0.042754,  0.046883,  0.042035,  0.031787,  0.019706,  0.008601,  0.000168, -0.005008, -0.007198, -0.007168,
-0.005847, -0.005872, -0.004085,  0.000137,  0.007016,  0.016075,  0.025929,  0.034289,  0.038243,  0.034875,  0.022107, -0.000400, -0.030581, -0.063816, -0.093470, -0.112161, -0.113515, -0.093997, -0.054284,  0.000294,  0.060485,  0.114971,  0.152907,  0.166496,  0.152907,  0.114971,  0.060485,  0.000294, -0.054284, -0.093997, -0.113515, -0.112161, -0.093470, -0.063816, -0.030581, -0.000400,  0.022107,  0.034875,  0.038243,  0.034289,  0.025929,  0.016075,  0.007016,  0.000137, -0.004085, -0.005872, -0.005847,
-0.004711, -0.004731, -0.003292,  0.000110,  0.005653,  0.012951,  0.020891,  0.027627,  0.030813,  0.028099,  0.017812, -0.000322, -0.024640, -0.051417, -0.075310, -0.090369, -0.091460, -0.075734, -0.043737,  0.000237,  0.048733,  0.092633,  0.123198,  0.134147,  0.123198,  0.092633,  0.048733,  0.000237, -0.043737, -0.075734, -0.091460, -0.090369, -0.075310, -0.051417, -0.024640, -0.000322,  0.017812,  0.028099,  0.030813,  0.027627,  0.020891,  0.012951,  0.005653,  0.000110, -0.003292, -0.004731, -0.004711,
-0.003749, -0.003765, -0.002620,  0.000088,  0.004499,  0.010307,  0.016626,  0.021986,  0.024522,  0.022362,  0.014175, -0.000256, -0.019609, -0.040919, -0.059933, -0.071918, -0.072786, -0.060271, -0.034807,  0.000188,  0.038783,  0.073719,  0.098044,  0.106757,  0.098044,  0.073720,  0.038783,  0.000188, -0.034807, -0.060271, -0.072786, -0.071918, -0.059933, -0.040919, -0.019609, -0.000256,  0.014175,  0.022362,  0.024522,  0.021986,  0.016626,  0.010307,  0.004499,  0.000088, -0.002620, -0.003765, -0.003749,
-0.002947, -0.002960, -0.002059,  0.000069,  0.003536,  0.008102,  0.013069,  0.017282,  0.019275,  0.017578,  0.011142, -0.000201, -0.015414, -0.032165, -0.047111, -0.056531, -0.057214, -0.047377, -0.027360,  0.000148,  0.030486,  0.057948,  0.077068,  0.083917,  0.077068,  0.057948,  0.030486,  0.000148, -0.027360, -0.047377, -0.057214, -0.056531, -0.047111, -0.032165, -0.015414, -0.000201,  0.011142,  0.017578,  0.019275,  0.017282,  0.013069,  0.008102,  0.003536,  0.000069, -0.002059, -0.002960, -0.002947,
-0.002288, -0.002298, -0.001599,  0.000054,  0.002746,  0.006290,  0.010147,  0.013418,  0.014966,  0.013647,  0.008651, -0.000156, -0.011967, -0.024973, -0.036578, -0.043892, -0.044422, -0.036784, -0.021243,  0.000115,  0.023669,  0.044991,  0.059837,  0.065155,  0.059837,  0.044992,  0.023669,  0.000115, -0.021243, -0.036784, -0.044422, -0.043892, -0.036578, -0.024973, -0.011967, -0.000156,  0.008651,  0.013647,  0.014966,  0.013418,  0.010147,  0.006290,  0.002746,  0.000054, -0.001599, -0.002298, -0.002288,
-0.001755, -0.001762, -0.001226,  0.000041,  0.002106,  0.004824,  0.007782,  0.010290,  0.011477,  0.010466,  0.006634, -0.000120, -0.009178, -0.019152, -0.028051, -0.033660, -0.034066, -0.028209, -0.016291,  0.000088,  0.018152,  0.034503,  0.045888,  0.049966,  0.045888,  0.034503,  0.018152,  0.000088, -0.016291, -0.028209, -0.034066, -0.033660, -0.028051, -0.019152, -0.009178, -0.000120,  0.006634,  0.010466,  0.011477,  0.010290,  0.007782,  0.004824,  0.002106,  0.000041, -0.001226, -0.001762, -0.001755,
-0.001329, -0.001335, -0.000929,  0.000031,  0.001595,  0.003654,  0.005894,  0.007795,  0.008694,  0.007928,  0.005025, -0.000091, -0.006952, -0.014507, -0.021248, -0.025497, -0.025805, -0.021368, -0.012340,  0.000067,  0.013750,  0.026136,  0.034759,  0.037848,  0.034759,  0.026136,  0.013750,  0.000067, -0.012340, -0.021368, -0.025805, -0.025497, -0.021248, -0.014507, -0.006952, -0.000091,  0.005025,  0.007928,  0.008694,  0.007795,  0.005894,  0.003654,  0.001595,  0.000031, -0.000929, -0.001335, -0.001329
};


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