<br>Hi Maarten,<br><br>Thanks for your detailed question.<br><br>I tried your code, and...<br><br> It works for me<br><br><br>See the full code below.<br><br>I compiled it in Ubuntu Linux 9.10,<br>with Gcc 4.4, and the CVS version of ITK.<br>
<br>It compiled without errors and ran without<br>crashing.<br><br>It is likely that you have a configuration<br>problem in your build.<br><br><br>Have you run the program in a debugger ?<br><br><br> Regards,<br><br><br>
Luis<br><br><br>--------------------------------------------------------<br>#include "itkImage.h"<br>#include "itkMatrix.h"<br><br>int main()<br>{<br> const unsigned int Dimension = 3;<br>
typedef itk::Matrix<double, 4, 4> MatrixType;<br> typedef itk::Image<MatrixType, Dimension> MatrixFieldType;<br> MatrixFieldType::Pointer matrixfield = MatrixFieldType::New();<br> MatrixFieldType::IndexType index;<br>
index.Fill( 0 );<br> MatrixFieldType::SizeType size;<br> size.Fill( 1 );<br> MatrixFieldType::RegionType region( index, size );<br> matrixfield->SetRegions( region );<br> matrixfield->Allocate();<br>
MatrixType matrix;<br> matrix.SetIdentity();<br> matrixfield->SetPixel( index, matrix );<br><br> return 0;<br>}<br><br>-------------------------------------------------------------------<br><br><div class="gmail_quote">
On Tue, May 25, 2010 at 12:53 PM, Maarten Beek <span dir="ltr"><<a href="mailto:beekmaarten@yahoo.com">beekmaarten@yahoo.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<div><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>
const unsigned int Dimension = 3;<br> typedef itk::Matrix<double, 4, 4> MatrixType;<br> typedef itk::Image<MatrixType, Dimension> MatrixFieldType;<br> MatrixFieldType::Pointer matrixfield = MatrixFieldType::New();<br>
MatrixFieldType::IndexType index;<br> index.Fill( 0 );<br> MatrixFieldType::SizeType size;<br> size.Fill( 1 );<br> MatrixFieldType::RegionType region( index, size );<br> matrixfield->SetRegions( region );<br>
matrixfield->Allocate();<br> MatrixType matrix;<br> matrix.SetIdentity();<br> matrixfield->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>
</div><br>_____________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at<br>
<a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Kitware offers ITK Training Courses, for more information visit:<br>
<a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at:<br>
<a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
<br></blockquote></div><br>