<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <title></title>
</head>
<body bgcolor="#ffffff">
 
<div></div>
 
<div
 style="font-family: arial; font-style: normal; font-variant: normal; font-weight: normal; font-size: 10pt; line-height: normal; font-stretch: normal; font-size-adjust: none;"></div>
 
<div><font face="Arial" size="2"></font><font face="Arial" size="2">Hi <font
 face="Times New Roman" size="3">Guorong,</font></font></div>
 
<div>&nbsp;</div>
 
<div><font face="Arial" size="2">QUESTION 2:</font></div>
 
<div><font face="Arial" size="2"></font>&nbsp;</div>
 
<div><font face="Arial" size="2">The SymmetricForcesDemonsRegistrationFilter
is  based on&nbsp;formula (5) out of Thirion's paper:</font></div>
 
<div><font face="Arial" size="2">Jean-Philippe Thirion: Fast Non-Rigid Matching
of  3D Medical Images; May 1995; Publication no. 2547 at INRIA (</font><font
 face="Arial" size="2"><a href="http://www.inria.fr/rrrt/rr-2547.html">http://www.inria.fr/rrrt/rr-2547.html</a>).</font></div>
 
<div><font face="Arial" size="2"></font>&nbsp;</div>
 
<div><font face="Arial" size="2">It is basiacally the same as the original
formula  (4) except that instead of just considering one image gradient (the
one of the  fixed image), the moving image gradient is included as well.
This leads to  better results where there is a large deformation.</font></div>
 
<div><font face="Arial" size="2"></font>&nbsp;</div>
 
<div><font face="Arial" size="2">The implementation of this formula can be
found in  the file SymmetricForcesDemonsRegistrationFunction on&nbsp;line  163-323.</font></div>
 
<div><font face="Arial" size="2"></font>&nbsp;</div>
 
<div><font face="Arial" size="2">Hope this answers your question,</font></div>
 
<div><font face="Arial" size="2">Corinne</font></div>
 
<div><font face="Arial" size="2"></font>&nbsp;</div>
 
<div><font face="Arial" size="2">___________________________________________________________</font></div>
 
<div><font face="Arial" size="2"></font><br>
Message: 2<br>
Date: Thu, 15 Jul 2004  15:12:55 +0800<br>
From: "Guorong Wu" <br>
Subject: [Insight-users] Question  about<br>
SymmetricForcesDemonsRegistrationFilter  and<br>
LevelSetMotionRegisrationFilter<br>
To: ITK Mailling List<br>
Message-ID:  <a class="moz-txt-link-rfc2396E"
 href="mailto:20040715071305.26CA27B09EE@mx1.sjtu.edu.cn">&lt;20040715071305.26CA27B09EE@mx1.sjtu.edu.cn&gt;</a><br>
Content-Type:  text/plain; charset="gb2312"<br>
<br>
Hello all:<br>
In the directory  /Insight/Examples/Registration there are three demons like
non-rigid  registration examples: DeformableRegistration2.cxx, DeformableRegistration3.cxx
 DeformableRegistration5.cxx respectively. DeformableRegistration2.cxx is
based  on Thirion's original demons algorithm, and DeformableRegistration3.cxx
use  symmetric forces which can avoid the non-uniformity of deformation field,
while  DeformRegistration5.cxx implements a deformable registration algorithm
that  align a fixed image and a moving image under level set motion. <br>
I'v tested  those three programs for sevel MR brain images finding that levelsetmotionfilter
 is the most efficient one. Symmetricforcefilter does great better result
than  Thirion's but the computation time is three times than demonsregistrationfilter.
 The match result of levelsetmotionfilter is about 10% better than  symmetricforcefilter
in every test while takes much less (about 30%) time than  demonsregistraionfilter.
<br>
I am so interested in  itkSymmetricForcesDemonsRegistrationFilter and  itkLevelSetMotionRegistrationFilter
that I look into the source code of these  algorithms, but I am got confused
in some implementations in the souce code. So  I will be very appreciated
if someone would give me some help. Thanks in  advance.<br>
<br>
QUESTION1:&nbsp; This is about LevelSetMotionRegistrationFilter.  I'v read the
paper mentioned in  /Code/Algorithms/itkLevelSetMotionRegistrationFilter.h.
In this paper, the  process of registration is controled by two PEDs (equation
23 and equation 27).  But in itkLevelSetMotionRegistrationFunction I could
only find the  implementation of equation 27 in function ComputeUpdate(...)
at line 399. Could  someone tell me where is the implementation of equation
23?<br>
QUESTION2:  SysmmetricForcesRegistrationFilter is the variant of triditional
demons which is  contributed by Corinne Mattmann. I could see the modification
through the source  code but cann't tell why does it take the effect. Would
you like to recommmend  some paper?<br>
QUESTION3: Most of the registraion methods is proposed in the  multi-resolution
framework. But the example programs does not use  multi-resolution. If I
try to add it to the examples what kind of work I should  do or what should
I take careful?<br>
<br>
Good Luck.  <br>
<br>
<br>
????????????????Guorong  Wu<br>
????????????????grwu@sjtu.edu.cn<br>
????????????????????2004-07-15<br>
<br>
=============================================================<br>
Guorong  Wu <br>
The Department of Computer Science and Technology<br>
Shanghai Jiao Tong  University<br>
1954 Huashan RD<br>
Shanghai China&nbsp;&nbsp;&nbsp;  <br>
200030<br>
Tel:&nbsp; +86 21 3226 1236&nbsp; , 6293  2089<br>
=============================================================<br>
--------------  next part --------------<br>
A non-text attachment was scrubbed...<br>
Name:  fox.gif<br>
Type: image/gif<br>
Size: 9519 bytes<br>
Desc: not available<br>
Url :  <a class="moz-txt-link-freetext"
 href="http://public.kitware.com/pipermail/insight-users/attachments/20040715/c431ca5c/fox-0001.gif">http://public.kitware.com/pipermail/insight-users/attachments/20040715/c431ca5c/fox-0001.gif</a><br>
<br>
<br>
<br>
</div>
</body>
</html>