diff options
Diffstat (limited to 'src/VolumeGeometry2D.cpp')
| -rw-r--r-- | src/VolumeGeometry2D.cpp | 282 |
1 files changed, 282 insertions, 0 deletions
diff --git a/src/VolumeGeometry2D.cpp b/src/VolumeGeometry2D.cpp new file mode 100644 index 0000000..01d9b9a --- /dev/null +++ b/src/VolumeGeometry2D.cpp @@ -0,0 +1,282 @@ +/* +----------------------------------------------------------------------- +Copyright 2012 iMinds-Vision Lab, University of Antwerp + +Contact: astra@ua.ac.be +Website: http://astra.ua.ac.be + + +This file is part of the +All Scale Tomographic Reconstruction Antwerp Toolbox ("ASTRA Toolbox"). + +The ASTRA Toolbox is free software: you can redistribute it and/or modify +it under the terms of the GNU General Public License as published by +the Free Software Foundation, either version 3 of the License, or +(at your option) any later version. + +The ASTRA Toolbox 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 General Public License for more details. + +You should have received a copy of the GNU General Public License +along with the ASTRA Toolbox. If not, see <http://www.gnu.org/licenses/>. + +----------------------------------------------------------------------- +$Id$ +*/ + +#include "astra/VolumeGeometry2D.h" + +#include <boost/lexical_cast.hpp> +#include <cmath> + +namespace astra +{ + +//---------------------------------------------------------------------------------------- +// Check all variable values +bool CVolumeGeometry2D::_check() +{ + ASTRA_CONFIG_CHECK(m_iGridColCount > 0, "VolumeGeometry2D", "GridColCount must be strictly positive."); + ASTRA_CONFIG_CHECK(m_iGridRowCount > 0, "VolumeGeometry2D", "GridRowCount must be strictly positive."); + ASTRA_CONFIG_CHECK(m_fWindowMinX < m_fWindowMaxX, "VolumeGeometry2D", "WindowMinX should be lower than WindowMaxX."); + ASTRA_CONFIG_CHECK(m_fWindowMinY < m_fWindowMaxY, "VolumeGeometry2D", "WindowMinY should be lower than WindowMaxY."); + + ASTRA_CONFIG_CHECK(m_iGridTotCount == (m_iGridColCount * m_iGridRowCount), "VolumeGeometry2D", "Internal configuration error."); + ASTRA_CONFIG_CHECK(m_fWindowLengthX == (m_fWindowMaxX - m_fWindowMinX), "VolumeGeometry2D", "Internal configuration error."); + ASTRA_CONFIG_CHECK(m_fWindowLengthY == (m_fWindowMaxY - m_fWindowMinY), "VolumeGeometry2D", "Internal configuration error."); + ASTRA_CONFIG_CHECK(m_fWindowArea == (m_fWindowLengthX * m_fWindowLengthY), "VolumeGeometry2D", "Internal configuration error."); + ASTRA_CONFIG_CHECK(m_fPixelLengthX == (m_fWindowLengthX / (float32)m_iGridColCount), "VolumeGeometry2D", "Internal configuration error."); + ASTRA_CONFIG_CHECK(m_fPixelLengthY == (m_fWindowLengthY / (float32)m_iGridRowCount), "VolumeGeometry2D", "Internal configuration error."); + + ASTRA_CONFIG_CHECK(m_fPixelArea == (m_fPixelLengthX * m_fPixelLengthY), "VolumeGeometry2D", "Internal configuration error."); + ASTRA_CONFIG_CHECK(fabsf(m_fDivPixelLengthX * m_fPixelLengthX - 1.0f) < eps, "VolumeGeometry2D", "Internal configuration error."); + ASTRA_CONFIG_CHECK(fabsf(m_fDivPixelLengthY * m_fPixelLengthY - 1.0f) < eps, "VolumeGeometry2D", "Internal configuration error."); + + return true; +} + +//---------------------------------------------------------------------------------------- +// Clear all member variables, setting all numeric variables to 0 and all pointers to NULL. +void CVolumeGeometry2D::clear() +{ + m_iGridColCount = 0; + m_iGridRowCount = 0; + m_iGridTotCount = 0; + + m_fWindowLengthX = 0.0f; + m_fWindowLengthY = 0.0f; + m_fWindowArea = 0.0f; + + m_fPixelLengthX = 0.0f; + m_fPixelLengthY = 0.0f; + m_fPixelArea = 0.0f; + + m_fDivPixelLengthX = 0.0f; + m_fDivPixelLengthY = 0.0f; + + m_fWindowMinX = 0.0f; + m_fWindowMinY = 0.0f; + m_fWindowMaxX = 0.0f; + m_fWindowMaxY = 0.0f; + + m_bInitialized = false; +} + +//---------------------------------------------------------------------------------------- +// Default constructor. +CVolumeGeometry2D::CVolumeGeometry2D() : configCheckData(0) +{ + clear(); +} + +//---------------------------------------------------------------------------------------- +// Default constructor +CVolumeGeometry2D::CVolumeGeometry2D(int _iGridColCount, int _iGridRowCount) + : configCheckData(0) +{ + clear(); + initialize(_iGridColCount, _iGridRowCount); +} + +//---------------------------------------------------------------------------------------- +// Constructor. +CVolumeGeometry2D::CVolumeGeometry2D(int _iGridColCount, + int _iGridRowCount, + float32 _fWindowMinX, + float32 _fWindowMinY, + float32 _fWindowMaxX, + float32 _fWindowMaxY) +{ + clear(); + initialize(_iGridColCount, + _iGridRowCount, + _fWindowMinX, + _fWindowMinY, + _fWindowMaxX, + _fWindowMaxY); +} + +//---------------------------------------------------------------------------------------- +// Destructor. +CVolumeGeometry2D::~CVolumeGeometry2D() +{ + if (m_bInitialized) { + clear(); + } +} + +//---------------------------------------------------------------------------------------- +// Clone +CVolumeGeometry2D* CVolumeGeometry2D::clone() +{ + CVolumeGeometry2D* res = new CVolumeGeometry2D(); + res->m_bInitialized = m_bInitialized; + res->m_iGridColCount = m_iGridColCount; + res->m_iGridRowCount = m_iGridRowCount; + res->m_iGridTotCount = m_iGridTotCount; + res->m_fWindowLengthX = m_fWindowLengthX; + res->m_fWindowLengthY = m_fWindowLengthY; + res->m_fWindowArea = m_fWindowArea; + res->m_fPixelLengthX = m_fPixelLengthX; + res->m_fPixelLengthY = m_fPixelLengthY; + res->m_fPixelArea = m_fPixelArea; + res->m_fDivPixelLengthX = m_fDivPixelLengthX; + res->m_fDivPixelLengthY = m_fDivPixelLengthY; + res->m_fWindowMinX = m_fWindowMinX; + res->m_fWindowMinY = m_fWindowMinY; + res->m_fWindowMaxX = m_fWindowMaxX; + res->m_fWindowMaxY = m_fWindowMaxY; + return res; +} + +//---------------------------------------------------------------------------------------- +// Initialization witha COnfig object +bool CVolumeGeometry2D::initialize(const Config& _cfg) +{ + ASTRA_ASSERT(_cfg.self); + ConfigStackCheck<CVolumeGeometry2D> CC("VolumeGeometry2D", this, _cfg); + + // uninitialize if the object was initialized before + if (m_bInitialized) { + clear(); + } + + // Required: GridColCount + XMLNode* node = _cfg.self->getSingleNode("GridColCount"); + ASTRA_CONFIG_CHECK(node, "ReconstructionGeometry2D", "No GridColCount tag specified."); + m_iGridColCount = boost::lexical_cast<int>(node->getContent()); + ASTRA_DELETE(node); + CC.markNodeParsed("GridColCount"); + + // Required: GridRowCount + node = _cfg.self->getSingleNode("GridRowCount"); + ASTRA_CONFIG_CHECK(node, "ReconstructionGeometry2D", "No GridRowCount tag specified."); + m_iGridRowCount = boost::lexical_cast<int>(node->getContent()); + ASTRA_DELETE(node); + CC.markNodeParsed("GridRowCount"); + + // Optional: Window minima and maxima + m_fWindowMinX = _cfg.self->getOptionNumerical("WindowMinX", -m_iGridColCount/2.0f); + m_fWindowMaxX = _cfg.self->getOptionNumerical("WindowMaxX", m_iGridColCount/2.0f); + m_fWindowMinY = _cfg.self->getOptionNumerical("WindowMinY", -m_iGridRowCount/2.0f); + m_fWindowMaxY = _cfg.self->getOptionNumerical("WindowMaxY", m_iGridRowCount/2.0f); + CC.markOptionParsed("WindowMinX"); + CC.markOptionParsed("WindowMaxX"); + CC.markOptionParsed("WindowMinY"); + CC.markOptionParsed("WindowMaxY"); + + _calculateDependents(); + + // success + m_bInitialized = _check(); + return m_bInitialized; +} + +//---------------------------------------------------------------------------------------- +// Initialization. +bool CVolumeGeometry2D::initialize(int _iGridColCount, int _iGridRowCount) +{ + return initialize(_iGridColCount, + _iGridRowCount, + -_iGridColCount/2.0f, + -_iGridRowCount/2.0f, + _iGridColCount/2.0f, + _iGridRowCount/2.0f); +} + +//---------------------------------------------------------------------------------------- +// Initialization. +bool CVolumeGeometry2D::initialize(int _iGridColCount, + int _iGridRowCount, + float32 _fWindowMinX, + float32 _fWindowMinY, + float32 _fWindowMaxX, + float32 _fWindowMaxY) +{ + if (m_bInitialized) { + clear(); + } + + m_iGridColCount = _iGridColCount; + m_iGridRowCount = _iGridRowCount; + + m_fWindowMinX = _fWindowMinX; + m_fWindowMinY = _fWindowMinY; + m_fWindowMaxX = _fWindowMaxX; + m_fWindowMaxY = _fWindowMaxY; + + _calculateDependents(); + + m_bInitialized = _check(); + return m_bInitialized; +} + +void CVolumeGeometry2D::_calculateDependents() +{ + m_iGridTotCount = (m_iGridColCount * m_iGridRowCount); + + m_fWindowLengthX = (m_fWindowMaxX - m_fWindowMinX); + m_fWindowLengthY = (m_fWindowMaxY - m_fWindowMinY); + m_fWindowArea = (m_fWindowLengthX * m_fWindowLengthY); + + m_fPixelLengthX = (m_fWindowLengthX / (float32)m_iGridColCount); + m_fPixelLengthY = (m_fWindowLengthY / (float32)m_iGridRowCount); + m_fPixelArea = (m_fPixelLengthX * m_fPixelLengthY); + + m_fDivPixelLengthX = ((float32)m_iGridColCount / m_fWindowLengthX); // == (1.0f / m_fPixelLengthX); + m_fDivPixelLengthY = ((float32)m_iGridRowCount / m_fWindowLengthY); // == (1.0f / m_fPixelLengthY); +} + +//---------------------------------------------------------------------------------------- +// is of type +bool CVolumeGeometry2D::isEqual(CVolumeGeometry2D* _pGeom2) const +{ + if (_pGeom2 == NULL) return false; + + // both objects must be initialized + if (!m_bInitialized || !_pGeom2->m_bInitialized) return false; + + // check all values + if (m_iGridColCount != _pGeom2->m_iGridColCount) return false; + if (m_iGridRowCount != _pGeom2->m_iGridRowCount) return false; + if (m_iGridTotCount != _pGeom2->m_iGridTotCount) return false; + if (m_fWindowLengthX != _pGeom2->m_fWindowLengthX) return false; + if (m_fWindowLengthY != _pGeom2->m_fWindowLengthY) return false; + if (m_fWindowArea != _pGeom2->m_fWindowArea) return false; + if (m_fPixelLengthX != _pGeom2->m_fPixelLengthX) return false; + if (m_fPixelLengthY != _pGeom2->m_fPixelLengthY) return false; + if (m_fPixelArea != _pGeom2->m_fPixelArea) return false; + if (m_fDivPixelLengthX != _pGeom2->m_fDivPixelLengthX) return false; + if (m_fDivPixelLengthY != _pGeom2->m_fDivPixelLengthY) return false; + if (m_fWindowMinX != _pGeom2->m_fWindowMinX) return false; + if (m_fWindowMinY != _pGeom2->m_fWindowMinY) return false; + if (m_fWindowMaxX != _pGeom2->m_fWindowMaxX) return false; + if (m_fWindowMaxY != _pGeom2->m_fWindowMaxY) return false; + + return true; +} +//---------------------------------------------------------------------------------------- + +} // namespace astra |