рдореЗрдбрд┐рдпрди рдлрд╝рд┐рд▓реНрдЯрд░ рд╕рд╛рдорд╛рдиреНрдпрдХрд░рдг

рдЕрдореВрд░реНрдд


рдпрд╣ рд▓реЗрдЦ рдПрдХ рдЕрджреНрд╡рд┐рддреАрдп рдлрд┐рд▓реНрдЯрд░ рдХреЗ рдмрд╛рд░реЗ рдореЗрдВ рдмрд╛рдд рдХрд░рддрд╛ рд╣реИ, рдЬрд┐рд╕рдХреЗ рдмрд╛рд░реЗ рдореЗрдВ рдПрдХ рд▓реЗрдЦ 1990 рдореЗрдВ рд╕рд╛рдордиреЗ рдЖрдпрд╛ рдерд╛: рдорд╛рд╕реНрд▓реЛрд╡ рдП.рдПрдо., рд╕рд░реНрдЧреЗрд╡ рд╡реА.рд╡реА. рд░реИрдВрдХ рд╕рд┐рдЧреНрдирд▓ рдкреНрд░реЛрд╕реЗрд╕рд┐рдВрдЧ // рдХрдВрдкреНрдпреВрдЯрд░ рдСрдкреНрдЯрд┐рдХреНрд╕ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░ рдПрдХ рд░реЗрдЦреАрдп рд╡рд┐рдХреГрдд рдкреНрд░рдгрд╛рд▓реА рдХреА рдкрд╣рдЪрд╛рдиред - рдПрдоред, 1990ред - рдЕрдВрдХ 6ред - S.97-102ред рдЗрд╕ рдПрд▓реНрдЧреЛрд░рд┐рдереНрдо рдХреЛ "рд░реИрдВрдХ рдкреНрд░реЛрд╕реЗрд╕рд┐рдВрдЧ рдПрд▓реНрдЧреЛрд░рд┐рдердо" рдХрд╣рд╛ рдЬрд╛рддрд╛ рд╣реИ рдФрд░ рд╡рд╛рд╕реНрддрд╡ рдореЗрдВ рдорд╛рдзреНрдпрд┐рдХрд╛ рдлрд┐рд▓реНрдЯрд░ рдХрд╛ рдПрдХ рд╕рд╛рдорд╛рдиреНрдпреАрдХрд░рдг рд╣реИред
рдЗрд╕ рдлрд╝рд┐рд▓реНрдЯрд░ рдХрд╛ рдЙрдкрдпреЛрдЧ рджреЛ рдорд╛рдорд▓реЛрдВ рдореЗрдВ рдЙрдЪрд┐рдд рд╣реИ - рд╢реЛрд░ рдХреЛ рджрдмрд╛рдиреЗ рдФрд░ рдЧреНрд░реАрд╕ рдХреЛ рдХрдо рдХрд░рдиреЗ рдХреЗ рд▓рд┐рдПред
рдЫрд╡рд┐
рдЪрд┐рддреНрд░рд╛ 1 - рдореВрд▓ рдЫрд╡рд┐, 2 - рдзреБрдВрдзрд▓реА рдФрд░ рд╢реЛрд░ рдирдордХред


рд╕рд┐рдЧреНрдирд▓ рдкреНрд░реЛрд╕реЗрд╕рд┐рдВрдЧ рдПрд▓реНрдЧреЛрд░рд┐рдердо - рдЪрд░рдо рдлрд╝рд┐рд▓реНрдЯрд░


