const unsigned int Dimension = 3;

			typedef short    InputPixelType;
			typedef short    OutputPixelType;
	
			typedef itk::Image< InputPixelType,  Dimension >   InputImageType;
			typedef itk::Image< OutputPixelType, Dimension >   OutputImageType;
	
			
	
			cout<<"1"<<endl;
			
	
	
			typedef itk::Image< unsigned char,  3 >  UnsignedCharImageType;
			typedef itk::Image< float,  3 >   FloatImageType;
	
	
			vtkImageData * imageVtk = image_view->GetInput();
		    
			typedef itk::VTKImageToImageFilter<InputImageType> vtk2itkType;
			vtk2itkType::Pointer vtk2itk = vtk2itkType::New();
			vtk2itk->SetInput(imageVtk);
			vtk2itk->Update();
		     
			cout<<"2"<<endl;
	
			
			InputImageType::Pointer image = vtk2itk->GetOutput();
			InputImageType::RegionType inputRegion = image->GetLargestPossibleRegion();
		
			int nx = inputRegion.GetSize(0);
			int ny = inputRegion.GetSize(1);
			int nz = inputRegion.GetSize(2);
			

			 typedef itk::CoherenceEnhancingDiffusionImageFilter<InputImageType> DiffusionFilterType;
			 DiffusionFilterType::Pointer diffusionFilter = DiffusionFilterType::New();
			 diffusionFilter->SetInput(image);
			 const double diffusionTime = 2;
			 cout<<"3"<<endl;
			 diffusionFilter->SetDiffusionTime(diffusionTime);
			 const double lambda = 0.0001;
			 diffusionFilter->SetLambda(lambda);
			 const char * enhancement = "Isotropic";
			 cout<<"4"<<endl;
             diffusionFilter->SetEnhancement(DiffusionFilterType::Isotropic);
			 const double noiseScale = 1;
			 diffusionFilter->SetNoiseScale(noiseScale);
			 const double featureScale = 2;
			 cout<<"5"<<endl;
			 diffusionFilter->SetFeatureScale(featureScale);
			 const double exponent = 0.8;
			 diffusionFilter->SetExponent(exponent);
			 diffusionFilter->Update();

			 cout<<"6"<<endl;
    


			prepareImageViews();
			renderViewportsItk(diffusionFilter->GetOutput());