template< class TImage >
void CreateImage(typename TImage::Pointer image)
{
typedef TImage ImageType;
typename ImageType::RegionType region(start,size);
image->SetRegions(region);
image->Allocate();
image->FillBuffer( -2 );
}
int main(int, char*[])
{
typedef double PixelType;
FilterType::Pointer filter = FilterType::New();
FilterType::Pointer filter2 =
dynamic_cast<FilterType*>(filter->CreateAnother().GetPointer());
ImageType::Pointer image = ImageType::New();
CreateImage< ImageType >(image);
filter2->SetInput(image);
filter2->Update();
std::cout << filter2->GetOutput()->GetPixel(index) << std::endl;
return EXIT_SUCCESS;
}