multimodality non-rigid demon algorithm image registration

There is a small bug that images with odd valued dimensions now gives an error. Will be updated tomorrow...

I try to complie it with other machine, thanks. Would like to know if you have version for color image~ Thanks~

thanks a lot. I know what you mean. But what I mean is change the Idiff as something like this:
Idiff_color=M_color-S_color;
Idiff = Idiff_color(:,1)+Idiff_color(:,2)+Idiff_color(:,3) ; Idiff = Idiff/3 ;
where M_color is color version of M previous obtained.

I test the code by adding the following:

% Register image 1 to image 2
[I3,Tx,Ty]=demonregistration(I1,I2,[],true);

[h,w] = size(Tx) ;
Ty = zeros(h,w) ;
Tx = 5*ones(h,w) ;

I3=movepixels(I1,Tx,Ty);

but the result is not what I expected.
The image is shifted up instead of left~~~
Indeed what is Tx,Ty doing? is it a deformation field? forward or backward warping?

*K K,
Thanks for your comments. The method uses backward warping. Today, I will look at my code if I accidentally swapped Tx and Ty. This week or the next I will upload a new version, which can be used with any Matlab optimizer.

* Franz Rudolf
This transformation algorithm is called: non-parametric, iconic feature based, using point features.

Please upgrade to version 10b if you experience Matlab crashes with affine transformation in version 10.

*Michel Kocher...
After this weekend i will add a (previous) single threaded version, which works with the Mac computer.

*yichen fan
I have never seen this error.
Please have a look at the help of Mex : http://www.mathworks.com/access/helpdesk/help/techdoc/index.html?/access/helpdesk/help/techdoc/ref/mex.htm

Maybee Mex or the LCC compiler is not installed with Matlab.
You can also trie the single threaded c-code files...

Can you add single thread codes of the latest version. I am using linux and it gives error.
In addition, as using in linux, compilation of mex files gives error just because of the comment lines. You can avoid such errors by replacing "//" with "/*" and of course add "*/" at the end of the line.

Hi, Sir, your codes are definitely good. When I use your 3D example in register_volumes, the results shown by showcs3
looks weird. Also the values range from -0.2 to 0.3 for example.
Can you explain why it. On the other hand, when I used Ireg2=movepixels(Istatic,Fx,Fy,Fz); showcs3(Ireg2);
Ireg2, however, turns out to be surprising good. thanks.

Good work. Thank you.

One question: When you determine the bin size in function "registration_error_mutual_info", you use:

bins=round(numel(V)^(1/ndims(V)));

May I know what this formulation is obtained?

*gx,

The equation "bins=round(numel(V)^(1/ndims(V)));" is not from literature. You can also use a fix bin size of 255, or so.

- Mutual information only works if there are overlaps in the histogram bins.
- More histogram bins will give better registration accuracy.
Thus this equation gives a good trade-off based on the number of available pixels ...

Hello, I know it may be a stupid, but I would like to find a function which could compute a pixel say, (x1, y1) in the moving image 's location in the static image, say (x2, y2). It seems to me that the (x1, y1) and (x2, y2) are related through the so-called transformation . Could you explain the detailed relationship among them and what exactly the direction of x, y and z in your code? (from left to right, up to down and top tp buttom?) Thanks

sorry i made a mistake .The result in multi-registration is also good