MATLAB: 1D median filtering

signal processing

how to apply 1D median filter with 100ms window length?

Best Answer

Assuming your data are stored as a row vector D, and the sampling interval is T ms, and you have the Image Processing Toolbox, you could do this:

windowlength = 100;
windowsize = ceil(windowlength/T);
D_filtered = medfilt2(D, [1 windowsize]);

If the data are stored as a column vector, use

D_filtered = medfilt2(D, [windowsize 1]);

instead of the final line.

Related Solutions

MATLAB: How to apply a 1-D median filter and a moving average filter to this set of data

Try medfilt1():

filteredData = medfilt1(data, windowWidth);

MATLAB: Median filter

To understand the medial filter, maybe it is better if we start in 1D.

Let's suppose the following array:

a = [3 2 1 6 5];

Also, let's suppose you define a window size of n = 3.

You place yourself in the first value of the array, that is 3. Since your window size is 3, you take the sample you are, one sample to the left and one sample to the right. Since this is the first value, there is no sample to the left, so we will just consider a zero:

0 3 2

From this window, you sort the values from minimum to maximum and get the median (the value in the middle)

0 2 3

In this case, the value in the middle is 2. This is the value you write in the filtered vector a_

a_ = [2 2 1 6 5];

Then you go to the next value of a, that is a 2. You take the left sample (which has already been modified) and the right sample

2 2 1

Your sort them and choose the median value (in this case, a 2 again)

1 2 2

So your filtered vector a_ is now

a_ = [2 2 1 6 5];

An so on. In the end, you will have a filtered vector a_

a_ = [2 2 2 5 5];

A median filter in images works the same way, only in 2D. So you take not only the values (pixels) that are left or right, but all the values that surround the sample (pixel) you are in.

In MATLAB, check medfilt1 and medfilt2 ;) .

doc medfilt1;
doc medfilt2;

Anyway, take a look at this useful tutorial http://www.librow.com/articles/article-1 by Chernenko. Even though it is not coded in MATLAB, the explanation is very, very clear.

Keep coding!

Best Answer

Related Solutions

MATLAB: How to apply a 1-D median filter and a moving average filter to this set of data

MATLAB: Median filter

Related Question