Authors belong to:
- University of LYON http://www.universite-lyon.fr/
- - Léon Bérard cancer center http://oncora1.lyon.fnclcc.fr
+ - Léon Bérard cancer center http://www.centreleonberard.fr
- CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr
This software is distributed WITHOUT ANY WARRANTY; without even
- BSD See included LICENSE.txt file
- CeCILL-B http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
- ======================================================================-====*/
+ ===========================================================================**/
#include "vvLandmarks.h"
#include <ios>
#include "vtkPointData.h"
#include "clitkCommon.h"
+//--------------------------------------------------------------------
vvLandmarks::vvLandmarks(int size)
{
mLandmarks.resize(0);
}
mPolyData = vtkPolyData::New();
mIds = vtkFloatArray::New();
+ mLabels = vtkStringArray::New();
+ mLabels->SetName("labels");
}
+//--------------------------------------------------------------------
+
+//--------------------------------------------------------------------
vvLandmarks::~vvLandmarks()
{
for (unsigned int i = 0; i < mPoints.size(); i++) {
mIds->Delete();
if (mPolyData)
mPolyData->Delete();
+ if (mLabels)
+ mLabels->Delete();
}
+//--------------------------------------------------------------------
+
+//--------------------------------------------------------------------
void vvLandmarks::AddLandmark(float x,float y,float z,float t,double value)
{
vvLandmark point;
+ vtkIdType idPoint;
point.coordinates[0] = x;
point.coordinates[1] = y;
point.coordinates[2] = z;
point.coordinates[3] = t;
point.pixel_value=value;
mLandmarks.push_back(point);
- mPoints[int(t)]->InsertNextPoint(x,y,z);
-
- /*std::stringstream numberVal;
- numberVal << (mLandmarks.size()-1);
- vvLandmarksGlyph *number = vvLandmarksGlyph::New();
- number->SetText(numberVal.str().c_str());
- number->BackingOff();
- mText.push_back(number);*/
+
+ idPoint = mPoints[int(t)]->InsertNextPoint(x,y,z);
+ std::string str_vtkIdType; // string which will contain the result
+ std::ostringstream convert; // stream used for the conversion
+ convert << idPoint; // insert the textual representation of 'idPoint' in the characters in the stream
+ str_vtkIdType = convert.str(); // set 'str_vtkIdType' to the contents of the stream
+ mLabels->InsertNextValue(str_vtkIdType.c_str());
+
+ std::stringstream numberVal;
+ numberVal << (mLandmarks.size()-1);
+ /*
+ vvLandmarksGlyph *number = vvLandmarksGlyph::New();
+ number->SetText(numberVal.str().c_str());
+ number->BackingOff();
+ DD(numberVal.str().c_str());
+ mText.push_back(number);
+ */
+
mIds->InsertNextTuple1(0.55);
//mIds->InsertTuple1(mLandmarks.size(),mLandmarks.size());
SetTime(int(t));
}
+//--------------------------------------------------------------------
+
+//--------------------------------------------------------------------
void vvLandmarks::RemoveLastLandmark()
{
mPoints[mLandmarks.back().coordinates[3]]->SetNumberOfPoints(
mPoints[mLandmarks.back().coordinates[3]]->GetNumberOfPoints()-1);
- mPolyData->Modified();
- //mText.pop_back();
+ // mText.pop_back();
mLandmarks.pop_back();
mIds->RemoveLastTuple();
+ mLabels->SetNumberOfValues(mLabels->GetNumberOfValues()-1);
+ mLabels->Modified();
+ mPolyData->Modified();
+}
+//--------------------------------------------------------------------
+
+
+//--------------------------------------------------------------------
+void vvLandmarks::RemoveLandmark(int index)
+{
+ // erase a vtkPoint by shifiting the array .
+ // not a problem here because there are no
+ // pologyons linking the points
+ int npoints = mPoints[mLandmarks[index].coordinates[3]]->GetNumberOfPoints();
+ int t = mLandmarks[index].coordinates[3];
+ for (int i = index; i < npoints - 1; i++) {
+ mPoints[t]->InsertPoint(i, mPoints[t]->GetPoint(i+1));
+ std::string str_i; // string which will contain the result
+ std::ostringstream convert; // stream used for the conversion
+ convert << i; // insert the textual representation of 'i' in the characters in the stream
+ str_i = convert.str(); // set 'str_i' to the contents of the stream
+ mLabels->SetValue(i,str_i.c_str());
+ }
+ mPoints[t]->SetNumberOfPoints(npoints-1);
+ mLabels->SetNumberOfValues(npoints-1);
+ mLabels->Modified();
+ mPolyData->Modified();
+
+ mLandmarks.erase(mLandmarks.begin() + index);
+ mIds->RemoveLastTuple();
}
+//--------------------------------------------------------------------
+
+//--------------------------------------------------------------------
void vvLandmarks::ChangeComments(int index, std::string comments)
{
mLandmarks[index].comments = comments;
}
+//--------------------------------------------------------------------
+
+//--------------------------------------------------------------------
double vvLandmarks::GetPixelValue(int index)
{
return mLandmarks[index].pixel_value;
}
+//--------------------------------------------------------------------
+
+//--------------------------------------------------------------------
float* vvLandmarks::GetCoordinates(int index)
{
return mLandmarks[index].coordinates;
}
+//--------------------------------------------------------------------
+
+//--------------------------------------------------------------------
std::string vvLandmarks::GetComments(int index)
{
return mLandmarks[index].comments;
}
+//--------------------------------------------------------------------
+
+//--------------------------------------------------------------------
void vvLandmarks::LoadFile(std::string filename)
{
std::ifstream fp(filename.c_str(), std::ios::in|std::ios::binary);
}
mFilename = filename;
mLandmarks.clear();
+ vtkIdType idPoint;
char line[255];
for (unsigned int i = 0; i < mPoints.size(); i++)
mPoints[i]->SetNumberOfPoints(0);
bool first_line=true;
while (fp.getline(line,255)) {
+ // DD(line);
std::string stringline = line;
if (first_line) {
first_line=false;
continue;
}
point.coordinates[0] = atof(replace_dots(stringline.substr(previousSpace,space - previousSpace)).c_str());
+ // DD(point.coordinates[0]);
previousSpace = space;
space = stringline.find(" ", previousSpace+1);
if (space < -1 || space > (int)stringline.size()) {
continue;
}
point.coordinates[1] = atof(replace_dots(stringline.substr(previousSpace,space - previousSpace)).c_str());
+ // DD(point.coordinates[1]);
previousSpace = space;
space = stringline.find(" ", previousSpace+1);
if (space < -1 || space > (int)stringline.size()) {
continue;
}
point.pixel_value = atof(replace_dots(stringline.substr(previousSpace,space - previousSpace)).c_str());
+ // DD(point.pixel_value);
} else {
point.pixel_value=0.; //Not in file
point.coordinates[3]=0.;
previousSpace = space;
//this is the maximum size of comments
space = (stringline.find("\n", previousSpace+1) < 254 ? stringline.find("\n", previousSpace+1) : 254);
- point.comments = stringline.substr(previousSpace,space - (previousSpace)).c_str();
+ if (previousSpace != -1) {
+ point.comments = stringline.substr(previousSpace,space - (previousSpace)).c_str();
+ }
+ // DD(point.comments);
mLandmarks.push_back(point);
mIds->InsertNextTuple1(0.55);
- mPoints[int(point.coordinates[3])]->InsertNextPoint(
+ idPoint = mPoints[int(point.coordinates[3])]->InsertNextPoint(
point.coordinates[0],point.coordinates[1],point.coordinates[2]);
+ std::string str_vtkIdType; // string which will contain the result
+ std::ostringstream convert; // stream used for the conversion
+ convert << idPoint; // insert the textual representation of 'idPoint' in the characters in the stream
+ str_vtkIdType = convert.str(); // set 'str_vtkIdType' to the contents of the stream
+ mLabels->InsertNextValue(str_vtkIdType.c_str());
}
}
SetTime(0);
}
+//--------------------------------------------------------------------
+
+//--------------------------------------------------------------------
bool vvLandmarks::ErrorMsg(int num,const char * text)
{
std::cerr << "error when loading point " << num << " at " << text << std::endl;
return false;
}
+//--------------------------------------------------------------------
+
+//--------------------------------------------------------------------
void vvLandmarks::SaveFile(std::string filename)
{
std::string fileContent = "LANDMARKS1\n"; //File format version identification
fp << fileContent.c_str()<< std::endl;
fp.close();
}
+//--------------------------------------------------------------------
+
+//--------------------------------------------------------------------
void vvLandmarks::SetTime(int time)
{
if (time >= 0 && time <= ((int)mPoints.size() -1)) {
mPolyData->SetPoints(mPoints[time]);
mPolyData->GetPointData()->SetScalars(mIds);
+ mPolyData->GetPointData()->AddArray(mLabels);
mPolyData->Modified();
mPolyData->Update();
}
}
+//--------------------------------------------------------------------
+
+//--------------------------------------------------------------------
std::string vvLandmarks::replace_dots(std::string input)
{
///Replaces the dots used in the file with the decimal separator in use on the platform
}
return input;
}
+//--------------------------------------------------------------------
git://git.creatis.insa-lyon.fr / clitk.git / blobdiff