<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman,new york,times,serif;font-size:12pt"><div>Hi all,<br><br>I am trying to create a one voxel image (that will grow depending on the number of steps in the algorithm) using to following code:<br>&nbsp;&nbsp;&nbsp; const unsigned int Dimension = 3;<br>&nbsp;&nbsp;&nbsp; typedef itk::Matrix&lt;double, 4, 4&gt; MatrixType;<br>&nbsp;&nbsp;&nbsp; typedef itk::Image&lt;MatrixType, Dimension&gt; MatrixFieldType;<br>&nbsp;&nbsp;&nbsp; MatrixFieldType::Pointer matrixfield = MatrixFieldType::New();<br>&nbsp;&nbsp;&nbsp; MatrixFieldType::IndexType index;<br>&nbsp;&nbsp;&nbsp; index.Fill( 0 );<br>&nbsp;&nbsp;&nbsp; MatrixFieldType::SizeType size;<br>&nbsp;&nbsp;&nbsp; size.Fill( 1 );<br>&nbsp;&nbsp;&nbsp; MatrixFieldType::RegionType region( index, size );<br>&nbsp;&nbsp;&nbsp; matrixfield-&gt;SetRegions( region );<br>&nbsp;&nbsp;&nbsp;
 matrixfield-&gt;Allocate();<br>&nbsp;&nbsp;&nbsp; MatrixType matrix;<br>&nbsp;&nbsp;&nbsp; matrix.SetIdentity();<br>&nbsp;&nbsp;&nbsp; matrixfield-&gt;SetPixel( index, matrix );<br>However this code crashes.<br>When I fill size with 2's, the code seems to work.<br><br>Any ideas? Maarten<br></div>
</div><br>

      </body></html>