--- /dev/null
+
+#include <cmath>
+
+#include "ModelShowNPoints.h"
+
+
+//----------------------------------------------------------------------
+ModelShowNPoints::ModelShowNPoints()
+{
+ firsttime=true;
+}
+//----------------------------------------------------------------------
+
+ModelShowNPoints::~ModelShowNPoints()
+{
+}
+
+//------------------------------------------------------------------------
+void ModelShowNPoints::SetRadio(double radio)
+{
+ mradio = radio;
+}
+
+
+//------------------------------------------------------------------------
+double ModelShowNPoints::GetRadio()
+{
+ return mradio;
+}
+
+
+//------------------------------------------------------------------------
+std::vector<int> ModelShowNPoints::GetLstPointsX()
+{
+ return lstPointsX;
+}
+
+//------------------------------------------------------------------------
+std::vector<int> ModelShowNPoints::GetLstPointsY()
+{
+ return lstPointsY;
+}
+
+//------------------------------------------------------------------------
+std::vector<int> ModelShowNPoints::GetLstPointsZ()
+{
+ return lstPointsZ;
+}
+
+//------------------------------------------------------------------------
+std::vector<std::string> ModelShowNPoints::GetLstLabels()
+{
+ return lstLabels;
+}
+
+//------------------------------------------------------------------------
+void ModelShowNPoints::SetReferencePoint(std::vector<int> ppoint)
+{
+ mReferencePoint = ppoint;
+}
+
+//------------------------------------------------------------------------
+std::vector<int> ModelShowNPoints::GetReferencePoint()
+{
+ return mReferencePoint;
+}
+
+//------------------------------------------------------------------------
+void ModelShowNPoints::SetImage(vtkImageData *image)
+{
+ this->mimage=image;
+}
+
+
+//------------------------------------------------------------------------
+std::string ModelShowNPoints::CleanSpaces(std::string ss)
+{
+ int i;
+ while( (i=ss.find(32))>=0 )
+ {
+ ss.replace(i,1,"_");
+ }
+ return ss;
+}
+
+//------------------------------------------------------------------------
+void ModelShowNPoints::GetIdPoint(int id, int *x, int *y, int *z)
+{
+ *x=lstPointsX[id];
+ *y=lstPointsY[id];
+ *z=lstPointsZ[id];
+}
+
+//------------------------------------------------------------------------
+std::string ModelShowNPoints::GetIdLabel(int id)
+{
+ return lstLabels[id];
+}
+
+//------------------------------------------------------------------------
+vtkImageData *ModelShowNPoints::GetImage()
+{
+ return mimage;
+}
+
+//------------------------------------------------------------------------
+void ModelShowNPoints::AddPoint(int x, int y, int z, std::string label)
+{
+ lstPointsX.push_back( x );
+ lstPointsY.push_back( y );
+ lstPointsZ.push_back( z );
+ std::string strLabel = CleanSpaces( label );
+ lstLabels.push_back( strLabel );
+}
+
+//------------------------------------------------------------------------
+double ModelShowNPoints::Distance(double dX0, double dY0, double dZ0, double dX1, double dY1, double dZ1)//CFT
+{
+ return sqrt((dX1 - dX0)*(dX1 - dX0) + (dY1 - dY0)*(dY1 - dY0) + (dZ1 - dZ0)*(dZ1 - dZ0));
+}
+
+//------------------------------------------------------------------------
+int ModelShowNPoints::InsertPoint(int x, int y, int z, std::string label)
+{
+ if(lstPointsX.size()>1)
+ {
+ std::vector<int> dTotal;
+ int pos = 1;
+ int a,b,res;
+
+ //Calcule distance for each pair of points
+ for(int i = 0; i<(int)lstPointsX.size()-1 ; i++)
+ {
+ a = Distance(x, y, z, lstPointsX[i], lstPointsY[i], lstPointsZ[i]);
+ b = Distance(x, y, z, lstPointsX[i+1], lstPointsY[i+1], lstPointsZ[i+1]);
+ res = a + b;
+ dTotal.push_back (res);
+ }
+ //Gets the smallest distance
+ int smallTMP = dTotal[0];
+ for (int j = 0; j < (int) dTotal.size(); j++)
+ {
+ if(dTotal[j]<smallTMP)
+ {
+ smallTMP=dTotal[j];
+ pos = j+1;
+ }
+ }
+
+ std::vector<int>::iterator it;
+ //Insert the point in the list of points
+ it = lstPointsX.begin();
+ lstPointsX.insert( it+pos, x );
+ it = lstPointsY.begin();
+ lstPointsY.insert( it+pos, y );
+ it = lstPointsZ.begin();
+ lstPointsZ.insert( it+pos, z );
+
+ std::string strLabel = CleanSpaces( label );
+
+ std::vector<std::string>::iterator itS;
+ itS = lstLabels.begin();
+ //Insert Label in list of labels
+ lstLabels.insert( itS+pos, strLabel );
+ return pos;
+ } else {
+ return -1;
+ }// if size lst X
+
+}
+
+
+//------------------------------------------------------------------------
+void ModelShowNPoints::SavePoints(std::string filename)
+{
+ std::string tmpLabel;
+ FILE *ff;
+ ff = fopen( filename.c_str() , "w+" );
+ if (ff!=NULL)
+ {
+ int i , size = (int) (lstPointsX.size());
+ fprintf(ff,"NumberOfPoints %d \n",size);
+ fprintf(ff," X\tY\tZ\tvalue\tLabel\n");
+ int x, y, z;
+ double value;
+ for (i=0; i<size; i++)
+ {
+ x=lstPointsX[i];
+ y=lstPointsY[i];
+ z=lstPointsZ[i];
+ value= mimage->GetScalarComponentAsDouble(x,y,z,0);
+ if (lstLabels[i]!="")
+ {
+ tmpLabel=lstLabels[i];
+ } else{
+ tmpLabel="<_VOID_>";
+ }
+ fprintf(ff,"%d\t%d\t%d\t%f\t%s\n", x , y , z , value , tmpLabel.c_str());
+ } // for
+ fclose(ff);
+ } else { // else ff
+ printf("ModelShowNPoints::SavePoints ...Error... creating file");
+ } //ff
+}
+
+//------------------------------------------------------------------------
+int ModelShowNPoints::ReadPoints(std::string filename)
+{
+ int i,size;
+ char chartmp[256];
+ FILE *ff;
+ ff = fopen( filename.c_str() , "r+" );
+ if (ff!=NULL)
+ {
+ fscanf(ff," %s %d",chartmp,&size);
+ fscanf(ff," %s %s %s %s %s",chartmp, chartmp,chartmp,chartmp,chartmp );
+
+ float value;
+ int x,y,z;
+ for (i=0; i<size; i++)
+ {
+ fscanf(ff,"%d%d%d%f%s",&x,&y,&z,&value,chartmp ); // x,y,z,value,label
+ if (strcmp(chartmp,"<_VOID_>")==0) { strcpy(chartmp,""); }
+ AddPoint(x,y,z, chartmp );
+ }
+ fclose(ff);
+ } else { // else ff
+ printf("ModelShowNPoints::LoadPoints ...Error... reading file");
+ } //ff
+ return size;
+}
+
+//------------------------------------------------------------------------
+int ModelShowNPoints::GetNearestPoint()
+{
+ int id=-1;
+ int i, size=(int)(lstPointsX.size());
+ double radioMin=10000000;
+ for ( i=0 ; i<size; i++ )
+ {
+ double rx = mReferencePoint[0] - lstPointsX [i];
+ double ry = mReferencePoint[1] - lstPointsY [i];
+ double rz = mReferencePoint[2] - lstPointsZ [i];
+ double radio = rx*rx + ry*ry + rz*rz;
+ if ( radio <= radioMin)
+ {
+ radioMin=radio;
+ id=i;
+ } // if
+ } // for
+ return id;
+}
+
+//------------------------------------------------------------------------
+int ModelShowNPoints::GetLstPointsSize()
+{
+ return lstPointsX.size();
+}
+
+//------------------------------------------------------------------------
+void ModelShowNPoints::SetPointId_mReferencePoint(int id)
+{
+ lstPointsX[id] = mReferencePoint[0];
+ lstPointsY[id] = mReferencePoint[1];
+ lstPointsZ[id] = mReferencePoint[2];
+}
+
+
+//------------------------------------------------------------------------
+int ModelShowNPoints::IdInsidePoint()
+{
+ int id=-1;
+ int i, size=(int)(lstPointsX.size());
+ double spc[3];
+ if(mimage ==NULL)
+ {
+ printf("WidgetShowNPoints::IdInsidePoint image not set\n");
+ return -1;
+ }else{
+ mimage->GetSpacing(spc);
+ for ( i=0 ; i<size; i++ )
+ {
+ double rx = spc[0]*(mReferencePoint[0] - lstPointsX [i]);
+ double ry = spc[1]*(mReferencePoint[1] - lstPointsY [i]);
+ double rz = spc[2]*(mReferencePoint[2] - lstPointsZ [i]);
+ if ( rx*rx + ry*ry + rz*rz <= mradio*mradio)
+ {
+ id=i;
+ } // if
+ } // for
+ return id;
+ } // if
+}
+
+//------------------------------------------------------------------------
+int ModelShowNPoints::RenamePoint(std::string label)
+{
+ int id=IdInsidePoint();
+ if (id>=0)
+ {
+ std::string strLabel = CleanSpaces( label );
+ lstLabels[id] = strLabel;
+ }
+ return id;
+}
+
+//----------------------------------------------------------------------
+void ModelShowNPoints::ErasePoint(int id)
+{
+ lstPointsX.erase( lstPointsX.begin()+id );
+ lstPointsY.erase( lstPointsY.begin()+id );
+ lstPointsZ.erase( lstPointsZ.begin()+id );
+ lstLabels.erase( lstLabels.begin()+id );
+}
+
+//----------------------------------------------------------------------
+void ModelShowNPoints::SetFirstTime(bool value)
+{
+ firsttime=value;
+}
+
+//----------------------------------------------------------------------
+bool ModelShowNPoints::GetFirstTime()
+{
+ return firsttime;
+}
+
+
git://git.creatis.insa-lyon.fr / creaMaracasVisu.git / blobdiff