<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
</style>
</head>
<body class='hmmessage'>
Hi, Luis,<div><br></div><div>Thank you for your reply.</div><div><br></div><div>I have 2GB RAM. Here is the error message if I run more than 62 slices:&nbsp;</div><div>Unhandled exception at 0x7c812afb in RegTest.exe: Microsoft C++ exception: itk::InvalidRequestedRegionError at memory location 0x0013fbbc..</div><div><br>My CT volumes are liver images from two different phases.&nbsp;</div><div><br></div><div>I also did another test on these two volumes by switch them (switch fixed and moving volumes). If I run 62 slices,, it would come with same memory problem (see above error message). This time if I run 60 slices, it would&nbsp;reach "CONVERGENCE: REL_REDUCTION_OF_F &lt;= FACTR*EPSMCH" after one iteration; but it seems 59 slices allow me to run 34 iterations and gave me a good result.</div><div><br></div><div>So, it seems the good results came with "reach&nbsp;Projected Gradient Tolerance". If it&nbsp;reaches "CONVERGENCE: REL_REDUCTION_OF_F &lt;= FACTR*EPSMCH" first, it come with a "bad" result. How can I control this situation?</div><div><br></div><div>Thanks.</div><div><br></div><div>Sean</div><div><br>&gt; Date: Wed, 8 Jul 2009 17:39:25 -0400<br>&gt; From: luis.ibanez@kitware.com<br>&gt; To: kevinseanlee@live.com<br>&gt; CC: insight-users@itk.org<br>&gt; Subject: Re: [Insight-users] problems on deformable registration 8<br>&gt; <br>&gt; <br>&gt; Hi Sean,<br>&gt; <br>&gt; This is very suspicious.<br>&gt; <br>&gt; <br>&gt; A) your images do not seem to be that big...<br>&gt; <br>&gt;     One image of:<br>&gt; <br>&gt;       64 slices of 256 x 256 in 16bits<br>&gt; <br>&gt;     will take 8 Mb<br>&gt; <br>&gt; <br>&gt;     How much RAM do you have ?<br>&gt; <br>&gt; <br>&gt; B) Just removing one slice (out of 60) from your<br>&gt;     dataset shouldn't lead the registration to<br>&gt;     behave differently.<br>&gt; <br>&gt; <br>&gt;     Unless...<br>&gt; <br>&gt;     There is some special content in that slice...<br>&gt; <br>&gt;     Are any objects in Slice 60 that have very<br>&gt;     different intensities to those objects visible<br>&gt;     in slice 59  ?<br>&gt; <br>&gt; <br>&gt; ...and...<br>&gt; <br>&gt; <br>&gt; C) How are you removing that slice ?<br>&gt; <br>&gt;     Are you using a region of interest extraction first ?<br>&gt; <br>&gt; <br>&gt; <br>&gt; <br>&gt; Please let us know,<br>&gt; <br>&gt; <br>&gt;      Thanks<br>&gt; <br>&gt; <br>&gt; <br>&gt;          Luis<br>&gt; <br>&gt; <br>&gt; <br>&gt; ---------------<br>&gt; Sean lee wrote:<br>&gt; &gt; Hello,<br>&gt; &gt; <br>&gt; &gt; I am trying to use the example of deformable registration 8 on two ct <br>&gt; &gt; image volumes (64 256X256 slices). The problem is If i use more than 62 <br>&gt; &gt; slices, I would have memory problem; I run 60 slice, I got a perfect <br>&gt; &gt; results with 20 iterations; but if I run 59 slices, it only run one <br>&gt; &gt; iteration and reaches "CONVERGENCE: REL_REDUCTION_OF_F &lt;= FACTR*EPSMCH". <br>&gt; &gt; I tried to reduce the Cost Function Convergence Factor to 1e+1 and it <br>&gt; &gt; still comes same results. What could be the problem and any suggested <br>&gt; &gt; solutions? Thanks. Sean<br>&gt; &gt; &gt; <br>&gt; &gt; <br>&gt; &gt; ------------------------------------------------------------------------<br>&gt; &gt; <br>&gt; &gt; _____________________________________<br>&gt; &gt; Powered by www.kitware.com<br>&gt; &gt; <br>&gt; &gt; Visit other Kitware open-source projects at<br>&gt; &gt; http://www.kitware.com/opensource/opensource.html<br>&gt; &gt; <br>&gt; &gt; Please keep messages on-topic and check the ITK FAQ at: http://www.itk.org/Wiki/ITK_FAQ<br>&gt; &gt; <br>&gt; &gt; Follow this link to subscribe/unsubscribe:<br>&gt; &gt; http://www.itk.org/mailman/listinfo/insight-users<br></div><br /><hr />Windows Live™ SkyDrive™: Get 25 GB of free online storage.   <a href='http://windowslive.com/online/skydrive?ocid=TXT_TAGLM_WL_SD_25GB_062009' target='_new'>Get it on your BlackBerry or iPhone.</a></body>
</html>