if, for loop, median

2011 4 11
1 个回答
回答已采纳
1 次查看（30 天）

0 个投票

hi, I am new to matlab programming. I wish to write a program involving image processing, which serves the following purpose: consider a 3x3 window in a normalised 512x512 image. Compare each value of the pixel with 0.3. If the value of image pixel is less than 0.3, sort the vector, and find the median of all the pixels with value less than 0.3. And if value is more than 0.3, neglect that pixel value, and find median of all remaining pixels which have values less than 0.3. All of this should result in a new image. Thanks and Regards Manpreet Kapoor

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

 采纳的回答

Walter Roberson
Walter Roberson 2011-4-11

1 个投票

Is this 3 x 3 window to be slid overlapping over the whole image, or is the image to be broken up in to 3 x 3 blocks? If it is to be broken up in to 3 x 3 blocks, there is a problem because 512 is not divisible by 3.
You can use sorting to find the median if you want, but there is a median() function. For example, if P are the pixels in the window, then
median(P(P<0.3))

10 个评论
显示 8更早的评论 隐藏 8更早的评论

Manpreet Kapoor
Manpreet Kapoor 2011-4-11
walter,
please dont concentrate on the size of the image, wat I want is the program to implement the logic given. I just want to see the way the functions are designed....
Walter Roberson
Walter Roberson 2011-4-11
The logic for overlapping blocks is different than for non-overlapping blocks.
If you are using non-overlapping blocks, see
http://www.mathworks.com/matlabcentral/answers/5163-image-block-block-dividing
Manpreet Kapoor
Manpreet Kapoor 2011-4-11
walter,
it is overlapping window, coz the median needs to be found for every pixel window, thus forming a new image. I know the function to be used is nlfilter, wat i want to know is the sequence of commands to perform the above said function. I basically need to learn the usage of functions in a program.
Walter Roberson
Walter Roberson 2011-4-11
I'm not familiar with nlfilter, but it does appear suitable
B = nlfilter(A,[3,3],@(x) median(x(x<0.3)));
The rest of the problem statement is noise, two equivalent conditions.
Manpreet Kapoor
Manpreet Kapoor 2011-4-11
hi
this command does not work, even wen I defined both A and B.... Wat to do?
Walter Roberson
Walter Roberson 2011-4-11
Please be more specific about what you mean by "does not work". Do you have the image processing toolbox?
You might want to try the example given in the documentation,
http://www.mathworks.com/help/toolbox/images/ref/nlfilter.html
Manpreet Kapoor
Manpreet Kapoor 2011-4-11
walter,
I have already gone through the above document earlier, wat I am facing trouble is in the the designing of 'fun' which performs the median of the pixel values less than 0.3.
I already mentioned I am new to matlab and programming watsoever, so cant understand the implementation of 'fun'. Please help.....
Walter Roberson
Walter Roberson 2011-4-11
Okay, so let's go back to you being more specific about the command not working. What exactly do you mean by that? If the cameraman example from the documentation worked for you, then my very small variation on it should work. The cameraman example had median(x(:)) and I have median(x(x<0.3)) . For an array, x(x<0.3) will return a vector suitable for input to median, and that vector will consist only of the pixels with value less than 0.3 .
Manpreet Kapoor
Manpreet Kapoor 2011-4-11
Sorry for the misunderstanding earlier, the logic u gave is absolutely correct, but the problem I have now is that nlfilter takes almost an hour to code one image, hence another command I found faster is colfilt.
can u offer any help in this....???
any help will be highly appreciated.
Thanks and Regards
Walter Roberson
Walter Roberson 2011-4-11
Sorry, I am not familiar with colfilt .
I would think it likely that blkproc() or blockproc() would be much faster than nlfilter.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Neighborhood and Block Processing 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by