рдПрд▓реНрдЧреЛрд░рд┐рдереНрдо рдореЗрдВ рдПрдХ рдирдИ рдЫрд╡рд┐ рдореЗрдВ рдкреНрд░рд╕рдВрд╕реНрдХрд░рдг рдХреЗ рдкрд░рд┐рдгрд╛рдо рдХреЛ рд░рд┐рдХреЙрд░реНрдб рдХрд░рдиреЗ рдХреЗ рд╕рд╛рде рдПрдХ рд╕реНрдерд╛рдиреАрдп рд╡рд┐рдВрдбреЛ рджреНрд╡рд╛рд░рд╛ рдЫрд╡рд┐ рдХреЛ рд╕рдВрд╕рд╛рдзрд┐рдд рдХрд░рдирд╛ рд╢рд╛рдорд┐рд▓ рд╣реИ:
  1. рд╣рдо рдирд┐рд░реНрджреЗрд╢рд╛рдВрдХ (I, J) рдХреЗ рд╕рд╛рде рдЫрд╡рд┐ рдмрд┐рдВрджреБ рдкрд░ рд╣реИрдВ - рд╡рд░реНрддрдорд╛рди рдЧрдгрдирд╛
  2. рд╡рд░реНрддрдорд╛рди рд╕рдВрджрд░реНрдн рдХреЗ рдЖрд╕рдкрд╛рд╕, рд╣рдо рдЖрдХрд╛рд░ NxN рдХреЗ рд╕реНрдерд╛рдиреАрдп рдкрдбрд╝реЛрд╕ рдкрд░ рд╡рд┐рдЪрд╛рд░ рдХрд░рддреЗ рд╣реИрдВ
  3. рдПрдХ рд╕реНрдерд╛рдиреАрдп рдкрдбрд╝реЛрд╕ рдореЗрдВ рдЕрдВрдХреЛрдВ рдХреЗ рдореВрд▓реНрдпреЛрдВ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддреЗ рд╣реБрдП, рдПрдХ рдкрд░рд┐рд╡рд░реНрддрдирд╢реАрд▓ рд╢реНрд░реГрдВрдЦрд▓рд╛ рдХрд╛ рдирд┐рд░реНрдорд╛рдг рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ, рдЬрд┐рд╕реЗ рд╣рдо рдкреА рджреНрд╡рд╛рд░рд╛ рджрд░реНрд╢рд╛рддреЗ рд╣реИрдВред рдЗрд╕ рдкрдВрдХреНрддрд┐ рдХрд╛ рдЖрдХрд╛рд░ N * N рд╣реИред
  4. рдкрд░рд┐рдгрд╛рдореА рдЫрд╡рд┐ рдореЗрдВ, рд╡рд░реНрддрдорд╛рди рдЧрдгрдирд╛ рдирд┐рдореНрдирд▓рд┐рдЦрд┐рдд рдирд┐рдпрдо рдкрд░ рдЪрд▓рддреА рд╣реИ:

рдЫрд╡рд┐
рдЬрд╣рд╛рдВ k рдПрд▓реНрдЧреЛрд░рд┐рдердо рдХрд╛ рдкреИрд░рд╛рдореАрдЯрд░ рд╣реИ, N рдПрдХ рд╡рд┐рд╖рдо рд╕рдВрдЦреНрдпрд╛ рд╣реИ, Im1 рдореВрд▓ рдЫрд╡рд┐ рд╣реИ, Im2 рдкрд░рд┐рдгрд╛рдореА рдЫрд╡рд┐ рд╣реИред
рдпрджрд┐ k = (N ^ 2 + 1) / 2 - рдЕрд░реНрдерд╛рдд, рд╡реИрдЪрд╛рд░рд┐рдХ рд╢реНрд░реГрдВрдЦрд▓рд╛ рдХрд╛ рдХреЗрдВрджреНрд░ - рдпрд╣ рдлрд╝рд┐рд▓реНрдЯрд░ рдПрдХ рдЬреНрдЮрд╛рдд рдорд╛рдзреНрдпрд┐рдХрд╛ рдлрд╝рд┐рд▓реНрдЯрд░ рдмрди рдЬрд╛рддрд╛ рд╣реИред рднрд╡рд┐рд╖реНрдп рдореЗрдВ, рдЗрд╕ рдкреИрд░рд╛рдореАрдЯрд░ рдХреЛ рдЗрдВрдбреЗрдВрдЯреЗрд╢рди рдХрд╣рд╛ рдЬрд╛рдПрдЧрд╛ред

рдЪрд░рдо рдлрд╝рд┐рд▓реНрдЯрд░ рдЧреБрдг


