My Project
3d/deformer.hh
Go to the documentation of this file.
1/* -*- mia-c++ -*-
2 *
3 * This file is part of MIA - a toolbox for medical image analysis
4 * Copyright (c) Leipzig, Madrid 1999-2017 Gert Wollny
5 *
6 * MIA is free software; you can redistribute it and/or modify
7 * it under the terms of the GNU General Public License as published by
8 * the Free Software Foundation; either version 3 of the License, or
9 * (at your option) any later version.
10 *
11 * This program is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14 * GNU General Public License for more details.
15 *
16 * You should have received a copy of the GNU General Public License
17 * along with MIA; if not, see <http://www.gnu.org/licenses/>.
18 *
19 */
20
21#ifndef reg3d_deformer_hh
22#define reg3d_deformer_hh
23
24#include <memory>
25#include <mia/core/parallel.hh>
26
27
28#include <mia/3d/image.hh>
29#include <mia/3d/filter.hh>
32
33
35
42struct FDeformer3D: public TFilter<P3DImage> {
44 m_vf(vf),
45 m_ipfac(ipfac)
46 {
47 }
48 template <typename T>
49 P3DImage operator () (const T3DImage<T>& image) const
50 {
51 T3DImage<T> *timage = new T3DImage<T>(m_vf.get_size(), image);
52 P3DImage result(timage);
53 this->operator()(image, *timage);
54 return result;
55 }
56
57 template <typename T>
58 P3DImage operator () (const T3DImage<T>& image, T3DImage<T>& result) const
59 {
60 assert(result.get_size() == m_vf.get_size());
61 std::unique_ptr<T3DConvoluteInterpolator<T>> interp(m_ipfac.create(image.data()));
62 const auto& rinterp = *interp;
63 auto callback = [this, &rinterp, &result](const C1DParallelRange & range) {
64 CThreadMsgStream thread_stream;
65 auto cache = rinterp.create_cache();
66
67 for (auto z = range.begin(); z != range.end(); ++z) {
68 auto r = result.begin_at(0, 0, z);
69 auto v = m_vf.begin_at(0, 0, z);
70
71 for (size_t y = 0; y < result.get_size().y; ++y)
72 for (size_t x = 0; x < result.get_size().x; ++x, ++r, ++v)
73 *r = rinterp(C3DFVector(x - v->x, y - v->y, z - v->z), cache);
74 }
75 };
76 pfor(C1DParallelRange(0, result.get_size().z, 1), callback);
77 return P3DImage();
78 }
79
80private:
81 C3DFVectorfield m_vf;
83};
84
85
87
88#endif
C3DImage::Pointer P3DImage
define a shortcut to the 3D image shared pointer.
Definition 3d/image.hh:135
T3DVector< float > C3DFVector
A float 3D Vector.
Definition 3d/vector.hh:409
a 3D field of floating point single accuracy 3D vectors
A factory to create interpolators of a given type by providing input data.
T3DConvoluteInterpolator< T > * create(const T3DDatafield< T > &src) const __attribute__((warn_unused_result))
This class is used to handle syncronizized output of logging output in a multi-threaded environment.
const C3DBounds & get_size() const
const_iterator begin_at(size_t x, size_t y, size_t z) const
Specific type of the 3D images that hold real pixel data.
Definition 3d/image.hh:150
const_iterator begin_at(size_t x, size_t y, size_t z) const
constant iterator starting at the given location
Definition 3d/image.hh:372
virtual const C3DBounds & get_size() const
const T3DDatafield< T > & data() const
read only access to the underlying data
T y
vector element
Definition 3d/vector.hh:54
T z
vector element
Definition 3d/vector.hh:56
T x
vector element
Definition 3d/vector.hh:52
#define NS_MIA_BEGIN
conveniance define to start the mia namespace
Definition defines.hh:33
#define NS_MIA_END
conveniance define to end the mia namespace
Definition defines.hh:36
void pfor(Range range, const Func &f)
A filter to transform an image.
P3DImage operator()(const T3DImage< T > &image) const
FDeformer3D(const C3DFVectorfield &vf, const C3DInterpolatorFactory &ipfac)
base class for all filer type functors.