manualViewPerpPlaneContour.cxx

Go to the documentation of this file.

#include "manualViewPerpPlaneContour.h"

#include "vtkPlane2DView.h"

manualViewPerpPlaneContour::manualViewPerpPlaneContour()
{
}
// ----------------------------------------------------------------------------
manualViewPerpPlaneContour::~manualViewPerpPlaneContour()
{
}


// ----------------------------------------------------------------------------
manualViewPerpPlaneContour * manualViewPerpPlaneContour :: Clone()
{
        manualViewPerpPlaneContour * clone = new manualViewPerpPlaneContour();
        CopyAttributesTo(clone);
        return clone;
}

// ---------------------------------------------------------------------------

void manualViewPerpPlaneContour::CopyAttributesTo( manualViewPerpPlaneContour * cloneObject)
{
        // Fathers object
        manualViewContour::CopyAttributesTo(cloneObject);
}




// ----------------------------------------------------------------------------
void manualViewPerpPlaneContour::UpdateViewPoint(int id)
{  // virtual
        double x,y,z;
        manualPoint *mp = _manContModel->GetManualPoint(id);
        x = mp->GetX();
        y = mp->GetY();
        z = mp->GetZ();
        FilterCordinateXYZ(x,y,z);

        _lstViewPoints[id]->SetPositionXY( x , y ,GetRange(), z );
        if ((z>=-1) && (z<=1))
        {
                _lstViewPoints[id]->GetVtkActor()->VisibilityOn();
        } else {
                _lstViewPoints[id]->GetVtkActor()->VisibilityOff();
        }
}
// ----------------------------------------------------------------------------
void manualViewPerpPlaneContour::FilterCordinateXYZ (double &x, double &y, double &z)
{
        vtkPlane2DView *vtkplane2Dview = (vtkPlane2DView*)( this->GetWxVtkBaseView() );
        vtkMPRBaseData *vtkmprbasedata = vtkplane2Dview->GetVtkmprbasedata();
        vtkTransform *transf1 = vtkTransform::New();
        vtkTransform *transf2 = vtkmprbasedata->GetTransformOrientation();
        transf1->SetMatrix( transf2->GetMatrix() ); 
        transf1->Inverse();

        double in[4], out[4];
        double center = vtkplane2Dview->GetImgSize() / 2;

        in[0] = x - vtkmprbasedata->GetX();
        in[1] = y - vtkmprbasedata->GetY();
        in[2] = z - vtkmprbasedata->GetZ();
        in[3] = 0;

        transf1->MultiplyPoint(in,out);
        z = out[0];
        x = out[1]+center;
        y = out[2]+center;

        double spc[3];
        this->GetWxVtkBaseView()->GetSpacing(spc);
        x = x / spc[0];
        y = y / spc[1];
        z = z / spc[2];

        transf1->Delete();
}

// ----------------------------------------------------------------------------
void manualViewPerpPlaneContour::TransfromCoordViewWorld( double &X,double &Y,double &Z,int type )
{
        Z = 0;
        vtkPlane2DView *vtkplane2Dview = (vtkPlane2DView*)( this->GetWxVtkBaseView() );
        vtkplane2Dview->TransfromCoordViewWorld2(X,Y,Z);
}


