<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META content="text/html; charset=iso-8859-1" http-equiv=Content-Type>
<META name=GENERATOR content="MSHTML 8.00.6001.19258"></HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT size=2 face="Courier New">Gareth,</FONT></DIV>
<DIV><FONT size=2 face="Courier New"></FONT> </DIV>
<DIV><FONT size=2 face="Courier New">Sorry for the problem you experienced. I am
a main developer for GPU demons in ITK4.2. Would you please send me the images
and the testing parameters you used? </FONT></DIV>
<DIV><FONT size=2 face="Courier New"></FONT> </DIV>
<DIV><FONT size=2 face="Courier New">Thank you for reporting this
issue!</FONT></DIV>
<DIV><FONT size=2 face="Courier New"></FONT> </DIV>
<DIV><FONT size=2 face=Arial>Bao</FONT></DIV>
<BLOCKQUOTE
style="BORDER-LEFT: #000000 2px solid; PADDING-LEFT: 5px; PADDING-RIGHT: 0px; MARGIN-LEFT: 5px; MARGIN-RIGHT: 0px"
dir=ltr>
<DIV style="FONT: 10pt arial">----- Original Message ----- </DIV>
<DIV
style="FONT: 10pt arial; BACKGROUND: #e4e4e4; font-color: black"><B>From:</B>
<A title=gee@mail.med.upenn.edu href="mailto:gee@mail.med.upenn.edu">James
Gee</A> </DIV>
<DIV style="FONT: 10pt arial"><B>To:</B> <A title=baohua@mail.med.upenn.edu
href="mailto:baohua@mail.med.upenn.edu">Baohua Wu</A> ; <A
title=stnava@gmail.com href="mailto:stnava@gmail.com">Avants Brian</A> </DIV>
<DIV style="FONT: 10pt arial"><B>Sent:</B> Thursday, July 19, 2012 9:30
AM</DIV>
<DIV style="FONT: 10pt arial"><B>Subject:</B> Fwd: [Insight-users] GPU Demons
registration in ITK 4.2</DIV>
<DIV><BR></DIV>
<DIV><BR><BR>Sent from my iPhone</DIV>
<DIV><BR>Begin forwarded message:<BR><BR></DIV>
<BLOCKQUOTE type="cite">
<DIV><B>From:</B> Gareth Price <<A
href="mailto:Gareth.Price@physics.cr.man.ac.uk">Gareth.Price@physics.cr.man.ac.uk</A>><BR><B>Date:</B>
July 19, 2012 8:24:40 AM EDT<BR><B>To:</B> <<A
href="mailto:insight-users@itk.org">insight-users@itk.org</A>><BR><B>Subject:</B>
<B>[Insight-users] GPU Demons registration in ITK
4.2</B><BR><BR></DIV></BLOCKQUOTE>
<DIV></DIV>
<BLOCKQUOTE type="cite">
<DIV>
<META name=Generator content="Microsoft Word 10 (filtered)">
<STYLE>@font-face {
font-family: Wingdings;
}
@font-face {
font-family: Calibri;
}
@page Section1 {size: 612.0pt 792.0pt; margin: 72.0pt 90.0pt 72.0pt 90.0pt; }
P.MsoNormal {
MARGIN: 0cm 0cm 0pt; FONT-FAMILY: "Times New Roman"; FONT-SIZE: 12pt
}
LI.MsoNormal {
MARGIN: 0cm 0cm 0pt; FONT-FAMILY: "Times New Roman"; FONT-SIZE: 12pt
}
DIV.MsoNormal {
MARGIN: 0cm 0cm 0pt; FONT-FAMILY: "Times New Roman"; FONT-SIZE: 12pt
}
H1 {
PAGE-BREAK-AFTER: avoid; MARGIN: 12pt 0cm 3pt; FONT-FAMILY: Arial; FONT-SIZE: 16pt
}
H2 {
PAGE-BREAK-AFTER: avoid; FONT-STYLE: italic; MARGIN: 12pt 0cm 3pt; FONT-FAMILY: Arial; FONT-SIZE: 14pt
}
H3 {
PAGE-BREAK-AFTER: avoid; MARGIN: 12pt 0cm 3pt; FONT-FAMILY: Arial; FONT-SIZE: 13pt
}
A:link {
COLOR: blue; TEXT-DECORATION: underline
}
SPAN.MsoHyperlink {
COLOR: blue; TEXT-DECORATION: underline
}
A:visited {
COLOR: purple; TEXT-DECORATION: underline
}
SPAN.MsoHyperlinkFollowed {
COLOR: purple; TEXT-DECORATION: underline
}
P.EPSRCHeading1 {
PAGE-BREAK-AFTER: avoid; MARGIN: 0cm 0cm 6pt; FONT-FAMILY: Calibri; FONT-SIZE: 14pt; FONT-WEIGHT: bold
}
LI.EPSRCHeading1 {
PAGE-BREAK-AFTER: avoid; MARGIN: 0cm 0cm 6pt; FONT-FAMILY: Calibri; FONT-SIZE: 14pt; FONT-WEIGHT: bold
}
DIV.EPSRCHeading1 {
PAGE-BREAK-AFTER: avoid; MARGIN: 0cm 0cm 6pt; FONT-FAMILY: Calibri; FONT-SIZE: 14pt; FONT-WEIGHT: bold
}
P.EPSRCHeading2 {
PAGE-BREAK-AFTER: avoid; TEXT-TRANSFORM: uppercase; FONT-STYLE: italic; MARGIN: 3pt 0cm 0pt; FONT-FAMILY: Calibri; FONT-SIZE: 12pt; FONT-WEIGHT: bold
}
LI.EPSRCHeading2 {
PAGE-BREAK-AFTER: avoid; TEXT-TRANSFORM: uppercase; FONT-STYLE: italic; MARGIN: 3pt 0cm 0pt; FONT-FAMILY: Calibri; FONT-SIZE: 12pt; FONT-WEIGHT: bold
}
DIV.EPSRCHeading2 {
PAGE-BREAK-AFTER: avoid; TEXT-TRANSFORM: uppercase; FONT-STYLE: italic; MARGIN: 3pt 0cm 0pt; FONT-FAMILY: Calibri; FONT-SIZE: 12pt; FONT-WEIGHT: bold
}
P.EPSRCHeading3 {
PAGE-BREAK-AFTER: avoid; FONT-STYLE: italic; MARGIN: 3pt 0cm 0pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt; FONT-WEIGHT: bold
}
LI.EPSRCHeading3 {
PAGE-BREAK-AFTER: avoid; FONT-STYLE: italic; MARGIN: 3pt 0cm 0pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt; FONT-WEIGHT: bold
}
DIV.EPSRCHeading3 {
PAGE-BREAK-AFTER: avoid; FONT-STYLE: italic; MARGIN: 3pt 0cm 0pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt; FONT-WEIGHT: bold
}
P.EPSRCNormal {
LINE-HEIGHT: 12pt; TEXT-INDENT: 17pt; MARGIN: 0cm 0cm 3pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
LI.EPSRCNormal {
LINE-HEIGHT: 12pt; TEXT-INDENT: 17pt; MARGIN: 0cm 0cm 3pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
DIV.EPSRCNormal {
LINE-HEIGHT: 12pt; TEXT-INDENT: 17pt; MARGIN: 0cm 0cm 3pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
P.EPSRCBullets {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 3pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
LI.EPSRCBullets {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 3pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
DIV.EPSRCBullets {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 3pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
P.EPSRCBulletsCxSpFirst {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 0pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
LI.EPSRCBulletsCxSpFirst {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 0pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
DIV.EPSRCBulletsCxSpFirst {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 0pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
P.EPSRCBulletsCxSpMiddle {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 0pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
LI.EPSRCBulletsCxSpMiddle {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 0pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
DIV.EPSRCBulletsCxSpMiddle {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 0pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
P.EPSRCBulletsCxSpLast {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 3pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
LI.EPSRCBulletsCxSpLast {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 3pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
DIV.EPSRCBulletsCxSpLast {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 3pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
P.EPSRCNumbering {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 3pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
LI.EPSRCNumbering {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 3pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
DIV.EPSRCNumbering {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 3pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
P.EPSRCNumberingCxSpFirst {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 0pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
LI.EPSRCNumberingCxSpFirst {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 0pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
DIV.EPSRCNumberingCxSpFirst {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 0pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
P.EPSRCNumberingCxSpMiddle {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 0pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
LI.EPSRCNumberingCxSpMiddle {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 0pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
DIV.EPSRCNumberingCxSpMiddle {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 0pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
P.EPSRCNumberingCxSpLast {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 3pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
LI.EPSRCNumberingCxSpLast {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 3pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
DIV.EPSRCNumberingCxSpLast {
LINE-HEIGHT: 12pt; TEXT-INDENT: -11.35pt; MARGIN: 0cm 0cm 3pt 19.85pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
P.StyleEPSRCNormalCentered {
TEXT-ALIGN: center; LINE-HEIGHT: 12pt; TEXT-INDENT: 17pt; MARGIN: 0cm 0cm 3pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
LI.StyleEPSRCNormalCentered {
TEXT-ALIGN: center; LINE-HEIGHT: 12pt; TEXT-INDENT: 17pt; MARGIN: 0cm 0cm 3pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
DIV.StyleEPSRCNormalCentered {
TEXT-ALIGN: center; LINE-HEIGHT: 12pt; TEXT-INDENT: 17pt; MARGIN: 0cm 0cm 3pt; FONT-FAMILY: Calibri; FONT-SIZE: 11pt
}
P.CVList {
TEXT-ALIGN: justify; TEXT-INDENT: -18pt; MARGIN: 0cm 0cm 6pt 36pt; FONT-FAMILY: "Times New Roman"; COLOR: black; FONT-SIZE: 12pt
}
LI.CVList {
TEXT-ALIGN: justify; TEXT-INDENT: -18pt; MARGIN: 0cm 0cm 6pt 36pt; FONT-FAMILY: "Times New Roman"; COLOR: black; FONT-SIZE: 12pt
}
DIV.CVList {
TEXT-ALIGN: justify; TEXT-INDENT: -18pt; MARGIN: 0cm 0cm 6pt 36pt; FONT-FAMILY: "Times New Roman"; COLOR: black; FONT-SIZE: 12pt
}
SPAN.EmailStyle25 {
FONT-FAMILY: Arial; COLOR: windowtext
}
DIV.Section1 {
page: Section1
}
OL {
MARGIN-BOTTOM: 0cm
}
UL {
MARGIN-BOTTOM: 0cm
}
</STYLE>
<DIV class=Section1>
<P class=MsoNormal><FONT size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt">Hi
there,</SPAN></FONT></P>
<P class=MsoNormal><FONT size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt"></SPAN></FONT> </P>
<P class=MsoNormal><FONT size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt">I was wondering if
someone might be able to offer advice regarding the new ITK 4.2
</SPAN></FONT><FONT size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt">GPU</SPAN></FONT><FONT
size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt"> implementation of the
demons algorithm.</SPAN></FONT></P>
<P class=MsoNormal><FONT size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt"></SPAN></FONT> </P>
<P class=MsoNormal><FONT size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt">I have an application
using the standard CPU implementation to register 3D CT images (512x512x~70
voxels). This behaves perfectly when registering synthetic deformations to
their original image. I typically use a smoothing kernel standard deviation
of 10.0 with no update field smoothing. With my images, the metric reduces
exponentially as expected from ~8000 to 350 over 50 iterations with the RMS
update field change showing a similar but slower trend from ~0.08 to 0.04.
After 50 iterations the solutions have made good progress towards the
original images. </SPAN></FONT></P>
<P class=MsoNormal><FONT size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt"></SPAN></FONT> </P>
<P class=MsoNormal><FONT size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt">Porting this code to the
</SPAN></FONT><FONT size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt">GPU</SPAN></FONT><FONT
size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt"> implementation
following the examples in the itkGPUDemonsRegistrationFilterTests (i.e.
swapping the internal image type for GPUImageType and
GPUDemonsRegistrationFilter for DemonsRegistrationFilter) provides an
immediate ~x15 increase in speed on my machine (GeForce GTX480). However
using the same input data and smoothing parameters, the metric and RMS
update field trends, although still exponential show changes from ~8000 to
1500 (metric) and ~0.3 to 0.2 (update field). The result is that the
deformation field produced is much too large. Changing smoothing parameters
does not change this behaviour, whilst enabling update field smoothing
crashes the process (in GPUDemonsRegistrationFilter::ApplyUpdate). It is
very strange that the similarity metric is continuing to reduce as the
deformation field pushes the moving image through and past its correct fixed
image alignment.</SPAN></FONT></P>
<P class=MsoNormal><FONT size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt"></SPAN></FONT> </P>
<P class=MsoNormal><FONT size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt">Would someone be able to
tell me if this is expected behaviour/are known bugs, or offer advice if I
am using the </SPAN></FONT><FONT size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt">GPU</SPAN></FONT><FONT
size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt"> registration filter
incorrectly?</SPAN></FONT></P>
<P class=MsoNormal><FONT size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt"></SPAN></FONT> </P>
<P class=MsoNormal><FONT size=2 face="Courier New"><SPAN
style="FONT-FAMILY: 'Courier New'; FONT-SIZE: 10pt">Many thanks,
Gareth</SPAN></FONT></P>
<P class=MsoNormal><FONT size=2 face=Arial><SPAN
style="FONT-FAMILY: Arial; FONT-SIZE: 10pt"></SPAN></FONT> </P></DIV>
<DIV>
<P></P>
<HR>
****************************************************************<BR>This
e-mail and any files transmitted with it are confidential<BR>and solely for
the use of the intended recipient. If you have<BR>received this e-mail in
error you should not disseminate,<BR>distribute or copy it. Please notify
the sender immediately and<BR>delete this e-mail from your
system.<BR>****************************************************************<BR>
<P></P></DIV></DIV></BLOCKQUOTE>
<BLOCKQUOTE type="cite">
<DIV><SPAN>_____________________________________</SPAN><BR><SPAN>Powered by
<A
href="http://www.kitware.com">www.kitware.com</A></SPAN><BR><SPAN></SPAN><BR><SPAN>Visit
other Kitware open-source projects at</SPAN><BR><SPAN><A
href="http://www.kitware.com/opensource/opensource.html">http://www.kitware.com/opensource/opensource.html</A></SPAN><BR><SPAN></SPAN><BR><SPAN>Kitware
offers ITK Training Courses, for more information visit:</SPAN><BR><SPAN><A
href="http://www.kitware.com/products/protraining.php">http://www.kitware.com/products/protraining.php</A></SPAN><BR><SPAN></SPAN><BR><SPAN>Please
keep messages on-topic and check the ITK FAQ at:</SPAN><BR><SPAN><A
href="http://www.itk.org/Wiki/ITK_FAQ">http://www.itk.org/Wiki/ITK_FAQ</A></SPAN><BR><SPAN></SPAN><BR><SPAN>Follow
this link to subscribe/unsubscribe:</SPAN><BR><SPAN><A
href="http://www.itk.org/mailman/listinfo/insight-users">http://www.itk.org/mailman/listinfo/insight-users</A></SPAN><BR></DIV></BLOCKQUOTE></BLOCKQUOTE></BODY></HTML>