1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkMetaEllipseConverter.txx
Language: C++
Date: $Date$
Version: $Revision$
Copyright (c) Insight Software Consortium. All rights reserved.
See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notices for more information.
=========================================================================*/
#ifndef __itkMetaEllipseConverter_txx
#define __itkMetaEllipseConverter_txx
#include "itkMetaEllipseConverter.h"
namespace itk
{
/** Constructor */
template <unsigned int NDimensions>
MetaEllipseConverter<NDimensions>
::MetaEllipseConverter()
{
}
/** Convert a metaEllipse into an ellipse SpatialObject */
template <unsigned int NDimensions>
typename MetaEllipseConverter<NDimensions>::SpatialObjectPointer
MetaEllipseConverter<NDimensions>
::MetaEllipseToEllipseSpatialObject(MetaEllipse * ellipse)
{
SpatialObjectPointer spatialObject = SpatialObjectType::New();
double spacing[NDimensions];
typename SpatialObjectType::ArrayType radius;
for(unsigned int i=0;i<NDimensions;i++)
{
radius[i]=ellipse->Radius()[i];
spacing[i]=ellipse->ElementSpacing()[i];
}
spatialObject->GetIndexToObjectTransform()->SetScaleComponent(spacing);
spatialObject->SetRadius(radius);
spatialObject->GetProperty()->SetName(ellipse->Name());
spatialObject->SetId(ellipse->ID());
spatialObject->SetParentId(ellipse->ParentID());
spatialObject->GetProperty()->SetRed(ellipse->Color()[0]);
spatialObject->GetProperty()->SetGreen(ellipse->Color()[1]);
spatialObject->GetProperty()->SetBlue(ellipse->Color()[2]);
spatialObject->GetProperty()->SetAlpha(ellipse->Color()[3]);
return spatialObject;
}
/** Convert an ellipse SpatialObject into a metaEllipse */
template <unsigned int NDimensions>
MetaEllipse*
MetaEllipseConverter<NDimensions>
::EllipseSpatialObjectToMetaEllipse(SpatialObjectType * spatialObject)
{
MetaEllipse* ellipse = new MetaEllipse(NDimensions);
float* radius = new float[NDimensions];
for(unsigned int i=0;i<NDimensions;i++)
{
radius[i] = spatialObject->GetRadius()[i];
}
if(spatialObject->GetParent())
{
ellipse->ParentID(spatialObject->GetParent()->GetId());
}
ellipse->Radius(radius);
ellipse->ID(spatialObject->GetId());
ellipse->Color(spatialObject->GetProperty()->GetRed(),
spatialObject->GetProperty()->GetGreen(),
spatialObject->GetProperty()->GetBlue(),
spatialObject->GetProperty()->GetAlpha());
for(unsigned int i=0;i<NDimensions;i++)
{
ellipse->ElementSpacing(i, spatialObject->GetIndexToObjectTransform()
->GetScaleComponent()[i]);
}
delete []radius;
return ellipse;
}
/** Read a meta file give the type */
template <unsigned int NDimensions>
typename MetaEllipseConverter<NDimensions>::SpatialObjectPointer
MetaEllipseConverter<NDimensions>
::ReadMeta(const char* name)
{
SpatialObjectPointer spatialObject;
MetaEllipse* ellipse = new MetaEllipse();
ellipse->Read(name);
spatialObject = MetaEllipseToEllipseSpatialObject(ellipse);
return spatialObject;
}
/** Write a meta ellipse file */
template <unsigned int NDimensions>
bool
MetaEllipseConverter<NDimensions>
::WriteMeta(SpatialObjectType* spatialObject,const char* name)
{
MetaEllipse* ellipse = EllipseSpatialObjectToMetaEllipse(spatialObject);
ellipse->Write(name);
return true;
}
} // end namespace itk
#endif
|
