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
|
#ifndef SEEN_TRANSF_MAT_3x4_H
#define SEEN_TRANSF_MAT_3x4_H
/*
* 3x4 transformation matrix to map points from projective 3-space into the projective plane
*
* Authors:
* Maximilian Albert <Anhalter42@gmx.de>
*
* Copyright (C) 2007 Authors
*
* Released under GNU GPL, read the file 'COPYING' for more information
*/
#include "proj_pt.h"
#include "axis-manip.h"
namespace Proj {
class TransfMat3x4 {
public:
TransfMat3x4();
TransfMat3x4(Pt2 vp_x, Pt2 vp_y, Pt2 vp_z, Pt2 origin);
TransfMat3x4(TransfMat3x4 const &rhs);
Pt2 column (Proj::Axis axis) const;
Pt2 image (Pt3 const &point);
Pt3 preimage (Geom::Point const &pt, double coord = 0, Axis = Z);
void set_image_pt (Proj::Axis axis, Proj::Pt2 const &pt);
void toggle_finite (Proj::Axis axis);
double get_infinite_angle (Proj::Axis axis) {
if (has_finite_image(axis)) {
return Geom::infinity();
}
Pt2 vp(column(axis));
return Geom::atan2(Geom::Point(vp[0], vp[1])) * 180.0/M_PI;
}
void set_infinite_direction (Proj::Axis axis, double angle) { // angle is in degrees
g_return_if_fail(tmat[2][axis] == 0); // don't set directions for finite VPs
double a = angle * M_PI/180;
Geom::Point pt(tmat[0][axis], tmat[1][axis]);
double rad = Geom::L2(pt);
set_image_pt(axis, Proj::Pt2(cos (a) * rad, sin (a) * rad, 0.0));
}
inline bool has_finite_image (Proj::Axis axis) { return (tmat[2][axis] != 0.0); }
gchar * pt_to_str (Proj::Axis axis);
bool operator==(const TransfMat3x4 &rhs) const;
TransfMat3x4 operator*(Geom::Affine const &A) const;
TransfMat3x4 &operator*=(Geom::Affine const &A);
void print() const;
void copy_tmat(double rhs[3][4]);
private:
// FIXME: Is changing a single column allowed when a projective coordinate system is specified!?!?!
void normalize_column (Proj::Axis axis);
inline void set_column (Proj::Axis axis, Proj::Pt2 pt) {
tmat[0][axis] = pt[0];
tmat[1][axis] = pt[1];
tmat[2][axis] = pt[2];
}
double tmat[3][4];
};
} // namespace Proj
#endif /* __TRANSF_MAT_3x4_H__ */
/*
Local Variables:
mode:c++
c-file-style:"stroustrup"
c-file-offsets:((innamespace . 0)(inline-open . 0)(case-label . +))
indent-tabs-mode:nil
fill-column:99
End:
*/
// vim: filetype=cpp:expandtab:shiftwidth=4:tabstop=8:softtabstop=4:fileencoding=utf-8:textwidth=99 :
|