orthanc-stone: UnitTestsSources/GeometryToolboxTests.cpp annotate

annotate UnitTestsSources/GeometryToolboxTests.cpp @ 1871:7053b8a0aaec

upgrade to year 2022

author	Sebastien Jodogne <s.jodogne@gmail.com>
date	Tue, 11 Jan 2022 11:18:48 +0100
parents	3889ae96d2e9
children

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
1871 7053b8a0aaec upgrade to year 2022 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1870 diff changeset	5 * Copyright (C) 2017-2022 Osimis S.A., Belgium
7053b8a0aaec upgrade to year 2022 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1870 diff changeset	6 * Copyright (C) 2021-2022 Sebastien Jodogne, ICTEAM UCLouvain, Belgium
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	7 *
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	8 * 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	9 * 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	10 * 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	11 * 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	12 *
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	13 * 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	14 * 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	15 * 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	16 * Affero General Public License for more details.
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	17 *
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	18 * 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	19 * 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	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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	23 #include <gtest/gtest.h>
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	24
1765 177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	25 #include "../OrthancStone/Sources/Toolbox/DicomInstanceParameters.h"
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	26 #include "../OrthancStone/Sources/Toolbox/FiniteProjectiveCamera.h"
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	27 #include "../OrthancStone/Sources/Toolbox/GenericToolbox.h"
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	28 #include "../OrthancStone/Sources/Toolbox/GeometryToolbox.h"
1765 177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	29 #include "../OrthancStone/Sources/Toolbox/SlicesSorter.h"
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	30
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	31 #include <Logging.h>
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	32 #include <OrthancException.h>
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	33
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	34 #include <boost/math/special_functions/round.hpp>
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	37 TEST(GeometryToolbox, Interpolation)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	38 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	39 using namespace OrthancStone::GeometryToolbox;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	40
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	41 // 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	42 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	43
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	44 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	45 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	46 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	47 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	48
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	49 ASSERT_DOUBLE_EQ(123.35f, ComputeTrilinearInterpolationUnitSquare
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	50 (0.5f, 0.2f, 0.7f,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	51 91, 210, 162, 95,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	52 51, 190, 80, 92));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	53
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	54 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	55 ComputeTrilinearInterpolationUnitSquare(0.5f, 0.2f, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	56 91, 210, 162, 95,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	57 51, 190, 80, 92));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	58
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	59 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	60 ComputeTrilinearInterpolationUnitSquare(0.5f, 0.2f, 1,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	61 91, 210, 162, 95,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	62 51, 190, 80, 92));
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	66 static bool CompareMatrix(const OrthancStone::Matrix& a,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	67 const OrthancStone::Matrix& b,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	68 double threshold = 0.00000001)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	69 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	70 if (a.size1() != b.size1() \|\|
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	71 a.size2() != b.size2())
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	72 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	73 return false;
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	76 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	77 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	78 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	79 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	80 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	81 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	82 LOG(ERROR) << "Too large difference in component ("
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	83 << i << "," << j << "): " << a(i,j) << " != " << b(i,j);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	84 return false;
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	89 return true;
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	93 static bool CompareVector(const OrthancStone::Vector& a,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	94 const OrthancStone::Vector& b,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	95 double threshold = 0.00000001)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	96 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	97 if (a.size() != b.size())
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	98 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	99 return false;
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	102 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	103 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	104 if (fabs(a(i) - b(i)) > threshold)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	105 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	106 LOG(ERROR) << "Too large difference in component "
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	107 << i << ": " << a(i) << " != " << b(i);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	108 return false;
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	112 return true;
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	117 TEST(FiniteProjectiveCamera, Decomposition1)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	118 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	119 // 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	120 // edition" (page 163)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	121 const double p[12] = {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	122 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	123 -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	124 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	125 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	126
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	127 OrthancStone::FiniteProjectiveCamera camera(p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	128 ASSERT_EQ(3u, camera.GetMatrix().size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	129 ASSERT_EQ(4u, camera.GetMatrix().size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	130 ASSERT_EQ(3u, camera.GetIntrinsicParameters().size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	131 ASSERT_EQ(3u, camera.GetIntrinsicParameters().size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	132 ASSERT_EQ(3u, camera.GetRotation().size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	133 ASSERT_EQ(3u, camera.GetRotation().size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	134 ASSERT_EQ(3u, camera.GetCenter().size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	135
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	136 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	137 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	138 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	139
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	140 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	141 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	142 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	143 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	144 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	145 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	146
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	147 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	148 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	149 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	150
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	151 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	152 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	153 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	154 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	155 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	156 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	157 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	158 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	159 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	160
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	161 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	162
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	163 OrthancStone::FiniteProjectiveCamera camera2(camera.GetIntrinsicParameters(),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	164 camera.GetRotation(),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	165 camera.GetCenter());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	166
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	167 ASSERT_TRUE(CompareMatrix(camera.GetMatrix(), camera2.GetMatrix()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	168 ASSERT_TRUE(CompareMatrix(camera.GetIntrinsicParameters(), camera2.GetIntrinsicParameters()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	169 ASSERT_TRUE(CompareMatrix(camera.GetRotation(), camera2.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	170 ASSERT_TRUE(CompareVector(camera.GetCenter(), camera2.GetCenter()));
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	174 TEST(FiniteProjectiveCamera, Decomposition2)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	175 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	176 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	177 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	178 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	179 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	180
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	181 OrthancStone::FiniteProjectiveCamera camera(p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	182 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	183
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	184 OrthancStone::FiniteProjectiveCamera camera2(k, r, c);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	185 ASSERT_TRUE(CompareMatrix(camera.GetMatrix(), camera2.GetMatrix(), 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	186 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	187 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	188 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	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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	192 TEST(FiniteProjectiveCamera, Decomposition3)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	193 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	194 const double p[] = { 10, 0, 0, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	195 0, 20, 0, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	196 0, 0, 30, 0 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	197
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	198 OrthancStone::FiniteProjectiveCamera camera(p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	199 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	200 ASSERT_DOUBLE_EQ(10, camera.GetIntrinsicParameters() (0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	201 ASSERT_DOUBLE_EQ(20, camera.GetIntrinsicParameters() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	202 ASSERT_DOUBLE_EQ(30, camera.GetIntrinsicParameters() (2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	203 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	204 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	205 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	206 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	207 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	208 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (2));
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	212 TEST(FiniteProjectiveCamera, Decomposition4)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	213 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	214 const double p[] = { 1, 0, 0, 10,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	215 0, 1, 0, 20,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	216 0, 0, 1, 30 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	217
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	218 OrthancStone::FiniteProjectiveCamera camera(p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	219 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	220 ASSERT_DOUBLE_EQ(1, camera.GetIntrinsicParameters() (0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	221 ASSERT_DOUBLE_EQ(1, camera.GetIntrinsicParameters() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	222 ASSERT_DOUBLE_EQ(1, camera.GetIntrinsicParameters() (2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	223 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	224 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	225 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	226 ASSERT_DOUBLE_EQ(-10, camera.GetCenter() (0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	227 ASSERT_DOUBLE_EQ(-20, camera.GetCenter() (1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	228 ASSERT_DOUBLE_EQ(-30, camera.GetCenter() (2));
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	232 TEST(FiniteProjectiveCamera, Decomposition5)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	233 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	234 const double p[] = { 0, 0, 10, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	235 0, 20, 0, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	236 30, 0, 0, 0 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	237
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	238 OrthancStone::FiniteProjectiveCamera camera(p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	239 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	240 ASSERT_DOUBLE_EQ(-10, camera.GetIntrinsicParameters() (0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	241 ASSERT_DOUBLE_EQ(20, camera.GetIntrinsicParameters() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	242 ASSERT_DOUBLE_EQ(30, camera.GetIntrinsicParameters() (2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	243 ASSERT_DOUBLE_EQ(-1, camera.GetRotation() (0, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	244 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	245 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (2, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	246 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	247 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	248 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	249
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	250 OrthancStone::FiniteProjectiveCamera camera2(camera.GetIntrinsicParameters(),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	251 camera.GetRotation(),
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	252 camera.GetCenter());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	253 ASSERT_TRUE(CompareMatrix(camera.GetMatrix(), camera2.GetMatrix()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	254 ASSERT_TRUE(CompareMatrix(camera.GetIntrinsicParameters(), camera2.GetIntrinsicParameters()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	255 ASSERT_TRUE(CompareMatrix(camera.GetRotation(), camera2.GetRotation()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	256 ASSERT_TRUE(CompareVector(camera.GetCenter(), camera2.GetCenter()));
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	260 static double GetCosAngle(const OrthancStone::Vector& a,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	261 const OrthancStone::Vector& b)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	262 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	263 // 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	264 // 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	265 return boost::numeric::ublas::inner_prod(a, b) /
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	266 (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	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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	270 TEST(FiniteProjectiveCamera, Ray)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	271 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	272 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	273 -2951.517707, -1501.019129, -285.785281, 637891.819097,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	274 0.008528, 0.003067, -0.999959, 2491.764918 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	275
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	276 const OrthancStone::FiniteProjectiveCamera camera(pp);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	277
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	278 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	279 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	280 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	281
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	282 // 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	283 // (256,256). The image has dimensions 512x512.
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	284 OrthancStone::Vector o =
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	285 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	286 OrthancStone::Vector ax =
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	287 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	288 OrthancStone::Vector ay =
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	289 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	290
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	291 OrthancStone::CoordinateSystem3D imagePlane(o, ax, ay);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	292
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	293 // Back-projection of the principal point
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	294 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	295 OrthancStone::Vector ray = camera.GetRayDirection(256, 256);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	296
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	297 // 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	298 // (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	299 // ("-1" corresponds to "cos(pi)").
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	300 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	301
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	302 // 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	303 double x, y;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	304 camera.ApplyFinite(x, y, camera.GetCenter() - ray);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	305
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	306 ASSERT_NEAR(256, x, 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	307 ASSERT_NEAR(256, y, 0.00001);
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	310 // 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	311 std::vector<double> cx, cy;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	312 cx.push_back(0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	313 cy.push_back(0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	314 cx.push_back(512);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	315 cy.push_back(0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	316 cx.push_back(512);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	317 cy.push_back(512);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	318 cx.push_back(0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	319 cy.push_back(512);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	320
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	321 bool first = true;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	322 double angle;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	323
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	324 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	325 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	326 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	327
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	328 // 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	329 // the 4 corners
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	330 double a = GetCosAngle(ray, imagePlane.GetNormal());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	331 if (first)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	332 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	333 first = false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	334 angle = a;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	335 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	336 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	337 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	338 ASSERT_NEAR(angle, a, 0.000001);
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	341 // 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	342 double x, y;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	343 camera.ApplyFinite(x, y, camera.GetCenter() - ray);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	344
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	345 ASSERT_NEAR(cx[i], x, 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	346 ASSERT_NEAR(cy[i], y, 0.00001);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	347
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	348 // Alternative construction, by computing the intersection of the
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	349 // ray with the image plane
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	350 OrthancStone::Vector p;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	351 ASSERT_TRUE(imagePlane.IntersectLine(p, camera.GetCenter(), -ray));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	352 imagePlane.ProjectPoint(x, y, p);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	353 ASSERT_NEAR(cx[i], x + 256, 0.01);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	354 ASSERT_NEAR(cy[i], y + 256, 0.01);
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	359 TEST(Matrix, Inverse1)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	360 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	361 OrthancStone::Matrix a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	362
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	363 a.resize(0, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	364 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	365 ASSERT_EQ(0u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	366 ASSERT_EQ(0u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	367
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	368 a.resize(2, 3);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	369 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	370
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	371 a.resize(1, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	372 a(0, 0) = 45.0;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	373
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	374 ASSERT_DOUBLE_EQ(45, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	375 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	376 ASSERT_EQ(1u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	377 ASSERT_EQ(1u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	378 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	379
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	380 a(0, 0) = 0;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	381 ASSERT_DOUBLE_EQ(0, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	382 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	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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	386 TEST(Matrix, Inverse2)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	387 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	388 OrthancStone::Matrix a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	389 a.resize(2, 2);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	390 a(0, 0) = 4;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	391 a(0, 1) = 3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	392 a(1, 0) = 3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	393 a(1, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	394
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	395 ASSERT_DOUBLE_EQ(-1, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	396 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	397 ASSERT_EQ(2u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	398 ASSERT_EQ(2u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	399
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	400 ASSERT_DOUBLE_EQ(-2, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	401 ASSERT_DOUBLE_EQ(3, b(0, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	402 ASSERT_DOUBLE_EQ(3, b(1, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	403 ASSERT_DOUBLE_EQ(-4, b(1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	404
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	405 a(0, 0) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	406 a(0, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	407 a(1, 0) = 3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	408 a(1, 1) = 4;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	409
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	410 ASSERT_DOUBLE_EQ(-2, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	411 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	412
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	413 ASSERT_DOUBLE_EQ(-2, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	414 ASSERT_DOUBLE_EQ(1, b(0, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	415 ASSERT_DOUBLE_EQ(1.5, b(1, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	416 ASSERT_DOUBLE_EQ(-0.5, b(1, 1));
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	420 TEST(Matrix, Inverse3)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	421 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	422 OrthancStone::Matrix a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	423 a.resize(3, 3);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	424 a(0, 0) = 7;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	425 a(0, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	426 a(0, 2) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	427 a(1, 0) = 0;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	428 a(1, 1) = 3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	429 a(1, 2) = -1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	430 a(2, 0) = -3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	431 a(2, 1) = 4;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	432 a(2, 2) = -2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	433
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	434 ASSERT_DOUBLE_EQ(1, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	435 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	436 ASSERT_EQ(3u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	437 ASSERT_EQ(3u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	438
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	439 ASSERT_DOUBLE_EQ(-2, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	440 ASSERT_DOUBLE_EQ(8, b(0, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	441 ASSERT_DOUBLE_EQ(-5, b(0, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	442 ASSERT_DOUBLE_EQ(3, b(1, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	443 ASSERT_DOUBLE_EQ(-11, b(1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	444 ASSERT_DOUBLE_EQ(7, b(1, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	445 ASSERT_DOUBLE_EQ(9, b(2, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	446 ASSERT_DOUBLE_EQ(-34, b(2, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	447 ASSERT_DOUBLE_EQ(21, b(2, 2));
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	450 a(0, 0) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	451 a(0, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	452 a(0, 2) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	453 a(1, 0) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	454 a(1, 1) = 0;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	455 a(1, 2) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	456 a(2, 0) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	457 a(2, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	458 a(2, 2) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	459
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	460 ASSERT_DOUBLE_EQ(2, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	461 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	462 ASSERT_EQ(3u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	463 ASSERT_EQ(3u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	464
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	465 ASSERT_DOUBLE_EQ(-1, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	466 ASSERT_DOUBLE_EQ(1, b(0, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	467 ASSERT_DOUBLE_EQ(1, b(0, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	468 ASSERT_DOUBLE_EQ(0, b(1, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	469 ASSERT_DOUBLE_EQ(-0.5, b(1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	470 ASSERT_DOUBLE_EQ(0.5, b(1, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	471 ASSERT_DOUBLE_EQ(1, b(2, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	472 ASSERT_DOUBLE_EQ(0, b(2, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	473 ASSERT_DOUBLE_EQ(-1, b(2, 2));
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	477 TEST(Matrix, Inverse4)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	478 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	479 OrthancStone::Matrix a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	480 a.resize(4, 4);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	481 a(0, 0) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	482 a(0, 1) = 1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	483 a(0, 2) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	484 a(0, 3) = -3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	485 a(1, 0) = -2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	486 a(1, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	487 a(1, 2) = -1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	488 a(1, 3) = -1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	489 a(2, 0) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	490 a(2, 1) = 2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	491 a(2, 2) = -3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	492 a(2, 3) = -1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	493 a(3, 0) = 3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	494 a(3, 1) = -2;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	495 a(3, 2) = -3;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	496 a(3, 3) = -1;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	497
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	498 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	499 ASSERT_EQ(4u, b.size1());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	500 ASSERT_EQ(4u, b.size2());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	501
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	502 b *= 134.0; // This is the determinant
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	503
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	504 ASSERT_DOUBLE_EQ(8, b(0, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	505 ASSERT_DOUBLE_EQ(-44, b(0, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	506 ASSERT_DOUBLE_EQ(30, b(0, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	507 ASSERT_DOUBLE_EQ(-10, b(0, 3));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	508 ASSERT_DOUBLE_EQ(2, b(1, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	509 ASSERT_DOUBLE_EQ(-11, b(1, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	510 ASSERT_DOUBLE_EQ(41, b(1, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	511 ASSERT_DOUBLE_EQ(-36, b(1, 3));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	512 ASSERT_DOUBLE_EQ(16, b(2, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	513 ASSERT_DOUBLE_EQ(-21, b(2, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	514 ASSERT_DOUBLE_EQ(-7, b(2, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	515 ASSERT_DOUBLE_EQ(-20, b(2, 3));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	516 ASSERT_DOUBLE_EQ(-28, b(3, 0));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	517 ASSERT_DOUBLE_EQ(-47, b(3, 1));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	518 ASSERT_DOUBLE_EQ(29, b(3, 2));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	519 ASSERT_DOUBLE_EQ(-32, b(3, 3));
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	523 TEST(FiniteProjectiveCamera, Calibration)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	524 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	525 unsigned int volumeWidth = 512;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	526 unsigned int volumeHeight = 512;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	527 unsigned int volumeDepth = 110;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	528
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	529 OrthancStone::Vector camera = OrthancStone::LinearAlgebra::CreateVector
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	530 (-1000, -5000, -static_cast<double>(volumeDepth) * 32);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	531
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	532 OrthancStone::Vector principalPoint = OrthancStone::LinearAlgebra::CreateVector
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	533 (volumeWidth/2, volumeHeight/2, volumeDepth * 2);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	534
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	535 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	536
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	537 double swapv[9] = { 1, 0, 0,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	538 0, -1, 512,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	539 0, 0, 1 };
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	540 OrthancStone::Matrix swap;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	541 OrthancStone::LinearAlgebra::FillMatrix(swap, 3, 3, swapv);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	542
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	543 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	544 p /= p(2,3);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	545
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	546 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	547 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	548 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	549 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	550 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	551 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	552 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	553 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	554 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	555 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	556 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	557 ASSERT_NEAR( 1, p(2,3), 0.0000001);
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	561 static bool IsEqualRotationVector(OrthancStone::Vector a,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	562 OrthancStone::Vector b)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	563 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	564 if (a.size() != b.size() \|\|
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	565 a.size() != 3)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	566 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	567 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	568 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	569 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	570 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	571 OrthancStone::LinearAlgebra::NormalizeVector(a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	572 OrthancStone::LinearAlgebra::NormalizeVector(b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	573 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	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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	578 TEST(GeometryToolbox, AlignVectorsWithRotation)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	579 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	580 OrthancStone::Vector a, b;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	581 OrthancStone::Matrix r;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	582
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	583 OrthancStone::LinearAlgebra::AssignVector(a, -200, 200, -846.63);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	584 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	585
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	586 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	587 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	588 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	589
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	590 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, b, a);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	591 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	592 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	593
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	594 OrthancStone::LinearAlgebra::AssignVector(a, 1, 0, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	595 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	596 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	597 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	598 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	599
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	600 OrthancStone::LinearAlgebra::AssignVector(a, 0, 1, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	601 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	602 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	603 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	604 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	605
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	606 OrthancStone::LinearAlgebra::AssignVector(a, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	607 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	608 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	609 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	610 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	611
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	612 OrthancStone::LinearAlgebra::AssignVector(a, 0, 0, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	613 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	614 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	615
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	616 // TODO: Deal with opposite vectors
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 /*
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	619 OrthancStone::LinearAlgebra::AssignVector(a, 0, 0, -1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	620 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	621 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	622 OrthancStone::LinearAlgebra::Print(r);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	623 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	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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	628 static bool IsEqualVectorL1(OrthancStone::Vector a,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	629 OrthancStone::Vector b)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	630 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	631 if (a.size() != b.size())
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	632 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	633 return false;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	634 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	635 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	636 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	637 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	638 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	639 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	640 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	641 return false;
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	645 return true;
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	650 TEST(VolumeImageGeometry, Basic)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	651 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	652 using namespace OrthancStone;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	653
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	654 VolumeImageGeometry g;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	655 g.SetSizeInVoxels(10, 20, 30);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	656 g.SetVoxelDimensions(1, 2, 3);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	657
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	658 Vector p = g.GetCoordinates(0, 0, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	659 ASSERT_EQ(3u, p.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	660 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	661 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	662 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	663
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	664 p = g.GetCoordinates(1, 1, 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	665 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	666 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	667 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	668
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	669 VolumeProjection proj;
1766 ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	670 bool isOpposite;
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	671 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	672 ASSERT_EQ(VolumeProjection_Axial, proj);
1766 ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	673 ASSERT_FALSE(isOpposite);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	674 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	675 ASSERT_EQ(VolumeProjection_Coronal, proj);
1766 ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	676 ASSERT_FALSE(isOpposite);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	677 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	678 ASSERT_EQ(VolumeProjection_Sagittal, proj);
1766 ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	679 ASSERT_FALSE(isOpposite);
1679 5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	680
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	681 ASSERT_EQ(10u, g.GetProjectionWidth(VolumeProjection_Axial));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	682 ASSERT_EQ(20u, g.GetProjectionHeight(VolumeProjection_Axial));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	683 ASSERT_EQ(30u, g.GetProjectionDepth(VolumeProjection_Axial));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	684 ASSERT_EQ(10u, g.GetProjectionWidth(VolumeProjection_Coronal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	685 ASSERT_EQ(30u, g.GetProjectionHeight(VolumeProjection_Coronal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	686 ASSERT_EQ(20u, g.GetProjectionDepth(VolumeProjection_Coronal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	687 ASSERT_EQ(20u, g.GetProjectionWidth(VolumeProjection_Sagittal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	688 ASSERT_EQ(30u, g.GetProjectionHeight(VolumeProjection_Sagittal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	689 ASSERT_EQ(10u, g.GetProjectionDepth(VolumeProjection_Sagittal));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	690
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	691 p = g.GetVoxelDimensions(VolumeProjection_Axial);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	692 ASSERT_EQ(3u, p.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	693 ASSERT_DOUBLE_EQ(1, p[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	694 ASSERT_DOUBLE_EQ(2, p[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	695 ASSERT_DOUBLE_EQ(3, p[2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	696 p = g.GetVoxelDimensions(VolumeProjection_Coronal);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	697 ASSERT_EQ(3u, p.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	698 ASSERT_DOUBLE_EQ(1, p[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	699 ASSERT_DOUBLE_EQ(3, p[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	700 ASSERT_DOUBLE_EQ(2, p[2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	701 p = g.GetVoxelDimensions(VolumeProjection_Sagittal);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	702 ASSERT_EQ(3u, p.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	703 ASSERT_DOUBLE_EQ(2, p[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	704 ASSERT_DOUBLE_EQ(3, p[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	705 ASSERT_DOUBLE_EQ(1, p[2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	706
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	707 // Loop over all the voxels of the volume
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	708 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	709 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	710 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	711
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	712 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	713 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	714 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	715
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	716 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	717 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	718 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	719
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	720 const float sx = 1.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	721 const float sy = 2.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	722 const float sz = 3.0f;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	723
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	724 Vector p = g.GetCoordinates(xx, yy, zz);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	725
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	726 Vector q = (g.GetAxialGeometry().MapSliceToWorldCoordinates(
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	727 static_cast<double>(x) * sx,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	728 static_cast<double>(y) * sy) +
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	729 z * sz * g.GetAxialGeometry().GetNormal());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	730 ASSERT_TRUE(IsEqualVectorL1(p, q));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	731
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	732 q = (g.GetCoronalGeometry().MapSliceToWorldCoordinates(
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	733 static_cast<double>(x) * sx,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	734 static_cast<double>(g.GetDepth() - 1 - z) * sz) +
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	735 y * sy * g.GetCoronalGeometry().GetNormal());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	736 ASSERT_TRUE(IsEqualVectorL1(p, q));
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 /**
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	739 * WARNING: In sagittal geometry, the normal points to
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	740 * 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	741 * the right-hand coordinate system. Hence the "-".
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	742 **/
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	743 q = (g.GetSagittalGeometry().MapSliceToWorldCoordinates(
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	744 static_cast<double>(y) * sy,
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	745 static_cast<double>(g.GetDepth() - 1 - z) * sz) +
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	746 x * sx * (-g.GetSagittalGeometry().GetNormal()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	747 ASSERT_TRUE(IsEqualVectorL1(p, q));
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	752 ASSERT_EQ(0, (int) VolumeProjection_Axial);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	753 ASSERT_EQ(1, (int) VolumeProjection_Coronal);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	754 ASSERT_EQ(2, (int) VolumeProjection_Sagittal);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	755
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	756 for (int p = 0; p < 3; p++)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	757 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	758 VolumeProjection projection = (VolumeProjection) p;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	759 const CoordinateSystem3D& s = g.GetProjectionGeometry(projection);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	760
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	761 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	762
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	763 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	764 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	765 CoordinateSystem3D plane = g.GetProjectionSlice(projection, i);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	766
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	767 if (projection == VolumeProjection_Sagittal)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	768 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	769 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	770 (-s.GetNormal()) * g.GetVoxelDimensions(projection)[2]));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	771 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	772 else
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	773 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	774 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	775 s.GetNormal() * g.GetVoxelDimensions(projection)[2]));
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	778 ASSERT_TRUE(IsEqualVectorL1(plane.GetAxisX(), s.GetAxisX()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	779 ASSERT_TRUE(IsEqualVectorL1(plane.GetAxisY(), s.GetAxisY()));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	780
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	781 unsigned int slice;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	782 VolumeProjection q;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	783 ASSERT_TRUE(g.DetectSlice(q, slice, plane));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	784 ASSERT_EQ(projection, q);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	785 ASSERT_EQ(i, slice);
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	791 TEST(LinearAlgebra, ParseVectorLocale)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	792 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	793 OrthancStone::Vector v;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	794
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	795 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "1.2"));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	796 ASSERT_EQ(1u, v.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	797 ASSERT_DOUBLE_EQ(1.2, v[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	798
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	799 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	800 ASSERT_EQ(1u, v.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	801 ASSERT_DOUBLE_EQ(-120.0, v[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	802
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	803 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	804 ASSERT_EQ(2u, v.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	805 ASSERT_DOUBLE_EQ(-0.01, v[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	806 ASSERT_DOUBLE_EQ(2.0, v[1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	807
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	808 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	809 ASSERT_EQ(2u, v.size());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	810 ASSERT_DOUBLE_EQ(1.3671875, v[0]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	811 ASSERT_DOUBLE_EQ(1.3671875, v[1]);
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	814 TEST(GenericToolbox, NormalizeUuid)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	815 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	816 std::string ref("44ca5051-14ef-4d2f-8bd7-db20bfb61fbb");
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 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	819 std::string u("44ca5051-14ef-4d2f-8bd7-db20bfb61fbb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	820 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	821 ASSERT_EQ(ref, u);
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	824 // space left
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	825 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	826 std::string u(" 44ca5051-14ef-4d2f-8bd7-db20bfb61fbb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	827 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	828 ASSERT_EQ(ref, u);
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	831 // space right
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	832 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	833 std::string u("44ca5051-14ef-4d2f-8bd7-db20bfb61fbb ");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	834 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	835 ASSERT_EQ(ref, u);
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	838 // space l & r
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	839 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	840 std::string u(" 44ca5051-14ef-4d2f-8bd7-db20bfb61fbb ");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	841 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	842 ASSERT_EQ(ref, u);
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	845 // space left + case
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	846 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	847 std::string u(" 44CA5051-14ef-4d2f-8bd7-dB20bfb61fbb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	848 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	849 ASSERT_EQ(ref, u);
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	852 // space right + case
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	853 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	854 std::string u("44ca5051-14EF-4D2f-8bd7-db20bfb61fbB ");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	855 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	856 ASSERT_EQ(ref, u);
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	859 // space l & r + case
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	860 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	861 std::string u(" 44cA5051-14Ef-4d2f-8bD7-db20bfb61fbb ");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	862 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	863 ASSERT_EQ(ref, u);
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	866 // no
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	867 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	868 std::string u(" 44ca5051-14ef-4d2f-8bd7- db20bfb61fbb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	869 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	870 ASSERT_NE(ref, u);
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	873 // no
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	874 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	875 std::string u("44ca5051-14ef-4d2f-8bd7-db20bfb61fb");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	876 OrthancStone::GenericToolbox::NormalizeUuid(u);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	877 ASSERT_NE(ref, u);
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
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	882 TEST(CoordinateSystem3D, Basic)
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	883 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	884 using namespace OrthancStone;
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 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	887 CoordinateSystem3D c;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	888 ASSERT_FALSE(c.IsValid());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	889 ASSERT_DOUBLE_EQ(c.GetNormal()[0], 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	890 ASSERT_DOUBLE_EQ(c.GetNormal()[1], 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	891 ASSERT_DOUBLE_EQ(c.GetNormal()[2], 1);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	892
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(0, 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(5, 0, 0)));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	895 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	896 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	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 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	900 CoordinateSystem3D c("nope1", "nope2");
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	901 ASSERT_FALSE(c.IsValid());
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	902 ASSERT_DOUBLE_EQ(c.GetNormal()[0], 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	903 ASSERT_DOUBLE_EQ(c.GetNormal()[1], 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	904 ASSERT_DOUBLE_EQ(c.GetNormal()[2], 1);
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 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	908 // 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	909 CoordinateSystem3D c = CoordinateSystem3D::CreateFromPlaneGeneralForm(2, 4, -4, -6);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	910 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	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 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	914 // https://mathinsight.org/distance_point_plane_examples
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	915 CoordinateSystem3D c = CoordinateSystem3D::CreateFromPlaneGeneralForm(2, -2, 5, 8);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	916 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	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 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	920 // 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	921 const Vector a = LinearAlgebra::CreateVector(3, 6, 9);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	922 const Vector b = LinearAlgebra::CreateVector(9, 6, 3);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	923 const Vector c = LinearAlgebra::CreateVector(6, -9, 9);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	924 CoordinateSystem3D d = CoordinateSystem3D::CreateFromThreePoints(a, b, c);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	925 ASSERT_DOUBLE_EQ(0, d.ComputeDistance(a));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	926 ASSERT_DOUBLE_EQ(0, d.ComputeDistance(b));
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	927 ASSERT_DOUBLE_EQ(0, d.ComputeDistance(c));
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 {
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	931 // 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	932 const Vector a = LinearAlgebra::CreateVector(1, -2, 0);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	933 const Vector b = LinearAlgebra::CreateVector(3, 1, 4);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	934 const Vector c = LinearAlgebra::CreateVector(0, -1, 2);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	935 CoordinateSystem3D d = CoordinateSystem3D::CreateFromThreePoints(a, b, c);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	936 double r = d.GetNormal() [0] / 2.0;
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	937 ASSERT_DOUBLE_EQ(-8 * r, d.GetNormal() [1]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	938 ASSERT_DOUBLE_EQ(5 * r, d.GetNormal() [2]);
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	939 }
5b8b88e5bfd6 successfully running unit tests in WebAssembly Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	940 }
1765 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
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	943 TEST(SlicesSorter, HFP)
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	944 {
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	945 // 2021-04-27-repro-bug-HFP-HFS-cartman
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 {
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	948 // 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	949 // 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	950 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	951
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	952 /**
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	953 * 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	954 * 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	955 * "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	956 * 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	957 * shifted by the depth of the volume.
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
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	960 const double spacingXY = 1.171875;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	961 const double width = 512;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	962 const double height = 512;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	963 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	964
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	965 OrthancStone::VolumeImageGeometry geometry;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	966 geometry.SetAxialGeometry(system);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	967 geometry.SetSizeInVoxels(width, height, depth);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	968 geometry.SetVoxelDimensions(spacingXY, spacingXY, 1 /* pixel spacing Z */);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	969
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	970 OrthancStone::Vector p, q;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	971
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	972 p = OrthancStone::LinearAlgebra::CreateVector(0.5 / double(width),
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	973 0.5 / double(height),
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	974 0.5 / double(depth), 1);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	975 q = OrthancStone::LinearAlgebra::Product(geometry.GetTransform(), p);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	976 ASSERT_FLOAT_EQ(300, q[0]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	977 ASSERT_FLOAT_EQ(302.5, q[1]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	978 ASSERT_FLOAT_EQ(323.11, q[2]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	979 ASSERT_FLOAT_EQ(1, q[3]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	980
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	981 p = OrthancStone::LinearAlgebra::CreateVector((width - 0.5) / double(width),
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	982 (height - 0.5) / double(height),
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	983 (depth - 0.5) / double(depth), 1);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	984 q = OrthancStone::LinearAlgebra::Product(geometry.GetTransform(), p);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	985
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	986 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	987 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	988 ASSERT_FLOAT_EQ(690.11, q[2]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	989 ASSERT_FLOAT_EQ(1, q[3]);
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 {
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	993 // 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	994 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	995 const double spacingXY = 3;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	996 const double width = 146; // Columns
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	997 const double height = 84; // Row
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	998 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	999
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1000 OrthancStone::VolumeImageGeometry geometry;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1001 geometry.SetAxialGeometry(system);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1002 geometry.SetSizeInVoxels(width, height, depth);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1003 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	1004
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1005 OrthancStone::Vector p, q;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1006
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1007 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	1008 q = OrthancStone::LinearAlgebra::Product(geometry.GetTransform(), p);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1009 ASSERT_FLOAT_EQ(-217.0492, q[0]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1010 ASSERT_FLOAT_EQ(-161.4141, q[1]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1011 ASSERT_FLOAT_EQ(376.61, q[2]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1012 ASSERT_FLOAT_EQ(1, q[3]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1013
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1014 p = OrthancStone::LinearAlgebra::CreateVector((width - 0.5) / double(width),
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1015 (height - 0.5) / double(height),
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1016 (depth - 0.5) / double(depth), 1);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1017 q = OrthancStone::LinearAlgebra::Product(geometry.GetTransform(), p);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1018
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1019 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	1020 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	1021 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	1022 ASSERT_FLOAT_EQ(1, q[3]);
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
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1025 for (unsigned int upward = 0; upward < 2; upward++)
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1026 {
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1027 OrthancStone::SlicesSorter slices;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1028
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1029 for (unsigned int i = 0; i < 368; i++)
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1030 {
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1031 unsigned int z = (upward ? 323 + i : 690 - i);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1032 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	1033 "-1\\0\\0\\0\\-1\\0");
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1034 slices.AddSlice(p, new Orthanc::SingleValueObject<unsigned int>(z));
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
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1037 slices.Sort();
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1038
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1039 double spacingZ;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1040 ASSERT_TRUE(slices.ComputeSpacingBetweenSlices(spacingZ));
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1041 ASSERT_FLOAT_EQ(1, spacingZ);
1766 ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1042 ASSERT_TRUE(slices.AreAllSlicesDistinct());
1765 177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1043
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1044 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	1045 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	1046 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	1047
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1048 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	1049 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	1050 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	1051
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1052 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	1053 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	1054
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1055 const unsigned int width = 512;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1056 const unsigned int height = 512;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1057 const unsigned int depth = 368;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1058 const double spacingXY = 1.171875;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1059
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1060 OrthancStone::VolumeImageGeometry geometry;
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1061 geometry.SetSizeInVoxels(width, height, depth);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1062 geometry.SetVoxelDimensions(spacingXY, spacingXY, 1);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1063 geometry.SetAxialGeometry(slices.GetSliceGeometry(0));
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1064
1766 ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1065 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	1066 ASSERT_EQ(3u, q.size());
1765 177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1067 ASSERT_FLOAT_EQ(300 + spacingXY / 2.0, q[0]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1068 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	1069 ASSERT_FLOAT_EQ(323.11 - 0.5, q[2]);
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1070
1766 ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1071 q = geometry.GetCoordinates(1, 1, 1);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1072 ASSERT_EQ(3u, q.size());
1765 177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1073 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	1074 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	1075 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	1076
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1077 OrthancStone::VolumeProjection projection;
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1078 unsigned int slice;
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1079 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	1080 ASSERT_EQ(OrthancStone::VolumeProjection_Axial, projection);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1081 ASSERT_EQ(367u, slice);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1082 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	1083 ASSERT_EQ(OrthancStone::VolumeProjection_Axial, projection);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1084 ASSERT_EQ(0u, slice);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1085
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1086 // 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	1087 ASSERT_TRUE(geometry.DetectSlice(projection, slice,
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1088 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	1089 ASSERT_EQ(OrthancStone::VolumeProjection_Axial, projection);
ae18b5cfbd3e VolumeImageGeometry::DetectProjection() reports if normal is opposite Sebastien Jodogne <s.jodogne@gmail.com> parents: 1765 diff changeset	1090 ASSERT_EQ(376u - 323u, slice);
1765 177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1091 }
177f93c8371c unit test SlicesSorter.HFP Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	1092 }

Mercurial > hg > orthanc-stone

annotate UnitTestsSources/GeometryToolboxTests.cpp @ 1871:7053b8a0aaec