register: Registration based on the NMI criteria

In coldfir3/surf: Surface and topological analysis for tribology

Description

bla bla

Usage

register(src, tar, method = "gen", par0 = NULL, lower = NULL,
  upper = NULL, ran = c(0.25, 0.25, 0), control = NULL, verbosity = 2)

Arguments

src, tar	cimg objects
method	character indicating the method to be used
par0	numeric vector for the initial guess for the registration parameters
lower, upper	numeric vector for the bounds of the optimization problem
ran	numeric vector for the range around par0 for lower and upper (when not provided)
control	list to control the parameters of the optimization routine
verbosity	Numeric indicating the level of verbosity is displayed

Value

bla bla

Examples

# simulated shif to check if the taylor approximation is correct
imA <- ground[[1]]
imB <- transform(imA, 0.51, -0.23, 0, method = 'taylor')
round(register(imA, imB, method = 'taylor'),2)

##################
## The following examples were exluded from the automatic
## run to save compilation time. However, they should all be working.
# Now lets check the differences between the algorithms for a real pair of images
## Not run: 
imB <- ground[[10]]
par.taylor <- register(imA, imB, method = 'taylor')
par.grad   <- register(imA, imB, method = 'grad')
par.gen    <- register(imA, imB, method = 'gen')
par.sa    <- register(imA, imB, method = 'sa')
imA.grad <- transform(imA, par.grad[1], par.grad[2], 0)
imA.gen  <- transform(imA, par.gen[1], par.gen[2], 0)
imA.sa   <- transform(imA, par.sa[1], par.sa[2], 0)
NMI(imA,imB)
NMI(imA.grad,imB)
NMI(imA.gen,imB)
NMI(imA.sa,imB)

## End(Not run)
# Now, lets see if we can visualize the motion of the sample
## Not run: 
pars.taylor <- t(sapply(ground, register, tar = ground[[1]], method = 'taylor'))
pars.grad   <- t(sapply(ground, register, tar = ground[[1]], method = 'grad'))
pars.gen    <- t(sapply(ground, register, tar = ground[[1]], method = 'gen'))
pars.sa     <- t(sapply(ground, register, tar = ground[[1]], method = 'sa'))
plot(0, pch = NA, asp = 1,
 xlim = range(rbind(pars.taylor, pars.grad, pars.gen, pars.sa)[,1]),
 ylim = range(rbind(pars.taylor, pars.grad, pars.gen, pars.sa)[,2]))
points(pars.taylor[,1:2], col = 'orange')
points(pars.grad[,1:2], col = 'skyblue')
points(pars.gen[,1:2], col = 'black')
points(pars.sa[,1:2], col = 'red')
lines(pars.taylor[,1:2], col = 'orange')
lines(pars.grad[,1:2], col = 'skyblue')
lines(pars.gen[,1:2], col = 'black')
lines(pars.sa[,1:2], col = 'red')
plot(0, pch = NA, xlim = c(2,80), ylim = c(0.5,0.7))
lines(sapply(2:80, function(i) NMI(ground[[1]], transform(ground[[i]],
 pars.taylor[i,1], pars.taylor[i,2], pars.taylor[i,3]))), col = 'orange')
lines(sapply(2:80, function(i) NMI(ground[[1]], transform(ground[[i]],
 pars.grad[i,1], pars.grad[i,2], pars.grad[i,3]))), col = 'skyblue')
lines(sapply(2:80, function(i) NMI(ground[[1]], transform(ground[[i]],
 pars.gen[i,1], pars.gen[i,2], pars.gen[i,3]))), col = 'black')
lines(sapply(2:80, function(i) NMI(ground[[1]], transform(ground[[i]],
 pars.sa[i,1], pars.sa[i,2], pars.sa[i,3]))), col = 'red')

## End(Not run)

coldfir3/surf documentation built on May 13, 2019, 8:49 p.m.