imager (version 0.14)

periodic.part: Compute the periodic part of an image, using the periodic/smooth decomposition of Moisan (2009)

Description

Moisan (2009) defines an additive image decomposition im = periodic + smooth where the periodic part shouldn't be too far from the original image. The periodic part can be used in frequency-domain analyses, to reduce the artifacts induced by non-periodicity.

Usage

periodic.part(im)

Arguments

im

an image

Value

• an image

References

L. Moisan, Periodic plus Smooth Image Decomposition,J. Math. Imaging Vision, vol. 39:2, pp. 161-179, 2011

Examples

imname <- system.file('extdata/parrots.png',package='imager')
im <- load.image(imname) %>% subim(x <= 512)
layout(t(1:3))
plot(im,main="Original image")
periodic.part(im) %>% plot(main="Periodic part")
#The smooth error is the difference between
#the original image and its periodic part
(im-periodic.part(im)) %>% plot(main="Smooth part")