orthanc-stone: UnitTestsSources/UnitTestsMain.cpp annotate

annotate UnitTestsSources/UnitTestsMain.cpp @ 1327:4f8db2d202c8 broker

OrthancSeriesProgressiveLoader now has two modes that can be selected at object creation : - progressive (will first load jpeg50, then jpeg90 then PAM) - non-progressive (will directly load PAM (uncompressed)) Please note that the slice loading order remains dynamic and depending upon the slice that the client code wishes to extract from the volume.

author	Benjamin Golinvaux <bgo@osimis.io>
date	Wed, 25 Mar 2020 14:34:27 +0100
parents	1f877e0846fe
children	5630c2fb7b0f

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.
28956ed68280 agpl license Sebastien Jodogne <s.jodogne@gmail.com> parents: 40 diff changeset	16 *
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
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	22 #include "gtest/gtest.h"
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	23
732 c35e98d22764 move Deprecated classes to a separate folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 689 diff changeset	24 #include "../Framework/Deprecated/Layers/FrameRenderer.h"
c35e98d22764 move Deprecated classes to a separate folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 689 diff changeset	25 #include "../Framework/Deprecated/Toolbox/DownloadStack.h"
792 4fe4b221a31f deprecating MessagingToolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 732 diff changeset	26 #include "../Framework/Deprecated/Toolbox/MessagingToolbox.h"
732 c35e98d22764 move Deprecated classes to a separate folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 689 diff changeset	27 #include "../Framework/Deprecated/Toolbox/OrthancSlicesLoader.h"
1167 ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	28 #include "../Framework/StoneInitialization.h"
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	29 #include "../Framework/Toolbox/FiniteProjectiveCamera.h"
732 c35e98d22764 move Deprecated classes to a separate folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 689 diff changeset	30 #include "../Framework/Toolbox/GeometryToolbox.h"
87 4a541cd4fa83 OrthancVolumeImageLoader Sebastien Jodogne <s.jodogne@gmail.com> parents: 82 diff changeset	31 #include "../Framework/Volumes/ImageBuffer3D.h"
113 2eca030792aa using the Orthanc Framework Sebastien Jodogne <s.jodogne@gmail.com> parents: 102 diff changeset	32 #include "../Platforms/Generic/OracleWebService.h"
2eca030792aa using the Orthanc Framework Sebastien Jodogne <s.jodogne@gmail.com> parents: 102 diff changeset	33
212 5412adf19980 resort to OrthancFramework Sebastien Jodogne <s.jodogne@gmail.com> parents: 201 diff changeset	34 #include <Core/HttpClient.h>
5412adf19980 resort to OrthancFramework Sebastien Jodogne <s.jodogne@gmail.com> parents: 201 diff changeset	35 #include <Core/Images/ImageProcessing.h>
5412adf19980 resort to OrthancFramework Sebastien Jodogne <s.jodogne@gmail.com> parents: 201 diff changeset	36 #include <Core/Logging.h>
5412adf19980 resort to OrthancFramework Sebastien Jodogne <s.jodogne@gmail.com> parents: 201 diff changeset	37 #include <Core/MultiThreading/SharedMessageQueue.h>
5412adf19980 resort to OrthancFramework Sebastien Jodogne <s.jodogne@gmail.com> parents: 201 diff changeset	38 #include <Core/OrthancException.h>
87 4a541cd4fa83 OrthancVolumeImageLoader Sebastien Jodogne <s.jodogne@gmail.com> parents: 82 diff changeset	39
68 1526d38ef6da SliceSorter Sebastien Jodogne <s.jodogne@gmail.com> parents: 67 diff changeset	40 #include <boost/lexical_cast.hpp>
70 f73aed014bde OrthancAsynchronousWebService Sebastien Jodogne <s.jodogne@gmail.com> parents: 69 diff changeset	41 #include <boost/date_time/posix_time/posix_time.hpp>
f73aed014bde OrthancAsynchronousWebService Sebastien Jodogne <s.jodogne@gmail.com> parents: 69 diff changeset	42 #include <boost/thread/thread.hpp>
102 fcec0ab44054 display volumes Sebastien Jodogne <s.jodogne@gmail.com> parents: 93 diff changeset	43 #include <boost/math/special_functions/round.hpp>
fcec0ab44054 display volumes Sebastien Jodogne <s.jodogne@gmail.com> parents: 93 diff changeset	44
68 1526d38ef6da SliceSorter Sebastien Jodogne <s.jodogne@gmail.com> parents: 67 diff changeset	45
144 9b83f30fc1c0 bilinear and trilinear interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 135 diff changeset	46 TEST(GeometryToolbox, Interpolation)
9b83f30fc1c0 bilinear and trilinear interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 135 diff changeset	47 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	48 using namespace OrthancStone::GeometryToolbox;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	49
144 9b83f30fc1c0 bilinear and trilinear interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 135 diff changeset	50 // 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	51 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	52
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	53 ASSERT_FLOAT_EQ(91, ComputeBilinearInterpolationUnitSquare(0, 0, 91, 210, 162, 95));
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	54 ASSERT_FLOAT_EQ(210, ComputeBilinearInterpolationUnitSquare(1, 0, 91, 210, 162, 95));
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	55 ASSERT_FLOAT_EQ(162, ComputeBilinearInterpolationUnitSquare(0, 1, 91, 210, 162, 95));
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	56 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	57
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	58 ASSERT_FLOAT_EQ(123.35f, ComputeTrilinearInterpolationUnitSquare
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	59 (0.5f, 0.2f, 0.7f,
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	60 91, 210, 162, 95,
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	61 51, 190, 80, 92));
176 ab9c799f5de1 fix interpolation with negative coordinates Sebastien Jodogne <s.jodogne@gmail.com> parents: 165 diff changeset	62
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	63 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	64 ComputeTrilinearInterpolationUnitSquare(0.5f, 0.2f, 0,
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));
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	67
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	68 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	69 ComputeTrilinearInterpolationUnitSquare(0.5f, 0.2f, 1,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	70 91, 210, 162, 95,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 176 diff changeset	71 51, 190, 80, 92));
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	72 }
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	73
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	74
162 77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	75 static bool CompareMatrix(const OrthancStone::Matrix& a,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	76 const OrthancStone::Matrix& b,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	77 double threshold = 0.00000001)
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 if (a.size1() != b.size1() \|\|
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	80 a.size2() != b.size2())
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 return false;
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
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	85 for (size_t i = 0; i < a.size1(); i++)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	86 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	87 for (size_t j = 0; j < a.size2(); j++)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	88 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	89 if (fabs(a(i, j) - b(i, j)) > threshold)
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 LOG(ERROR) << "Too large difference in component ("
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	92 << i << "," << j << "): " << a(i,j) << " != " << b(i,j);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	93 return false;
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
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	98 return true;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	99 }
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
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	102 static bool CompareVector(const OrthancStone::Vector& a,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	103 const OrthancStone::Vector& b,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	104 double threshold = 0.00000001)
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 if (a.size() != b.size())
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 return false;
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
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	111 for (size_t i = 0; i < a.size(); i++)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	112 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	113 if (fabs(a(i) - b(i)) > threshold)
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 LOG(ERROR) << "Too large difference in component "
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	116 << i << ": " << a(i) << " != " << b(i);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	117 return false;
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 return true;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	122 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	123
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	124
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	125
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	126 TEST(FiniteProjectiveCamera, Decomposition1)
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	127 {
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	128 // Example 6.2 of "Multiple View Geometry in Computer Vision - 2nd
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	129 // edition" (page 163)
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	130 const double p[12] = {
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	131 3.53553e+2, 3.39645e+2, 2.77744e+2, -1.44946e+6,
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	132 -1.03528e+2, 2.33212e+1, 4.59607e+2, -6.32525e+5,
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	133 7.07107e-1, -3.53553e-1, 6.12372e-1, -9.18559e+2
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	134 };
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	135
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	136 OrthancStone::FiniteProjectiveCamera camera(p);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	137 ASSERT_EQ(3u, camera.GetMatrix().size1());
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	138 ASSERT_EQ(4u, camera.GetMatrix().size2());
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	139 ASSERT_EQ(3u, camera.GetIntrinsicParameters().size1());
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	140 ASSERT_EQ(3u, camera.GetIntrinsicParameters().size2());
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	141 ASSERT_EQ(3u, camera.GetRotation().size1());
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	142 ASSERT_EQ(3u, camera.GetRotation().size2());
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	143 ASSERT_EQ(3u, camera.GetCenter().size());
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	144
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	145 ASSERT_NEAR(1000.0, camera.GetCenter()[0], 0.01);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	146 ASSERT_NEAR(2000.0, camera.GetCenter()[1], 0.01);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	147 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	148
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	149 ASSERT_NEAR(468.2, camera.GetIntrinsicParameters() (0, 0), 0.1);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	150 ASSERT_NEAR(91.2, camera.GetIntrinsicParameters() (0, 1), 0.1);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	151 ASSERT_NEAR(300.0, camera.GetIntrinsicParameters() (0, 2), 0.1);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	152 ASSERT_NEAR(427.2, camera.GetIntrinsicParameters() (1, 1), 0.1);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	153 ASSERT_NEAR(200.0, camera.GetIntrinsicParameters() (1, 2), 0.1);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	154 ASSERT_NEAR(1.0, camera.GetIntrinsicParameters() (2, 2), 0.1);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	155
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	156 ASSERT_NEAR(0, camera.GetIntrinsicParameters() (1, 0), 0.0000001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	157 ASSERT_NEAR(0, camera.GetIntrinsicParameters() (2, 0), 0.0000001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	158 ASSERT_NEAR(0, camera.GetIntrinsicParameters() (2, 1), 0.0000001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	159
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	160 ASSERT_NEAR(0.41380, camera.GetRotation() (0, 0), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	161 ASSERT_NEAR(0.90915, camera.GetRotation() (0, 1), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	162 ASSERT_NEAR(0.04708, camera.GetRotation() (0, 2), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	163 ASSERT_NEAR(-0.57338, camera.GetRotation() (1, 0), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	164 ASSERT_NEAR(0.22011, camera.GetRotation() (1, 1), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	165 ASSERT_NEAR(0.78917, camera.GetRotation() (1, 2), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	166 ASSERT_NEAR(0.70711, camera.GetRotation() (2, 0), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	167 ASSERT_NEAR(-0.35355, camera.GetRotation() (2, 1), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	168 ASSERT_NEAR(0.61237, camera.GetRotation() (2, 2), 0.00001);
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	169
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	170 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	171
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	172 OrthancStone::FiniteProjectiveCamera camera2(camera.GetIntrinsicParameters(),
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	173 camera.GetRotation(),
197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	174 camera.GetCenter());
162 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 ASSERT_TRUE(CompareMatrix(camera.GetMatrix(), camera2.GetMatrix()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	177 ASSERT_TRUE(CompareMatrix(camera.GetIntrinsicParameters(), camera2.GetIntrinsicParameters()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	178 ASSERT_TRUE(CompareMatrix(camera.GetRotation(), camera2.GetRotation()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	179 ASSERT_TRUE(CompareVector(camera.GetCenter(), camera2.GetCenter()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	180 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	181
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	182
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	183 TEST(FiniteProjectiveCamera, Decomposition2)
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 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	186 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	187 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	188 const double c[] = { 243.558936, -145.863085, 411.585964 };
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	189
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	190 OrthancStone::FiniteProjectiveCamera camera(p);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	191 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	192
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	193 OrthancStone::FiniteProjectiveCamera camera2(k, r, c);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	194 ASSERT_TRUE(CompareMatrix(camera.GetMatrix(), camera2.GetMatrix(), 1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	195 ASSERT_TRUE(CompareMatrix(camera.GetIntrinsicParameters(), camera2.GetIntrinsicParameters(), 0.001));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	196 ASSERT_TRUE(CompareMatrix(camera.GetRotation(), camera2.GetRotation(), 0.000001));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	197 ASSERT_TRUE(CompareVector(camera.GetCenter(), camera2.GetCenter(), 0.0001));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	198 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	199
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	200
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	201 TEST(FiniteProjectiveCamera, Decomposition3)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	202 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	203 const double p[] = { 10, 0, 0, 0,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	204 0, 20, 0, 0,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	205 0, 0, 30, 0 };
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	206
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	207 OrthancStone::FiniteProjectiveCamera camera(p);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	208 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	209 ASSERT_DOUBLE_EQ(10, camera.GetIntrinsicParameters() (0, 0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	210 ASSERT_DOUBLE_EQ(20, camera.GetIntrinsicParameters() (1, 1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	211 ASSERT_DOUBLE_EQ(30, camera.GetIntrinsicParameters() (2, 2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	212 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (0, 0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	213 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (1, 1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	214 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (2, 2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	215 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	216 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	217 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	218 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	219
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	220
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	221 TEST(FiniteProjectiveCamera, Decomposition4)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	222 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	223 const double p[] = { 1, 0, 0, 10,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	224 0, 1, 0, 20,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	225 0, 0, 1, 30 };
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	226
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	227 OrthancStone::FiniteProjectiveCamera camera(p);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	228 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	229 ASSERT_DOUBLE_EQ(1, camera.GetIntrinsicParameters() (0, 0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	230 ASSERT_DOUBLE_EQ(1, camera.GetIntrinsicParameters() (1, 1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	231 ASSERT_DOUBLE_EQ(1, camera.GetIntrinsicParameters() (2, 2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	232 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (0, 0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	233 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (1, 1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	234 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (2, 2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	235 ASSERT_DOUBLE_EQ(-10, camera.GetCenter() (0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	236 ASSERT_DOUBLE_EQ(-20, camera.GetCenter() (1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	237 ASSERT_DOUBLE_EQ(-30, camera.GetCenter() (2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	238 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	239
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	240
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	241 TEST(FiniteProjectiveCamera, Decomposition5)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	242 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	243 const double p[] = { 0, 0, 10, 0,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	244 0, 20, 0, 0,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	245 30, 0, 0, 0 };
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	246
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	247 OrthancStone::FiniteProjectiveCamera camera(p);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	248 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(camera.GetRotation()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	249 ASSERT_DOUBLE_EQ(-10, camera.GetIntrinsicParameters() (0, 0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	250 ASSERT_DOUBLE_EQ(20, camera.GetIntrinsicParameters() (1, 1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	251 ASSERT_DOUBLE_EQ(30, camera.GetIntrinsicParameters() (2, 2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	252 ASSERT_DOUBLE_EQ(-1, camera.GetRotation() (0, 2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	253 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (1, 1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	254 ASSERT_DOUBLE_EQ(1, camera.GetRotation() (2, 0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	255 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (0));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	256 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (1));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	257 ASSERT_DOUBLE_EQ(0, camera.GetCenter() (2));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	258
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	259 OrthancStone::FiniteProjectiveCamera camera2(camera.GetIntrinsicParameters(),
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	260 camera.GetRotation(),
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	261 camera.GetCenter());
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	262 ASSERT_TRUE(CompareMatrix(camera.GetMatrix(), camera2.GetMatrix()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	263 ASSERT_TRUE(CompareMatrix(camera.GetIntrinsicParameters(), camera2.GetIntrinsicParameters()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	264 ASSERT_TRUE(CompareMatrix(camera.GetRotation(), camera2.GetRotation()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	265 ASSERT_TRUE(CompareVector(camera.GetCenter(), camera2.GetCenter()));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	266 }
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	267
162 77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	268
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	269 static double GetCosAngle(const OrthancStone::Vector& a,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	270 const OrthancStone::Vector& 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 // Returns the cosine of the angle between two vectors
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	273 // https://en.wikipedia.org/wiki/Dot_product#Geometric_definition
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	274 return boost::numeric::ublas::inner_prod(a, b) /
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	275 (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	276 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	277
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	278
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	279 TEST(FiniteProjectiveCamera, Ray)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	280 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	281 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	282 -2951.517707, -1501.019129, -285.785281, 637891.819097,
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	283 0.008528, 0.003067, -0.999959, 2491.764918 };
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	284
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	285 const OrthancStone::FiniteProjectiveCamera camera(pp);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	286
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	287 ASSERT_NEAR(-21.2492, camera.GetCenter() (0), 0.0001);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	288 ASSERT_NEAR(-7.64234, camera.GetCenter() (1), 0.00001);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	289 ASSERT_NEAR(2491.66, camera.GetCenter() (2), 0.01);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	290
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	291 // 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	292 // (256,256). The image has dimensions 512x512.
165 8d50e6be565d LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 164 diff changeset	293 OrthancStone::Vector o =
8d50e6be565d LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 164 diff changeset	294 OrthancStone::LinearAlgebra::CreateVector(7.009620, 2.521030, -821.942000);
8d50e6be565d LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 164 diff changeset	295 OrthancStone::Vector ax =
8d50e6be565d LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 164 diff changeset	296 OrthancStone::LinearAlgebra::CreateVector(-0.453219, 0.891399, -0.001131);
8d50e6be565d LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 164 diff changeset	297 OrthancStone::Vector ay =
8d50e6be565d LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 164 diff changeset	298 OrthancStone::LinearAlgebra::CreateVector(-0.891359, -0.453210, -0.008992);
8d50e6be565d LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 164 diff changeset	299
162 77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	300 OrthancStone::CoordinateSystem3D imagePlane(o, ax, ay);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	301
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	302 // Back-projection of the principal point
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	303 {
162 77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	304 OrthancStone::Vector ray = camera.GetRayDirection(256, 256);
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 // The principal axis vector is orthogonal to the image plane
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	307 // (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	308 // ("-1" corresponds to "cos(pi)").
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	309 ASSERT_NEAR(-1, GetCosAngle(ray, imagePlane.GetNormal()), 0.0000001);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	310
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	311 // Forward projection of principal axis, resulting in the principal point
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	312 double x, y;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	313 camera.ApplyFinite(x, y, camera.GetCenter() - ray);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	314
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	315 ASSERT_NEAR(256, x, 0.00001);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	316 ASSERT_NEAR(256, y, 0.00001);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	317 }
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	318
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	319 // Back-projection of the 4 corners of the image
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	320 std::vector<double> cx, cy;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	321 cx.push_back(0);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	322 cy.push_back(0);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	323 cx.push_back(512);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	324 cy.push_back(0);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	325 cx.push_back(512);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	326 cy.push_back(512);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	327 cx.push_back(0);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	328 cy.push_back(512);
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 bool first = true;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	331 double angle;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	332
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	333 for (size_t i = 0; i < cx.size(); i++)
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	334 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	335 OrthancStone::Vector ray = camera.GetRayDirection(cx[i], cy[i]);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	336
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	337 // 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	338 // the 4 corners
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	339 double a = GetCosAngle(ray, imagePlane.GetNormal());
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	340 if (first)
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	341 {
162 77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	342 first = false;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	343 angle = a;
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 else
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	346 {
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	347 ASSERT_NEAR(angle, a, 0.000001);
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	348 }
162 77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	349
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	350 // 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	351 double x, y;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	352 camera.ApplyFinite(x, y, camera.GetCenter() - ray);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	353
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	354 ASSERT_NEAR(cx[i], x, 0.00001);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	355 ASSERT_NEAR(cy[i], y, 0.00001);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	356
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	357 // Alternative construction, by computing the intersection of the
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	358 // ray with the image plane
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	359 OrthancStone::Vector p;
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	360 ASSERT_TRUE(imagePlane.IntersectLine(p, camera.GetCenter(), -ray));
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	361 imagePlane.ProjectPoint(x, y, p);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	362 ASSERT_NEAR(cx[i], x + 256, 0.01);
77715c340767 unit testing FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 161 diff changeset	363 ASSERT_NEAR(cy[i], y + 256, 0.01);
161 197a5ddaf68c FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 144 diff changeset	364 }
144 9b83f30fc1c0 bilinear and trilinear interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 135 diff changeset	365 }
87 4a541cd4fa83 OrthancVolumeImageLoader Sebastien Jodogne <s.jodogne@gmail.com> parents: 82 diff changeset	366
20 946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	367
163 8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	368 TEST(Matrix, Inverse1)
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	369 {
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	370 OrthancStone::Matrix a, b;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	371
164 432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	372 a.resize(0, 0);
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	373 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	374 ASSERT_EQ(0u, b.size1());
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	375 ASSERT_EQ(0u, b.size2());
163 8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	376
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	377 a.resize(2, 3);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	378 ASSERT_THROW(OrthancStone::LinearAlgebra::InvertMatrix(b, a), Orthanc::OrthancException);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	379
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	380 a.resize(1, 1);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	381 a(0, 0) = 45.0;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	382
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	383 ASSERT_DOUBLE_EQ(45, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	384 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	385 ASSERT_EQ(1u, b.size1());
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	386 ASSERT_EQ(1u, b.size2());
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	387 ASSERT_DOUBLE_EQ(1.0 / 45.0, b(0, 0));
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 a(0, 0) = 0;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	390 ASSERT_DOUBLE_EQ(0, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	391 ASSERT_THROW(OrthancStone::LinearAlgebra::InvertMatrix(b, a), Orthanc::OrthancException);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	392 }
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	393
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	394
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	395 TEST(Matrix, Inverse2)
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	396 {
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	397 OrthancStone::Matrix a, b;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	398 a.resize(2, 2);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	399 a(0, 0) = 4;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	400 a(0, 1) = 3;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	401 a(1, 0) = 3;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	402 a(1, 1) = 2;
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(-1, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	405 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	406 ASSERT_EQ(2u, b.size1());
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	407 ASSERT_EQ(2u, b.size2());
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 ASSERT_DOUBLE_EQ(-2, b(0, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	410 ASSERT_DOUBLE_EQ(3, b(0, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	411 ASSERT_DOUBLE_EQ(3, b(1, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	412 ASSERT_DOUBLE_EQ(-4, b(1, 1));
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 a(0, 0) = 1;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	415 a(0, 1) = 2;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	416 a(1, 0) = 3;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	417 a(1, 1) = 4;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	418
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	419 ASSERT_DOUBLE_EQ(-2, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	420 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
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 ASSERT_DOUBLE_EQ(-2, b(0, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	423 ASSERT_DOUBLE_EQ(1, b(0, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	424 ASSERT_DOUBLE_EQ(1.5, b(1, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	425 ASSERT_DOUBLE_EQ(-0.5, b(1, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	426 }
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	427
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	428
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	429 TEST(Matrix, Inverse3)
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	430 {
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	431 OrthancStone::Matrix a, b;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	432 a.resize(3, 3);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	433 a(0, 0) = 7;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	434 a(0, 1) = 2;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	435 a(0, 2) = 1;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	436 a(1, 0) = 0;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	437 a(1, 1) = 3;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	438 a(1, 2) = -1;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	439 a(2, 0) = -3;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	440 a(2, 1) = 4;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	441 a(2, 2) = -2;
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(1, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	444 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	445 ASSERT_EQ(3u, b.size1());
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	446 ASSERT_EQ(3u, b.size2());
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	447
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	448 ASSERT_DOUBLE_EQ(-2, b(0, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	449 ASSERT_DOUBLE_EQ(8, b(0, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	450 ASSERT_DOUBLE_EQ(-5, b(0, 2));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	451 ASSERT_DOUBLE_EQ(3, b(1, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	452 ASSERT_DOUBLE_EQ(-11, b(1, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	453 ASSERT_DOUBLE_EQ(7, b(1, 2));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	454 ASSERT_DOUBLE_EQ(9, b(2, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	455 ASSERT_DOUBLE_EQ(-34, b(2, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	456 ASSERT_DOUBLE_EQ(21, b(2, 2));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	457
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	458
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	459 a(0, 0) = 1;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	460 a(0, 1) = 2;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	461 a(0, 2) = 2;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	462 a(1, 0) = 1;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	463 a(1, 1) = 0;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	464 a(1, 2) = 1;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	465 a(2, 0) = 1;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	466 a(2, 1) = 2;
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	467 a(2, 2) = 1;
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(2, OrthancStone::LinearAlgebra::ComputeDeterminant(a));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	470 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	471 ASSERT_EQ(3u, b.size1());
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	472 ASSERT_EQ(3u, b.size2());
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	473
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	474 ASSERT_DOUBLE_EQ(-1, b(0, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	475 ASSERT_DOUBLE_EQ(1, b(0, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	476 ASSERT_DOUBLE_EQ(1, b(0, 2));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	477 ASSERT_DOUBLE_EQ(0, b(1, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	478 ASSERT_DOUBLE_EQ(-0.5, b(1, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	479 ASSERT_DOUBLE_EQ(0.5, b(1, 2));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	480 ASSERT_DOUBLE_EQ(1, b(2, 0));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	481 ASSERT_DOUBLE_EQ(0, b(2, 1));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	482 ASSERT_DOUBLE_EQ(-1, b(2, 2));
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	483 }
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	484
8c5b24892ed2 LinearAlgebra::InvertMatrix Sebastien Jodogne <s.jodogne@gmail.com> parents: 162 diff changeset	485
164 432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	486 TEST(Matrix, Inverse4)
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	487 {
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	488 OrthancStone::Matrix a, b;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	489 a.resize(4, 4);
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	490 a(0, 0) = 2;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	491 a(0, 1) = 1;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	492 a(0, 2) = 2;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	493 a(0, 3) = -3;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	494 a(1, 0) = -2;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	495 a(1, 1) = 2;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	496 a(1, 2) = -1;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	497 a(1, 3) = -1;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	498 a(2, 0) = 2;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	499 a(2, 1) = 2;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	500 a(2, 2) = -3;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	501 a(2, 3) = -1;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	502 a(3, 0) = 3;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	503 a(3, 1) = -2;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	504 a(3, 2) = -3;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	505 a(3, 3) = -1;
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	506
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	507 OrthancStone::LinearAlgebra::InvertMatrix(b, a);
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	508 ASSERT_EQ(4u, b.size1());
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	509 ASSERT_EQ(4u, b.size2());
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	510
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	511 b *= 134.0; // This is the determinant
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	512
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	513 ASSERT_DOUBLE_EQ(8, b(0, 0));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	514 ASSERT_DOUBLE_EQ(-44, b(0, 1));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	515 ASSERT_DOUBLE_EQ(30, b(0, 2));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	516 ASSERT_DOUBLE_EQ(-10, b(0, 3));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	517 ASSERT_DOUBLE_EQ(2, b(1, 0));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	518 ASSERT_DOUBLE_EQ(-11, b(1, 1));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	519 ASSERT_DOUBLE_EQ(41, b(1, 2));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	520 ASSERT_DOUBLE_EQ(-36, b(1, 3));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	521 ASSERT_DOUBLE_EQ(16, b(2, 0));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	522 ASSERT_DOUBLE_EQ(-21, b(2, 1));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	523 ASSERT_DOUBLE_EQ(-7, b(2, 2));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	524 ASSERT_DOUBLE_EQ(-20, b(2, 3));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	525 ASSERT_DOUBLE_EQ(-28, b(3, 0));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	526 ASSERT_DOUBLE_EQ(-47, b(3, 1));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	527 ASSERT_DOUBLE_EQ(29, b(3, 2));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	528 ASSERT_DOUBLE_EQ(-32, b(3, 3));
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	529 }
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	530
432b1f812d14 inversion of general matrices Sebastien Jodogne <s.jodogne@gmail.com> parents: 163 diff changeset	531
188 45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	532 TEST(FiniteProjectiveCamera, Calibration)
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	533 {
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	534 unsigned int volumeWidth = 512;
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	535 unsigned int volumeHeight = 512;
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	536 unsigned int volumeDepth = 110;
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	537
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	538 OrthancStone::Vector camera = OrthancStone::LinearAlgebra::CreateVector
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	539 (-1000, -5000, -static_cast<double>(volumeDepth) * 32);
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 OrthancStone::Vector principalPoint = OrthancStone::LinearAlgebra::CreateVector
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	542 (volumeWidth/2, volumeHeight/2, volumeDepth * 2);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	543
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	544 OrthancStone::FiniteProjectiveCamera c(camera, principalPoint, 0, 512, 512, 1, 1);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	545
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	546 double swapv[9] = { 1, 0, 0,
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	547 0, -1, 512,
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	548 0, 0, 1 };
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	549 OrthancStone::Matrix swap;
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	550 OrthancStone::LinearAlgebra::FillMatrix(swap, 3, 3, swapv);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	551
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	552 OrthancStone::Matrix p = OrthancStone::LinearAlgebra::Product(swap, c.GetMatrix());
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	553 p /= p(2,3);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	554
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	555 ASSERT_NEAR( 1.04437, p(0,0), 0.00001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	556 ASSERT_NEAR(-0.0703111, p(0,1), 0.00000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	557 ASSERT_NEAR(-0.179283, p(0,2), 0.000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	558 ASSERT_NEAR( 61.7431, p(0,3), 0.0001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	559 ASSERT_NEAR( 0.11127, p(1,0), 0.000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	560 ASSERT_NEAR(-0.595541, p(1,1), 0.000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	561 ASSERT_NEAR( 0.872211, p(1,2), 0.000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	562 ASSERT_NEAR( 203.748, p(1,3), 0.001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	563 ASSERT_NEAR( 3.08593e-05, p(2,0), 0.0000000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	564 ASSERT_NEAR( 0.000129138, p(2,1), 0.000000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	565 ASSERT_NEAR( 9.18901e-05, p(2,2), 0.0000000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	566 ASSERT_NEAR( 1, p(2,3), 0.0000001);
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	567 }
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	568
45b03b04a777 calibration of FiniteProjectiveCamera Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	569
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	570 static bool IsEqualRotationVector(OrthancStone::Vector a,
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	571 OrthancStone::Vector b)
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	572 {
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	573 if (a.size() != b.size() \|\|
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	574 a.size() != 3)
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	575 {
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	576 return false;
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	577 }
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	578 else
189 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 OrthancStone::LinearAlgebra::NormalizeVector(a);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	581 OrthancStone::LinearAlgebra::NormalizeVector(b);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	582 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	583 }
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	584 }
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	585
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 TEST(GeometryToolbox, AlignVectorsWithRotation)
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	588 {
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	589 OrthancStone::Vector a, b;
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	590 OrthancStone::Matrix r;
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	591
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	592 OrthancStone::LinearAlgebra::AssignVector(a, -200, 200, -846.63);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	593 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	594
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	595 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	596 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	597 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	598
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	599 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, b, a);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	600 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	601 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	602
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	603 OrthancStone::LinearAlgebra::AssignVector(a, 1, 0, 0);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	604 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	605 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	606 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	607 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	608
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	609 OrthancStone::LinearAlgebra::AssignVector(a, 0, 1, 0);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	610 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	611 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	612 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	613 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	614
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	615 OrthancStone::LinearAlgebra::AssignVector(a, 0, 0, 1);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	616 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	617 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	618 ASSERT_TRUE(OrthancStone::LinearAlgebra::IsRotationMatrix(r));
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	619 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	620
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	621 OrthancStone::LinearAlgebra::AssignVector(a, 0, 0, 0);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	622 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	623 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	624
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	625 // TODO: Deal with opposite vectors
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	626
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	627 /*
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	628 OrthancStone::LinearAlgebra::AssignVector(a, 0, 0, -1);
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	629 OrthancStone::LinearAlgebra::AssignVector(b, 0, 0, 1);
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	630 OrthancStone::GeometryToolbox::AlignVectorsWithRotation(r, a, b);
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	631 OrthancStone::LinearAlgebra::Print(r);
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	632 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	633 */
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	634 }
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	635
307 be2660b6e40a wip: commands + status update am@osimis.io parents: 257 diff changeset	636 TEST(MessagingToolbox, ParseJson)
be2660b6e40a wip: commands + status update am@osimis.io parents: 257 diff changeset	637 {
be2660b6e40a wip: commands + status update am@osimis.io parents: 257 diff changeset	638 Json::Value response;
be2660b6e40a wip: commands + status update am@osimis.io parents: 257 diff changeset	639 std::string source = "{\"command\":\"panel:takeDarkImage\",\"commandType\":\"simple\",\"args\":{}}";
792 4fe4b221a31f deprecating MessagingToolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 732 diff changeset	640 ASSERT_TRUE(Deprecated::MessagingToolbox::ParseJson(response, source.c_str(), source.size()));
307 be2660b6e40a wip: commands + status update am@osimis.io parents: 257 diff changeset	641 }
20 946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	642
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	643
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	644
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	645 static bool IsEqualVectorL1(OrthancStone::Vector a,
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	646 OrthancStone::Vector b)
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	647 {
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	648 if (a.size() != b.size())
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	649 {
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	650 return false;
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 else
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	653 {
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	654 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	655 {
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	656 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	657 {
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	658 return false;
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	659 }
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	660 }
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	661
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	662 return true;
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	663 }
189 964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	664 }
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	665
964118e7e6de unit test for AlignVectorsWithRotation Sebastien Jodogne <s.jodogne@gmail.com> parents: 188 diff changeset	666
685 ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	667 TEST(VolumeImageGeometry, Basic)
ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	668 {
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	669 using namespace OrthancStone;
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	670
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	671 VolumeImageGeometry g;
953 118fc5c85d07 Fixed call to renamed method SetSize --> SetSizeInVoxels Benjamin Golinvaux <bgo@osimis.io> parents: 809 diff changeset	672 g.SetSizeInVoxels(10, 20, 30);
685 ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	673 g.SetVoxelDimensions(1, 2, 3);
ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	674
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	675 Vector p = g.GetCoordinates(0, 0, 0);
685 ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	676 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	677 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	678 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	679 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	680
ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	681 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	682 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	683 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	684 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	685
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	686 VolumeProjection proj;
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	687 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	688 ASSERT_EQ(VolumeProjection_Axial, proj);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	689 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	690 ASSERT_EQ(VolumeProjection_Coronal, proj);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	691 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	692 ASSERT_EQ(VolumeProjection_Sagittal, proj);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	693
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	694 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	695 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	696 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	697 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	698 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	699 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	700 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	701 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	702 ASSERT_EQ(10u, g.GetProjectionDepth(VolumeProjection_Sagittal));
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	703
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	704 p = g.GetVoxelDimensions(VolumeProjection_Axial);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	705 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	706 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	707 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	708 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	709 p = g.GetVoxelDimensions(VolumeProjection_Coronal);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	710 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	711 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	712 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	713 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	714 p = g.GetVoxelDimensions(VolumeProjection_Sagittal);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	715 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	716 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	717 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	718 ASSERT_DOUBLE_EQ(1, p[2]);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	719
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	720 // 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	721 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	722 {
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	723 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	724
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	725 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	726 {
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	727 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	728
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	729 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	730 {
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	731 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	732
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	733 const float sx = 1.0f;
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	734 const float sy = 2.0f;
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	735 const float sz = 3.0f;
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	736
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	737 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	738
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	739 Vector q = (g.GetAxialGeometry().MapSliceToWorldCoordinates(
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	740 static_cast<double>(x) * sx,
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	741 static_cast<double>(y) * sy) +
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	742 z * sz * g.GetAxialGeometry().GetNormal());
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	743 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	744
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	745 q = (g.GetCoronalGeometry().MapSliceToWorldCoordinates(
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	746 static_cast<double>(x) * sx,
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	747 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	748 y * sy * g.GetCoronalGeometry().GetNormal());
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	749 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	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 * 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	753 * 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	754 * 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	755 **/
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	756 q = (g.GetSagittalGeometry().MapSliceToWorldCoordinates(
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	757 static_cast<double>(y) * sy,
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	758 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	759 x * sx * (-g.GetSagittalGeometry().GetNormal()));
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	760 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	761 }
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	762 }
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	763 }
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	764
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	765 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	766 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	767 ASSERT_EQ(2, (int) VolumeProjection_Sagittal);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	768
93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	769 for (int p = 0; p < 3; p++)
93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	770 {
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	771 VolumeProjection projection = (VolumeProjection) p;
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	772 const CoordinateSystem3D& s = g.GetProjectionGeometry(projection);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	773
809 e8fdf29cd0ca VolumeImageGeometry::GetProjectionSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 792 diff changeset	774 ASSERT_THROW(g.GetProjectionSlice(projection, g.GetProjectionDepth(projection)), Orthanc::OrthancException);
e8fdf29cd0ca VolumeImageGeometry::GetProjectionSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 792 diff changeset	775
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	776 for (unsigned int i = 0; i < g.GetProjectionDepth(projection); i++)
93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	777 {
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	778 CoordinateSystem3D plane = g.GetProjectionSlice(projection, i);
809 e8fdf29cd0ca VolumeImageGeometry::GetProjectionSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 792 diff changeset	779
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	780 if (projection == VolumeProjection_Sagittal)
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	781 {
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	782 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	783 (-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	784 }
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	785 else
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	786 {
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	787 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	788 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	789 }
1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	790
1172 6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	791 ASSERT_TRUE(IsEqualVectorL1(plane.GetAxisX(), s.GetAxisX()));
6def5bfba867 fix unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 1170 diff changeset	792 ASSERT_TRUE(IsEqualVectorL1(plane.GetAxisY(), s.GetAxisY()));
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	793
93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	794 unsigned int slice;
1170 1644de437a7b fixes related to swapped normal in sagittal geometry Sebastien Jodogne <s.jodogne@gmail.com> parents: 1167 diff changeset	795 VolumeProjection q;
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	796 ASSERT_TRUE(g.DetectSlice(q, slice, plane));
93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	797 ASSERT_EQ(projection, q);
1167 ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	798 ASSERT_EQ(i, slice);
689 93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	799 }
93a8949a1ef7 VolumeImageGeometry::DetectSlice() Sebastien Jodogne <s.jodogne@gmail.com> parents: 685 diff changeset	800 }
685 ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	801 }
ea1a5b963798 unit test: VolumeImageGeometry.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	802
1167 ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	803
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	804 TEST(LinearAlgebra, ParseVectorLocale)
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	805 {
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	806 OrthancStone::Vector v;
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.2"));
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	809 ASSERT_EQ(1u, 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.2, v[0]);
1167 ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	811
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	812 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "-1.2e+2"));
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	813 ASSERT_EQ(1u, v.size());
1309 1f877e0846fe Fixed some warnings + fixed build for C++03 Benjamin Golinvaux <bgo@osimis.io> parents: 1279 diff changeset	814 ASSERT_DOUBLE_EQ(-120.0, v[0]);
1167 ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	815
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	816 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "-1e-2\\2"));
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	817 ASSERT_EQ(2u, v.size());
1309 1f877e0846fe Fixed some warnings + fixed build for C++03 Benjamin Golinvaux <bgo@osimis.io> parents: 1279 diff changeset	818 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	819 ASSERT_DOUBLE_EQ(2.0, v[1]);
1167 ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	820
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	821 ASSERT_TRUE(OrthancStone::LinearAlgebra::ParseVector(v, "1.3671875\\1.3671875"));
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	822 ASSERT_EQ(2u, v.size());
1309 1f877e0846fe Fixed some warnings + fixed build for C++03 Benjamin Golinvaux <bgo@osimis.io> parents: 1279 diff changeset	823 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	824 ASSERT_DOUBLE_EQ(1.3671875, v[1]);
1167 ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	825 }
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	826
ad4e21df4e40 enriching OrthancStone::StoneInitialize() Sebastien Jodogne <s.jodogne@gmail.com> parents: 999 diff changeset	827
20 946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	828 int main(int argc, char **argv)
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	829 {
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	830 Orthanc::Logging::Initialize();
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	831 Orthanc::Logging::EnableInfoLevel(true);
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	832
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	833 ::testing::InitGoogleTest(&argc, argv);
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	834 int result = RUN_ALL_TESTS();
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	835
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	836 Orthanc::Logging::Finalize();
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	837
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	838 return result;
946377d1c992 skeleton for unit tests Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	839 }

Mercurial > hg > orthanc-stone

annotate UnitTestsSources/UnitTestsMain.cpp @ 1327:4f8db2d202c8 broker