imageDomainToSpatialMatrix.RdTransform a d-dimensional image into a d by p matrix that encodes the spatial coordinates of the image domain.
imageDomainToSpatialMatrix(image, mask)matrix is output
img <- antsImageRead(getANTsRData("r16"))
spmat <- imageDomainToSpatialMatrix(img)