mask6odd.c

一个有关eigenface的matlab 例程

float mask6odd[529] = {
-0.001848, -0.004201, -0.004574,  0.000319,  0.010702,  0.020283,  0.018431, -0.000520, -0.028195, -0.044568, -0.033795,  0.000000,  0.033795,  0.044568,  0.028195,  0.000520, -0.018431, -0.020283, -0.010702, -0.000319,  0.004574,  0.004201,  0.001848,
-0.003098, -0.007042, -0.007667,  0.000534,  0.017939,  0.033999,  0.030895, -0.000872, -0.047262, -0.074708, -0.056649,  0.000000,  0.056649,  0.074708,  0.047262,  0.000872, -0.030895, -0.033999, -0.017939, -0.000534,  0.007667,  0.007042,  0.003098,
-0.004943, -0.011237, -0.012234,  0.000853,  0.028627,  0.054256,  0.049302, -0.001392, -0.075421, -0.119219, -0.090401,  0.000000,  0.090401,  0.119219,  0.075421,  0.001392, -0.049302, -0.054256, -0.028627, -0.000853,  0.012234,  0.011237,  0.004943,
-0.007510, -0.017071, -0.018586,  0.001296,  0.043491,  0.082425,  0.074899, -0.002115, -0.114578, -0.181116, -0.137336,  0.000000,  0.137336,  0.181116,  0.114578,  0.002115, -0.074899, -0.082425, -0.043491, -0.001296,  0.018586,  0.017071,  0.007510,
-0.010861, -0.024689, -0.026881,  0.001874,  0.062898,  0.119208,  0.108323, -0.003059, -0.165709, -0.261940, -0.198624,  0.000000,  0.198624,  0.261940,  0.165709,  0.003059, -0.108323, -0.119208, -0.062898, -0.001874,  0.026881,  0.024689,  0.010861,
-0.014954, -0.033993, -0.037010,  0.002580,  0.086600,  0.164129,  0.149142, -0.004211, -0.228153, -0.360645, -0.273470,  0.000000,  0.273470,  0.360645,  0.228153,  0.004211, -0.149142, -0.164129, -0.086600, -0.002580,  0.037010,  0.033993,  0.014954,
-0.019601, -0.044555, -0.048510,  0.003382,  0.113509,  0.215128,  0.195484, -0.005520, -0.299046, -0.472707, -0.358444,  0.000000,  0.358445,  0.472707,  0.299046,  0.005520, -0.195484, -0.215128, -0.113509, -0.003382,  0.048510,  0.044555,  0.019601,
-0.024458, -0.055596, -0.060531,  0.004220,  0.141637,  0.268438,  0.243926, -0.006888, -0.373151, -0.589846, -0.447269,  0.000000,  0.447269,  0.589846,  0.373151,  0.006888, -0.243926, -0.268438, -0.141637, -0.004220,  0.060531,  0.055596,  0.024458,
-0.029053, -0.066043, -0.071905,  0.005013,  0.168251,  0.318878,  0.289760, -0.008182, -0.443266, -0.700679, -0.531311,  0.000000,  0.531311,  0.700679,  0.443266,  0.008182, -0.289760, -0.318878, -0.168251, -0.005013,  0.071905,  0.066043,  0.029053,
-0.032856, -0.074686, -0.081315,  0.005669,  0.190271,  0.360611,  0.327682, -0.009253, -0.501279, -0.792379, -0.600846,  0.000000,  0.600846,  0.792379,  0.501278,  0.009253, -0.327682, -0.360611, -0.190271, -0.005669,  0.081315,  0.074686,  0.032856,
-0.035372, -0.080406, -0.087543,  0.006103,  0.204843,  0.388228,  0.352778, -0.009962, -0.539669, -0.853064, -0.646862,  0.000000,  0.646862,  0.853064,  0.539669,  0.009962, -0.352778, -0.388228, -0.204843, -0.006103,  0.087543,  0.080406,  0.035372,
-0.036253, -0.082409, -0.089723,  0.006255,  0.209944,  0.397896,  0.361563, -0.010210, -0.553109, -0.874309, -0.662971, -0.000000,  0.662971,  0.874309,  0.553109,  0.010210, -0.361563, -0.397896, -0.209944, -0.006255,  0.089723,  0.082409,  0.036253,
-0.035372, -0.080406, -0.087543,  0.006103,  0.204843,  0.388228,  0.352778, -0.009962, -0.539669, -0.853065, -0.646862, -0.000000,  0.646862,  0.853065,  0.539669,  0.009962, -0.352778, -0.388228, -0.204843, -0.006103,  0.087543,  0.080406,  0.035372,
-0.032856, -0.074686, -0.081315,  0.005669,  0.190271,  0.360611,  0.327682, -0.009253, -0.501279, -0.792379, -0.600846, -0.000000,  0.600846,  0.792379,  0.501279,  0.009253, -0.327682, -0.360611, -0.190271, -0.005669,  0.081315,  0.074686,  0.032856,
-0.029053, -0.066043, -0.071905,  0.005013,  0.168251,  0.318878,  0.289760, -0.008182, -0.443266, -0.700679, -0.531311, -0.000000,  0.531311,  0.700679,  0.443266,  0.008182, -0.289760, -0.318878, -0.168251, -0.005013,  0.071905,  0.066043,  0.029053,
-0.024458, -0.055596, -0.060531,  0.004220,  0.141637,  0.268438,  0.243926, -0.006888, -0.373151, -0.589846, -0.447268, -0.000000,  0.447268,  0.589846,  0.373151,  0.006888, -0.243926, -0.268438, -0.141637, -0.004220,  0.060531,  0.055596,  0.024458,
-0.019601, -0.044555, -0.048510,  0.003382,  0.113509,  0.215128,  0.195484, -0.005520, -0.299046, -0.472707, -0.358444, -0.000000,  0.358444,  0.472707,  0.299046,  0.005520, -0.195484, -0.215128, -0.113509, -0.003382,  0.048510,  0.044555,  0.019601,
-0.014954, -0.033993, -0.037010,  0.002580,  0.086600,  0.164129,  0.149142, -0.004211, -0.228153, -0.360645, -0.273470, -0.000000,  0.273470,  0.360645,  0.228153,  0.004211, -0.149142, -0.164129, -0.086600, -0.002580,  0.037010,  0.033993,  0.014954,
-0.010861, -0.024689, -0.026881,  0.001874,  0.062898,  0.119208,  0.108323, -0.003059, -0.165709, -0.261940, -0.198624, -0.000000,  0.198624,  0.261940,  0.165709,  0.003059, -0.108323, -0.119208, -0.062899, -0.001874,  0.026881,  0.024689,  0.010861,
-0.007510, -0.017071, -0.018586,  0.001296,  0.043491,  0.082425,  0.074899, -0.002115, -0.114578, -0.181116, -0.137336, -0.000000,  0.137336,  0.181116,  0.114578,  0.002115, -0.074899, -0.082425, -0.043491, -0.001296,  0.018586,  0.017071,  0.007510,
-0.004943, -0.011237, -0.012234,  0.000853,  0.028627,  0.054256,  0.049302, -0.001392, -0.075421, -0.119219, -0.090401, -0.000000,  0.090401,  0.119219,  0.075421,  0.001392, -0.049302, -0.054256, -0.028627, -0.000853,  0.012234,  0.011237,  0.004943,
-0.003098, -0.007042, -0.007667,  0.000534,  0.017939,  0.033999,  0.030895, -0.000872, -0.047262, -0.074708, -0.056649, -0.000000,  0.056649,  0.074708,  0.047262,  0.000872, -0.030895, -0.033999, -0.017939, -0.000534,  0.007667,  0.007042,  0.003098,
-0.001848, -0.004201, -0.004574,  0.000319,  0.010702,  0.020283,  0.018431, -0.000520, -0.028195, -0.044568, -0.033795, -0.000000,  0.033795,  0.044568,  0.028195,  0.000520, -0.018431, -0.020283, -0.010702, -0.000319,  0.004574,  0.004201,  0.001848
};


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