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
|
/************************************************************************
*
* Copyright (C) 2017-2024 IRCAD France
* Copyright (C) 2017-2020 IHU Strasbourg
*
* This file is part of Sight.
*
* Sight is free software: you can redistribute it and/or modify it under
* the terms of the GNU Lesser General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* Sight is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* License along with Sight. If not, see <https://www.gnu.org/licenses/>.
*
***********************************************************************/
#include "image_center.hpp"
#include <data/helper/medical_image.hpp>
#include <geometry/data/matrix4.hpp>
#include <service/macros.hpp>
namespace sight::module::filter::image
{
//------------------------------------------------------------------------------
image_center::image_center() :
filter(m_signals)
{
}
//------------------------------------------------------------------------------
void image_center::configuring()
{
}
//------------------------------------------------------------------------------
void image_center::starting()
{
}
//------------------------------------------------------------------------------
void image_center::updating()
{
const auto image = m_image.const_lock();
SIGHT_ASSERT("Missing image '" << IMAGE_IN << "'", image);
const bool image_validity = data::helper::medical_image::check_image_validity(image.get_shared());
if(!image_validity)
{
SIGHT_WARN("Can not compute center of a invalid image.");
return;
}
auto matrix = m_transform.lock();
SIGHT_ASSERT("Missing matrix '" << TRANSFORM_INOUT << "'", matrix);
geometry::data::identity(*matrix);
const auto& origin = image->origin();
const auto& orientation = image->orientation();
// Get World transform
const glm::dmat4 world_transform {
orientation[0], orientation[3], orientation[6], 0,
orientation[1], orientation[4], orientation[7], 0,
orientation[2], orientation[5], orientation[8], 0,
origin[0], origin[1], origin[2], 1
};
// Get image center in mm
const auto& size = image->size();
const auto& spacing = image->spacing();
const glm::dvec4 image_center {
double(size[0]) * spacing[0] / 2.0,
double(size[1]) * spacing[1] / 2.0,
double(size[2]) * spacing[2] / 2.0,
1
};
// Compute world center
const auto world_center = world_transform * image_center;
// Update the output matrix
matrix->set_orientation(orientation);
matrix->set_position(world_center);
matrix->signal<data::matrix4::modified_signal_t>(data::matrix4::MODIFIED_SIG)->async_emit();
this->async_emit(signals::COMPUTED);
}
//------------------------------------------------------------------------------
void image_center::stopping()
{
}
//------------------------------------------------------------------------------
service::connections_t image_center::auto_connections() const
{
return {{IMAGE_IN, data::image::MODIFIED_SIG, service::slots::UPDATE}};
}
//------------------------------------------------------------------------------
} // namespace sight::module::filter::image
|
