orthanc-stone: UnitTestsSources/GeometryToolboxTests.cpp annotate

annotate UnitTestsSources/GeometryToolboxTests.cpp @ 1864:35ac2be493e2

todo

author	Sebastien Jodogne <s.jodogne@gmail.com>
date	Wed, 24 Nov 2021 12:25:19 +0100
parents	ae18b5cfbd3e
children	3889ae96d2e9

rev	line source
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1 /**
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	2 * Stone of Orthanc
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	4 * Department, University Hospital of Liege, Belgium
1739 9ac2a65d4172 upgrade to year 2021 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1679 diff changeset	5 * Copyright (C) 2017-2021 Osimis S.A., Belgium
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	6 *
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	7 * This program is free software: you can redistribute it and/or
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	8 * modify it under the terms of the GNU Affero General Public License
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	9 * as published by the Free Software Foundation, either version 3 of
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	10 * the License, or (at your option) any later version.
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	11 *
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	12 * This program is distributed in the hope that it will be useful, but
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	13 * WITHOUT ANY WARRANTY; without even the implied warranty of
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	15 * Affero General Public License for more details.
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	16 *
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	17 * You should have received a copy of the GNU Affero General Public License
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	18 * along with this program. If not, see <http://www.gnu.org/licenses/>.
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	19 **/
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	20
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	21
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	22 #include <gtest/gtest.h>
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	23
1765 177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	24 #include "../OrthancStone/Sources/Toolbox/DicomInstanceParameters.h"
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	25 #include "../OrthancStone/Sources/Toolbox/FiniteProjectiveCamera.h"
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	26 #include "../OrthancStone/Sources/Toolbox/GenericToolbox.h"
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	27 #include "../OrthancStone/Sources/Toolbox/GeometryToolbox.h"
1765 177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	28 #include "../OrthancStone/Sources/Toolbox/SlicesSorter.h"
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	29
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	30 #include <Logging.h>
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	31 #include <OrthancException.h>
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	32
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	33 #include <boost/math/special_functions/round.hpp>
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	34
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	35
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	36 TEST(GeometryToolbox, Interpolation)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	37 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	38 using namespace OrthancStone::GeometryToolbox;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	39
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	40 // https://en.wikipedia.org/wiki/Bilinear_interpolation#Application_in_image_processing
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	41 ASSERT_DOUBLE_EQ(146.1f, ComputeBilinearInterpolationUnitSquare(0.5f, 0.2f, 91, 210, 162, 95));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	42
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	43 ASSERT_DOUBLE_EQ(91, ComputeBilinearInterpolationUnitSquare(0, 0, 91, 210, 162, 95));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	44 ASSERT_DOUBLE_EQ(210, ComputeBilinearInterpolationUnitSquare(1, 0, 91, 210, 162, 95));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	45 ASSERT_DOUBLE_EQ(162, ComputeBilinearInterpolationUnitSquare(0, 1, 91, 210, 162, 95));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	46 ASSERT_DOUBLE_EQ(95, ComputeBilinearInterpolationUnitSquare(1, 1, 91, 210, 162, 95));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	47
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	48 ASSERT_DOUBLE_EQ(123.35f, ComputeTrilinearInterpolationUnitSquare
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	49 (0.5f, 0.2f, 0.7f,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	50 91, 210, 162, 95,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	51 51, 190, 80, 92));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	52
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	53 ASSERT_DOUBLE_EQ(ComputeBilinearInterpolationUnitSquare(0.5f, 0.2f, 91, 210, 162, 95),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	54 ComputeTrilinearInterpolationUnitSquare(0.5f, 0.2f, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	55 91, 210, 162, 95,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	56 51, 190, 80, 92));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	57
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	58 ASSERT_DOUBLE_EQ(ComputeBilinearInterpolationUnitSquare(0.5f, 0.2f, 51, 190, 80, 92),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	59 ComputeTrilinearInterpolationUnitSquare(0.5f, 0.2f, 1,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	60 91, 210, 162, 95,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	61 51, 190, 80, 92));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	62 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	63
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	64
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	65 static bool CompareMatrix(const OrthancStone::Matrix& a,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	66 const OrthancStone::Matrix& b,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	67 double threshold = 0.00000001)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	68 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	69 if (a.size1() != b.size1() \|\|
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	70 a.size2() != b.size2())
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	71 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	72 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	73 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	74
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	75 for (size_t i = 0; i < a.size1(); i++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	76 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	77 for (size_t j = 0; j < a.size2(); j++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	78 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	79 if (fabs(a(i, j) - b(i, j)) > threshold)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	80 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	81 LOG(ERROR) << "Too large difference in component ("
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	82 << i << "," << j << "): " << a(i,j) << " != " << b(i,j);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	83 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	84 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	85 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	86 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	87
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	88 return true;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	89 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	90
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	91
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	92 static bool CompareVector(const OrthancStone::Vector& a,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	93 const OrthancStone::Vector& b,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	94 double threshold = 0.00000001)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	95 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	96 if (a.size() != b.size())
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	97 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	98 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	99 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	100
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	101 for (size_t i = 0; i < a.size(); i++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	102 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	103 if (fabs(a(i) - b(i)) > threshold)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	104 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	105 LOG(ERROR) << "Too large difference in component "
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	106 << i << ": " << a(i) << " != " << b(i);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	107 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	108 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	109 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	110
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	111 return true;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	112 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	113
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	114
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	115
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	116 TEST(FiniteProjectiveCamera, Decomposition1)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	117 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	118 // Example 6.2 of "Multiple View Geometry in Computer Vision - 2nd
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	119 // edition" (page 163)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	120 const double p[12] = {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	121 3.53553e+2, 3.39645e+2, 2.77744e+2, -1.44946e+6,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	122 -1.03528e+2, 2.33212e+1, 4.59607e+2, -6.32525e+5,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	123 7.07107e-1, -3.53553e-1, 6.12372e-1, -9.18559e+2
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	124 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	125
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	126 OrthancStone::FiniteProjectiveCamera camera(p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	127 ASSERT_EQ(3u, camera.GetMatrix().size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	128 ASSERT_EQ(4u, camera.GetMatrix().size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	129 ASSERT_EQ(3u, camera.GetIntrinsicParameters().size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	130 ASSERT_EQ(3u, camera.GetIntrinsicParameters().size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	131 ASSERT_EQ(3u, camera.GetRotation().size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	132 ASSERT_EQ(3u, camera.GetRotation().size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	133 ASSERT_EQ(3u, camera.GetCenter().size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	134
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	135 ASSERT_NEAR(1000.0, camera.GetCenter()[0], 0.01);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	136 ASSERT_NEAR(2000.0, camera.GetCenter()[1], 0.01);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	137 ASSERT_NEAR(1500.0, camera.GetCenter()[2], 0.01);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	138
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	139 ASSERT_NEAR(468.2, camera.GetIntrinsicParameters() (0, 0), 0.1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	140 ASSERT_NEAR(91.2, camera.GetIntrinsicParameters() (0, 1), 0.1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	141 ASSERT_NEAR(300.0, camera.GetIntrinsicParameters() (0, 2), 0.1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	142 ASSERT_NEAR(427.2, camera.GetIntrinsicParameters() (1, 1), 0.1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	143 ASSERT_NEAR(200.0, camera.GetIntrinsicParameters() (1, 2), 0.1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	144 ASSERT_NEAR(1.0, camera.GetIntrinsicParameters() (2, 2), 0.1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	145
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	146 ASSERT_NEAR(0, camera.GetIntrinsicParameters() (1, 0), 0.0000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	147 ASSERT_NEAR(0, camera.GetIntrinsicParameters() (2, 0), 0.0000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	148 ASSERT_NEAR(0, camera.GetIntrinsicParameters() (2, 1), 0.0000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	149
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	150 ASSERT_NEAR(0.41380, camera.GetRotation() (0, 0), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	151 ASSERT_NEAR(0.90915, camera.GetRotation() (0, 1), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	152 ASSERT_NEAR(0.04708, camera.GetRotation() (0, 2), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	153 ASSERT_NEAR(-0.57338, camera.GetRotation() (1, 0), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	154 ASSERT_NEAR(0.22011, camera.GetRotation() (1, 1), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	155 ASSERT_NEAR(0.78917, camera.GetRotation() (1, 2), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	156 ASSERT_NEAR(0.70711, camera.GetRotation() (2, 0), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	157 ASSERT_NEAR(-0.35355, camera.GetRotation() (2, 1), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	158 ASSERT_NEAR(0.61237, camera.GetRotation() (2, 2), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	159
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	160 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	161
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	162 OrthancStone::FiniteProjectiveCamera camera2(camera.GetIntrinsicParameters(),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	163 camera.GetRotation(),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	164 camera.GetCenter());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	165
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	166 ASSERT_TRUE(CompareMatrix(camera.GetMatrix(), camera2.GetMatrix()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	167 ASSERT_TRUE(CompareMatrix(camera.GetIntrinsicParameters(), camera2.GetIntrinsicParameters()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	168 ASSERT_TRUE(CompareMatrix(camera.GetRotation(), camera2.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	169 ASSERT_TRUE(CompareVector(camera.GetCenter(), camera2.GetCenter()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	170 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	171
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	172
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	173 TEST(FiniteProjectiveCamera, Decomposition2)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	174 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	175 const double p[] = { 1188.111986, 580.205341, -808.445330, 128000.000000, -366.466264, 1446.510501, 418.499736, 128000.000000, -0.487118, 0.291726, -0.823172, 500.000000 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	176 const double k[] = { -1528.494743, 0.000000, 256.000000, 0.000000, 1528.494743, 256.000000, 0.000000, 0.000000, 1.000000 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	177 const double r[] = { -0.858893, -0.330733, 0.391047, -0.158171, 0.897503, 0.411668, -0.487118, 0.291726, -0.823172 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	178 const double c[] = { 243.558936, -145.863085, 411.585964 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	179
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	180 OrthancStone::FiniteProjectiveCamera camera(p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	181 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	182
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	183 OrthancStone::FiniteProjectiveCamera camera2(k, r, c);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	184 ASSERT_TRUE(CompareMatrix(camera.GetMatrix(), camera2.GetMatrix(), 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	185 ASSERT_TRUE(CompareMatrix(camera.GetIntrinsicParameters(), camera2.GetIntrinsicParameters(), 0.001));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	186 ASSERT_TRUE(CompareMatrix(camera.GetRotation(), camera2.GetRotation(), 0.000001));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	187 ASSERT_TRUE(CompareVector(camera.GetCenter(), camera2.GetCenter(), 0.0001));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	188 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	189
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	190
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	191 TEST(FiniteProjectiveCamera, Decomposition3)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	192 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	193 const double p[] = { 10, 0, 0, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	194 0, 20, 0, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	195 0, 0, 30, 0 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	196
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	197 OrthancStone::FiniteProjectiveCamera camera(p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	198 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	199 ASSERT_DOUBLE_EQ(10, camera.GetIntrinsicParameters() (0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	200 ASSERT_DOUBLE_EQ(20, camera.GetIntrinsicParameters() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	201 ASSERT_DOUBLE_EQ(30, camera.GetIntrinsicParameters() (2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	202 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	203 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	204 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	205 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	206 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	207 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	208 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	209
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	210
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	211 TEST(FiniteProjectiveCamera, Decomposition4)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	212 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	213 const double p[] = { 1, 0, 0, 10,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	214 0, 1, 0, 20,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	215 0, 0, 1, 30 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	216
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	217 OrthancStone::FiniteProjectiveCamera camera(p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	218 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	219 ASSERT_DOUBLE_EQ(1, camera.GetIntrinsicParameters() (0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	220 ASSERT_DOUBLE_EQ(1, camera.GetIntrinsicParameters() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	221 ASSERT_DOUBLE_EQ(1, camera.GetIntrinsicParameters() (2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	222 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	223 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	224 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	225 ASSERT_DOUBLE_EQ(-10, camera.GetCenter() (0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	226 ASSERT_DOUBLE_EQ(-20, camera.GetCenter() (1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	227 ASSERT_DOUBLE_EQ(-30, camera.GetCenter() (2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	228 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	229
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	230
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	231 TEST(FiniteProjectiveCamera, Decomposition5)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	232 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	233 const double p[] = { 0, 0, 10, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	234 0, 20, 0, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	235 30, 0, 0, 0 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	236
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	237 OrthancStone::FiniteProjectiveCamera camera(p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	238 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	239 ASSERT_DOUBLE_EQ(-10, camera.GetIntrinsicParameters() (0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	240 ASSERT_DOUBLE_EQ(20, camera.GetIntrinsicParameters() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	241 ASSERT_DOUBLE_EQ(30, camera.GetIntrinsicParameters() (2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	242 ASSERT_DOUBLE_EQ(-1, camera.GetRotation() (0, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	243 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	244 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (2, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	245 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	246 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	247 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	248
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	249 OrthancStone::FiniteProjectiveCamera camera2(camera.GetIntrinsicParameters(),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	250 camera.GetRotation(),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	251 camera.GetCenter());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	252 ASSERT_TRUE(CompareMatrix(camera.GetMatrix(), camera2.GetMatrix()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	253 ASSERT_TRUE(CompareMatrix(camera.GetIntrinsicParameters(), camera2.GetIntrinsicParameters()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	254 ASSERT_TRUE(CompareMatrix(camera.GetRotation(), camera2.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	255 ASSERT_TRUE(CompareVector(camera.GetCenter(), camera2.GetCenter()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	256 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	257
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	258
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	259 static double GetCosAngle(const OrthancStone::Vector& a,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	260 const OrthancStone::Vector& b)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	261 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	262 // Returns the cosine of the angle between two vectors
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	263 // https://en.wikipedia.org/wiki/Dot_product#Geometric_definition
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	264 return boost::numeric::ublas::inner_prod(a, b) /
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	265 (boost::numeric::ublas::norm_2(a) * boost::numeric::ublas::norm_2(b));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	266 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	267
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	268
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	269 TEST(FiniteProjectiveCamera, Ray)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	270 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	271 const double pp[] = { -1499.650894, 2954.618773, -259.737419, 637891.819097,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	272 -2951.517707, -1501.019129, -285.785281, 637891.819097,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	273 0.008528, 0.003067, -0.999959, 2491.764918 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	274
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	275 const OrthancStone::FiniteProjectiveCamera camera(pp);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	276
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	277 ASSERT_NEAR(-21.2492, camera.GetCenter() (0), 0.0001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	278 ASSERT_NEAR(-7.64234, camera.GetCenter() (1), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	279 ASSERT_NEAR(2491.66, camera.GetCenter() (2), 0.01);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	280
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	281 // Image plane that led to these parameters, with principal point at
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	282 // (256,256). The image has dimensions 512x512.
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	283 OrthancStone::Vector o =
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	284 OrthancStone::LinearAlgebra::CreateVector(7.009620, 2.521030, -821.942000);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	285 OrthancStone::Vector ax =
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	286 OrthancStone::LinearAlgebra::CreateVector(-0.453219, 0.891399, -0.001131);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	287 OrthancStone::Vector ay =
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	288 OrthancStone::LinearAlgebra::CreateVector(-0.891359, -0.453210, -0.008992);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	289
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	290 OrthancStone::CoordinateSystem3D imagePlane(o, ax, ay);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	291
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	292 // Back-projection of the principal point
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	293 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	294 OrthancStone::Vector ray = camera.GetRayDirection(256, 256);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	295
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	296 // The principal axis vector is orthogonal to the image plane
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	297 // (i.e. parallel to the plane normal), in the opposite direction
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	298 // ("-1" corresponds to "cos(pi)").
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	299 ASSERT_NEAR(-1, GetCosAngle(ray, imagePlane.GetNormal()), 0.0000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	300
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	301 // Forward projection of principal axis, resulting in the principal point
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	302 double x, y;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	303 camera.ApplyFinite(x, y, camera.GetCenter() - ray);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	304
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	305 ASSERT_NEAR(256, x, 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	306 ASSERT_NEAR(256, y, 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	307 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	308
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	309 // Back-projection of the 4 corners of the image
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	310 std::vector<double> cx, cy;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	311 cx.push_back(0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	312 cy.push_back(0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	313 cx.push_back(512);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	314 cy.push_back(0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	315 cx.push_back(512);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	316 cy.push_back(512);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	317 cx.push_back(0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	318 cy.push_back(512);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	319
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	320 bool first = true;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	321 double angle;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	322
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	323 for (size_t i = 0; i < cx.size(); i++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	324 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	325 OrthancStone::Vector ray = camera.GetRayDirection(cx[i], cy[i]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	326
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	327 // Check that the angle wrt. principal axis is the same for all
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	328 // the 4 corners
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	329 double a = GetCosAngle(ray, imagePlane.GetNormal());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	330 if (first)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	331 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	332 first = false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	333 angle = a;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	334 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	335 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	336 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	337 ASSERT_NEAR(angle, a, 0.000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	338 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	339
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	340 // Forward projection of the ray, going back to the original point
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	341 double x, y;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	342 camera.ApplyFinite(x, y, camera.GetCenter() - ray);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	343
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	344 ASSERT_NEAR(cx[i], x, 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	345 ASSERT_NEAR(cy[i], y, 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	346
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	347 // Alternative construction, by computing the intersection of the
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	348 // ray with the image plane
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	349 OrthancStone::Vector p;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	350 ASSERT_TRUE(imagePlane.IntersectLine(p, camera.GetCenter(), -ray));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	351 imagePlane.ProjectPoint(x, y, p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	352 ASSERT_NEAR(cx[i], x + 256, 0.01);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	353 ASSERT_NEAR(cy[i], y + 256, 0.01);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	354 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	355 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	356
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	357
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	358 TEST(Matrix, Inverse1)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	359 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	360 OrthancStone::Matrix a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	361
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	362 a.resize(0, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	363 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	364 ASSERT_EQ(0u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	365 ASSERT_EQ(0u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	366
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	367 a.resize(2, 3);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	368 ASSERT_THROW(OrthancStone::LinearAlgebra::InvertMatrix(b, a), Orthanc::OrthancException);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	369
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	370 a.resize(1, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	371 a(0, 0) = 45.0;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	372
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	373 ASSERT_DOUBLE_EQ(45, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	374 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	375 ASSERT_EQ(1u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	376 ASSERT_EQ(1u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	377 ASSERT_DOUBLE_EQ(1.0 / 45.0, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	378
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	379 a(0, 0) = 0;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	380 ASSERT_DOUBLE_EQ(0, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	381 ASSERT_THROW(OrthancStone::LinearAlgebra::InvertMatrix(b, a), Orthanc::OrthancException);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	382 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	383
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	384
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	385 TEST(Matrix, Inverse2)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	386 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	387 OrthancStone::Matrix a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	388 a.resize(2, 2);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	389 a(0, 0) = 4;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	390 a(0, 1) = 3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	391 a(1, 0) = 3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	392 a(1, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	393
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	394 ASSERT_DOUBLE_EQ(-1, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	395 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	396 ASSERT_EQ(2u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	397 ASSERT_EQ(2u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	398
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	399 ASSERT_DOUBLE_EQ(-2, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	400 ASSERT_DOUBLE_EQ(3, b(0, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	401 ASSERT_DOUBLE_EQ(3, b(1, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	402 ASSERT_DOUBLE_EQ(-4, b(1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	403
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	404 a(0, 0) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	405 a(0, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	406 a(1, 0) = 3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	407 a(1, 1) = 4;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	408
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	409 ASSERT_DOUBLE_EQ(-2, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	410 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	411
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	412 ASSERT_DOUBLE_EQ(-2, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	413 ASSERT_DOUBLE_EQ(1, b(0, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	414 ASSERT_DOUBLE_EQ(1.5, b(1, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	415 ASSERT_DOUBLE_EQ(-0.5, b(1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	416 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	417
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	418
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	419 TEST(Matrix, Inverse3)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	420 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	421 OrthancStone::Matrix a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	422 a.resize(3, 3);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	423 a(0, 0) = 7;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	424 a(0, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	425 a(0, 2) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	426 a(1, 0) = 0;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	427 a(1, 1) = 3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	428 a(1, 2) = -1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	429 a(2, 0) = -3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	430 a(2, 1) = 4;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	431 a(2, 2) = -2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	432
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	433 ASSERT_DOUBLE_EQ(1, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	434 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	435 ASSERT_EQ(3u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	436 ASSERT_EQ(3u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	437
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	438 ASSERT_DOUBLE_EQ(-2, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	439 ASSERT_DOUBLE_EQ(8, b(0, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	440 ASSERT_DOUBLE_EQ(-5, b(0, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	441 ASSERT_DOUBLE_EQ(3, b(1, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	442 ASSERT_DOUBLE_EQ(-11, b(1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	443 ASSERT_DOUBLE_EQ(7, b(1, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	444 ASSERT_DOUBLE_EQ(9, b(2, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	445 ASSERT_DOUBLE_EQ(-34, b(2, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	446 ASSERT_DOUBLE_EQ(21, b(2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	447
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	448
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	449 a(0, 0) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	450 a(0, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	451 a(0, 2) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	452 a(1, 0) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	453 a(1, 1) = 0;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	454 a(1, 2) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	455 a(2, 0) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	456 a(2, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	457 a(2, 2) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	458
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	459 ASSERT_DOUBLE_EQ(2, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	460 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	461 ASSERT_EQ(3u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	462 ASSERT_EQ(3u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	463
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	464 ASSERT_DOUBLE_EQ(-1, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	465 ASSERT_DOUBLE_EQ(1, b(0, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	466 ASSERT_DOUBLE_EQ(1, b(0, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	467 ASSERT_DOUBLE_EQ(0, b(1, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	468 ASSERT_DOUBLE_EQ(-0.5, b(1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	469 ASSERT_DOUBLE_EQ(0.5, b(1, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	470 ASSERT_DOUBLE_EQ(1, b(2, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	471 ASSERT_DOUBLE_EQ(0, b(2, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	472 ASSERT_DOUBLE_EQ(-1, b(2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	473 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	474
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	475
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	476 TEST(Matrix, Inverse4)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	477 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	478 OrthancStone::Matrix a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	479 a.resize(4, 4);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	480 a(0, 0) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	481 a(0, 1) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	482 a(0, 2) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	483 a(0, 3) = -3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	484 a(1, 0) = -2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	485 a(1, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	486 a(1, 2) = -1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	487 a(1, 3) = -1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	488 a(2, 0) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	489 a(2, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	490 a(2, 2) = -3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	491 a(2, 3) = -1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	492 a(3, 0) = 3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	493 a(3, 1) = -2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	494 a(3, 2) = -3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	495 a(3, 3) = -1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	496
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	497 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	498 ASSERT_EQ(4u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	499 ASSERT_EQ(4u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	500
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	501 b *= 134.0; // This is the determinant
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	502
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	503 ASSERT_DOUBLE_EQ(8, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	504 ASSERT_DOUBLE_EQ(-44, b(0, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	505 ASSERT_DOUBLE_EQ(30, b(0, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	506 ASSERT_DOUBLE_EQ(-10, b(0, 3));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	507 ASSERT_DOUBLE_EQ(2, b(1, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	508 ASSERT_DOUBLE_EQ(-11, b(1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	509 ASSERT_DOUBLE_EQ(41, b(1, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	510 ASSERT_DOUBLE_EQ(-36, b(1, 3));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	511 ASSERT_DOUBLE_EQ(16, b(2, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	512 ASSERT_DOUBLE_EQ(-21, b(2, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	513 ASSERT_DOUBLE_EQ(-7, b(2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	514 ASSERT_DOUBLE_EQ(-20, b(2, 3));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	515 ASSERT_DOUBLE_EQ(-28, b(3, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	516 ASSERT_DOUBLE_EQ(-47, b(3, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	517 ASSERT_DOUBLE_EQ(29, b(3, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	518 ASSERT_DOUBLE_EQ(-32, b(3, 3));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	519 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	520
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	521
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	522 TEST(FiniteProjectiveCamera, Calibration)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	523 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	524 unsigned int volumeWidth = 512;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	525 unsigned int volumeHeight = 512;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	526 unsigned int volumeDepth = 110;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	527
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	528 OrthancStone::Vector camera = OrthancStone::LinearAlgebra::CreateVector
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	529 (-1000, -5000, -static_cast<double>(volumeDepth) * 32);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	530
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	531 OrthancStone::Vector principalPoint = OrthancStone::LinearAlgebra::CreateVector
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	532 (volumeWidth/2, volumeHeight/2, volumeDepth * 2);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	533
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	534 OrthancStone::FiniteProjectiveCamera c(camera, principalPoint, 0, 512, 512, 1, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	535
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	536 double swapv[9] = { 1, 0, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	537 0, -1, 512,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	538 0, 0, 1 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	539 OrthancStone::Matrix swap;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	540 OrthancStone::LinearAlgebra::FillMatrix(swap, 3, 3, swapv);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	541
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	542 OrthancStone::Matrix p = OrthancStone::LinearAlgebra::Product(swap, c.GetMatrix());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	543 p /= p(2,3);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	544
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	545 ASSERT_NEAR( 1.04437, p(0,0), 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	546 ASSERT_NEAR(-0.0703111, p(0,1), 0.00000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	547 ASSERT_NEAR(-0.179283, p(0,2), 0.000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	548 ASSERT_NEAR( 61.7431, p(0,3), 0.0001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	549 ASSERT_NEAR( 0.11127, p(1,0), 0.000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	550 ASSERT_NEAR(-0.595541, p(1,1), 0.000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	551 ASSERT_NEAR( 0.872211, p(1,2), 0.000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	552 ASSERT_NEAR( 203.748, p(1,3), 0.001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	553 ASSERT_NEAR( 3.08593e-05, p(2,0), 0.0000000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	554 ASSERT_NEAR( 0.000129138, p(2,1), 0.000000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	555 ASSERT_NEAR( 9.18901e-05, p(2,2), 0.0000000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	556 ASSERT_NEAR( 1, p(2,3), 0.0000001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	557 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	558
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	559
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	560 static bool IsEqualRotationVector(OrthancStone::Vector a,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	561 OrthancStone::Vector b)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	562 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	563 if (a.size() != b.size() \|\|
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	564 a.size() != 3)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	565 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	566 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	567 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	568 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	569 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	570 OrthancStone::LinearAlgebra::NormalizeVector(a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	571 OrthancStone::LinearAlgebra::NormalizeVector(b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	572 return OrthancStone::LinearAlgebra::IsCloseToZero(boost::numeric::ublas::norm_2(a - b));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	573 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	574 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	575
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	576
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	577 TEST(GeometryToolbox, AlignVectorsWithRotation)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	578 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	579 OrthancStone::Vector a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	580 OrthancStone::Matrix r;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	581
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	582 OrthancStone::LinearAlgebra::AssignVector(a, -200, 200, -846.63);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	583 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	584
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	585 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	586 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	587 ASSERT_TRUE(IsEqualRotationVector(OrthancStone::LinearAlgebra::Product(r, a), b));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	588
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	589 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	590 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	591 ASSERT_TRUE(IsEqualRotationVector(OrthancStone::LinearAlgebra::Product(r, b), a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	592
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	593 OrthancStone::LinearAlgebra::AssignVector(a, 1, 0, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	594 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	595 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	596 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	597 ASSERT_TRUE(IsEqualRotationVector(OrthancStone::LinearAlgebra::Product(r, a), b));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	598
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	599 OrthancStone::LinearAlgebra::AssignVector(a, 0, 1, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	600 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	601 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	602 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	603 ASSERT_TRUE(IsEqualRotationVector(OrthancStone::LinearAlgebra::Product(r, a), b));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	604
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	605 OrthancStone::LinearAlgebra::AssignVector(a, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	606 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	607 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	608 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	609 ASSERT_TRUE(IsEqualRotationVector(OrthancStone::LinearAlgebra::Product(r, a), b));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	610
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	611 OrthancStone::LinearAlgebra::AssignVector(a, 0, 0, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	612 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	613 ASSERT_THROW(OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b), Orthanc::OrthancException);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	614
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	615 // TODO: Deal with opposite vectors
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	616
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	617 /*
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	618 OrthancStone::LinearAlgebra::AssignVector(a, 0, 0, -1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	619 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	620 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	621 OrthancStone::LinearAlgebra::Print(r);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	622 OrthancStone::LinearAlgebra::Print(boost::numeric::ublas::prod(r, a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	623 */
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	624 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	625
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	626
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	627 static bool IsEqualVectorL1(OrthancStone::Vector a,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	628 OrthancStone::Vector b)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	629 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	630 if (a.size() != b.size())
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	631 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	632 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	633 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	634 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	635 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	636 for (size_t i = 0; i < a.size(); i++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	637 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	638 if (!OrthancStone::LinearAlgebra::IsNear(a[i], b[i], 0.0001))
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	639 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	640 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	641 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	642 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	643
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	644 return true;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	645 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	646 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	647
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	648
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	649 TEST(VolumeImageGeometry, Basic)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	650 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	651 using namespace OrthancStone;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	652
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	653 VolumeImageGeometry g;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	654 g.SetSizeInVoxels(10, 20, 30);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	655 g.SetVoxelDimensions(1, 2, 3);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	656
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	657 Vector p = g.GetCoordinates(0, 0, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	658 ASSERT_EQ(3u, p.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	659 ASSERT_DOUBLE_EQ(-1.0 / 2.0, p[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	660 ASSERT_DOUBLE_EQ(-2.0 / 2.0, p[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	661 ASSERT_DOUBLE_EQ(-3.0 / 2.0, p[2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	662
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	663 p = g.GetCoordinates(1, 1, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	664 ASSERT_DOUBLE_EQ(-1.0 / 2.0 + 10.0 * 1.0, p[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	665 ASSERT_DOUBLE_EQ(-2.0 / 2.0 + 20.0 * 2.0, p[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	666 ASSERT_DOUBLE_EQ(-3.0 / 2.0 + 30.0 * 3.0, p[2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	667
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	668 VolumeProjection proj;
1766 ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	669 bool isOpposite;
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	670 ASSERT_TRUE(g.DetectProjection(proj, isOpposite, g.GetAxialGeometry().GetNormal()));
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	671 ASSERT_EQ(VolumeProjection_Axial, proj);
1766 ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	672 ASSERT_FALSE(isOpposite);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	673 ASSERT_TRUE(g.DetectProjection(proj, isOpposite, g.GetCoronalGeometry().GetNormal()));
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	674 ASSERT_EQ(VolumeProjection_Coronal, proj);
1766 ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	675 ASSERT_FALSE(isOpposite);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	676 ASSERT_TRUE(g.DetectProjection(proj, isOpposite, g.GetSagittalGeometry().GetNormal()));
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	677 ASSERT_EQ(VolumeProjection_Sagittal, proj);
1766 ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	678 ASSERT_FALSE(isOpposite);
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	679
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	680 ASSERT_EQ(10u, g.GetProjectionWidth(VolumeProjection_Axial));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	681 ASSERT_EQ(20u, g.GetProjectionHeight(VolumeProjection_Axial));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	682 ASSERT_EQ(30u, g.GetProjectionDepth(VolumeProjection_Axial));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	683 ASSERT_EQ(10u, g.GetProjectionWidth(VolumeProjection_Coronal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	684 ASSERT_EQ(30u, g.GetProjectionHeight(VolumeProjection_Coronal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	685 ASSERT_EQ(20u, g.GetProjectionDepth(VolumeProjection_Coronal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	686 ASSERT_EQ(20u, g.GetProjectionWidth(VolumeProjection_Sagittal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	687 ASSERT_EQ(30u, g.GetProjectionHeight(VolumeProjection_Sagittal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	688 ASSERT_EQ(10u, g.GetProjectionDepth(VolumeProjection_Sagittal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	689
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	690 p = g.GetVoxelDimensions(VolumeProjection_Axial);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	691 ASSERT_EQ(3u, p.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	692 ASSERT_DOUBLE_EQ(1, p[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	693 ASSERT_DOUBLE_EQ(2, p[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	694 ASSERT_DOUBLE_EQ(3, p[2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	695 p = g.GetVoxelDimensions(VolumeProjection_Coronal);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	696 ASSERT_EQ(3u, p.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	697 ASSERT_DOUBLE_EQ(1, p[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	698 ASSERT_DOUBLE_EQ(3, p[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	699 ASSERT_DOUBLE_EQ(2, p[2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	700 p = g.GetVoxelDimensions(VolumeProjection_Sagittal);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	701 ASSERT_EQ(3u, p.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	702 ASSERT_DOUBLE_EQ(2, p[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	703 ASSERT_DOUBLE_EQ(3, p[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	704 ASSERT_DOUBLE_EQ(1, p[2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	705
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	706 // Loop over all the voxels of the volume
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	707 for (unsigned int z = 0; z < g.GetDepth(); z++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	708 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	709 const float zz = (0.5f + static_cast<float>(z)) / static_cast<float>(g.GetDepth()); // Z-center of the voxel
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	710
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	711 for (unsigned int y = 0; y < g.GetHeight(); y++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	712 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	713 const float yy = (0.5f + static_cast<float>(y)) / static_cast<float>(g.GetHeight()); // Y-center of the voxel
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	714
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	715 for (unsigned int x = 0; x < g.GetWidth(); x++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	716 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	717 const float xx = (0.5f + static_cast<float>(x)) / static_cast<float>(g.GetWidth()); // X-center of the voxel
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	718
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	719 const float sx = 1.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	720 const float sy = 2.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	721 const float sz = 3.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	722
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	723 Vector p = g.GetCoordinates(xx, yy, zz);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	724
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	725 Vector q = (g.GetAxialGeometry().MapSliceToWorldCoordinates(
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	726 static_cast<double>(x) * sx,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	727 static_cast<double>(y) * sy) +
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	728 z * sz * g.GetAxialGeometry().GetNormal());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	729 ASSERT_TRUE(IsEqualVectorL1(p, q));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	730
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	731 q = (g.GetCoronalGeometry().MapSliceToWorldCoordinates(
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	732 static_cast<double>(x) * sx,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	733 static_cast<double>(g.GetDepth() - 1 - z) * sz) +
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	734 y * sy * g.GetCoronalGeometry().GetNormal());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	735 ASSERT_TRUE(IsEqualVectorL1(p, q));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	736
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	737 /**
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	738 * WARNING: In sagittal geometry, the normal points to
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	739 * REDUCING X-axis in the 3D world. This is necessary to keep
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	740 * the right-hand coordinate system. Hence the "-".
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	741 **/
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	742 q = (g.GetSagittalGeometry().MapSliceToWorldCoordinates(
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	743 static_cast<double>(y) * sy,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	744 static_cast<double>(g.GetDepth() - 1 - z) * sz) +
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	745 x * sx * (-g.GetSagittalGeometry().GetNormal()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	746 ASSERT_TRUE(IsEqualVectorL1(p, q));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	747 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	748 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	749 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	750
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	751 ASSERT_EQ(0, (int) VolumeProjection_Axial);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	752 ASSERT_EQ(1, (int) VolumeProjection_Coronal);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	753 ASSERT_EQ(2, (int) VolumeProjection_Sagittal);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	754
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	755 for (int p = 0; p < 3; p++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	756 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	757 VolumeProjection projection = (VolumeProjection) p;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	758 const CoordinateSystem3D& s = g.GetProjectionGeometry(projection);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	759
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	760 ASSERT_THROW(g.GetProjectionSlice(projection, g.GetProjectionDepth(projection)), Orthanc::OrthancException);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	761
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	762 for (unsigned int i = 0; i < g.GetProjectionDepth(projection); i++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	763 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	764 CoordinateSystem3D plane = g.GetProjectionSlice(projection, i);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	765
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	766 if (projection == VolumeProjection_Sagittal)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	767 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	768 ASSERT_TRUE(IsEqualVectorL1(plane.GetOrigin(), s.GetOrigin() + static_cast<double>(i) *
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	769 (-s.GetNormal()) * g.GetVoxelDimensions(projection)[2]));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	770 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	771 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	772 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	773 ASSERT_TRUE(IsEqualVectorL1(plane.GetOrigin(), s.GetOrigin() + static_cast<double>(i) *
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	774 s.GetNormal() * g.GetVoxelDimensions(projection)[2]));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	775 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	776
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	777 ASSERT_TRUE(IsEqualVectorL1(plane.GetAxisX(), s.GetAxisX()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	778 ASSERT_TRUE(IsEqualVectorL1(plane.GetAxisY(), s.GetAxisY()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	779
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	780 unsigned int slice;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	781 VolumeProjection q;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	782 ASSERT_TRUE(g.DetectSlice(q, slice, plane));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	783 ASSERT_EQ(projection, q);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	784 ASSERT_EQ(i, slice);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	785 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	786 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	787 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	788
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	789
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	790 TEST(LinearAlgebra, ParseVectorLocale)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	791 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	792 OrthancStone::Vector v;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	793
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	794 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "1.2"));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	795 ASSERT_EQ(1u, v.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	796 ASSERT_DOUBLE_EQ(1.2, v[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	797
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	798 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "-1.2e+2"));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	799 ASSERT_EQ(1u, v.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	800 ASSERT_DOUBLE_EQ(-120.0, v[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	801
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	802 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "-1e-2\\2"));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	803 ASSERT_EQ(2u, v.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	804 ASSERT_DOUBLE_EQ(-0.01, v[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	805 ASSERT_DOUBLE_EQ(2.0, v[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	806
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	807 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "1.3671875\\1.3671875"));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	808 ASSERT_EQ(2u, v.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	809 ASSERT_DOUBLE_EQ(1.3671875, v[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	810 ASSERT_DOUBLE_EQ(1.3671875, v[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	811 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	812
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	813 TEST(GenericToolbox, NormalizeUuid)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	814 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	815 std::string ref("44ca5051-14ef-4d2f-8bd7-db20bfb61fbb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	816
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	817 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	818 std::string u("44ca5051-14ef-4d2f-8bd7-db20bfb61fbb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	819 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	820 ASSERT_EQ(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	821 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	822
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	823 // space left
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	824 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	825 std::string u(" 44ca5051-14ef-4d2f-8bd7-db20bfb61fbb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	826 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	827 ASSERT_EQ(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	828 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	829
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	830 // space right
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	831 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	832 std::string u("44ca5051-14ef-4d2f-8bd7-db20bfb61fbb ");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	833 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	834 ASSERT_EQ(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	835 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	836
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	837 // space l & r
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	838 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	839 std::string u(" 44ca5051-14ef-4d2f-8bd7-db20bfb61fbb ");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	840 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	841 ASSERT_EQ(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	842 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	843
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	844 // space left + case
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	845 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	846 std::string u(" 44CA5051-14ef-4d2f-8bd7-dB20bfb61fbb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	847 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	848 ASSERT_EQ(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	849 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	850
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	851 // space right + case
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	852 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	853 std::string u("44ca5051-14EF-4D2f-8bd7-db20bfb61fbB ");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	854 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	855 ASSERT_EQ(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	856 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	857
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	858 // space l & r + case
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	859 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	860 std::string u(" 44cA5051-14Ef-4d2f-8bD7-db20bfb61fbb ");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	861 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	862 ASSERT_EQ(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	863 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	864
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	865 // no
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	866 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	867 std::string u(" 44ca5051-14ef-4d2f-8bd7- db20bfb61fbb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	868 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	869 ASSERT_NE(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	870 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	871
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	872 // no
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	873 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	874 std::string u("44ca5051-14ef-4d2f-8bd7-db20bfb61fb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	875 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	876 ASSERT_NE(ref, u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	877 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	878 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	879
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	880
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	881 TEST(CoordinateSystem3D, Basic)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	882 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	883 using namespace OrthancStone;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	884
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	885 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	886 CoordinateSystem3D c;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	887 ASSERT_FALSE(c.IsValid());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	888 ASSERT_DOUBLE_EQ(c.GetNormal()[0], 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	889 ASSERT_DOUBLE_EQ(c.GetNormal()[1], 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	890 ASSERT_DOUBLE_EQ(c.GetNormal()[2], 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	891
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	892 ASSERT_DOUBLE_EQ(0, c.ComputeDistance(LinearAlgebra::CreateVector(0, 0, 0)));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	893 ASSERT_DOUBLE_EQ(0, c.ComputeDistance(LinearAlgebra::CreateVector(5, 0, 0)));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	894 ASSERT_DOUBLE_EQ(0, c.ComputeDistance(LinearAlgebra::CreateVector(0, 5, 0)));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	895 ASSERT_DOUBLE_EQ(5, c.ComputeDistance(LinearAlgebra::CreateVector(0, 0, 5)));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	896 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	897
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	898 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	899 CoordinateSystem3D c("nope1", "nope2");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	900 ASSERT_FALSE(c.IsValid());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	901 ASSERT_DOUBLE_EQ(c.GetNormal()[0], 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	902 ASSERT_DOUBLE_EQ(c.GetNormal()[1], 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	903 ASSERT_DOUBLE_EQ(c.GetNormal()[2], 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	904 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	905
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	906 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	907 // https://www.vedantu.com/maths/perpendicular-distance-of-a-point-from-a-plane
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	908 CoordinateSystem3D c = CoordinateSystem3D::CreateFromPlaneGeneralForm(2, 4, -4, -6);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	909 ASSERT_DOUBLE_EQ(3, c.ComputeDistance(LinearAlgebra::CreateVector(0, 3, 6)));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	910 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	911
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	912 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	913 // https://mathinsight.org/distance_point_plane_examples
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	914 CoordinateSystem3D c = CoordinateSystem3D::CreateFromPlaneGeneralForm(2, -2, 5, 8);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	915 ASSERT_DOUBLE_EQ(39.0 / sqrt(33.0), c.ComputeDistance(LinearAlgebra::CreateVector(4, -4, 3)));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	916 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	917
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	918 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	919 // https://www.ck12.org/calculus/distance-between-a-point-and-a-plane/lesson/Distance-Between-a-Point-and-a-Plane-MAT-ALY/
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	920 const Vector a = LinearAlgebra::CreateVector(3, 6, 9);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	921 const Vector b = LinearAlgebra::CreateVector(9, 6, 3);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	922 const Vector c = LinearAlgebra::CreateVector(6, -9, 9);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	923 CoordinateSystem3D d = CoordinateSystem3D::CreateFromThreePoints(a, b, c);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	924 ASSERT_DOUBLE_EQ(0, d.ComputeDistance(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	925 ASSERT_DOUBLE_EQ(0, d.ComputeDistance(b));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	926 ASSERT_DOUBLE_EQ(0, d.ComputeDistance(c));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	927 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	928
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	929 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	930 // https://tutorial.math.lamar.edu/classes/calcii/eqnsofplanes.aspx
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	931 const Vector a = LinearAlgebra::CreateVector(1, -2, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	932 const Vector b = LinearAlgebra::CreateVector(3, 1, 4);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	933 const Vector c = LinearAlgebra::CreateVector(0, -1, 2);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	934 CoordinateSystem3D d = CoordinateSystem3D::CreateFromThreePoints(a, b, c);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	935 double r = d.GetNormal() [0] / 2.0;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	936 ASSERT_DOUBLE_EQ(-8 * r, d.GetNormal() [1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	937 ASSERT_DOUBLE_EQ(5 * r, d.GetNormal() [2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	938 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	939 }
1765 177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	940
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	941
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	942 TEST(SlicesSorter, HFP)
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	943 {
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	944 // 2021-04-27-repro-bug-HFP-HFS-cartman
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	945
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	946 {
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	947 // This is the last instance in the CT series ("InstanceNumber" is 368):
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	948 // CT1.2.752.243.1.1.20210202150623868.3730.61448.dcm
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	949 const OrthancStone::CoordinateSystem3D system("300\\302.5\\323.11", "-1\\0\\0\\0\\-1\\0");
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	950
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	951 /**
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	952 * The first instance in the series ("InstanceNumber" is 1) is
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	953 * CT1.2.752.243.1.1.20210202150623381.2000.76318.dcm, and its
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	954 * "ImagePositionPatient" is "300\\302.5\\690.11". It cannot be
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	955 * taken as the origin of the volume, otherwise the Z axis is
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	956 * shifted by the depth of the volume.
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	957 **/
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	958
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	959 const double spacingXY = 1.171875;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	960 const double width = 512;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	961 const double height = 512;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	962 const double depth = 368; // Number of instances in the series
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	963
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	964 OrthancStone::VolumeImageGeometry geometry;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	965 geometry.SetAxialGeometry(system);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	966 geometry.SetSizeInVoxels(width, height, depth);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	967 geometry.SetVoxelDimensions(spacingXY, spacingXY, 1 /* pixel spacing Z */);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	968
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	969 OrthancStone::Vector p, q;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	970
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	971 p = OrthancStone::LinearAlgebra::CreateVector(0.5 / double(width),
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	972 0.5 / double(height),
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	973 0.5 / double(depth), 1);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	974 q = OrthancStone::LinearAlgebra::Product(geometry.GetTransform(), p);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	975 ASSERT_FLOAT_EQ(300, q[0]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	976 ASSERT_FLOAT_EQ(302.5, q[1]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	977 ASSERT_FLOAT_EQ(323.11, q[2]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	978 ASSERT_FLOAT_EQ(1, q[3]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	979
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	980 p = OrthancStone::LinearAlgebra::CreateVector((width - 0.5) / double(width),
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	981 (height - 0.5) / double(height),
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	982 (depth - 0.5) / double(depth), 1);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	983 q = OrthancStone::LinearAlgebra::Product(geometry.GetTransform(), p);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	984
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	985 ASSERT_FLOAT_EQ(300.0 - (width - 1.0) * spacingXY, q[0]); // "X" is swapped
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	986 ASSERT_FLOAT_EQ(302.5 - (height - 1.0) * spacingXY, q[1]); // "Y" is swapped
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	987 ASSERT_FLOAT_EQ(690.11, q[2]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	988 ASSERT_FLOAT_EQ(1, q[3]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	989 }
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	990
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	991 {
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	992 // DOSE instance: RD1.2.752.243.1.1.20210202150624529.3790.85357_DoseTPS.dcm
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	993 OrthancStone::CoordinateSystem3D system("-217.0492\\-161.4141\\376.61", "1\\0\\0\\0\\1\\0");
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	994 const double spacingXY = 3;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	995 const double width = 146; // Columns
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	996 const double height = 84; // Row
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	997 const double depth = 86; // Number of frames, same as the length of "GridFrameOffsetVector"
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	998
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	999 OrthancStone::VolumeImageGeometry geometry;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1000 geometry.SetAxialGeometry(system);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1001 geometry.SetSizeInVoxels(width, height, depth);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1002 geometry.SetVoxelDimensions(spacingXY, spacingXY, 3 /* pixel spacing Z, cf. "GridFrameOffsetVector" */);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1003
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1004 OrthancStone::Vector p, q;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1005
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1006 p = OrthancStone::LinearAlgebra::CreateVector(0.5 / double(width), 0.5 / double(height), 0.5 / double(depth), 1);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1007 q = OrthancStone::LinearAlgebra::Product(geometry.GetTransform(), p);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1008 ASSERT_FLOAT_EQ(-217.0492, q[0]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1009 ASSERT_FLOAT_EQ(-161.4141, q[1]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1010 ASSERT_FLOAT_EQ(376.61, q[2]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1011 ASSERT_FLOAT_EQ(1, q[3]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1012
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1013 p = OrthancStone::LinearAlgebra::CreateVector((width - 0.5) / double(width),
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1014 (height - 0.5) / double(height),
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1015 (depth - 0.5) / double(depth), 1);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1016 q = OrthancStone::LinearAlgebra::Product(geometry.GetTransform(), p);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1017
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1018 ASSERT_FLOAT_EQ(-217.0492 + (width - 1.0) * spacingXY, q[0]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1019 ASSERT_FLOAT_EQ(-161.4141 + (height - 1.0) * spacingXY, q[1]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1020 ASSERT_FLOAT_EQ(376.61 + 255.0 /* last item in "GridFrameOffsetVector" */, q[2]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1021 ASSERT_FLOAT_EQ(1, q[3]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1022 }
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1023
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1024 for (unsigned int upward = 0; upward < 2; upward++)
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1025 {
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1026 OrthancStone::SlicesSorter slices;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1027
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1028 for (unsigned int i = 0; i < 368; i++)
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1029 {
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1030 unsigned int z = (upward ? 323 + i : 690 - i);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1031 OrthancStone::CoordinateSystem3D p("300\\302.5\\" + boost::lexical_cast<std::string>(z) + ".11",
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1032 "-1\\0\\0\\0\\-1\\0");
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1033 slices.AddSlice(p, new Orthanc::SingleValueObject<unsigned int>(z));
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1034 }
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1035
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1036 slices.Sort();
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1037
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1038 double spacingZ;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1039 ASSERT_TRUE(slices.ComputeSpacingBetweenSlices(spacingZ));
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1040 ASSERT_FLOAT_EQ(1, spacingZ);
1766 ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1041 ASSERT_TRUE(slices.AreAllSlicesDistinct());
1765 177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1042
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1043 ASSERT_FLOAT_EQ(300.0, slices.GetSliceGeometry(0).GetOrigin() [0]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1044 ASSERT_FLOAT_EQ(302.5, slices.GetSliceGeometry(0).GetOrigin() [1]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1045 ASSERT_FLOAT_EQ(323.11, slices.GetSliceGeometry(0).GetOrigin() [2]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1046
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1047 ASSERT_FLOAT_EQ(300.0, slices.GetSliceGeometry(367).GetOrigin() [0]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1048 ASSERT_FLOAT_EQ(302.5, slices.GetSliceGeometry(367).GetOrigin() [1]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1049 ASSERT_FLOAT_EQ(690.11, slices.GetSliceGeometry(367).GetOrigin() [2]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1050
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1051 ASSERT_EQ(323u, dynamic_cast<const Orthanc::SingleValueObject<unsigned int>&>(slices.GetSlicePayload(0)).GetValue());
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1052 ASSERT_EQ(690u, dynamic_cast<const Orthanc::SingleValueObject<unsigned int>&>(slices.GetSlicePayload(367)).GetValue());
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1053
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1054 const unsigned int width = 512;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1055 const unsigned int height = 512;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1056 const unsigned int depth = 368;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1057 const double spacingXY = 1.171875;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1058
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1059 OrthancStone::VolumeImageGeometry geometry;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1060 geometry.SetSizeInVoxels(width, height, depth);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1061 geometry.SetVoxelDimensions(spacingXY, spacingXY, 1);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1062 geometry.SetAxialGeometry(slices.GetSliceGeometry(0));
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1063
1766 ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1064 OrthancStone::Vector q = geometry.GetCoordinates(0, 0, 0);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1065 ASSERT_EQ(3u, q.size());
1765 177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1066 ASSERT_FLOAT_EQ(300 + spacingXY / 2.0, q[0]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1067 ASSERT_FLOAT_EQ(302.5 + spacingXY / 2.0, q[1]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1068 ASSERT_FLOAT_EQ(323.11 - 0.5, q[2]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1069
1766 ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1070 q = geometry.GetCoordinates(1, 1, 1);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1071 ASSERT_EQ(3u, q.size());
1765 177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1072 ASSERT_FLOAT_EQ(300 - double(width) * spacingXY + spacingXY / 2.0, q[0]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1073 ASSERT_FLOAT_EQ(302.5 - double(height) * spacingXY + spacingXY / 2.0, q[1]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1074 ASSERT_FLOAT_EQ(323.11 + 368.0 - 0.5, q[2]);
1766 ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1075
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1076 OrthancStone::VolumeProjection projection;
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1077 unsigned int slice;
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1078 ASSERT_TRUE(geometry.DetectSlice(projection, slice, OrthancStone::CoordinateSystem3D("300\\302.5\\690.11", "-1\\0\\0\\0\\-1\\0")));
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1079 ASSERT_EQ(OrthancStone::VolumeProjection_Axial, projection);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1080 ASSERT_EQ(367u, slice);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1081 ASSERT_TRUE(geometry.DetectSlice(projection, slice, OrthancStone::CoordinateSystem3D("300\\302.5\\323.11", "-1\\0\\0\\0\\-1\\0")));
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1082 ASSERT_EQ(OrthancStone::VolumeProjection_Axial, projection);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1083 ASSERT_EQ(0u, slice);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1084
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1085 // DOSE instance: RD1.2.752.243.1.1.20210202150624529.3790.85357_DoseTPS.dcm
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1086 ASSERT_TRUE(geometry.DetectSlice(projection, slice,
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1087 OrthancStone::CoordinateSystem3D("-217.0492\\-161.4141\\376.61", "1\\0\\0\\0\\1\\0")));
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1088 ASSERT_EQ(OrthancStone::VolumeProjection_Axial, projection);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1089 ASSERT_EQ(376u - 323u, slice);
1765 177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1090 }
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1091 }

Mercurial > hg > orthanc-stone

annotate UnitTestsSources/GeometryToolboxTests.cpp @ 1864:35ac2be493e2