robust nonrigid point set registration

MATLAB implementation of two nonrigid point set registration algorithms

Description

A unified MATLAB implementation of the following two papers:
(1) Haili Chui, Anand Rangarajan: A New Algorithm for Non-Rigid Point Matching. CVPR 2000: 2044-2051 http://www.cis.ufl.edu/~anand/students/chui/research.html
(2) Andriy Myronenko, Xubo B. Song, Miguel A. Carreira-Perpinan: Non-rigid point set registration: Coherent Point Drift. NIPS 2006: 1009-1016 http://www.csee.ogi.edu/~myron/matlab/cpd/

• MATLAB release: MATLAB 7 (R14)

[Bryn Lloyd]

The example demo seems to work. But I couldn't get the 3d data set (face.ini) to work.

I modified the code to:
dim=3;
config = gmmreg_load_config('./face.ini');
config.motion = 'grbf'; %'tps'
config.init_param = zeros(25,dim);
[fp,fm] = gmmreg_cpd(config);
DisplayPoints(fm,config.scene,dim);

But there is a dimension mismatch:

??? Error using ==> mtimes
Inner matrix dimensions must agree.

Error in ==> gmmreg_cpd at 78
model = model + basis*param;

[Bing Jian]

Bryn, please replace the number 25 in 'config.init_param = zeros(25,dim) ' by the actual number of control points specified in the ini file.

[Ruben]

Hi Bing Jian,

I'm trying to match two volumes defined by a number of contour points. I'd then like to apply the deformable transformation to points within the volume, so see how those have moved. However, no file is written like specified in the .ini file: final_params = ./output/final_params.txt

The matched set of points do get displayed in the figure though and the parameters are stored in the fp variable. Could you advise me on how to apply the transformation (this fp matrix) to other points to get information on deformations / displacements?

Thanks,
Ruben