- /* TODO
- std::cout << "6" << std::endl;
- { // begin
- TOutputImage* output = this->GetOutput( );
- output->SetBufferedRegion( output->GetRequestedRegion( ) );
- output->Allocate( );
- output->FillBuffer( -std::numeric_limits< _TScalar >::max( ) );
-
- itk::ImageRegionConstIterator< TOutputImage > rIt(
- output_roi, output_roi->GetRequestedRegion( )
- );
- itk::ImageRegionIterator< TOutputImage > oIt( output, roi );
- rIt.GoToBegin( );
- oIt.GoToBegin( );
- for( ; !rIt.IsAtEnd( ); ++rIt, ++oIt )
- oIt.Set( rIt.Get( ) );
-
- TPath* output_axis = this->GetOutputAxis( );
- output_axis->SetReferenceImage( output );
- for( unsigned long i = 0; i < output_axis_roi->GetSize( ); ++i )
- {
- TIndex v = output_axis_roi->GetVertex( i );
- for( unsigned int d = 0; d < TInputImage::ImageDimension; ++d )
- v[ d ] += roi.GetIndex( )[ d ];
- output_axis->AddVertex( v );
-
- } // rof
-
- } // end
- std::cout << "7" << std::endl;
- */
+ std::cout << "output" << std::endl;
+ { // begin
+ TOutputImage* output = this->GetOutput( );
+ output->SetBufferedRegion( output->GetRequestedRegion( ) );
+ output->Allocate( );
+ output->FillBuffer( -std::numeric_limits< _TScalar >::max( ) );
+
+ itk::ImageRegionConstIterator< TOutputImage > rIt(
+ out_dist, out_dist->GetRequestedRegion( )
+ );
+ itk::ImageRegionIterator< TOutputImage > oIt( output, roi );
+ rIt.GoToBegin( );
+ oIt.GoToBegin( );
+ for( ; !rIt.IsAtEnd( ); ++rIt, ++oIt )
+ oIt.Set( rIt.Get( ) );
+
+ TPath* output_axis = this->GetOutputAxis( );
+ output_axis->SetReferenceImage( output );
+ for( unsigned long i = 0; i < out_axis->GetSize( ); ++i )
+ {
+ TIndex v = out_axis->GetVertex( i );
+ for( unsigned int d = 0; d < TInputImage::ImageDimension; ++d )
+ v[ d ] += roi.GetIndex( )[ d ];
+ output_axis->AddVertex( v );
+
+ } // rof
+
+ } // end