рдЗрд╕ рдлрд╝рд┐рд▓реНрдЯрд░ рдХреЗ рдЧреБрдг рд╡реНрдпрд╡рд╣рд╛рд░ рдореЗрдВ рдмрд╣реБрдд рдЙрдкрдпреЛрдЧреА рд╣реЛрддреЗ рд╣реИрдВ, рдХреНрдпреЛрдВрдХрд┐ рдлрд╝рд┐рд▓реНрдЯрд░ рдЖрдкрдХреЛ рди рдХреЗрд╡рд▓ рд╢реЛрд░ рдХреА рднрд░рдкрд╛рдИ рдХрд░рдиреЗ рдХреА рдЕрдиреБрдорддрд┐ рджреЗрддрд╛ рд╣реИ, рдмрд▓реНрдХрд┐ (рдЖрдВрд╢рд┐рдХ рд░реВрдк рд╕реЗ) рдЫрд╡рд┐ рдзрдмреНрдмрд╛ рдХреЗ рдкреНрд░рднрд╛рд╡ рдХреЛ рд╕рдорд╛рдкреНрдд рдХрд░рддрд╛ рд╣реИред K = (N ^ 2 + 1) / 2 рдкрд░ рдЗрд╕ рдлрд╝рд┐рд▓реНрдЯрд░ рдХрд╛ рд╕реАрдорд┐рдд рдорд╛рдорд▓рд╛ рдПрдХ рдордзреНрдпрдпреБрдЧреАрди рдлрд╝рд┐рд▓реНрдЯрд░ рд╣реИ рдЬреЛ рдХреЗрд╡рд▓ рд╢реЛрд░ рдХреЛ рд╕рдорд╛рдкреНрдд рдХрд░рддрд╛ рд╣реИ, рд▓реЗрдХрд┐рди рд╕реАрдорд╛ рдХреЛ рдирд╣реАрдВ рдЫреВрддрд╛ рд╣реИ, рдФрд░ рдпрджрд┐ рдЫрд╡рд┐ рдзреБрдВрдзрд▓реА рд╣реИ, рддреЛ рдзрдмреНрдмрд╛ рдмрдирд╛ рд░рд╣реЗрдЧрд╛ред

K <(N ^ 2 + 1) / 2 рдкрд░, рд╢реЛрд░ рдереЛрдбрд╝рд╛ рдЦрд░рд╛рдм рдлрд╝рд┐рд▓реНрдЯрд░ рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИ, рд▓реЗрдХрд┐рди рдЫрд╡рд┐ рдХреА рддреАрд╡реНрд░рддрд╛ рдмрдврд╝ рдЬрд╛рддреА рд╣реИ, рдФрд░ k = 0 рдкрд░, рд╢реЛрд░ рдмрд┐рд▓реНрдХреБрд▓ рднреА рдлрд╝рд┐рд▓реНрдЯрд░реНрдб рдирд╣реАрдВ рд╣реЛрддрд╛ рд╣реИ, рд▓реЗрдХрд┐рди рд╕рдмрд╕реЗ рдордЬрдмреВрдд рддрд░реАрдХреЗ рд╕реЗ рдЧреНрд░реАрд╕ рд╕рдорд╛рдкреНрдд рд╣реЛ рдЬрд╛рдПрдЧрд╛ред

