Home > track > compute_radial_mean.m

compute_radial_mean

PURPOSE ^

SYNOPSIS ^

function radial_mean = compute_radial_mean(im_sub,rm_k)

DESCRIPTION ^

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function radial_mean = compute_radial_mean(im_sub,rm_k)
0002 
0003 rm_layer = max(rm_k(:));
0004 
0005  % compute radial means of 3 layers
0006 
0007 for i = 0:rm_layer 
0008     % current calculating region
0009     region_cur = rm_k == i;
0010     % number of pixel in that region
0011     sum_pix = sum ( region_cur(:) );
0012     try
0013     % sum of intensity of all the pixel in that region
0014     sum_i = sum ( im_sub (region_cur) );
0015     catch ME
0016         pause();
0017     end
0018     % the mean of the intensity of that region
0019     mean_rad = sum_i / sum_pix;                
0020     % radial mean vector for all 3 layers
0021     radial_mean(i+1) = mean_rad;
0022 end
0023 
0024 % normalize the vector (sum will = 1)
0025 radial_mean = radial_mean / sum(radial_mean);

Generated on Thu 17-Mar-2011 14:45:51 by m2html © 2005