Hi all,<br><br>I&#39;m trying to create an array of itkImages.<br><br>This is what i got so far:<br><br><br><br>int numLabels = 2;<br><br>  ImageType::Pointer* imageMatrix = new ImageType::Pointer[numLabels];<br> <br>  for( unsigned int i=0;i&lt;2;i++){<br>

    imageMatrix[i]= ImageType::New();<br>    imageMatrix[i]-&gt;SetSpacing(sp);<br>    imageMatrix[i]-&gt;SetOrigin(origin);<br>    imageMatrix[i]-&gt;SetDirection(direction);<br>    imageMatrix[i]-&gt;SetLargestPossibleRegion(image-&gt;GetLargestPossibleRegion());<br>

  }<br><br>  ImageType::SizeType sizeOfImage2 = imageMatrix[0]-&gt;GetLargestPossibleRegion().GetSize();<br>  std::cout &lt;&lt; &quot;Size2 &quot;;<br>  std::cout &lt;&lt; sizeOfImage2[0] &lt;&lt; &quot;, &quot; &lt;&lt; sizeOfImage2[1] &lt;&lt; &quot;, &quot; &lt;&lt; sizeOfImage2[2] &lt;&lt; std::endl &lt;&lt;std::endl;<br>

<br> // snipped a bit here that gets an indexValue of [0, 0, 0 ]<br><br>  typedef itk::Index&lt;3&gt; indexType;<br>  indexType blah2;<br>  blah2 = inputIt.GetIndex();<br>  std::cout &lt;&lt; &quot;index : &quot; &lt;&lt; blah2 &lt;&lt; std::endl;<br>

<br>  image-&gt;SetPixel(blah2,128);<br>  imageMatrix[0]-&gt;SetPixel(blah2,128);<br>  std::cout &lt;&lt; &quot;pixel: &quot; &lt;&lt; image-&gt;GetPixel(blah2) &lt;&lt; std::endl;<br><br><br><br>This is my output:<br><br>

Size2 512, 512, 52<br><br>here 2<br>index : [0, 0, 0]<br>Segmentation fault<br><br><br><br>So basically, when I&#39;m accessing the image volume, the program crashes.<br><br>Can anybody help me with this?<br><br>cheers, Michael<br>