рдЗрд╕ рдПрд▓реНрдЧреЛрд░рд┐рдереНрдо рдХреЗ рдПрдХ рдкреНрд░рд╛рдердорд┐рдХ рдХрд╛рд░реНрдпрд╛рдиреНрд╡рдпрди рдХреЗ рд╕рд╛рде рдкрд╛рдардХреЛрдВ рдХреЛ рдкрд░реЗрд╢рд╛рди рдирд╣реАрдВ рдХрд░рдиреЗ рдХреЗ рд▓рд┐рдП, рдореИрдВ рдпрд╣рд╛рдВ рдордЯрд▓рдм рдореЗрдВ рдХреЛрдб рдкреНрд░рджрд╛рди рдХрд░реВрдВрдЧрд╛ред
рдкрд░реАрдХреНрд╖рдг рд╕реНрдХреНрд░рд┐рдкреНрдЯ рдЫрд╡рд┐ рдХреЛ рдкрдврд╝реЗрдЧреА, рдЙрд╕реЗ рдзреБрдВрдзрд▓рд╛ рдХрд░реЗрдЧреА, рдФрд░ рд╢реЛрд░ рдЬреЛрдбрд╝реЗрдЧреАред рдлрд┐рд░ рдЗрди рдЪрд░рдо рдлрд╝рд┐рд▓реНрдЯрд░ рджреНрд╡рд╛рд░рд╛ рдЫрд╡рд┐ рдХреЛ рдкреБрдирд░реНрд╕реНрдерд╛рдкрд┐рдд рдХрд┐рдпрд╛ рдЬрд╛рддрд╛ рд╣реИред
I1 = imread( 'coins.png' ); h = ones(3,3) / 9; I2 = imfilter(I1,h) ; I3 = imnoise(I2, 'salt & pepper' ,0.02); I4 = im2_rang_filter (I3, 1, 2); figure; imagesc(I1); colormap gray; figure; imagesc(I3); colormap gray; figure; imagesc(I4) colormap gray; * This source code was highlighted with Source Code Highlighter .
  1. I1 = imread( 'coins.png' ); h = ones(3,3) / 9; I2 = imfilter(I1,h) ; I3 = imnoise(I2, 'salt & pepper' ,0.02); I4 = im2_rang_filter (I3, 1, 2); figure; imagesc(I1); colormap gray; figure; imagesc(I3); colormap gray; figure; imagesc(I4) colormap gray; * This source code was highlighted with Source Code Highlighter .
  2. I1 = imread( 'coins.png' ); h = ones(3,3) / 9; I2 = imfilter(I1,h) ; I3 = imnoise(I2, 'salt & pepper' ,0.02); I4 = im2_rang_filter (I3, 1, 2); figure; imagesc(I1); colormap gray; figure; imagesc(I3); colormap gray; figure; imagesc(I4) colormap gray; * This source code was highlighted with Source Code Highlighter .
  3. I1 = imread( 'coins.png' ); h = ones(3,3) / 9; I2 = imfilter(I1,h) ; I3 = imnoise(I2, 'salt & pepper' ,0.02); I4 = im2_rang_filter (I3, 1, 2); figure; imagesc(I1); colormap gray; figure; imagesc(I3); colormap gray; figure; imagesc(I4) colormap gray; * This source code was highlighted with Source Code Highlighter .
  4. I1 = imread( 'coins.png' ); h = ones(3,3) / 9; I2 = imfilter(I1,h) ; I3 = imnoise(I2, 'salt & pepper' ,0.02); I4 = im2_rang_filter (I3, 1, 2); figure; imagesc(I1); colormap gray; figure; imagesc(I3); colormap gray; figure; imagesc(I4) colormap gray; * This source code was highlighted with Source Code Highlighter .
  5. I1 = imread( 'coins.png' ); h = ones(3,3) / 9; I2 = imfilter(I1,h) ; I3 = imnoise(I2, 'salt & pepper' ,0.02); I4 = im2_rang_filter (I3, 1, 2); figure; imagesc(I1); colormap gray; figure; imagesc(I3); colormap gray; figure; imagesc(I4) colormap gray; * This source code was highlighted with Source Code Highlighter .
  6. I1 = imread( 'coins.png' ); h = ones(3,3) / 9; I2 = imfilter(I1,h) ; I3 = imnoise(I2, 'salt & pepper' ,0.02); I4 = im2_rang_filter (I3, 1, 2); figure; imagesc(I1); colormap gray; figure; imagesc(I3); colormap gray; figure; imagesc(I4) colormap gray; * This source code was highlighted with Source Code Highlighter .
  7. I1 = imread( 'coins.png' ); h = ones(3,3) / 9; I2 = imfilter(I1,h) ; I3 = imnoise(I2, 'salt & pepper' ,0.02); I4 = im2_rang_filter (I3, 1, 2); figure; imagesc(I1); colormap gray; figure; imagesc(I3); colormap gray; figure; imagesc(I4) colormap gray; * This source code was highlighted with Source Code Highlighter .
  8. I1 = imread( 'coins.png' ); h = ones(3,3) / 9; I2 = imfilter(I1,h) ; I3 = imnoise(I2, 'salt & pepper' ,0.02); I4 = im2_rang_filter (I3, 1, 2); figure; imagesc(I1); colormap gray; figure; imagesc(I3); colormap gray; figure; imagesc(I4) colormap gray; * This source code was highlighted with Source Code Highlighter .
  9. I1 = imread( 'coins.png' ); h = ones(3,3) / 9; I2 = imfilter(I1,h) ; I3 = imnoise(I2, 'salt & pepper' ,0.02); I4 = im2_rang_filter (I3, 1, 2); figure; imagesc(I1); colormap gray; figure; imagesc(I3); colormap gray; figure; imagesc(I4) colormap gray; * This source code was highlighted with Source Code Highlighter .
  10. I1 = imread( 'coins.png' ); h = ones(3,3) / 9; I2 = imfilter(I1,h) ; I3 = imnoise(I2, 'salt & pepper' ,0.02); I4 = im2_rang_filter (I3, 1, 2); figure; imagesc(I1); colormap gray; figure; imagesc(I3); colormap gray; figure; imagesc(I4) colormap gray; * This source code was highlighted with Source Code Highlighter .
  11. I1 = imread( 'coins.png' ); h = ones(3,3) / 9; I2 = imfilter(I1,h) ; I3 = imnoise(I2, 'salt & pepper' ,0.02); I4 = im2_rang_filter (I3, 1, 2); figure; imagesc(I1); colormap gray; figure; imagesc(I3); colormap gray; figure; imagesc(I4) colormap gray; * This source code was highlighted with Source Code Highlighter .
