Code covered by the BSD License  

Highlights from
Diffeomorphic Log Demons Image Registration

5.0
5.0 | 5 ratings Rate this file 42 Downloads (last 30 days) File Size: 195 KB File ID: #39194 Version: 1.7
image thumbnail

Diffeomorphic Log Demons Image Registration

by

 

27 Nov 2012 (Updated )

Simple matlab code for 2D and 3D image registration using the diffeomorphic log-demons algorithm

| Watch this File

File Information
Description

Simple matlab code for 2D and 3D image registration using the diffeomorphic log-demons algorithm
Code is provided in order to help the understanding of the Demons algorithm - Any comment or improvement is welcome

MATLAB release MATLAB 7.10 (R2010a)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (11)
05 May 2015 zhang le

I am confused that I can not get the same result showing in Fig. 6 ( Classical Circle to C registration example)

12 Jan 2015 Nikita

Nikita (view profile)

Please let me know what changes should be incorporated to work on RGB images.

Comment only
10 Jul 2014 Nicolas Yu

Some details should be changed for my purpose, and thanks.

28 Apr 2014 Xun Zhang

I'm running the 2D demons. I'm confused about findupdate.m. In findupdate.m, your code computes the gradient of the moving image or float image. As far as I have learnt from literature, gradient of fixed image are used in the expression for u. I'm not quite clear about this.

Comment only
18 Apr 2014 Herve Lombaert

@tianyu - You are raising an interesting point - There is a difference between Displacement field != Transformation (zero displacement leads to an identity matrix)

Typically, the identity matrix is added when computing the Jacobian determinant

See also this thread (from the itk code): http://www.cmake.org/Bug/bug_relationship_graph.php?bug_id=7327&graph=dependency&orientation=vertical (code is now including additional comments, perhaps this would help others)

Comment only
17 Apr 2014 tianyu fu

@Herve LomBaert: And I don't think the code is wrong.But I can't get the same result showing in the paper.This makes me very upset.

Comment only
17 Apr 2014 tianyu fu

@Herve LomBaert: Thank you for your reply.I have read many papers about the diffeomorphic Log-Demons.As you know,if the registration is diffeomorphic,so the Jacobian determinant of deformation must be positive.But,running your code,I find the determinant sometimes is negative.

Comment only
17 Apr 2014 Herve Lombaert

@tianyu - Indeed, this is the Log-domain Diffeomorphic version, the deformation field is defined as the exponential map of a velocity field via expfield() - (velocity field != deformation field - More details on the associated literature on the Diffeomorphic Log-Demons)

Comment only
17 Apr 2014 tianyu fu

Is it Diffeomorphic Image Registration?I do not think the Jacobi determinant of the transformation field is positive.

04 Mar 2014 Ernst Schwartz

@Hg: in image registration, you have to know how to deal with the boundaries of the images. A common choice is to have them fixed. So, in order not to over-regularize nor violate the diffeomorphism, a band of zeros is added during optimization.

09 Oct 2013 Hg

Hg (view profile)

May I know the purpose of the piggyback function? Where can I find the explanation for the algorithm used in the code other than the paper included in the package?

Updates
26 Dec 2012 1.2

cleanup in sample files

29 Jan 2013 1.5

Log domain

18 Apr 2014 1.6

Fix in 2D (adding comments regarding Jacobian computation)

19 Aug 2014 1.7

comment updates

Contact us