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
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkMetaGaussianConverter.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 __itkMetaGaussianConverter_txx
#define __itkMetaGaussianConverter_txx
#include "itkMetaGaussianConverter.h"
namespace itk
{
/** Constructor */
template <unsigned int NDimensions>
MetaGaussianConverter<NDimensions>
::MetaGaussianConverter()
{
}
/** Convert a metaGaussian into a gaussian SpatialObject */
template <unsigned int NDimensions>
typename MetaGaussianConverter<NDimensions>::SpatialObjectPointer
MetaGaussianConverter<NDimensions>
::MetaGaussianToGaussianSpatialObject(MetaGaussian * gaussian)
{
SpatialObjectPointer spatialObject = SpatialObjectType::New();
spatialObject->SetMaximum( gaussian->Maximum() );
spatialObject->SetRadius( gaussian->Radius() );
spatialObject->GetProperty()->SetName(gaussian->Name());
spatialObject->SetId(gaussian->ID());
spatialObject->SetParentId(gaussian->ParentID());
return spatialObject;
}
/** Convert a gaussian SpatialObject into a metaGaussian */
template <unsigned int NDimensions>
MetaGaussian*
MetaGaussianConverter<NDimensions>
::GaussianSpatialObjectToMetaGaussian(SpatialObjectType * spatialObject)
{
MetaGaussian* gaussian = new MetaGaussian(NDimensions);
if(spatialObject->GetParent())
{
gaussian->ParentID(spatialObject->GetParent()->GetId());
}
gaussian->Maximum( spatialObject->GetMaximum() );
gaussian->Radius( spatialObject->GetRadius() );
gaussian->ID(spatialObject->GetId());
return gaussian;
}
/** Read a meta file give the type */
template <unsigned int NDimensions>
typename MetaGaussianConverter<NDimensions>::SpatialObjectPointer
MetaGaussianConverter<NDimensions>
::ReadMeta(const char* name)
{
SpatialObjectPointer spatialObject;
MetaGaussian* gaussian = new MetaGaussian();
gaussian->Read(name);
spatialObject = MetaGaussianToGaussianSpatialObject(gaussian);
return spatialObject;
}
/** Write a meta gaussian file */
template <unsigned int NDimensions>
bool
MetaGaussianConverter<NDimensions>
::WriteMeta(SpatialObjectType* spatialObject,const char* name)
{
MetaGaussian* gaussian = GaussianSpatialObjectToMetaGaussian(spatialObject);
gaussian->Write(name);
return true;
}
} // end namespace itk
#endif
|
