X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=lib%2FmaracasVisuLib%2Fsrc%2Finterface%2FwxWindows%2Fwidgets%2FmanualContour%2FmanualContourModel.cpp;h=d5100a044ba3b880aa5d436f334bb74b12b982aa;hb=15112ae15da4222831b0e8e7c9be87daf9e99a1c;hp=f8369751bf63a9c61c893790852361d8776a308e;hpb=c6481486861a91e8ead2033566097c502488670e;p=creaMaracasVisu.git diff --git a/lib/maracasVisuLib/src/interface/wxWindows/widgets/manualContour/manualContourModel.cpp b/lib/maracasVisuLib/src/interface/wxWindows/widgets/manualContour/manualContourModel.cpp index f836975..d5100a0 100644 --- a/lib/maracasVisuLib/src/interface/wxWindows/widgets/manualContour/manualContourModel.cpp +++ b/lib/maracasVisuLib/src/interface/wxWindows/widgets/manualContour/manualContourModel.cpp @@ -1,3 +1,28 @@ +/*# --------------------------------------------------------------------- +# +# 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 "manualContourModel.h" // ---------------------------------------------------------------------------- @@ -9,53 +34,43 @@ manualContourModel::manualContourModel() _cntSplineX = vtkKochanekSpline::New( ); _cntSplineY = vtkKochanekSpline::New( ); _cntSplineZ = vtkKochanekSpline::New( ); - this->SetCloseContour(true); - _cntSplineX->SetDefaultTension( 0 ); _cntSplineX->SetDefaultBias( 0 ); _cntSplineX->SetDefaultContinuity( 0 ); - _cntSplineY->SetDefaultTension( 0 ); _cntSplineY->SetDefaultBias( 0 ); _cntSplineY->SetDefaultContinuity( 0 ); - _cntSplineZ->SetDefaultTension( 0 ); _cntSplineZ->SetDefaultBias( 0 ); _cntSplineZ->SetDefaultContinuity( 0 ); - //JSTG 25-02-08 ------------------------------------------------------------------------------------------------- - // this parameter is reset in the VIRTUAL manualContourBaseControler::Configure - _sizePointsContour = 100; //JSTG 25-02-08 The change in the inisialization of these variable is critical. - + _sizePointsContour = 500; //JSTG 25-02-08 The change in the inisialization of these variable is critical. _delta_JSTG = 0.0; //--------------------------------------------------------------------------------------------------------------- - //CMRU 17-08-09 ------------------------------------------------------------------------------------------------- _realSize = 0.0; - _label = ""; //--------------------------------------------------------------------------------------------------------------- - } // ---------------------------------------------------------------------------- manualContourModel::~manualContourModel() { int i,size=_lstPoints.size(); - for (i=0;iRemoveAllPoints(); _cntSplineY->RemoveAllPoints(); _cntSplineZ->RemoveAllPoints(); - _cntSplineX->Delete(); _cntSplineY->Delete(); _cntSplineZ->Delete(); } + // ---------------------------------------------------------------------------- int manualContourModel::AddPoint(double x,double y,double z) { @@ -63,7 +78,6 @@ int manualContourModel::AddPoint(double x,double y,double z) mp->SetPoint(x,y,z); AddManualPoint(mp); //UpdateSpline(); - return _lstPoints.size()-1; } // ---------------------------------------------------------------------------- @@ -162,6 +176,7 @@ void manualContourModel::DeleteAllPoints() void manualContourModel::MovePoint(int i,double dx,double dy,double dz) { +printf("EED manualContourModel::MovePoint %p \n", this); manualPoint *mp=_lstPoints[i]; double x=mp->GetX()+dx; double y=mp->GetY()+dy; @@ -312,7 +327,8 @@ void manualContourModel::UpdateSpline() // virtual _delta_JSTG = (double) (np) / double (_sizePointsContour - 1); //Without the -1 the curve is not close } else { // _delta_JSTG = (double) (np-1) / double (_sizePointsContour ); //Without the -1 the curve is not close - _delta_JSTG = (double) (np) / double (_sizePointsContour-1 ); //Without the -1 the curve is not close +//EED 9/7/2016 _delta_JSTG = (double) (np) / double (_sizePointsContour-1 ); //Without the -1 the curve is not close + _delta_JSTG = (double) (np-1) / double (_sizePointsContour-1 ); //Without the -1 the curve is not close } //----------------------------------------------------------------------------------------------------------- } @@ -339,7 +355,7 @@ void manualContourModel::GetSpline_i_Point(int i, double *x, double *y, double * //JSTG 25-02-08 --------------------------------------------------------------- void manualContourModel::GetSpline_t_Point(double t, double *x, double *y, double *z) { - if (_lstPoints.size()==0) + if (_lstPoints.size()==0) { *x = 0; *y = 0; @@ -362,8 +378,6 @@ void manualContourModel::GetSpline_t_Point(double t, double *x, double *y, doubl } // ---------------------------------------------------------------------------- - - /*void manualContourModel::GetSplinePoint(double t, double &x, double &y, double &z) { if (_lstPoints.size()==0) @@ -387,32 +401,30 @@ void manualContourModel::GetSpline_t_Point(double t, double *x, double *y, doubl z = _cntSplineZ->Evaluate(t); } }*/ + // ---------------------------------------------------------------------------- -double manualContourModel::GetPathSize() +double manualContourModel::GetPathSize(double *spc) { double result = 0; double x1,y1,z1; double x2,y2,z2; - // JSTG 25-02-08 ----------------------------- //double t,delta; //int i,np,nps; int i; //-------------------------------------------- - if (_lstPoints.size()==2) { - x1=_lstPoints[0]->GetX(); - y1=_lstPoints[0]->GetY(); - z1=_lstPoints[0]->GetZ(); - x2=_lstPoints[1]->GetX(); - y2=_lstPoints[1]->GetY(); - z2=_lstPoints[1]->GetZ(); + x1 = spc[0] * _lstPoints[0]->GetX(); + y1 = spc[1] * _lstPoints[0]->GetY(); + z1 = spc[2] * _lstPoints[0]->GetZ(); + x2 = spc[0] * _lstPoints[1]->GetX(); + y2 = spc[1] * _lstPoints[1]->GetY(); + z2 = spc[2] * _lstPoints[1]->GetZ(); result = sqrt( (x2-x1)*(x2-x1) + (y2-y1)*(y2-y1) + (z2-z1)*(z2-z1) ); } if (_lstPoints.size()>2) { - // JSTG 25-02-08 ------------------------------------------ //np = _lstPoints.size( ); //nps = 200; @@ -420,13 +432,18 @@ double manualContourModel::GetPathSize() UpdateSpline(); //GetSplinePoint(0,x1,y1,z1); GetSpline_i_Point(0,&x1,&y1,&z1); - + x1 = spc[0] * x1; + y1 = spc[1] * y1; + z1 = spc[2] * z1; //for( i = 1; i < nps; i++ ) for( i = 1; i < GetNumberOfPointsSpline(); i++ ) { //t = delta * (double)i; //GetSplinePoint(t,x2,y2,z2); GetSpline_i_Point(i,&x2,&y2,&z2); + x2 = spc[0] * x2; + y2 = spc[1] * y2; + z2 = spc[2] * z2; //--------------------------------------------------------- result=result + sqrt( (x2-x1)*(x2-x1) + (y2-y1)*(y2-y1) + (z2-z1)*(z2-z1) ); x1=x2; @@ -434,11 +451,11 @@ double manualContourModel::GetPathSize() z1=z2; }// for } - return result; } + // ---------------------------------------------------------------------------- -double manualContourModel::GetPathArea() +double manualContourModel::GetPathArea(double *spc) { double result = 555; if ((_lstPoints.size()>=3) && IfCloseContour()==true ) @@ -451,7 +468,6 @@ double manualContourModel::GetPathArea() double x2,y2,z2; bool okArea=true; int i, j; - // This uses Green's theorem: // A = 1/2 * sum( xiyi+1 - xi+1yi); pO == pN // A < 0 -> A = |A| (a negative value could raise because points are @@ -473,9 +489,13 @@ double manualContourModel::GetPathArea() GetSpline_i_Point(i,&x1,&y1,&z1); GetSpline_i_Point(j,&x2,&y2,&z2); //---------------------------------------------------------------- - area += - (x1 * y2 ) - - ( x2 * y1 ); + x1=spc[0]*x1; + y1=spc[1]*y1; + z1=spc[2]*z1; + x2=spc[0]*x2; + y2=spc[1]*y2; + z2=spc[2]*z2; + area += (x1 * y2 ) - ( x2 * y1 ); if (z1!=z2) { okArea=false; @@ -518,7 +538,6 @@ double manualContourModel::GetPathArea() // p[x,y,z] : data in // rp[x,y,z] : data out result point // rn[x,y,z] : data out result normal - void manualContourModel::GetNearestPointAndNormal(double *p, double *rp, double *rn) { double distMin=999999999; @@ -564,12 +583,11 @@ void manualContourModel::GetNearestPointAndNormal(double *p, double *rp, double rn[0] = x2-x1; rn[1] = y2-y1; rn[2] = z2-z1; - } - x2=x1; - y2=y1; - z2=z1; + } // if + x2 = x1; + y2 = y1; + z2 = z1; }// for - // JSTG 25-02-08 ------------------------------------------ //if (tback==0) if (iback==0) @@ -578,24 +596,22 @@ void manualContourModel::GetNearestPointAndNormal(double *p, double *rp, double //GetSplinePoint(t,x1,y1,z1); GetSpline_i_Point(i,&x1,&y1,&z1); //---------------------------------------------------------- - rn[0]=rp[0]-x1; - rn[1]=rp[1]-y1; - rn[2]=rp[2]-z1; - } - } - else - { + rn[0] = rp[0]-x1; + rn[1] = rp[1]-y1; + rn[2] = rp[2]-z1; + } // if iback + } else { rp[0] = 0; rp[1] = 0; rp[2] = 0; rn[0] = -1; rn[1] = 0; rn[2] = 0; - } + } // if np } // ---------------------------------------------------------------------------- -manualContourModel * manualContourModel :: Clone() // virtual +manualContourModel * manualContourModel::Clone() // virtual { manualContourModel * clone = new manualContourModel(); CopyAttributesTo(clone); @@ -655,9 +671,11 @@ void manualContourModel::SaveData(FILE *ff) { std::string etiqueta = GetLabel(); if(etiqueta.empty()) + { fprintf(ff,"Label: NO_LABEL\n"); - else + } else { fprintf(ff,"Label: %s\n",etiqueta.c_str()); + } // if etiqueta fprintf(ff,"Real_Size: %f\n",GetRealSize()); } @@ -670,8 +688,9 @@ void manualContourModel::OpenData(FILE *ff) std::cout<