orthanc-stone: UnitTestsSources/UnitTestsMain.cpp annotate

annotate UnitTestsSources/UnitTestsMain.cpp @ 1603:595c0952ef7e

focusing on osirix annotations in Stone Web viewer

author	Sebastien Jodogne <s.jodogne@gmail.com>
date	Wed, 28 Oct 2020 18:49:01 +0100
parents	4fb8fdf03314
children	78509230f0d7

rev	line source
20 946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1 /**
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	2 * Stone of Orthanc
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	4 * Department, University Hospital of Liege, Belgium
1270 2d8ab34c8c91 upgrade to year 2020 Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	5 * Copyright (C) 2017-2020 Osimis S.A., Belgium
20 946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	6 *
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	7 * This program is free software: you can redistribute it and/or
47 28956ed68280 agpl license Sebastien Jodogne <s.jodogne@gmail.com> parents: 40 diff changeset	8 * modify it under the terms of the GNU Affero General Public License
28956ed68280 agpl license Sebastien Jodogne <s.jodogne@gmail.com> parents: 40 diff changeset	9 * as published by the Free Software Foundation, either version 3 of
28956ed68280 agpl license Sebastien Jodogne <s.jodogne@gmail.com> parents: 40 diff changeset	10 * the License, or (at your option) any later version.
20 946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	11 *
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	12 * This program is distributed in the hope that it will be useful, but
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	13 * WITHOUT ANY WARRANTY; without even the implied warranty of
47 28956ed68280 agpl license Sebastien Jodogne <s.jodogne@gmail.com> parents: 40 diff changeset	14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
28956ed68280 agpl license Sebastien Jodogne <s.jodogne@gmail.com> parents: 40 diff changeset	15 * Affero General Public License for more details.
1596 4fb8fdf03314 removed annoying whitespace Sebastien Jodogne <s.jodogne@gmail.com> parents: 1591 diff changeset	16 *
47 28956ed68280 agpl license Sebastien Jodogne <s.jodogne@gmail.com> parents: 40 diff changeset	17 * You should have received a copy of the GNU Affero General Public License
20 946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	18 * along with this program. If not, see <http://www.gnu.org/licenses/>.
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	19 **/
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	20
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	21
1455 30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1413 diff changeset	22 #include <gtest/gtest.h>
20 946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	23
1587 a1405ab3a91c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: 1541 diff changeset	24 #include "../OrthancStone/Sources/StoneInitialization.h"
a1405ab3a91c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: 1541 diff changeset	25 #include "../OrthancStone/Sources/Toolbox/FiniteProjectiveCamera.h"
1591 5887a4f8594b moving platform-specific files out of the "OrthancStone" folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 1587 diff changeset	26 #include "../OrthancStone/Sources/Toolbox/GenericToolbox.h"
1587 a1405ab3a91c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: 1541 diff changeset	27 #include "../OrthancStone/Sources/Toolbox/GeometryToolbox.h"
a1405ab3a91c reorganization Sebastien Jodogne <s.jodogne@gmail.com> parents: 1541 diff changeset	28 #include "../OrthancStone/Sources/Volumes/ImageBuffer3D.h"
113 2eca030792aa using the Orthanc Framework Sebastien Jodogne <s.jodogne@gmail.com> parents: 102 diff changeset	29
1455 30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1413 diff changeset	30 #include <Images/ImageProcessing.h>
30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1413 diff changeset	31 #include <Logging.h>
30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1413 diff changeset	32 #include <MultiThreading/SharedMessageQueue.h>
30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1413 diff changeset	33 #include <OrthancException.h>
87 4a541cd4fa83 OrthancVolumeImageLoader Sebastien Jodogne <s.jodogne@gmail.com> parents: 82 diff changeset	34
68 1526d38ef6da SliceSorter Sebastien Jodogne <s.jodogne@gmail.com> parents: 67 diff changeset	35 #include <boost/lexical_cast.hpp>
70 f73aed014bde OrthancAsynchronousWebService Sebastien Jodogne <s.jodogne@gmail.com> parents: 69 diff changeset	36 #include <boost/date_time/posix_time/posix_time.hpp>
f73aed014bde OrthancAsynchronousWebService Sebastien Jodogne <s.jodogne@gmail.com> parents: 69 diff changeset	37 #include <boost/thread/thread.hpp>
102 fcec0ab44054 display volumes Sebastien Jodogne <s.jodogne@gmail.com> parents: 93 diff changeset	38 #include <boost/math/special_functions/round.hpp>
fcec0ab44054 display volumes Sebastien Jodogne <s.jodogne@gmail.com> parents: 93 diff changeset	39
68 1526d38ef6da SliceSorter Sebastien Jodogne <s.jodogne@gmail.com> parents: 67 diff changeset	40
144 9b83f30fc1c0 bilinear and trilinear interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 135 diff changeset	41 TEST(GeometryToolbox, Interpolation)
9b83f30fc1c0 bilinear and trilinear interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 135 diff changeset	42 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	43 using namespace OrthancStone::GeometryToolbox;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	44
144 9b83f30fc1c0 bilinear and trilinear interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 135 diff changeset	45 // https://en.wikipedia.org/wiki/Bilinear_interpolation#Application_in_image_processing
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	46 ASSERT_FLOAT_EQ(146.1f, ComputeBilinearInterpolationUnitSquare(0.5f, 0.2f, 91, 210, 162, 95));
176 ab9c799f5de1 fix interpolation with negative coordinates Sebastien Jodogne <s.jodogne@gmail.com> parents: 165 diff changeset	47
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	48 ASSERT_FLOAT_EQ(91, ComputeBilinearInterpolationUnitSquare(0, 0, 91, 210, 162, 95));
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	49 ASSERT_FLOAT_EQ(210, ComputeBilinearInterpolationUnitSquare(1, 0, 91, 210, 162, 95));
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	50 ASSERT_FLOAT_EQ(162, ComputeBilinearInterpolationUnitSquare(0, 1, 91, 210, 162, 95));
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	51 ASSERT_FLOAT_EQ(95, ComputeBilinearInterpolationUnitSquare(1, 1, 91, 210, 162, 95));
176 ab9c799f5de1 fix interpolation with negative coordinates Sebastien Jodogne <s.jodogne@gmail.com> parents: 165 diff changeset	52
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	53 ASSERT_FLOAT_EQ(123.35f, ComputeTrilinearInterpolationUnitSquare
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	54 (0.5f, 0.2f, 0.7f,
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	55 91, 210, 162, 95,
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	56 51, 190, 80, 92));
176 ab9c799f5de1 fix interpolation with negative coordinates Sebastien Jodogne <s.jodogne@gmail.com> parents: 165 diff changeset	57
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	58 ASSERT_FLOAT_EQ(ComputeBilinearInterpolationUnitSquare(0.5f, 0.2f, 91, 210, 162, 95),
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	59 ComputeTrilinearInterpolationUnitSquare(0.5f, 0.2f, 0,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	60 91, 210, 162, 95,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	61 51, 190, 80, 92));
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	62
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	63 ASSERT_FLOAT_EQ(ComputeBilinearInterpolationUnitSquare(0.5f, 0.2f, 51, 190, 80, 92),
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	64 ComputeTrilinearInterpolationUnitSquare(0.5f, 0.2f, 1,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	65 91, 210, 162, 95,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	66 51, 190, 80, 92));
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	67 }
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	68
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	69
162 77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	70 static bool CompareMatrix(const OrthancStone::Matrix& a,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	71 const OrthancStone::Matrix& b,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	72 double threshold = 0.00000001)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	73 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	74 if (a.size1() != b.size1() \|\|
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	75 a.size2() != b.size2())
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	76 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	77 return false;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	78 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	79
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	80 for (size_t i = 0; i < a.size1(); i++)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	81 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	82 for (size_t j = 0; j < a.size2(); j++)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	83 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	84 if (fabs(a(i, j) - b(i, j)) > threshold)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	85 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	86 LOG(ERROR) << "Too large difference in component ("
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	87 << i << "," << j << "): " << a(i,j) << " != " << b(i,j);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	88 return false;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	89 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	90 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	91 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	92
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	93 return true;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	94 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	95
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	96
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	97 static bool CompareVector(const OrthancStone::Vector& a,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	98 const OrthancStone::Vector& b,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	99 double threshold = 0.00000001)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	100 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	101 if (a.size() != b.size())
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	102 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	103 return false;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	104 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	105
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	106 for (size_t i = 0; i < a.size(); i++)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	107 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	108 if (fabs(a(i) - b(i)) > threshold)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	109 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	110 LOG(ERROR) << "Too large difference in component "
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	111 << i << ": " << a(i) << " != " << b(i);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	112 return false;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	113 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	114 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	115
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	116 return true;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	117 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	118
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	119
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	120
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	121 TEST(FiniteProjectiveCamera, Decomposition1)
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	122 {
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	123 // Example 6.2 of "Multiple View Geometry in Computer Vision - 2nd
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	124 // edition" (page 163)
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	125 const double p[12] = {
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	126 3.53553e+2, 3.39645e+2, 2.77744e+2, -1.44946e+6,
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	127 -1.03528e+2, 2.33212e+1, 4.59607e+2, -6.32525e+5,
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	128 7.07107e-1, -3.53553e-1, 6.12372e-1, -9.18559e+2
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	129 };
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	130
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	131 OrthancStone::FiniteProjectiveCamera camera(p);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	132 ASSERT_EQ(3u, camera.GetMatrix().size1());
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	133 ASSERT_EQ(4u, camera.GetMatrix().size2());
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	134 ASSERT_EQ(3u, camera.GetIntrinsicParameters().size1());
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	135 ASSERT_EQ(3u, camera.GetIntrinsicParameters().size2());
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	136 ASSERT_EQ(3u, camera.GetRotation().size1());
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	137 ASSERT_EQ(3u, camera.GetRotation().size2());
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	138 ASSERT_EQ(3u, camera.GetCenter().size());
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	139
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	140 ASSERT_NEAR(1000.0, camera.GetCenter()[0], 0.01);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	141 ASSERT_NEAR(2000.0, camera.GetCenter()[1], 0.01);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	142 ASSERT_NEAR(1500.0, camera.GetCenter()[2], 0.01);
144 9b83f30fc1c0 bilinear and trilinear interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 135 diff changeset	143
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	144 ASSERT_NEAR(468.2, camera.GetIntrinsicParameters() (0, 0), 0.1);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	145 ASSERT_NEAR(91.2, camera.GetIntrinsicParameters() (0, 1), 0.1);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	146 ASSERT_NEAR(300.0, camera.GetIntrinsicParameters() (0, 2), 0.1);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	147 ASSERT_NEAR(427.2, camera.GetIntrinsicParameters() (1, 1), 0.1);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	148 ASSERT_NEAR(200.0, camera.GetIntrinsicParameters() (1, 2), 0.1);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	149 ASSERT_NEAR(1.0, camera.GetIntrinsicParameters() (2, 2), 0.1);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	150
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	151 ASSERT_NEAR(0, camera.GetIntrinsicParameters() (1, 0), 0.0000001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	152 ASSERT_NEAR(0, camera.GetIntrinsicParameters() (2, 0), 0.0000001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	153 ASSERT_NEAR(0, camera.GetIntrinsicParameters() (2, 1), 0.0000001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	154
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	155 ASSERT_NEAR(0.41380, camera.GetRotation() (0, 0), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	156 ASSERT_NEAR(0.90915, camera.GetRotation() (0, 1), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	157 ASSERT_NEAR(0.04708, camera.GetRotation() (0, 2), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	158 ASSERT_NEAR(-0.57338, camera.GetRotation() (1, 0), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	159 ASSERT_NEAR(0.22011, camera.GetRotation() (1, 1), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	160 ASSERT_NEAR(0.78917, camera.GetRotation() (1, 2), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	161 ASSERT_NEAR(0.70711, camera.GetRotation() (2, 0), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	162 ASSERT_NEAR(-0.35355, camera.GetRotation() (2, 1), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	163 ASSERT_NEAR(0.61237, camera.GetRotation() (2, 2), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	164
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	165 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	166
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	167 OrthancStone::FiniteProjectiveCamera camera2(camera.GetIntrinsicParameters(),
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	168 camera.GetRotation(),
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	169 camera.GetCenter());
162 77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	170
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	171 ASSERT_TRUE(CompareMatrix(camera.GetMatrix(), camera2.GetMatrix()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	172 ASSERT_TRUE(CompareMatrix(camera.GetIntrinsicParameters(), camera2.GetIntrinsicParameters()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	173 ASSERT_TRUE(CompareMatrix(camera.GetRotation(), camera2.GetRotation()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	174 ASSERT_TRUE(CompareVector(camera.GetCenter(), camera2.GetCenter()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	175 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	176
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	177
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	178 TEST(FiniteProjectiveCamera, Decomposition2)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	179 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	180 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 };
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	181 const double k[] = { -1528.494743, 0.000000, 256.000000, 0.000000, 1528.494743, 256.000000, 0.000000, 0.000000, 1.000000 };
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	182 const double r[] = { -0.858893, -0.330733, 0.391047, -0.158171, 0.897503, 0.411668, -0.487118, 0.291726, -0.823172 };
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	183 const double c[] = { 243.558936, -145.863085, 411.585964 };
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	184
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	185 OrthancStone::FiniteProjectiveCamera camera(p);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	186 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	187
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	188 OrthancStone::FiniteProjectiveCamera camera2(k, r, c);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	189 ASSERT_TRUE(CompareMatrix(camera.GetMatrix(), camera2.GetMatrix(), 1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	190 ASSERT_TRUE(CompareMatrix(camera.GetIntrinsicParameters(), camera2.GetIntrinsicParameters(), 0.001));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	191 ASSERT_TRUE(CompareMatrix(camera.GetRotation(), camera2.GetRotation(), 0.000001));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	192 ASSERT_TRUE(CompareVector(camera.GetCenter(), camera2.GetCenter(), 0.0001));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	193 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	194
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	195
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	196 TEST(FiniteProjectiveCamera, Decomposition3)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	197 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	198 const double p[] = { 10, 0, 0, 0,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	199 0, 20, 0, 0,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	200 0, 0, 30, 0 };
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	201
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	202 OrthancStone::FiniteProjectiveCamera camera(p);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	203 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	204 ASSERT_DOUBLE_EQ(10, camera.GetIntrinsicParameters() (0, 0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	205 ASSERT_DOUBLE_EQ(20, camera.GetIntrinsicParameters() (1, 1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	206 ASSERT_DOUBLE_EQ(30, camera.GetIntrinsicParameters() (2, 2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	207 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (0, 0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	208 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (1, 1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	209 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (2, 2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	210 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	211 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	212 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	213 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	214
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	215
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	216 TEST(FiniteProjectiveCamera, Decomposition4)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	217 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	218 const double p[] = { 1, 0, 0, 10,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	219 0, 1, 0, 20,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	220 0, 0, 1, 30 };
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	221
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	222 OrthancStone::FiniteProjectiveCamera camera(p);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	223 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	224 ASSERT_DOUBLE_EQ(1, camera.GetIntrinsicParameters() (0, 0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	225 ASSERT_DOUBLE_EQ(1, camera.GetIntrinsicParameters() (1, 1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	226 ASSERT_DOUBLE_EQ(1, camera.GetIntrinsicParameters() (2, 2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	227 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (0, 0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	228 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (1, 1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	229 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (2, 2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	230 ASSERT_DOUBLE_EQ(-10, camera.GetCenter() (0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	231 ASSERT_DOUBLE_EQ(-20, camera.GetCenter() (1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	232 ASSERT_DOUBLE_EQ(-30, camera.GetCenter() (2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	233 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	234
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	235
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	236 TEST(FiniteProjectiveCamera, Decomposition5)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	237 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	238 const double p[] = { 0, 0, 10, 0,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	239 0, 20, 0, 0,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	240 30, 0, 0, 0 };
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	241
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	242 OrthancStone::FiniteProjectiveCamera camera(p);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	243 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	244 ASSERT_DOUBLE_EQ(-10, camera.GetIntrinsicParameters() (0, 0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	245 ASSERT_DOUBLE_EQ(20, camera.GetIntrinsicParameters() (1, 1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	246 ASSERT_DOUBLE_EQ(30, camera.GetIntrinsicParameters() (2, 2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	247 ASSERT_DOUBLE_EQ(-1, camera.GetRotation() (0, 2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	248 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (1, 1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	249 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (2, 0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	250 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	251 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	252 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	253
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	254 OrthancStone::FiniteProjectiveCamera camera2(camera.GetIntrinsicParameters(),
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	255 camera.GetRotation(),
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	256 camera.GetCenter());
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	257 ASSERT_TRUE(CompareMatrix(camera.GetMatrix(), camera2.GetMatrix()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	258 ASSERT_TRUE(CompareMatrix(camera.GetIntrinsicParameters(), camera2.GetIntrinsicParameters()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	259 ASSERT_TRUE(CompareMatrix(camera.GetRotation(), camera2.GetRotation()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	260 ASSERT_TRUE(CompareVector(camera.GetCenter(), camera2.GetCenter()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	261 }
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	262
162 77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	263
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	264 static double GetCosAngle(const OrthancStone::Vector& a,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	265 const OrthancStone::Vector& b)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	266 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	267 // Returns the cosine of the angle between two vectors
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	268 // https://en.wikipedia.org/wiki/Dot_product#Geometric_definition
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	269 return boost::numeric::ublas::inner_prod(a, b) /
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	270 (boost::numeric::ublas::norm_2(a) * boost::numeric::ublas::norm_2(b));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	271 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	272
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	273
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	274 TEST(FiniteProjectiveCamera, Ray)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	275 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	276 const double pp[] = { -1499.650894, 2954.618773, -259.737419, 637891.819097,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	277 -2951.517707, -1501.019129, -285.785281, 637891.819097,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	278 0.008528, 0.003067, -0.999959, 2491.764918 };
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	279
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	280 const OrthancStone::FiniteProjectiveCamera camera(pp);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	281
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	282 ASSERT_NEAR(-21.2492, camera.GetCenter() (0), 0.0001);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	283 ASSERT_NEAR(-7.64234, camera.GetCenter() (1), 0.00001);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	284 ASSERT_NEAR(2491.66, camera.GetCenter() (2), 0.01);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	285
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	286 // Image plane that led to these parameters, with principal point at
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	287 // (256,256). The image has dimensions 512x512.
165 8d50e6be565d LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 164 diff changeset	288 OrthancStone::Vector o =
8d50e6be565d LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 164 diff changeset	289 OrthancStone::LinearAlgebra::CreateVector(7.009620, 2.521030, -821.942000);
8d50e6be565d LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 164 diff changeset	290 OrthancStone::Vector ax =
8d50e6be565d LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 164 diff changeset	291 OrthancStone::LinearAlgebra::CreateVector(-0.453219, 0.891399, -0.001131);
8d50e6be565d LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 164 diff changeset	292 OrthancStone::Vector ay =
8d50e6be565d LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 164 diff changeset	293 OrthancStone::LinearAlgebra::CreateVector(-0.891359, -0.453210, -0.008992);
8d50e6be565d LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 164 diff changeset	294
162 77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	295 OrthancStone::CoordinateSystem3D imagePlane(o, ax, ay);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	296
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	297 // Back-projection of the principal point
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	298 {
162 77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	299 OrthancStone::Vector ray = camera.GetRayDirection(256, 256);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	300
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	301 // The principal axis vector is orthogonal to the image plane
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	302 // (i.e. parallel to the plane normal), in the opposite direction
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	303 // ("-1" corresponds to "cos(pi)").
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	304 ASSERT_NEAR(-1, GetCosAngle(ray, imagePlane.GetNormal()), 0.0000001);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	305
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	306 // Forward projection of principal axis, resulting in the principal point
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	307 double x, y;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	308 camera.ApplyFinite(x, y, camera.GetCenter() - ray);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	309
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	310 ASSERT_NEAR(256, x, 0.00001);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	311 ASSERT_NEAR(256, y, 0.00001);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	312 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	313
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	314 // Back-projection of the 4 corners of the image
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	315 std::vector<double> cx, cy;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	316 cx.push_back(0);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	317 cy.push_back(0);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	318 cx.push_back(512);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	319 cy.push_back(0);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	320 cx.push_back(512);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	321 cy.push_back(512);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	322 cx.push_back(0);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	323 cy.push_back(512);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	324
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	325 bool first = true;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	326 double angle;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	327
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	328 for (size_t i = 0; i < cx.size(); i++)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	329 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	330 OrthancStone::Vector ray = camera.GetRayDirection(cx[i], cy[i]);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	331
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	332 // Check that the angle wrt. principal axis is the same for all
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	333 // the 4 corners
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	334 double a = GetCosAngle(ray, imagePlane.GetNormal());
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	335 if (first)
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	336 {
162 77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	337 first = false;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	338 angle = a;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	339 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	340 else
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	341 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	342 ASSERT_NEAR(angle, a, 0.000001);
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	343 }
162 77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	344
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	345 // Forward projection of the ray, going back to the original point
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	346 double x, y;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	347 camera.ApplyFinite(x, y, camera.GetCenter() - ray);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	348
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	349 ASSERT_NEAR(cx[i], x, 0.00001);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	350 ASSERT_NEAR(cy[i], y, 0.00001);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	351
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	352 // Alternative construction, by computing the intersection of the
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	353 // ray with the image plane
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	354 OrthancStone::Vector p;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	355 ASSERT_TRUE(imagePlane.IntersectLine(p, camera.GetCenter(), -ray));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	356 imagePlane.ProjectPoint(x, y, p);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	357 ASSERT_NEAR(cx[i], x + 256, 0.01);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	358 ASSERT_NEAR(cy[i], y + 256, 0.01);
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	359 }
144 9b83f30fc1c0 bilinear and trilinear interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 135 diff changeset	360 }
87 4a541cd4fa83 OrthancVolumeImageLoader Sebastien Jodogne <s.jodogne@gmail.com> parents: 82 diff changeset	361
20 946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	362
163 8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	363 TEST(Matrix, Inverse1)
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	364 {
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	365 OrthancStone::Matrix a, b;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	366
164 432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	367 a.resize(0, 0);
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	368 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	369 ASSERT_EQ(0u, b.size1());
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	370 ASSERT_EQ(0u, b.size2());
163 8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	371
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	372 a.resize(2, 3);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	373 ASSERT_THROW(OrthancStone::LinearAlgebra::InvertMatrix(b, a), Orthanc::OrthancException);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	374
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	375 a.resize(1, 1);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	376 a(0, 0) = 45.0;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	377
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	378 ASSERT_DOUBLE_EQ(45, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	379 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	380 ASSERT_EQ(1u, b.size1());
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	381 ASSERT_EQ(1u, b.size2());
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	382 ASSERT_DOUBLE_EQ(1.0 / 45.0, b(0, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	383
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	384 a(0, 0) = 0;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	385 ASSERT_DOUBLE_EQ(0, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	386 ASSERT_THROW(OrthancStone::LinearAlgebra::InvertMatrix(b, a), Orthanc::OrthancException);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	387 }
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	388
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	389
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	390 TEST(Matrix, Inverse2)
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	391 {
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	392 OrthancStone::Matrix a, b;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	393 a.resize(2, 2);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	394 a(0, 0) = 4;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	395 a(0, 1) = 3;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	396 a(1, 0) = 3;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	397 a(1, 1) = 2;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	398
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	399 ASSERT_DOUBLE_EQ(-1, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	400 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	401 ASSERT_EQ(2u, b.size1());
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	402 ASSERT_EQ(2u, b.size2());
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	403
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	404 ASSERT_DOUBLE_EQ(-2, b(0, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	405 ASSERT_DOUBLE_EQ(3, b(0, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	406 ASSERT_DOUBLE_EQ(3, b(1, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	407 ASSERT_DOUBLE_EQ(-4, b(1, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	408
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	409 a(0, 0) = 1;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	410 a(0, 1) = 2;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	411 a(1, 0) = 3;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	412 a(1, 1) = 4;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	413
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	414 ASSERT_DOUBLE_EQ(-2, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	415 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	416
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	417 ASSERT_DOUBLE_EQ(-2, b(0, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	418 ASSERT_DOUBLE_EQ(1, b(0, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	419 ASSERT_DOUBLE_EQ(1.5, b(1, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	420 ASSERT_DOUBLE_EQ(-0.5, b(1, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	421 }
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	422
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	423
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	424 TEST(Matrix, Inverse3)
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	425 {
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	426 OrthancStone::Matrix a, b;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	427 a.resize(3, 3);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	428 a(0, 0) = 7;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	429 a(0, 1) = 2;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	430 a(0, 2) = 1;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	431 a(1, 0) = 0;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	432 a(1, 1) = 3;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	433 a(1, 2) = -1;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	434 a(2, 0) = -3;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	435 a(2, 1) = 4;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	436 a(2, 2) = -2;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	437
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	438 ASSERT_DOUBLE_EQ(1, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	439 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	440 ASSERT_EQ(3u, b.size1());
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	441 ASSERT_EQ(3u, b.size2());
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	442
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	443 ASSERT_DOUBLE_EQ(-2, b(0, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	444 ASSERT_DOUBLE_EQ(8, b(0, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	445 ASSERT_DOUBLE_EQ(-5, b(0, 2));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	446 ASSERT_DOUBLE_EQ(3, b(1, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	447 ASSERT_DOUBLE_EQ(-11, b(1, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	448 ASSERT_DOUBLE_EQ(7, b(1, 2));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	449 ASSERT_DOUBLE_EQ(9, b(2, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	450 ASSERT_DOUBLE_EQ(-34, b(2, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	451 ASSERT_DOUBLE_EQ(21, b(2, 2));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	452
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	453
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	454 a(0, 0) = 1;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	455 a(0, 1) = 2;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	456 a(0, 2) = 2;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	457 a(1, 0) = 1;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	458 a(1, 1) = 0;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	459 a(1, 2) = 1;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	460 a(2, 0) = 1;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	461 a(2, 1) = 2;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	462 a(2, 2) = 1;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	463
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	464 ASSERT_DOUBLE_EQ(2, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	465 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	466 ASSERT_EQ(3u, b.size1());
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	467 ASSERT_EQ(3u, b.size2());
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	468
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	469 ASSERT_DOUBLE_EQ(-1, b(0, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	470 ASSERT_DOUBLE_EQ(1, b(0, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	471 ASSERT_DOUBLE_EQ(1, b(0, 2));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	472 ASSERT_DOUBLE_EQ(0, b(1, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	473 ASSERT_DOUBLE_EQ(-0.5, b(1, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	474 ASSERT_DOUBLE_EQ(0.5, b(1, 2));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	475 ASSERT_DOUBLE_EQ(1, b(2, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	476 ASSERT_DOUBLE_EQ(0, b(2, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	477 ASSERT_DOUBLE_EQ(-1, b(2, 2));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	478 }
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	479
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	480
164 432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	481 TEST(Matrix, Inverse4)
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	482 {
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	483 OrthancStone::Matrix a, b;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	484 a.resize(4, 4);
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	485 a(0, 0) = 2;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	486 a(0, 1) = 1;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	487 a(0, 2) = 2;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	488 a(0, 3) = -3;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	489 a(1, 0) = -2;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	490 a(1, 1) = 2;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	491 a(1, 2) = -1;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	492 a(1, 3) = -1;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	493 a(2, 0) = 2;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	494 a(2, 1) = 2;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	495 a(2, 2) = -3;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	496 a(2, 3) = -1;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	497 a(3, 0) = 3;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	498 a(3, 1) = -2;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	499 a(3, 2) = -3;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	500 a(3, 3) = -1;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	501
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	502 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	503 ASSERT_EQ(4u, b.size1());
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	504 ASSERT_EQ(4u, b.size2());
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	505
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	506 b *= 134.0; // This is the determinant
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	507
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	508 ASSERT_DOUBLE_EQ(8, b(0, 0));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	509 ASSERT_DOUBLE_EQ(-44, b(0, 1));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	510 ASSERT_DOUBLE_EQ(30, b(0, 2));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	511 ASSERT_DOUBLE_EQ(-10, b(0, 3));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	512 ASSERT_DOUBLE_EQ(2, b(1, 0));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	513 ASSERT_DOUBLE_EQ(-11, b(1, 1));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	514 ASSERT_DOUBLE_EQ(41, b(1, 2));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	515 ASSERT_DOUBLE_EQ(-36, b(1, 3));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	516 ASSERT_DOUBLE_EQ(16, b(2, 0));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	517 ASSERT_DOUBLE_EQ(-21, b(2, 1));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	518 ASSERT_DOUBLE_EQ(-7, b(2, 2));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	519 ASSERT_DOUBLE_EQ(-20, b(2, 3));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	520 ASSERT_DOUBLE_EQ(-28, b(3, 0));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	521 ASSERT_DOUBLE_EQ(-47, b(3, 1));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	522 ASSERT_DOUBLE_EQ(29, b(3, 2));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	523 ASSERT_DOUBLE_EQ(-32, b(3, 3));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	524 }
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	525
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	526
188 45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	527 TEST(FiniteProjectiveCamera, Calibration)
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	528 {
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	529 unsigned int volumeWidth = 512;
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	530 unsigned int volumeHeight = 512;
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	531 unsigned int volumeDepth = 110;
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	532
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	533 OrthancStone::Vector camera = OrthancStone::LinearAlgebra::CreateVector
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	534 (-1000, -5000, -static_cast<double>(volumeDepth) * 32);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	535
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	536 OrthancStone::Vector principalPoint = OrthancStone::LinearAlgebra::CreateVector
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	537 (volumeWidth/2, volumeHeight/2, volumeDepth * 2);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	538
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	539 OrthancStone::FiniteProjectiveCamera c(camera, principalPoint, 0, 512, 512, 1, 1);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	540
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	541 double swapv[9] = { 1, 0, 0,
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	542 0, -1, 512,
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	543 0, 0, 1 };
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	544 OrthancStone::Matrix swap;
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	545 OrthancStone::LinearAlgebra::FillMatrix(swap, 3, 3, swapv);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	546
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	547 OrthancStone::Matrix p = OrthancStone::LinearAlgebra::Product(swap, c.GetMatrix());
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	548 p /= p(2,3);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	549
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	550 ASSERT_NEAR( 1.04437, p(0,0), 0.00001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	551 ASSERT_NEAR(-0.0703111, p(0,1), 0.00000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	552 ASSERT_NEAR(-0.179283, p(0,2), 0.000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	553 ASSERT_NEAR( 61.7431, p(0,3), 0.0001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	554 ASSERT_NEAR( 0.11127, p(1,0), 0.000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	555 ASSERT_NEAR(-0.595541, p(1,1), 0.000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	556 ASSERT_NEAR( 0.872211, p(1,2), 0.000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	557 ASSERT_NEAR( 203.748, p(1,3), 0.001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	558 ASSERT_NEAR( 3.08593e-05, p(2,0), 0.0000000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	559 ASSERT_NEAR( 0.000129138, p(2,1), 0.000000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	560 ASSERT_NEAR( 9.18901e-05, p(2,2), 0.0000000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	561 ASSERT_NEAR( 1, p(2,3), 0.0000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	562 }
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	563
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	564
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	565 static bool IsEqualRotationVector(OrthancStone::Vector a,
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	566 OrthancStone::Vector b)
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	567 {
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	568 if (a.size() != b.size() \|\|
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	569 a.size() != 3)
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	570 {
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	571 return false;
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	572 }
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	573 else
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	574 {
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	575 OrthancStone::LinearAlgebra::NormalizeVector(a);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	576 OrthancStone::LinearAlgebra::NormalizeVector(b);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	577 return OrthancStone::LinearAlgebra::IsCloseToZero(boost::numeric::ublas::norm_2(a - b));
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	578 }
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	579 }
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	580
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	581
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	582 TEST(GeometryToolbox, AlignVectorsWithRotation)
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	583 {
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	584 OrthancStone::Vector a, b;
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	585 OrthancStone::Matrix r;
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	586
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	587 OrthancStone::LinearAlgebra::AssignVector(a, -200, 200, -846.63);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	588 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	589
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	590 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	591 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	592 ASSERT_TRUE(IsEqualRotationVector(OrthancStone::LinearAlgebra::Product(r, a), b));
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	593
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	594 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, b, a);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	595 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	596 ASSERT_TRUE(IsEqualRotationVector(OrthancStone::LinearAlgebra::Product(r, b), a));
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	597
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	598 OrthancStone::LinearAlgebra::AssignVector(a, 1, 0, 0);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	599 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	600 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	601 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	602 ASSERT_TRUE(IsEqualRotationVector(OrthancStone::LinearAlgebra::Product(r, a), b));
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	603
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	604 OrthancStone::LinearAlgebra::AssignVector(a, 0, 1, 0);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	605 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	606 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	607 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	608 ASSERT_TRUE(IsEqualRotationVector(OrthancStone::LinearAlgebra::Product(r, a), b));
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	609
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	610 OrthancStone::LinearAlgebra::AssignVector(a, 0, 0, 1);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	611 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	612 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	613 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	614 ASSERT_TRUE(IsEqualRotationVector(OrthancStone::LinearAlgebra::Product(r, a), b));
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	615
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	616 OrthancStone::LinearAlgebra::AssignVector(a, 0, 0, 0);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	617 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	618 ASSERT_THROW(OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b), Orthanc::OrthancException);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	619
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	620 // TODO: Deal with opposite vectors
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	621
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	622 /*
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	623 OrthancStone::LinearAlgebra::AssignVector(a, 0, 0, -1);
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	624 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	625 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	626 OrthancStone::LinearAlgebra::Print(r);
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	627 OrthancStone::LinearAlgebra::Print(boost::numeric::ublas::prod(r, a));
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	628 */
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	629 }
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	630
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	631
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	632 static bool IsEqualVectorL1(OrthancStone::Vector a,
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	633 OrthancStone::Vector b)
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	634 {
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	635 if (a.size() != b.size())
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	636 {
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	637 return false;
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	638 }
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	639 else
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	640 {
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	641 for (size_t i = 0; i < a.size(); i++)
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	642 {
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	643 if (!OrthancStone::LinearAlgebra::IsNear(a[i], b[i], 0.0001))
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	644 {
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	645 return false;
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	646 }
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	647 }
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	648
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	649 return true;
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	650 }
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	651 }
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	652
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	653
685 ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	654 TEST(VolumeImageGeometry, Basic)
ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	655 {
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	656 using namespace OrthancStone;
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	657
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	658 VolumeImageGeometry g;
953 118fc5c85d07 Fixed call to renamed method SetSize --> SetSizeInVoxels Benjamin Golinvaux <bgo@osimis.io> parents: 809 diff changeset	659 g.SetSizeInVoxels(10, 20, 30);
685 ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	660 g.SetVoxelDimensions(1, 2, 3);
ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	661
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	662 Vector p = g.GetCoordinates(0, 0, 0);
685 ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	663 ASSERT_EQ(3u, p.size());
999 2d69b8bee484 Added tests for Dicom structure set classes (loaders and utils) Benjamin Golinvaux <bgo@osimis.io> parents: 953 diff changeset	664 ASSERT_DOUBLE_EQ(-1.0 / 2.0, p[0]);
2d69b8bee484 Added tests for Dicom structure set classes (loaders and utils) Benjamin Golinvaux <bgo@osimis.io> parents: 953 diff changeset	665 ASSERT_DOUBLE_EQ(-2.0 / 2.0, p[1]);
2d69b8bee484 Added tests for Dicom structure set classes (loaders and utils) Benjamin Golinvaux <bgo@osimis.io> parents: 953 diff changeset	666 ASSERT_DOUBLE_EQ(-3.0 / 2.0, p[2]);
685 ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	667
ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	668 p = g.GetCoordinates(1, 1, 1);
999 2d69b8bee484 Added tests for Dicom structure set classes (loaders and utils) Benjamin Golinvaux <bgo@osimis.io> parents: 953 diff changeset	669 ASSERT_DOUBLE_EQ(-1.0 / 2.0 + 10.0 * 1.0, p[0]);
2d69b8bee484 Added tests for Dicom structure set classes (loaders and utils) Benjamin Golinvaux <bgo@osimis.io> parents: 953 diff changeset	670 ASSERT_DOUBLE_EQ(-2.0 / 2.0 + 20.0 * 2.0, p[1]);
2d69b8bee484 Added tests for Dicom structure set classes (loaders and utils) Benjamin Golinvaux <bgo@osimis.io> parents: 953 diff changeset	671 ASSERT_DOUBLE_EQ(-3.0 / 2.0 + 30.0 * 3.0, p[2]);
685 ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	672
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	673 VolumeProjection proj;
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	674 ASSERT_TRUE(g.DetectProjection(proj, g.GetAxialGeometry().GetNormal()));
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	675 ASSERT_EQ(VolumeProjection_Axial, proj);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	676 ASSERT_TRUE(g.DetectProjection(proj, g.GetCoronalGeometry().GetNormal()));
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	677 ASSERT_EQ(VolumeProjection_Coronal, proj);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	678 ASSERT_TRUE(g.DetectProjection(proj, g.GetSagittalGeometry().GetNormal()));
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	679 ASSERT_EQ(VolumeProjection_Sagittal, proj);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	680
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	681 ASSERT_EQ(10u, g.GetProjectionWidth(VolumeProjection_Axial));
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	682 ASSERT_EQ(20u, g.GetProjectionHeight(VolumeProjection_Axial));
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	683 ASSERT_EQ(30u, g.GetProjectionDepth(VolumeProjection_Axial));
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	684 ASSERT_EQ(10u, g.GetProjectionWidth(VolumeProjection_Coronal));
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	685 ASSERT_EQ(30u, g.GetProjectionHeight(VolumeProjection_Coronal));
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	686 ASSERT_EQ(20u, g.GetProjectionDepth(VolumeProjection_Coronal));
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	687 ASSERT_EQ(20u, g.GetProjectionWidth(VolumeProjection_Sagittal));
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	688 ASSERT_EQ(30u, g.GetProjectionHeight(VolumeProjection_Sagittal));
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	689 ASSERT_EQ(10u, g.GetProjectionDepth(VolumeProjection_Sagittal));
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	690
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	691 p = g.GetVoxelDimensions(VolumeProjection_Axial);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	692 ASSERT_EQ(3u, p.size());
999 2d69b8bee484 Added tests for Dicom structure set classes (loaders and utils) Benjamin Golinvaux <bgo@osimis.io> parents: 953 diff changeset	693 ASSERT_DOUBLE_EQ(1, p[0]);
2d69b8bee484 Added tests for Dicom structure set classes (loaders and utils) Benjamin Golinvaux <bgo@osimis.io> parents: 953 diff changeset	694 ASSERT_DOUBLE_EQ(2, p[1]);
2d69b8bee484 Added tests for Dicom structure set classes (loaders and utils) Benjamin Golinvaux <bgo@osimis.io> parents: 953 diff changeset	695 ASSERT_DOUBLE_EQ(3, p[2]);
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	696 p = g.GetVoxelDimensions(VolumeProjection_Coronal);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	697 ASSERT_EQ(3u, p.size());
999 2d69b8bee484 Added tests for Dicom structure set classes (loaders and utils) Benjamin Golinvaux <bgo@osimis.io> parents: 953 diff changeset	698 ASSERT_DOUBLE_EQ(1, p[0]);
2d69b8bee484 Added tests for Dicom structure set classes (loaders and utils) Benjamin Golinvaux <bgo@osimis.io> parents: 953 diff changeset	699 ASSERT_DOUBLE_EQ(3, p[1]);
2d69b8bee484 Added tests for Dicom structure set classes (loaders and utils) Benjamin Golinvaux <bgo@osimis.io> parents: 953 diff changeset	700 ASSERT_DOUBLE_EQ(2, p[2]);
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	701 p = g.GetVoxelDimensions(VolumeProjection_Sagittal);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	702 ASSERT_EQ(3u, p.size());
999 2d69b8bee484 Added tests for Dicom structure set classes (loaders and utils) Benjamin Golinvaux <bgo@osimis.io> parents: 953 diff changeset	703 ASSERT_DOUBLE_EQ(2, p[0]);
2d69b8bee484 Added tests for Dicom structure set classes (loaders and utils) Benjamin Golinvaux <bgo@osimis.io> parents: 953 diff changeset	704 ASSERT_DOUBLE_EQ(3, p[1]);
2d69b8bee484 Added tests for Dicom structure set classes (loaders and utils) Benjamin Golinvaux <bgo@osimis.io> parents: 953 diff changeset	705 ASSERT_DOUBLE_EQ(1, p[2]);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	706
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	707 // Loop over all the voxels of the volume
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	708 for (unsigned int z = 0; z < g.GetDepth(); z++)
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	709 {
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	710 const float zz = (0.5f + static_cast<float>(z)) / static_cast<float>(g.GetDepth()); // Z-center of the voxel
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	711
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	712 for (unsigned int y = 0; y < g.GetHeight(); y++)
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	713 {
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	714 const float yy = (0.5f + static_cast<float>(y)) / static_cast<float>(g.GetHeight()); // Y-center of the voxel
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	715
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	716 for (unsigned int x = 0; x < g.GetWidth(); x++)
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	717 {
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	718 const float xx = (0.5f + static_cast<float>(x)) / static_cast<float>(g.GetWidth()); // X-center of the voxel
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	719
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	720 const float sx = 1.0f;
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	721 const float sy = 2.0f;
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	722 const float sz = 3.0f;
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	723
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	724 Vector p = g.GetCoordinates(xx, yy, zz);
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	725
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	726 Vector q = (g.GetAxialGeometry().MapSliceToWorldCoordinates(
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	727 static_cast<double>(x) * sx,
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	728 static_cast<double>(y) * sy) +
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	729 z * sz * g.GetAxialGeometry().GetNormal());
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	730 ASSERT_TRUE(IsEqualVectorL1(p, q));
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	731
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	732 q = (g.GetCoronalGeometry().MapSliceToWorldCoordinates(
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	733 static_cast<double>(x) * sx,
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	734 static_cast<double>(g.GetDepth() - 1 - z) * sz) +
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	735 y * sy * g.GetCoronalGeometry().GetNormal());
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	736 ASSERT_TRUE(IsEqualVectorL1(p, q));
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	737
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	738 /**
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	739 * WARNING: In sagittal geometry, the normal points to
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	740 * REDUCING X-axis in the 3D world. This is necessary to keep
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	741 * the right-hand coordinate system. Hence the "-".
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	742 **/
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	743 q = (g.GetSagittalGeometry().MapSliceToWorldCoordinates(
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	744 static_cast<double>(y) * sy,
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	745 static_cast<double>(g.GetDepth() - 1 - z) * sz) +
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	746 x * sx * (-g.GetSagittalGeometry().GetNormal()));
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	747 ASSERT_TRUE(IsEqualVectorL1(p, q));
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	748 }
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	749 }
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	750 }
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	751
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	752 ASSERT_EQ(0, (int) VolumeProjection_Axial);
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	753 ASSERT_EQ(1, (int) VolumeProjection_Coronal);
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	754 ASSERT_EQ(2, (int) VolumeProjection_Sagittal);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	755
93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	756 for (int p = 0; p < 3; p++)
93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	757 {
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	758 VolumeProjection projection = (VolumeProjection) p;
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	759 const CoordinateSystem3D& s = g.GetProjectionGeometry(projection);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	760
809 e8fdf29cd0ca VolumeImageGeometry::GetProjectionSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 792 diff changeset	761 ASSERT_THROW(g.GetProjectionSlice(projection, g.GetProjectionDepth(projection)), Orthanc::OrthancException);
e8fdf29cd0ca VolumeImageGeometry::GetProjectionSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 792 diff changeset	762
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	763 for (unsigned int i = 0; i < g.GetProjectionDepth(projection); i++)
93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	764 {
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	765 CoordinateSystem3D plane = g.GetProjectionSlice(projection, i);
809 e8fdf29cd0ca VolumeImageGeometry::GetProjectionSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 792 diff changeset	766
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	767 if (projection == VolumeProjection_Sagittal)
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	768 {
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	769 ASSERT_TRUE(IsEqualVectorL1(plane.GetOrigin(), s.GetOrigin() + static_cast<double>(i) *
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	770 (-s.GetNormal()) * g.GetVoxelDimensions(projection)[2]));
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	771 }
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	772 else
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	773 {
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	774 ASSERT_TRUE(IsEqualVectorL1(plane.GetOrigin(), s.GetOrigin() + static_cast<double>(i) *
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	775 s.GetNormal() * g.GetVoxelDimensions(projection)[2]));
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	776 }
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	777
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	778 ASSERT_TRUE(IsEqualVectorL1(plane.GetAxisX(), s.GetAxisX()));
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	779 ASSERT_TRUE(IsEqualVectorL1(plane.GetAxisY(), s.GetAxisY()));
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	780
93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	781 unsigned int slice;
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	782 VolumeProjection q;
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	783 ASSERT_TRUE(g.DetectSlice(q, slice, plane));
93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	784 ASSERT_EQ(projection, q);
1167 ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	785 ASSERT_EQ(i, slice);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	786 }
93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	787 }
685 ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	788 }
ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	789
1167 ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	790
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	791 TEST(LinearAlgebra, ParseVectorLocale)
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	792 {
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	793 OrthancStone::Vector v;
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	794
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	795 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "1.2"));
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	796 ASSERT_EQ(1u, v.size());
1309 1f877e0846fe Fixed some warnings + fixed build for C++03 Benjamin Golinvaux <bgo@osimis.io> parents: 1279 diff changeset	797 ASSERT_DOUBLE_EQ(1.2, v[0]);
1167 ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	798
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	799 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "-1.2e+2"));
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	800 ASSERT_EQ(1u, v.size());
1309 1f877e0846fe Fixed some warnings + fixed build for C++03 Benjamin Golinvaux <bgo@osimis.io> parents: 1279 diff changeset	801 ASSERT_DOUBLE_EQ(-120.0, v[0]);
1167 ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	802
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	803 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "-1e-2\\2"));
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	804 ASSERT_EQ(2u, v.size());
1309 1f877e0846fe Fixed some warnings + fixed build for C++03 Benjamin Golinvaux <bgo@osimis.io> parents: 1279 diff changeset	805 ASSERT_DOUBLE_EQ(-0.01, v[0]);
1f877e0846fe Fixed some warnings + fixed build for C++03 Benjamin Golinvaux <bgo@osimis.io> parents: 1279 diff changeset	806 ASSERT_DOUBLE_EQ(2.0, v[1]);
1167 ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	807
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	808 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "1.3671875\\1.3671875"));
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	809 ASSERT_EQ(2u, v.size());
1309 1f877e0846fe Fixed some warnings + fixed build for C++03 Benjamin Golinvaux <bgo@osimis.io> parents: 1279 diff changeset	810 ASSERT_DOUBLE_EQ(1.3671875, v[0]);
1f877e0846fe Fixed some warnings + fixed build for C++03 Benjamin Golinvaux <bgo@osimis.io> parents: 1279 diff changeset	811 ASSERT_DOUBLE_EQ(1.3671875, v[1]);
1167 ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	812 }
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	813
1591 5887a4f8594b moving platform-specific files out of the "OrthancStone" folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 1587 diff changeset	814 TEST(GenericToolbox, NormalizeUuid)
1413 cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	815 {
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	816 std::string ref("44ca5051-14ef-4d2f-8bd7-db20bfb61fbb");
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	817
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	818 {
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	819 std::string u("44ca5051-14ef-4d2f-8bd7-db20bfb61fbb");
1591 5887a4f8594b moving platform-specific files out of the "OrthancStone" folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 1587 diff changeset	820 OrthancStone::GenericToolbox::NormalizeUuid(u);
1413 cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	821 ASSERT_EQ(ref, u);
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	822 }
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	823
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	824 // space left
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	825 {
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	826 std::string u(" 44ca5051-14ef-4d2f-8bd7-db20bfb61fbb");
1591 5887a4f8594b moving platform-specific files out of the "OrthancStone" folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 1587 diff changeset	827 OrthancStone::GenericToolbox::NormalizeUuid(u);
1413 cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	828 ASSERT_EQ(ref, u);
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	829 }
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	830
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	831 // space right
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	832 {
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	833 std::string u("44ca5051-14ef-4d2f-8bd7-db20bfb61fbb ");
1591 5887a4f8594b moving platform-specific files out of the "OrthancStone" folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 1587 diff changeset	834 OrthancStone::GenericToolbox::NormalizeUuid(u);
1413 cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	835 ASSERT_EQ(ref, u);
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	836 }
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	837
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	838 // space l & r
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	839 {
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	840 std::string u(" 44ca5051-14ef-4d2f-8bd7-db20bfb61fbb ");
1591 5887a4f8594b moving platform-specific files out of the "OrthancStone" folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 1587 diff changeset	841 OrthancStone::GenericToolbox::NormalizeUuid(u);
1413 cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	842 ASSERT_EQ(ref, u);
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	843 }
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	844
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	845 // space left + case
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	846 {
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	847 std::string u(" 44CA5051-14ef-4d2f-8bd7-dB20bfb61fbb");
1591 5887a4f8594b moving platform-specific files out of the "OrthancStone" folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 1587 diff changeset	848 OrthancStone::GenericToolbox::NormalizeUuid(u);
1413 cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	849 ASSERT_EQ(ref, u);
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	850 }
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	851
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	852 // space right + case
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	853 {
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	854 std::string u("44ca5051-14EF-4D2f-8bd7-db20bfb61fbB ");
1591 5887a4f8594b moving platform-specific files out of the "OrthancStone" folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 1587 diff changeset	855 OrthancStone::GenericToolbox::NormalizeUuid(u);
1413 cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	856 ASSERT_EQ(ref, u);
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	857 }
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	858
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	859 // space l & r + case
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	860 {
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	861 std::string u(" 44cA5051-14Ef-4d2f-8bD7-db20bfb61fbb ");
1591 5887a4f8594b moving platform-specific files out of the "OrthancStone" folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 1587 diff changeset	862 OrthancStone::GenericToolbox::NormalizeUuid(u);
1413 cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	863 ASSERT_EQ(ref, u);
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	864 }
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	865
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	866 // no
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	867 {
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	868 std::string u(" 44ca5051-14ef-4d2f-8bd7- db20bfb61fbb");
1591 5887a4f8594b moving platform-specific files out of the "OrthancStone" folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 1587 diff changeset	869 OrthancStone::GenericToolbox::NormalizeUuid(u);
1413 cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	870 ASSERT_NE(ref, u);
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	871 }
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	872
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	873 // no
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	874 {
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	875 std::string u("44ca5051-14ef-4d2f-8bd7-db20bfb61fb");
1591 5887a4f8594b moving platform-specific files out of the "OrthancStone" folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 1587 diff changeset	876 OrthancStone::GenericToolbox::NormalizeUuid(u);
1413 cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	877 ASSERT_NE(ref, u);
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	878 }
cde379b9d1d2 LoaderCache is now extensible + factored the uuid normalizing Benjamin Golinvaux <bgo@osimis.io> parents: 1403 diff changeset	879 }
1167 ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	880
1477 5732edec7cbd sorting frames in 3D Sebastien Jodogne <s.jodogne@gmail.com> parents: 1455 diff changeset	881
20 946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	882 int main(int argc, char **argv)
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	883 {
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	884 Orthanc::Logging::Initialize();
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	885 Orthanc::Logging::EnableInfoLevel(true);
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	886
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	887 ::testing::InitGoogleTest(&argc, argv);
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	888 int result = RUN_ALL_TESTS();
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	889
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	890 Orthanc::Logging::Finalize();
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	891
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	892 return result;
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	893 }

Mercurial > hg > orthanc-stone

annotate UnitTestsSources/UnitTestsMain.cpp @ 1603:595c0952ef7e