median filter on a region of interest

13 次查看（过去 30 天）
David Levi
David Levi 2022-4-26
编辑： DGM 2022-4-26
In image processing toolbox, is there a way to apply median filter only on specific region of interest (and not on the whole image)?.
If there is, can you give me an example written in code?
Thanks you all !
  2 个评论
Matt J
Matt J 2022-4-26
Would it be a rectangular region? If not, how would the filter handle pixel neighborhoods that straddle the ROI boundary?
David Levi
David Levi 2022-4-26
The roi is not a rectangular region.
Like you suggested , ignoring all pixels outside the ROI and set them to NaN might be a good idea.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（2 个）

Matt J
Matt J 2022-4-26
编辑：Matt J 2022-4-26
Apparently there's a way to interface with Python's median filter, which can ignore NaNs.
https://www.mathworks.com/matlabcentral/answers/598651-problem-incorporating-a-python-image-filter-in-matlab?s_tid=srchtitle
Presumably therefore, if you set all pixels outside the ROI to NaN, this would do what you want.
DGM
DGM 2022-4-26
编辑：DGM 2022-4-26
Ran in:
Just filter the whole image and then compose the two using a mask that defines the ROI.
A = imread('cameraman.tif');
% make a mask
imshow(A) % set up the axes
R = images.roi.Polygon(gca); % create ROI object
R.Position = [66 56;112 19;168 29;198 86;150 163;64 143;83 95];
mask = createMask(R); % create binary mask
% make a filtered copy
Amed = medfilt2(A,[5 5]);
% insert filtered region into original image
B = A;
B(mask) = Amed(mask);
imshow(B)
In this example, I used an ROI object to create a mask. You can use any method you want to create such a mask. The rest is just basic logical indexing.
  3 个评论
显示 1更早的评论
DGM
DGM 2022-4-26
编辑：DGM 2022-4-26
That is correct. A compromise may be that you could select a rectangular region which circumscribes the non-rectangular ROI (plus padding determined by the window size). You could then filter that region of the image and compose as before. The relative benefit depends how big the ROI is compared to the image. At least that would let you get away with using standard tools.
DGM
DGM 2022-4-26
编辑：DGM 2022-4-26
Ran in:
Actually, I didn't know this, but apparently you can manage to do this with roifilt2().
A = imread('cameraman.tif'); % same picture as before
mask = imread('mk.png'); % same mask as before
F = @(x) medfilt2(x,[5 5]);
B = roifilt2(A,mask,F);
imshow(B)
I don't know which is faster, but doing it by composition is at least more flexible.
EDIT: at least with cursory testing, using roifilt2() on such an image/mask isn't any faster than the logical composition example. Upon opening roifilt2(), it appears to be doing exactly what I described above -- filtering a circumscribed rectangular region and then doing the composition via simple logical indexing. With that in mind, there may be relative speed benefits if the ROI is much smaller than it is in these examples.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 ROI-Based Processing 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by