0001 function imConfidence = radial_mean(im,kernelSize)
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020
0021
0022 rm_k_max = 2;
0023 kh = floor(kernelSize/2);
0024 temp = zeros (kernelSize,kernelSize);
0025 temp(kh+1,kh+1) = 1;
0026 rm_k = floor(bwdist(temp));
0027 max_dist = max ( rm_k(:) );
0028 f = double(rm_k_max+1) / double ( floor(max_dist)+1 );
0029 rm_k = floor( rm_k * f );
0030
0031
0032
0033 kernelInnerRing = double(rm_k == 0);
0034 kernelMiddleRing = double(rm_k == 1);
0035 kernelOutterRing = double(rm_k == 2);
0036
0037
0038
0039 im = double(im);
0040 r1 = imfilter(im,kernelInnerRing,'replicate');
0041 r2 = imfilter(im,kernelMiddleRing,'replicate');
0042 r3 = imfilter(im,kernelOutterRing,'replicate');
0043
0044
0045
0046 imConfidence = r1 ./ (r1 + r2 + r3);
0047