Difference between revisions of "Gaussian and sieve filters"

From BanghamLab
Jump to navigation Jump to search
Line 7: Line 7:
 
If, on the other hand the one behind is the value that has just been computed using the previous window then the recursive median is computed. It is similar to an 'irr' filter as opposed to a 'fir' filter. Stage (scale) 1 window of 3. The filter removes extrema of length (scale) 1. Stage 2 also has a window of 3 but is computed on a basis of 5 values. This filter removes extrema of scale 2. And so forth. Each stage simplifies the signal until, at the bottom, it is entirely removed. It turns out that extrema at each scale are removed in a nice regular way <span style="color:red;">it does preserve [http://en.wikipedia.org/wiki/Scale_space scale-space] and that is '''''GOOD'''''</span>.
 
If, on the other hand the one behind is the value that has just been computed using the previous window then the recursive median is computed. It is similar to an 'irr' filter as opposed to a 'fir' filter. Stage (scale) 1 window of 3. The filter removes extrema of length (scale) 1. Stage 2 also has a window of 3 but is computed on a basis of 5 values. This filter removes extrema of scale 2. And so forth. Each stage simplifies the signal until, at the bottom, it is entirely removed. It turns out that extrema at each scale are removed in a nice regular way <span style="color:red;">it does preserve [http://en.wikipedia.org/wiki/Scale_space scale-space] and that is '''''GOOD'''''</span>.
 
|}
 
|}
The output depends on the computation kernel the recursive median is one. Others include connected set openings ('o'), closings ('c') and composites of these 'M' ('o' followed by 'c' at each stage) or 'N' ('c' followed by 'o'). These last two yield outputs that are almost indistinguishable from the recursive median. (I have just realised how confusing it is referring to the recursive mean as an 'm' sieve, I shall start using the inebriation 'v' as much as possible instead.)<br><br>
+
The output depends on the computation kernel the recursive median is one. Others include connected set openings ('o'), closings ('c') and composites of these 'M' ('o' followed by 'c' at each stage) or 'N' ('c' followed by 'o'). These last two yield outputs that are almost indistinguishable from the recursive median. (I have just realised how confusing it is referring to the recursive mean as an 'm' sieve, I shall start using the inebriation 'v' as much as possible instead.) This is dull stuff.<br><br>
 
+
#FFBF00
This is dull stuff compared to something that has been thought about much less, let alone exploited - <span style="color:red;">a PhD project waiting to happen.</span>  Sieves (and the recursive median filter is what I call one of the sieves) <span style="color:red;">are idempotent. In other words having made one pass through the data at any particular scale, making another pass through the result changes nothing.</span> <br><br>
+
This is dull stuff <span style="color:#C46210;"> compared to something that has been thought about much less, let alone exploited</span>- <span style="color:red;">a PhD project waiting to happen.</span>  Sieves (and the recursive median filter is what I call one of the sieves) <span style="color:red;">are idempotent. In other words having made one pass through the data at any particular scale, making another pass through the result changes nothing.</span> <br><br>
This is not like a linear (diffusion) filter where repeated passes at the same scale simply smooth the signal away. It means that one could build entirely new filtering schema, here are some [[new filtering schema| that are (for me) last minute.]] They begin to look a little biological.
+
This is not like a linear (diffusion) filter where repeated passes at the same scale simply smooth the signal away. It means that one <span style="color:#C46210;">could build entirely new filtering schema</span>, here are some [[new filtering schema| ideas.]] They begin to look a little biological.

Revision as of 20:47, 6 August 2014

Siv1-gaussiansieve.png</wikiflv>


OK its a another crap Figure. I spent ages getting the trace above each new one to show the moving window used to compute the output of each stage.

Left panel: Gaussian filter. Stage (scale) <math>sigma</math>=1. Stage 2 <math>sigma</math>=2 on the previous stage. And so forth. Each stage simplifies the signal until, at the bottom, it fades away. It does preserve scale-space and that is GOOD.

Right panel: recursive median. Whatever the scale each computation requires three numbers, the middle one, one ahead and one behind. If the one behind comes directly from the previous stage the median is found.

If, on the other hand the one behind is the value that has just been computed using the previous window then the recursive median is computed. It is similar to an 'irr' filter as opposed to a 'fir' filter. Stage (scale) 1 window of 3. The filter removes extrema of length (scale) 1. Stage 2 also has a window of 3 but is computed on a basis of 5 values. This filter removes extrema of scale 2. And so forth. Each stage simplifies the signal until, at the bottom, it is entirely removed. It turns out that extrema at each scale are removed in a nice regular way it does preserve scale-space and that is GOOD.

The output depends on the computation kernel the recursive median is one. Others include connected set openings ('o'), closings ('c') and composites of these 'M' ('o' followed by 'c' at each stage) or 'N' ('c' followed by 'o'). These last two yield outputs that are almost indistinguishable from the recursive median. (I have just realised how confusing it is referring to the recursive mean as an 'm' sieve, I shall start using the inebriation 'v' as much as possible instead.) This is dull stuff.

  1. FFBF00

This is dull stuff compared to something that has been thought about much less, let alone exploited- a PhD project waiting to happen. Sieves (and the recursive median filter is what I call one of the sieves) are idempotent. In other words having made one pass through the data at any particular scale, making another pass through the result changes nothing.

This is not like a linear (diffusion) filter where repeated passes at the same scale simply smooth the signal away. It means that one could build entirely new filtering schema, here are some ideas. They begin to look a little biological.