Toolbox ia870 | List of Figures | Fig. 6.24 | Fig. 6.29

Figure 6.26 - Regional maxima of synthetic noise image

Description

Image whit the three Gaussian peaks with some additive noise. This figure illustrates the regional maxima detection of synthetic noise image.

The matter is illustrated in Fig. 6.26: (a) synthetic image with added noise; (b) regional maxima appearing randomly over the image; (c) reconstructive opening by a small disk resulting in the regional maximum operator detecting just the three main domes; (c) opening by a larger disk resulting in selection of just the wider dome; (d) -maxima resulting in just the tallest dome; and (f) the top-left dome resulting from -maxima filtering.

Demo Script

 1 import ia870 as MT
 2 import numpy as np
 3 
 4 f = adreadgray('MVBook/ggray.png')
 5 
 6 fn = MT.iaintersec(MT.ianeg(MT.iaframe(f)),MT.iaaddm(f, uint8(20 * np.random.rand(shape(f)[0],shape(f)[1]))));
 7 f1 = MT.iaregmax(fn, MT.iasecross())
 8 f2 = MT.iaregmax(MT.iaopenrec(fn, MT.iasedisk(3), MT.iasecross()), MT.iasecross())
 9 f3 = MT.iaregmax(MT.iaopen(fn, MT.iasedisk(15)), MT.iasecross())
10 f4 = MT.iaregmax(MT.iahmax(fn,150, MT.iasecross()), MT.iasecross())
11 f5 = MT.iaregmax (mmvmax(fn, 23000)) # Criar a função iavmax e iavmin
12 
13 adshow(f,'a')
14 adshow(f1,'b')
15 adshow(f2,'c')
16 adshow(f3, 'd')
17 adshow(f4,'e')
18 adshow(f5,'f')

a

b

c

d

e

f