/*# ---------------------------------------------------------------------
#
# Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image
#                        pour la Sant�)
# Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
# Previous Authors : Laurent Guigues, Jean-Pierre Roux
# CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil
#
#  This software is governed by the CeCILL-B license under French law and
#  abiding by the rules of distribution of free software. You can  use,
#  modify and/ or redistribute the software under the terms of the CeCILL-B
#  license as circulated by CEA, CNRS and INRIA at the following URL
#  http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
#  or in the file LICENSE.txt.
#
#  As a counterpart to the access to the source code and  rights to copy,
#  modify and redistribute granted by the license, users are provided only
#  with a limited warranty  and the software's author,  the holder of the
#  economic rights,  and the successive licensors  have only  limited
#  liability.
#
#  The fact that you are presently reading this means that you have had
#  knowledge of the CeCILL-B license and that you accept its terms.
# ------------------------------------------------------------------------ */

#include "wxMaracasSurfaceRenderingManagerDataMhd.h"

#include "vtkStripper.h"


wxMaracasSurfaceRenderingManagerDataMhd::wxMaracasSurfaceRenderingManagerDataMhd(vtkImageData* imagedata, std::string dataname, vtkRenderWindowInteractor*  interactor)
: wxMaracasSurfaceRenderingManagerData(NULL, dataname)
{

	this->setVtkImageData(imagedata);
	_dataname = dataname;
	
	_maxgreylevel = getMaxLevel(imagedata);
	_prop3D=NULL;			

        _tresholdFilter = vtkImageThreshold::New();
        _tresholdFilter->SetInput(this->_imagedata);
        _tresholdFilter->SetInValue(255);        
        _tresholdFilter->SetOutValue(0);


	_cubesFilter = vtkMarchingCubes::New();	
        _cubesFilter->SetInput(_tresholdFilter->GetOutput());
        //_cubesFilter->ComputeGradientsOn();
        _cubesFilter->ComputeScalarsOn();
        _cubesFilter->ComputeNormalsOn();
        _cubesFilter->SetNumberOfContours( 1 );
	_cleanFilter = vtkCleanPolyData::New();		
	_cleanFilter->SetInput ( _cubesFilter->GetOutput() );

        _smooth = vtkSmoothPolyDataFilter::New();
        _smooth->SetInput(_cleanFilter->GetOutput());
        _smooth->SetNumberOfIterations(6);
        _smooth->SetRelaxationFactor(0.3);
        //_smooth->FeatureEdgeSmoothingOff();

	_dataMapper = vtkPolyDataMapper::New( );
	_dataMapper->ScalarVisibilityOff( );
	_dataMapper->ImmediateModeRenderingOn();
	vtkActor* dataActor = vtkActor::New();

	//if(_boxWidgetS1){
	if(interactor){
            _boxWidgetS1 = vtkBoxWidget::New();
            _boxWidgetS1->SetInteractor( interactor );
            _boxWidgetS1->SetPlaceFactor(1.25);

            _boxWidgetS1->SetInput( this->_imagedata );
            _boxWidgetS1->PlaceWidget();
            boxSurfaceObserver* observer = boxSurfaceObserver::New();

            vtkStripper* striper = vtkStripper::New();
            striper->SetInput( _smooth->GetOutput() );
            //striper->SetInput( _cleanFilter->GetOutput() );

            striper->Update();
            _boxWidgetS1->SetInput(striper->GetOutput());
            //_boxWidgetS1->PlaceWidget();


            _tissuePlanes  = vtkPlanes::New();

            int x1,x2,y1,y2,z1,z2;
            this->_imagedata->GetExtent(x1,x2,y1,y2,z1,z2);
            _tissuePlanes->SetBounds  (x1,x2,y1,y2,z1,z2);



            _boxWidgetS1->GetPlanes( _tissuePlanes );

            _tissueClipper = vtkClipPolyData::New();
            _tissueClipper->SetInput( striper->GetOutput() );
            _tissueClipper->SetClipFunction( _tissuePlanes );
            _tissueClipper->InsideOutOn( );
            _dataMapper->SetInput( _tissueClipper->GetOutput() );
            observer->SetPlanes( _tissuePlanes );
            observer->SetActor( dataActor );
            _boxWidgetS1->AddObserver( vtkCommand::InteractionEvent		 , observer );


            _boxWidgetS1->HandlesOn ();
            _boxWidgetS1->EnabledOff();
	}else{		
                //_dataMapper->SetInput(_cleanFilter->GetOutput());
            _dataMapper->SetInput(_smooth->GetOutput());
	}
	
	dataActor->SetMapper(_dataMapper);	
        float cr=1,cg=0.5,cb=0.5;
        dataActor->GetProperty()->SetDiffuseColor(1,0.5,0.5   );
        dataActor->GetProperty()->SetSpecular(.3);
        dataActor->GetProperty()->SetSpecularPower(20);

        this->_prop3D = dataActor;
	

	this->changeIsoValue(this->_maxgreylevel);	


}

void wxMaracasSurfaceRenderingManagerDataMhd::enableBoxWidget(bool enable){

    if(_boxWidgetS1){
        if(enable){
            _boxWidgetS1->EnabledOn();
        }else{
            _boxWidgetS1->EnabledOff();
        }
    }else{
        cout<<"box widget not initialized!"<<endl;
    }
}


wxMaracasSurfaceRenderingManagerDataMhd::~wxMaracasSurfaceRenderingManagerDataMhd()
{
	_cubesFilter->Delete();
	_cleanFilter->Delete();
	_dataMapper->Delete();
}
	
void wxMaracasSurfaceRenderingManagerDataMhd::UpdateSurface()
{
	_cubesFilter->Update();    
	_cleanFilter->Update();
	_dataMapper->Update();	
}
/**
**	changes the isovalue in a prop3D
**/
void wxMaracasSurfaceRenderingManagerDataMhd::changeIsoValue(double value){	
		
   
    _tresholdFilter->ThresholdByLower(value);
    _tresholdFilter->Update();

    _cubesFilter->SetValue(0,255);
	_cubesFilter->Update();    
	_cleanFilter->Update();
	_dataMapper->Update();	
	
	
	
}
/**
**	changes the isovalue in a prop3D
**/
void wxMaracasSurfaceRenderingManagerDataMhd::changeIsoValue(double min, double max){

    _tresholdFilter->ThresholdBetween(min, max);
    _tresholdFilter->Update();
    _cubesFilter->SetValue(0,255);
        _cubesFilter->Update();
        _cleanFilter->Update();
        _dataMapper->Update();



}
int wxMaracasSurfaceRenderingManagerDataMhd::getMaxGreyLevel(){
	return _maxgreylevel;
}

/**
	** Get's the max grey level of the image
	**/
int wxMaracasSurfaceRenderingManagerDataMhd::getMaxLevel(vtkImageData* img){

	int ext[6], i, j, k,max=0;
	img->GetExtent(ext);

	for(i = ext[0]; i < ext[1];i++){
		for(j = ext[2]; j < ext[3];j++){
			for(k = ext[4]; k < ext[5];k++){
                unsigned short* ptr = (unsigned short*)img->GetScalarPointer(i,j,k);
				int temp = (int)*ptr;
				if(temp > max){
                    max = temp;
				}
			}
		}
	}
	return max;


}

/**
**	Sets the VTK image data
**/

void wxMaracasSurfaceRenderingManagerDataMhd::setVtkImageData(vtkImageData* imagedata){
	_imagedata = imagedata;
}