points2LabelImage: create a label image from landmark coordinates
In zarquon42b/RANTs: Playing with ANTs(R)

Description Usage Arguments Value Author(s) Examples

create a label image from landmark coordinates

1
2
3

points2LabelImage(pts, neighbours = 16, spacing = rep(1, ncol(pts)),
  margin = 0.1, labelcolors = seq_len(nrow(pts)), RAS2IJK = NULL,
  refimage = NULL)

`pts`	k x 2 or k x 3 matrix containing landmark coordinates
`neighbours`	integer: number of closest points in the image grid to label
`spacing`	spacing of the output image
`margin`	margin around the landmarks (determined by the extend of the bounding box) to be part of the image
`labelcolors`	vector of length nrow(pts) with a color for each landmark label
`RAS2IJK`	RAS2IJK 4x4 (3D) or 3x3 (2D) matrix with transform from RAS to IJK space (e.g. when landmarks are placed with Slicer).
`refimage`	optional: project the point label image into the image domain defined by refimage.

returns an antsImage

Schlager S

load(system.file("extdata/landmarks.RData",package="ANTsR"))
fi <- antsImageRead(getANTsRData("r16") ,2)
li <- points2LabelImage(fiLM,refimage=fi,RAS2IJK=diag(c(-1,-1,1)))
plot(fi,list(li),alpha=0.3,add=T)