deformGrid2d: visualise differences between two superimposed sets of 2D landmarks

Description

visualise differences between two superimposed sets of 2D landmarks by deforming a square grid based on a thin-plate spline interpolation

Usage

deformGrid2d(
  matrix,
  tarmatrix,
  ngrid = 0,
  lwd = 1,
  show = c(1:2),
  lines = TRUE,
  lcol = 1,
  lty = 2,
  col1 = 2,
  col2 = 3,
  pcaxis = FALSE,
  add = FALSE,
  wireframe = NULL,
  margin = 0.2,
  gridcol = "grey",
  gridlty = 1,
  cex1 = 1,
  cex2 = 1,
  ...
)

Value

if ngrid > 1 the coordinates of the displaced grid knots are returned.

Arguments

matrix: reference matrix containing 2D landmark coordinates or mesh of class "mesh3d"
tarmatrix: target matrix containing 2D landmark coordinates or mesh of class "mesh3d"
ngrid: number of grid lines to be plotted; ngrid=0 suppresses grid creation.
lwd: width of lines connecting landmarks.
show: integer (vector): if c(1:2) both configs will be plotted, show = 1 only plots the reference and show = 2 the target. plotted. Options are combinations of 1,2 and 3.
lines: logical: if TRUE, lines between landmarks will be plotted.
lcol: color of lines
lty: line type
col1: color of "matrix"
col2: color of "tarmat"
pcaxis: logical: align grid by shape's principal axes.
add: logical: if TRUE, output will be drawn on existing plot.
wireframe: list/vector containing row indices to be plotted as wireframe (see lineplot.)
margin: margin around the bounding box to draw the grid
gridcol: color of the grid
gridlty: linetype for grid
cex1: control size of points belonging to matrix
cex2: control size of points belonging to tarmatrix
...: additional parameters passed to plot

Author

Stefan Schlager

Examples

Run this code

if (require(shapes)) {
proc <- procSym(gorf.dat)
deformGrid2d(proc$mshape,proc$rotated[,,1],ngrid=5,pch=19)
}