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 // unsigned int npts = at.GetNumberOfValues() / 3;
148 assert(at.GetNumberOfValues() == static_cast<unsigned int>(mNbOfPoints)*3);
151 mData = vtkSmartPointer<vtkPoints>::New();
152 mData->SetDataTypeToDouble();
153 mData->SetNumberOfPoints(mNbOfPoints);
154 for(unsigned int i=0; i<mNbOfPoints; i++) {
157 p[1] = points[i*3+1];
158 p[2] = points[i*3+2];
159 mData->SetPoint(i, p);
160 if (mZ == -1) mZ = p[2];
165 std::cout << "ERROR ! contour not in the same slice" << std::endl;
174 bool clitk::DicomRT_Contour::Read(gdcm::SQItem * item)
177 // Contour type [Contour Geometric Type]
178 mType = item->GetEntryValue(0x3006,0x0042);
180 if (mType != "CLOSED_PLANAR " && mType != "POINT ") { ///WARNING to the space after the name ...
181 //std::cerr << "Skip this contour : type=" << mType << std::endl;
184 if (mType == "POINT ") {
185 std::cerr << "Warning: POINT type not fully supported. (don't use GetMesh() with this!)"
189 // Number of points [Number of Contour Points]
190 mNbOfPoints = parse_value<int>(item->GetEntryValue(0x3006,0x0046));
193 // Read values [Contour Data]
194 std::vector<float> points = parse_string<float>(item->GetEntryValue(0x3006,0x0050),'\\');
195 assert(points.size() == static_cast<unsigned int>(mNbOfPoints)*3);
198 mData = vtkSmartPointer<vtkPoints>::New();
199 mData->SetDataTypeToDouble();
200 mData->SetNumberOfPoints(mNbOfPoints);
201 for(unsigned int i=0; i<mNbOfPoints; i++) {
204 p[1] = points[i*3+1];
205 p[2] = points[i*3+2];
206 mData->SetPoint(i, p);
207 if (mZ == -1) mZ = p[2];
212 std::cout << "ERROR ! contour not in the same slice" << std::endl;
220 //--------------------------------------------------------------------
223 //--------------------------------------------------------------------
224 vtkPolyData * clitk::DicomRT_Contour::GetMesh()
226 if (!mMeshIsUpToDate) {
227 ComputeMeshFromDataPoints();
231 //--------------------------------------------------------------------
234 //--------------------------------------------------------------------
235 void clitk::DicomRT_Contour::SetMesh(vtkPolyData * mesh)
239 //--------------------------------------------------------------------
242 //--------------------------------------------------------------------
243 void clitk::DicomRT_Contour::ComputeMeshFromDataPoints()
245 // DD("ComputeMesh Contour");
246 mMesh = vtkSmartPointer<vtkPolyData>::New();
247 mMesh->Allocate(); //for cell structures
248 mPoints = vtkSmartPointer<vtkPoints>::New();
249 mMesh->SetPoints(mPoints);
251 for (unsigned int idx=0 ; idx<mNbOfPoints ; idx++) {
252 mMesh->GetPoints()->InsertNextPoint(mData->GetPoint(idx)[0],
253 mData->GetPoint(idx)[1],
254 mData->GetPoint(idx)[2]);
256 ids[1]=(ids[0]+1) % mNbOfPoints; //0-1,1-2,...,n-1-0
257 mMesh->GetLines()->InsertNextCell(2,ids);
259 mMeshIsUpToDate = true;
261 //--------------------------------------------------------------------
264 //--------------------------------------------------------------------
265 void clitk::DicomRT_Contour::ComputeDataPointsFromMesh()
267 DD("ComputeDataPointsFromMesh");
272 mMesh = vtkSmartPointer<vtkPolyData>::New();
273 mMesh->Allocate(); //for cell structures
274 mPoints = vtkSmartPointer<vtkPoints>::New();
275 mMesh->SetPoints(mPoints);
277 for (unsigned int idx=0 ; idx<mNbOfPoints ; idx++) {
278 mMesh->GetPoints()->InsertNextPoint(mData->GetPoint(idx)[0],
279 mData->GetPoint(idx)[1],
280 mData->GetPoint(idx)[2]);
282 ids[1]=(ids[0]+1) % mNbOfPoints; //0-1,1-2,...,n-1-0
283 mMesh->GetLines()->InsertNextCell(2,ids);
285 mMeshIsUpToDate = true;
288 //--------------------------------------------------------------------