1 /*=========================================================================
2 Program: vv http://www.creatis.insa-lyon.fr/rio/vv
3 Main authors : XX XX XX
6 - University of LYON http://www.universite-lyon.fr/
7 - Léon Bérard cancer center http://www.centreleonberard.fr
8 - CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr
10 This software is distributed WITHOUT ANY WARRANTY; without even
11 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
12 PURPOSE. See the copyright notices for more information.
14 It is distributed under dual licence
15 - BSD http://www.opensource.org/licenses/bsd-license.php
16 - CeCILL-B http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
18 =========================================================================*/
20 #include "clitkDicomRT_Contour.h"
21 #include <vtkCellArray.h>
23 #if GDCM_MAJOR_VERSION == 2
24 #include "gdcmAttribute.h"
28 //--------------------------------------------------------------------
29 clitk::DicomRT_Contour::DicomRT_Contour()
31 mMeshIsUpToDate = false;
35 //--------------------------------------------------------------------
38 //--------------------------------------------------------------------
39 clitk::DicomRT_Contour::~DicomRT_Contour()
43 //--------------------------------------------------------------------
47 //--------------------------------------------------------------------
48 void clitk::DicomRT_Contour::Print(std::ostream & os) const
50 DD("TODO : print Contours");
52 //--------------------------------------------------------------------
55 //--------------------------------------------------------------------
56 void clitk::DicomRT_Contour::UpdateDicomItem()
58 DD("DicomRT_Contour::UpdateDicomItem");
60 gdcm::DataSet & nestedds = mItem->GetNestedDataSet();
62 //NON CONSIDER CONTOUR ITEM NOT SEQ ITEM ?
64 // Contour type [Contour Geometric Type]
65 gdcm::Attribute<0x3006,0x0042> contgeotype;
66 contgeotype.SetFromDataSet( nestedds );
68 // Number of points [Number of Contour Points]
69 gdcm::Attribute<0x3006,0x0046> numcontpoints;
70 numcontpoints.SetFromDataSet( nestedds );
72 mNbOfPoints = numcontpoints.GetValue();
75 // Contour dicom values from DataPoints
76 //ComputeDataPointsFromMesh();
77 uint nb = mData->GetNumberOfPoints();
78 std::vector<double> points;
79 points.resize(mNbOfPoints*3);
80 for(unsigned int i=0; i<nb; i++) {
81 double * p = mData->GetPoint(i);
88 gdcm::Attribute<0x3006,0x0050> at;
89 gdcm::Tag tcontourdata(0x3006,0x0050);
90 gdcm::DataElement contourdata = nestedds.GetDataElement( tcontourdata );
91 at.SetFromDataElement( contourdata );
94 at.SetValues(&points[0], points.size(), false);
95 DD(at.GetValues()[0]);
98 nestedds.Replace(at.GetAsDataElement());
101 // Change Number of points [Number of Contour Points]
102 numcontpoints.SetValue(nb);
103 nestedds.Replace(numcontpoints.GetAsDataElement());
106 gdcm::DataElement aa = nestedds.GetDataElement( tcontourdata );
107 at.SetFromDataElement( aa );
108 const double* bb = at.GetValues();
112 //--------------------------------------------------------------------
115 //--------------------------------------------------------------------
116 #if GDCM_MAJOR_VERSION == 2
117 bool clitk::DicomRT_Contour::Read(gdcm::Item * item)
121 const gdcm::DataSet& nestedds2 = item->GetNestedDataSet();
123 // Contour type [Contour Geometric Type]
124 gdcm::Attribute<0x3006,0x0042> contgeotype;
125 contgeotype.SetFromDataSet( nestedds2 );
127 if (contgeotype.GetValue() != "CLOSED_PLANAR " && contgeotype.GetValue() != "POINT ") { ///WARNING to the space after the name ...
128 //std::cerr << "Skip this contour : type=" << mType << std::endl;
131 if (contgeotype.GetValue() == "POINT ") {
132 std::cerr << "Warning: POINT type not fully supported. (don't use GetMesh() with this!)"
136 gdcm::Attribute<0x3006,0x0046> numcontpoints;
137 numcontpoints.SetFromDataSet( nestedds2 );
138 // Number of points [Number of Contour Points]
139 mNbOfPoints = numcontpoints.GetValue();
142 gdcm::Attribute<0x3006,0x0050> at;
143 gdcm::Tag tcontourdata(0x3006,0x0050);
144 const gdcm::DataElement & contourdata = nestedds2.GetDataElement( tcontourdata );
145 at.SetFromDataElement( contourdata );
146 const double* points = at.GetValues();
147 <<<<<<< Updated upstream
150 // unsigned int npts = at.GetNumberOfValues() / 3;
153 // unsigned int npts = at.GetNumberOfValues() / 3;
154 >>>>>>> Stashed changes
156 assert(at.GetNumberOfValues() == static_cast<unsigned int>(mNbOfPoints)*3);
159 mData = vtkSmartPointer<vtkPoints>::New();
160 mData->SetDataTypeToDouble();
161 mData->SetNumberOfPoints(mNbOfPoints);
162 for(unsigned int i=0; i<mNbOfPoints; i++) {
165 p[1] = points[i*3+1];
166 p[2] = points[i*3+2];
167 mData->SetPoint(i, p);
168 if (mZ == -1) mZ = p[2];
173 std::cout << "ERROR ! contour not in the same slice" << std::endl;
182 bool clitk::DicomRT_Contour::Read(gdcm::SQItem * item)
185 // Contour type [Contour Geometric Type]
186 mType = item->GetEntryValue(0x3006,0x0042);
188 if (mType != "CLOSED_PLANAR " && mType != "POINT ") { ///WARNING to the space after the name ...
189 //std::cerr << "Skip this contour : type=" << mType << std::endl;
192 if (mType == "POINT ") {
193 std::cerr << "Warning: POINT type not fully supported. (don't use GetMesh() with this!)"
197 // Number of points [Number of Contour Points]
198 mNbOfPoints = parse_value<int>(item->GetEntryValue(0x3006,0x0046));
201 // Read values [Contour Data]
202 std::vector<float> points = parse_string<float>(item->GetEntryValue(0x3006,0x0050),'\\');
203 assert(points.size() == static_cast<unsigned int>(mNbOfPoints)*3);
206 mData = vtkSmartPointer<vtkPoints>::New();
207 mData->SetDataTypeToDouble();
208 mData->SetNumberOfPoints(mNbOfPoints);
209 for(unsigned int i=0; i<mNbOfPoints; i++) {
212 p[1] = points[i*3+1];
213 p[2] = points[i*3+2];
214 mData->SetPoint(i, p);
215 if (mZ == -1) mZ = p[2];
220 std::cout << "ERROR ! contour not in the same slice" << std::endl;
228 //--------------------------------------------------------------------
231 //--------------------------------------------------------------------
232 vtkPolyData * clitk::DicomRT_Contour::GetMesh()
234 if (!mMeshIsUpToDate) {
235 ComputeMeshFromDataPoints();
239 //--------------------------------------------------------------------
242 //--------------------------------------------------------------------
243 void clitk::DicomRT_Contour::SetMesh(vtkPolyData * mesh)
247 //--------------------------------------------------------------------
250 //--------------------------------------------------------------------
251 void clitk::DicomRT_Contour::ComputeMeshFromDataPoints()
253 // DD("ComputeMesh Contour");
254 mMesh = vtkSmartPointer<vtkPolyData>::New();
255 mMesh->Allocate(); //for cell structures
256 mPoints = vtkSmartPointer<vtkPoints>::New();
257 mMesh->SetPoints(mPoints);
259 for (unsigned int idx=0 ; idx<mNbOfPoints ; idx++) {
260 mMesh->GetPoints()->InsertNextPoint(mData->GetPoint(idx)[0],
261 mData->GetPoint(idx)[1],
262 mData->GetPoint(idx)[2]);
264 ids[1]=(ids[0]+1) % mNbOfPoints; //0-1,1-2,...,n-1-0
265 mMesh->GetLines()->InsertNextCell(2,ids);
267 mMeshIsUpToDate = true;
269 //--------------------------------------------------------------------
272 //--------------------------------------------------------------------
273 void clitk::DicomRT_Contour::ComputeDataPointsFromMesh()
275 DD("ComputeDataPointsFromMesh");
280 mMesh = vtkSmartPointer<vtkPolyData>::New();
281 mMesh->Allocate(); //for cell structures
282 mPoints = vtkSmartPointer<vtkPoints>::New();
283 mMesh->SetPoints(mPoints);
285 for (unsigned int idx=0 ; idx<mNbOfPoints ; idx++) {
286 mMesh->GetPoints()->InsertNextPoint(mData->GetPoint(idx)[0],
287 mData->GetPoint(idx)[1],
288 mData->GetPoint(idx)[2]);
290 ids[1]=(ids[0]+1) % mNbOfPoints; //0-1,1-2,...,n-1-0
291 mMesh->GetLines()->InsertNextCell(2,ids);
293 mMeshIsUpToDate = true;
296 //--------------------------------------------------------------------