<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'>


<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
Thanks Sergiy!!<br><br>The problem was the manually created image Dsk; I change the lines as you said me and now it runs perfectly.<br><br>Regards,<br>Daniela<br><br>&gt; Date: Fri, 6 Aug 2010 22:16:48 -0700<br>&gt; From: volkovych@gmail.com<br>&gt; To: insight-users@itk.org<br>&gt; Subject: Re: [Insight-users] Iterators: PROBLEM<br>&gt; <br>&gt; <br>&gt; I think the problem appears when the program tries to increment one of the<br>&gt; iterators.<br>&gt; Try to replace the "for" line to this:<br>&gt; <br>&gt; for (DskIt.GoToBegin(), DFltIt.GoToBegin(),<br>&gt; DSubIt.GoToBegin();!DskIt.IsAtEnd() &amp;&amp; !DFltIt.IsAtEnd() &amp;&amp;<br>&gt; !DSubIt.IsAtEnd(); ++DskIt, ++DFltIt, ++DSubIt)<br>&gt; <br>&gt; If it doesn't go to the internal part of the loop at all, it means that one<br>&gt; of your images is empty (the region wasn't allocated properly). You may<br>&gt; identify which iterator creates a problem by:<br>&gt; <br>&gt; DskIt.GoToBegin();<br>&gt; if (DskIt.IsAtEnd())<br>&gt;   std::cerr &lt;&lt; "Not ok ;-("&lt;&lt; std::endl;<br>&gt; <br>&gt; If the problem appeared due to the manually created image Dsk, try to<br>&gt; allocate it with this lines:<br>&gt; <br>&gt; ImageType::Pointer Dsk = ImageType::New(); <br>&gt; Dsk-&gt;SetRegions(DFlt-&gt;GetLargestPossibleRegion()); <br>&gt; Dsk-&gt;Allocate();<br>&gt; <br>&gt; But be aware that DFlt and DSub images should already have been filled by<br>&gt; the values before doing it.<br>&gt; <br>&gt; Also you may initialize iterators based on the LargestPossibleRegion:<br>&gt; <br>&gt; IteratorType DskIt (Dsk, Dsk-&gt;GetLargestPossibleRegion()); <br>&gt; ConstIteratorType DFltIt (DFlt, DFlt-&gt;GetLargestPossibleRegion()); <br>&gt; ConstIteratorType DSubIt (DSub, DSub-&gt;GetLargestPossibleRegion()); <br>&gt; <br>&gt; If it still doesn't work, please, place the larger part of code (I want to<br>&gt; see how all images are obtained to be sure that you haven't errors in the<br>&gt; other part of code).<br>&gt; <br>&gt; Sergiy.<br>&gt; -- <br>&gt; View this message in context: http://itk-insight-users.2283740.n2.nabble.com/Iterators-PROBLEM-tp5381369p5383008.html<br>&gt; Sent from the ITK Insight Users mailing list archive at Nabble.com.<br>&gt; _____________________________________<br>&gt; Powered by www.kitware.com<br>&gt; <br>&gt; Visit other Kitware open-source projects at<br>&gt; http://www.kitware.com/opensource/opensource.html<br>&gt; <br>&gt; Kitware offers ITK Training Courses, for more information visit:<br>&gt; http://www.kitware.com/products/protraining.html<br>&gt; <br>&gt; Please keep messages on-topic and check the ITK FAQ at:<br>&gt; http://www.itk.org/Wiki/ITK_FAQ<br>&gt; <br>&gt; Follow this link to subscribe/unsubscribe:<br>&gt; http://www.itk.org/mailman/listinfo/insight-users<br>                                               </body>
</html>