1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
|
/*LICENSE_START*/
/*
* Copyright 1995-2002 Washington University School of Medicine
*
* http://brainmap.wustl.edu
*
* This file is part of CARET.
*
* CARET is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* CARET is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with CARET; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*
*/
/*LICENSE_END*/
#include "BorderProjectionUnprojector.h"
/**
* Unproject a border projection point
*/
void
BorderProjectionUnprojector::unprojectBorderProjectionLink(const BorderProjectionLink* bl,
const CoordinateFile& cf,
float xyz[3], int& section, float& radius)
{
section = 0;
xyz[0] = 0.0;
xyz[1] = 0.0;
xyz[2] = 0.0;
radius = 0.0;
int vertices[3];
float areas[3];
bl->getData(section, vertices, areas, radius);
const float totalArea = areas[0] + areas[1] + areas[2];
if (totalArea > 0.0) {
const float* v1 = cf.getCoordinate(vertices[0]);
const float* v2 = cf.getCoordinate(vertices[1]);
const float* v3 = cf.getCoordinate(vertices[2]);
xyz[0] = (v1[0] * areas[1] + v2[0] * areas[2] + v3[0] * areas[0])
/ totalArea;
xyz[1] = (v1[1] * areas[1] + v2[1] * areas[2] + v3[1] * areas[0])
/ totalArea;
xyz[2] = (v1[2] * areas[1] + v2[2] * areas[2] + v3[2] * areas[0])
/ totalArea;
/*
xyz[0] = (v1[0] * areas[2] + v2[0] * areas[0] + v3[0] * areas[1])
/ totalArea;
xyz[1] = (v1[1] * areas[2] + v2[1] * areas[0] + v3[1] * areas[1])
/ totalArea;
xyz[2] = (v1[2] * areas[2] + v2[2] * areas[0] + v3[2] * areas[1])
/ totalArea;
*/
}
}
/**
* Constructor
*/
BorderProjectionUnprojector::BorderProjectionUnprojector()
{
}
/**
* Destructor
*/
BorderProjectionUnprojector::~BorderProjectionUnprojector()
{
}
/**
* Unproject a border projection file into surfaces projected borders
*/
void
BorderProjectionUnprojector::unprojectBorderProjections(const CoordinateFile& cf,
const BorderProjectionFile& bpf,
BorderFile& bf,
const int startAtProjection)
{
const int numProj = bpf.getNumberOfBorderProjections();
for (int i = startAtProjection; i < numProj; i++) {
const BorderProjection* bp = bpf.getBorderProjection(i);
const int numLinks = bp->getNumberOfLinks();
QString name;
float center[3];
float sampling, variance, topography, uncertainty;
bp->getData(name, center, sampling, variance, topography, uncertainty);
Border b(name, center, sampling, variance, topography, uncertainty);
b.setBorderColorIndex(bp->getBorderColorIndex());
for (int j = 0; j < numLinks; j++) {
const BorderProjectionLink* bpl = bp->getBorderProjectionLink(j);
int section;
float xyz[3];
float radius;
unprojectBorderProjectionLink(bpl, cf, xyz, section, radius);
b.addBorderLink(xyz, section, radius);
}
b.setBorderProjectionID(bp->getUniqueID());
b.setBorderColorIndex(bp->getBorderColorIndex());
bf.addBorder(b);
}
}
|