// ----------------------------------------------------------------------------
void manualViewPerpPlaneContour::RefreshContour() // virtual
{
//      manualViewContour::RefreshContour();

        double pp1[3];
        double pp2[3];
        double u;
                
// JSTG 25-02-08 ----------------------------------------------------------
//              tt;
//--------------------------------------------------------

        double pp[3];
        double ppB[3];
        double ppC[3];
        ppB[0] = 999999;
        double dist,distMin = 99999999;

        vtkPlane2DView *vtkplane2Dview = (vtkPlane2DView*)( this->GetWxVtkBaseView() );
        double center = vtkplane2Dview->GetImgSize() / 2;

        int i,np,nps;
        np              = GetNumberOfPoints( );

//JSTG 25-02-08 ------------------------------------------
        //double t,delta;
    //nps               = GetNumberOfPointsSpline(); 
        nps = _manContModel->GetNumberOfPointsSpline();
        //delta = ( double ) ( np  ) / ( double ) ( nps-1  );
        _manContModel->UpdateSpline();
//--------------------------------------------------------

        if ( np >= 2 )
        {
                for( i = 0; i < nps; i++ ) 
                {
// JSTG 25-02-08 ----------------------------------------------------------
                        //t     = delta * (double)i ;
                        //tt    = delta * (double)(i+1) ;
                        //_manContModel->GetSplinePoint(t ,pp1[0],pp1[1],pp1[2]);
                        //_manContModel->GetSplinePoint(tt,pp2[0],pp2[1],pp2[2]);
                        _manContModel->GetSpline_i_Point(i ,&pp1[0],&pp1[1],&pp1[2]);
                        _manContModel->GetSpline_i_Point(i+1,&pp2[0],&pp2[1],&pp2[2]);
//--------------------------------------------------------------------------
                        FilterCordinateXYZ(pp1[0],pp1[1],pp1[2]);
                        FilterCordinateXYZ(pp2[0],pp2[1],pp2[2]);
                        if (pp2[2]*pp1[2]<=0)
                        {
                                if (pp1[2]-pp2[2]!=0) { 
                                        u = -pp2[2] / (pp1[2]-pp2[2]);
                                } else  {
                                        u=9999999;
                                }
                                pp[0]    = ( pp1[0]-pp2[0] )*u + pp2[0];
                                pp[1]    = ( pp1[1]-pp2[1] )*u + pp2[1];
                                pp[2]    = ( pp1[2]-pp2[2] )*u + pp2[2];
                                ppC[0] = pp[0] - center;
                                ppC[1] = pp[1] - center;
                                ppC[2] = pp[2] ;
                                dist = sqrt( ppC[0]*ppC[0] + ppC[1]*ppC[1] + ppC[2]*ppC[2] );
                                if (dist<distMin)
                                {
                                        distMin=dist;
                                        ppB[0] = pp[0];
                                        ppB[1] = pp[1];
                                        ppB[2] = pp[2];
                                }
                        }

                }// for 
        } else {
                        _pts->SetPoint(0, 0 , 0 , 0);   
                        _pts->SetPoint(1, 0 , 0 , 0);   
        } // if

//EED 27 sep 2006
        ppB[0]=ppB[0]*_spc[0];
        ppB[1]=ppB[1]*_spc[1];

        _pts->SetPoint( 0, ppB[0]   , ppB[1]+1 , 1 );   
        _pts->SetPoint( 1, ppB[0]   , ppB[1]-1 , 1 );   
        _pts->SetPoint( 2, ppB[0]   , ppB[1]   , 1 );   
        _pts->SetPoint( 3, ppB[0]-1 , ppB[1]   , 1 );   
        _pts->SetPoint( 4, ppB[0]+1 , ppB[1]   , 1 );   


        for( i = 5; i < nps; i++ ) 
        {
                _pts->SetPoint( i, ppB[0] , ppB[1] , 1 );       
        }

}

//---------------------------------------------------------------------------
bool manualViewPerpPlaneContour::ifTouchContour( int x, int y, int z)
{
        bool ok=false;
        vtkPlane2DView *vtkplane2Dview = (vtkPlane2DView*)( this->GetWxVtkBaseView() );
        double X=x;
        double Y=y;
        double Z=z;
        vtkplane2Dview->TransFromCoordScreenToWorld(X,Y,Z);
        
        double ppA[3];
        _pts->GetPoint(0, ppA);
        if (sqrt( (ppA[0]-X)*(ppA[0]-X) + (ppA[1]-Y)*(ppA[1]-Y) ) <=2)
        {
                ok = true;
        }
        return ok;
}