Merge pull request #53 from blowekamp/FixVariableLengthVectorInitialization

Fix variable length vector initialization

Author: thewtex

include/itkCoocurrenceTextureFeaturesImageFilter.hxx

@@ -140,6 +140,7 @@ CoocurrenceTextureFeaturesImageFilter<TInputImage, TOutputImage, TMaskImage>
 
   // Creation of the output pixel type
   typename TOutputImage::PixelType outputPixel;
+  NumericTraits<typename TOutputImage::PixelType>::SetLength(outputPixel, outputPtr->GetNumberOfComponentsPerPixel());
 
   // Separation of the non-boundary region that will be processed in a different way
   NeighborhoodAlgorithm::ImageBoundaryFacesCalculator< DigitizedImageType > boundaryFacesCalculator;
@@ -149,12 +150,6 @@ CoocurrenceTextureFeaturesImageFilter<TInputImage, TOutputImage, TMaskImage>
 
   // Declaration of the variables useful to iterate over the all image region
   bool isInImage;
-  IndexType firstIndex;
-  for ( unsigned int i = 0; i < this->m_NeighborhoodRadius.Dimension; ++i )
-    {
-    firstIndex[i] = 0;
-    }
-  outputPixel = outputPtr->GetPixel(firstIndex);
   typename OffsetVector::ConstIterator offsets;
 
   // Declaration of the variables useful to iterate over the all the offsets

include/itkRunLengthTextureFeaturesImageFilter.hxx

@@ -136,15 +136,15 @@ RunLengthTextureFeaturesImageFilter<TInputImage, TOutputImage, TMaskImage>
                        ThreadIdType threadId)
 {
   // Get the inputs/outputs
-  typename TOutputImage::Pointer outputPtr = TOutputImage::New();
-  outputPtr = this->GetOutput();
+  TOutputImage * outputPtr = this->GetOutput();
 
   ProgressReporter progress( this,
                              threadId,
                              outputRegionForThread.GetNumberOfPixels() );
 
   // Creation of the output pixel type
   typename TOutputImage::PixelType outputPixel;
+  NumericTraits<typename TOutputImage::PixelType>::SetLength(outputPixel, outputPtr->GetNumberOfComponentsPerPixel());
 
   // Creation of a region with the same size as the neighborhood. This region
   // will be used to check if each voxel has already been visited.
@@ -174,7 +174,6 @@ RunLengthTextureFeaturesImageFilter<TInputImage, TOutputImage, TMaskImage>
 
   // Declaration of the variables useful to iterate over the all image region
   bool isInImage;
-  outputPixel = outputPtr->GetPixel(boolCurentInNeighborhoodIndex);
   typename OffsetVector::ConstIterator offsets;
 
   // Declaration of the variables useful to iterate over the all the offsets