Mercurial > hg > orthanc-stone
view OrthancStone/Sources/Toolbox/CoordinateSystem3D.h @ 2006:862c572021de
back to mainline
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Wed, 02 Nov 2022 17:10:04 +0100 |
parents | cbf54cd28d59 |
children | 07964689cb0b |
line wrap: on
line source
/** * Stone of Orthanc * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics * Department, University Hospital of Liege, Belgium * Copyright (C) 2017-2022 Osimis S.A., Belgium * Copyright (C) 2021-2022 Sebastien Jodogne, ICTEAM UCLouvain, Belgium * * This program 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. * * This program 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 this program. If not, see * <http://www.gnu.org/licenses/>. **/ #pragma once #include "../Scene2D/ScenePoint2D.h" #include "LinearAlgebra.h" #include "OrthancDatasets/IDicomDataset.h" #include <iosfwd> namespace OrthancStone { // Geometry of a 3D plane class CoordinateSystem3D { private: bool valid_; Vector origin_; Vector normal_; Vector axisX_; Vector axisY_; double d_; void CheckAndComputeNormal(); void Setup(const std::string& imagePositionPatient, const std::string& imageOrientationPatient); void SetupCanonical(); double GetOffset() const; public: CoordinateSystem3D() { SetupCanonical(); } friend std::ostream& operator<< (std::ostream& s, const CoordinateSystem3D& that); CoordinateSystem3D(const Vector& origin, const Vector& axisX, const Vector& axisY); CoordinateSystem3D(const std::string& imagePositionPatient, const std::string& imageOrientationPatient) { Setup(imagePositionPatient, imageOrientationPatient); } explicit CoordinateSystem3D(const IDicomDataset& dicom); explicit CoordinateSystem3D(const Orthanc::DicomMap& dicom); bool IsValid() const { return valid_; } const Vector& GetNormal() const { return normal_; } const Vector& GetOrigin() const // This is the "Image Position Patient" tag { return origin_; } const Vector& GetAxisX() const { return axisX_; } const Vector& GetAxisY() const { return axisY_; } void SetOrigin(const Vector& origin); Vector MapSliceToWorldCoordinates(double x, double y) const; Vector MapSliceToWorldCoordinates(const ScenePoint2D& p) const { return MapSliceToWorldCoordinates(p.GetX(), p.GetY()); } double ProjectAlongNormal(const Vector& point) const; void ProjectPoint(double& offsetX, double& offsetY, const Vector& point) const; ScenePoint2D ProjectPoint(const Vector& point) const { double x, y; ProjectPoint(x, y, point); return ScenePoint2D(x, y); } /* Alternated faster implementation (untested yet) */ void ProjectPoint2(double& offsetX, double& offsetY, const Vector& point) const; bool IntersectSegment(Vector& p, const Vector& edgeFrom, const Vector& edgeTo) const; bool IntersectLine(Vector& p, const Vector& origin, const Vector& direction) const; // Point-to-plane distance double ComputeDistance(const Vector& p) const; // Returns "false" is the two planes are not parallel static bool ComputeDistance(double& distance, const CoordinateSystem3D& a, const CoordinateSystem3D& b); // Normalize a cutting plane so that the origin (0,0,0) of the 3D // world is mapped to the origin of its (x,y) coordinate system static CoordinateSystem3D NormalizeCuttingPlane(const CoordinateSystem3D& plane); // Construct one possible coordinate system from the general form // of the equation of a plane "a*x+b*y+c*z+d=0". Note that the // axes are not determined in this case, and so they are chosen // arbitrarily. static CoordinateSystem3D CreateFromPlaneGeneralForm(double a, double b, double c, double d); static CoordinateSystem3D CreateFromThreePoints(const Vector& a, const Vector& b, const Vector& c); void GetOrientationMarkers(std::string& top /* out */, std::string& bottom /* out */, std::string& left /* out */, std::string& right /* out */) const; }; }