I1 = imread( 'coins.png' ); h = ones(3,3) / 9; I2 = imfilter(I1,h) ; I3 = imnoise(I2, 'salt & pepper' ,0.02); I4 = im2_rang_filter (I3, 1, 2); figure; imagesc(I1); colormap gray; figure; imagesc(I3); colormap gray; figure; imagesc(I4) colormap gray; * This source code was highlighted with Source Code Highlighter .


рдлрд╝рд┐рд▓реНрдЯрд░ рдлрд╝рдВрдХреНрд╢рди рдХреЛрдб im2_rang_filter :
function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  1. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  2. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  3. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  4. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  5. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  6. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  7. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  8. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  9. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  10. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  11. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  12. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  13. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  14. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  15. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
  16. function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .
function outImage= im2_rang_filter (aImage, aHalfWindowSize, aOtsup) [ver,hor] = size(aImage); wsize = (aHalfWindowSize*2+1)^2; result = zeros(ver,hor); for i = aHalfWindowSize+1 : (ver - aHalfWindowSize) for j = aHalfWindowSize+1 : (hor - aHalfWindowSize) wind = aImage((i-aHalfWindowSize) : (i + aHalfWindowSize), (j-aHalfWindowSize) : (j + aHalfWindowSize)); vec = reshape(wind,1,wsize); vec = sort(vec); if (abs(vec(aOtsup+1) - aImage(i,j)) < abs(vec(wsize - aOtsup) - aImage(i,j)) ) result(i,j) = vec(aOtsup+1); else result(i,j) = vec(wsize - aOtsup); end ; end ; end ; outImage = result; * This source code was highlighted with Source Code Highlighter .

рд╡рд┐рдЬреНрдЮрд╛рди рдкреНрд░рдпреЛрдЧ


рдиреАрдЪреЗ рджрд┐рдП рдЧрдП рдкреИрд░рд╛рдореАрдЯрд░ k = 1.2.5 рдХреЛ рдЕрд▓рдЧ рдХрд░рдиреЗ рдкрд░ 3 рдмрд╛рдп 3 рд╡рд┐рдВрдбреЛ рдХреЗ рд╕рд╛рде рдлрд╝рд┐рд▓реНрдЯрд░рд┐рдВрдЧ рдХреЗ рдкрд░рд┐рдгрд╛рдо рд╣реИрдВред рдкрд╣рд▓реА рддрд╕реНрд╡реАрд░ рдореВрд▓ рдЫрд╡рд┐ рд╣реИ, рджреВрд╕рд░реА рддрд╕реНрд╡реАрд░ рд╡рд┐рдХреГрдд рд╣реИ, рдлрд┐рд░ k = 1, 2, 5 рдХреЗ рд▓рд┐рдП рдлрд╝рд┐рд▓реНрдЯрд░реНрдб рдлрд╝рд┐рд▓реНрдЯрд░ рдПрдХ рджреВрд╕рд░реЗ рдХрд╛ рдЕрдиреБрд╕рд░рдг рдХрд░рддреЗ рд╣реИрдВ:
рдкреНрд░рд╛рд░рдВрднрд┐рдХ:
рдЫрд╡рд┐
рд╡рд┐рд░реВрдкрдг:
рдЫрд╡рд┐
K = 1 рдкрд░ рдкреБрдирд░реНрдирд┐рд░реНрдорд╛рдг рдХрд┐рдпрд╛ рдЧрдпрд╛:
рдЫрд╡рд┐
K = 2 рдкрд░ рдкреБрдирд░реНрдкреНрд░рд╛рдкреНрдд рдХрд┐рдпрд╛ рдЧрдпрд╛:
рдЫрд╡рд┐
K = 5 рдкрд░ рдкреБрдирд░реНрдирд┐рд░реНрдорд╛рдг рдХрд┐рдпрд╛ рдЧрдпрд╛:
рдЫрд╡рд┐

рдореИрдВ рд╡рд╛рд╕реНрддрд╡ рдореЗрдВ рдЗрд╕ рдлрд┐рд▓реНрдЯрд░ рдХреЛ рдкрд╕рдВрдж рдХрд░рддрд╛ рд╣реВрдВ рдФрд░ рдореИрдВ рд╕рдХреНрд░рд┐рдп рд░реВрдк рд╕реЗ рдЗрд╕рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддрд╛ рд╣реВрдВред рдкрд╛рдардХ рдЗрд╕ рдлрд┐рд▓реНрдЯрд░ рдореЗрдВ рд░реБрдЪрд┐ рд░рдЦрддреЗ рд╣реИрдВред

Source: https://habr.com/ru/post/In114551/


All Articles