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
|
/****************************************************************************
**
** Copyright (C) 2012 Nokia Corporation and/or its subsidiary(-ies).
** All rights reserved.
** Contact: http://www.qt-project.org/
**
** This file is part of the QtLocation module of the Qt Toolkit.
**
** $QT_BEGIN_LICENSE:LGPL$
** GNU Lesser General Public License Usage
** This file may be used under the terms of the GNU Lesser General Public
** License version 2.1 as published by the Free Software Foundation and
** appearing in the file LICENSE.LGPL included in the packaging of this
** file. Please review the following information to ensure the GNU Lesser
** General Public License version 2.1 requirements will be met:
** http://www.gnu.org/licenses/old-licenses/lgpl-2.1.html.
**
** In addition, as a special exception, Nokia gives you certain additional
** rights. These rights are described in the Nokia Qt LGPL Exception
** version 1.1, included in the file LGPL_EXCEPTION.txt in this package.
**
** GNU General Public License Usage
** Alternatively, this file may be used under the terms of the GNU General
** Public License version 3.0 as published by the Free Software Foundation
** and appearing in the file LICENSE.GPL included in the packaging of this
** file. Please review the following information to ensure the GNU General
** Public License version 3.0 requirements will be met:
** http://www.gnu.org/copyleft/gpl.html.
**
** Other Usage
** Alternatively, this file may be used in accordance with the terms and
** conditions contained in a signed written agreement between you and Nokia.
**
**
**
**
**
** $QT_END_LICENSE$
**
****************************************************************************/
#include "qgeoprojection_p.h"
#include "qgeocoordinate.h"
#include <QMatrix4x4>
#include <qnumeric.h>
#include <cmath>
#include "qdoublevector2d_p.h"
#include "qdoublevector3d_p.h"
QT_BEGIN_NAMESPACE
QGeoProjection::QGeoProjection() {}
QGeoProjection::~QGeoProjection() {}
QDoubleVector3D QGeoProjection::mercatorToPoint(const QDoubleVector2D &mercator) const
{
return this->coordToPoint(mercatorToCoord(mercator));
}
QDoubleVector2D QGeoProjection::pointToMercator(const QDoubleVector3D &point) const
{
return coordToMercator(this->pointToCoord(point));
}
QDoubleVector2D QGeoProjection::coordToMercator(const QGeoCoordinate &coord) const
{
const double pi = M_PI;
double lon = coord.longitude() / 360.0 + 0.5;
double lat = coord.latitude();
lat = 0.5 - (log(tan((pi / 4.0) + (pi / 2.0) * lat / 180.0)) / pi) / 2.0;
lat = qMax(0.0, lat);
lat = qMin(1.0, lat);
return QDoubleVector2D(lon, lat);
}
double QGeoProjection::realmod(const double a, const double b)
{
quint64 div = static_cast<quint64>(a / b);
return a - static_cast<double>(div) * b;
}
QGeoCoordinate QGeoProjection::mercatorToCoord(const QDoubleVector2D &mercator) const
{
const double pi = M_PI;
double fx = mercator.x();
double fy = mercator.y();
if (fy < 0.0)
fy = 0.0;
else if (fy > 1.0)
fy = 1.0;
double lat;
if (fy == 0.0)
lat = 90.0;
else if (fy == 1.0)
lat = -90.0;
else
lat = (180.0 / pi) * (2.0 * atan(exp(pi * (1.0 - 2.0 * fy))) - (pi / 2.0));
double lng;
if (fx >= 0) {
lng = realmod(fx, 1.0);
} else {
lng = realmod(1.0 - realmod(-1.0 * fx, 1.0), 1.0);
}
lng = lng * 360.0 - 180.0;
return QGeoCoordinate(lat, lng, 0.0);
}
QT_END_NAMESPACE
|
