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);