Subversion Repositories Projects

Rev

Rev 674 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 674 Rev 801
Line 22... Line 22...
22 
* Program URL   : http://qmapcontrol.sourceforge.net/
 22 
* Program URL   : http://qmapcontrol.sourceforge.net/
23 
*
 23 
*
24 
*/
 24 
*/
Line 25... Line 25...
25 
 
 25 
 
26 
#include "googlesatmapadapter.h"
 - 
 
27 
 
 - 
 
28 
#include <math.h>
 26 
#include "googlesatmapadapter.h"
29 
namespace qmapcontrol
 27 
namespace qmapcontrol
30 
{
 28 
{
31 
    GoogleSatMapAdapter::GoogleSatMapAdapter()
 29 
    GoogleSatMapAdapter::GoogleSatMapAdapter()
32 
            : TileMapAdapter("kh.google.com", "/kh?n=404&v=8&t=trtqtt", 256, 0, 19)
 30 
        : TileMapAdapter("khm1.google.com", "/kh/v=57&hl=de&x=%2&y=%3&z=%1&s=", 256, 0, 20)
33 
    {
 - 
 
34 
        // name = "googlesat";
 - 
 
35 
 
 - 
 
36 
        numberOfTiles = pow(2, current_zoom+0.0);
 - 
 
37 
        coord_per_x_tile = 360. / numberOfTiles;
 - 
 
38 
        coord_per_y_tile = 180. / numberOfTiles;
 31 
    {
Line 39... Line 32...
39 
    }
 32 
    }
40 
 
 33 
 
41 
    GoogleSatMapAdapter::~GoogleSatMapAdapter()
 34 
    GoogleSatMapAdapter::~GoogleSatMapAdapter()
42 
    {
 - 
 
43 
    }
 - 
 
44 
 
 - 
 
45 
    QString GoogleSatMapAdapter::getHost() const
 - 
 
46 
    {
 - 
 
47 
        int random = qrand() % 4;
 - 
 
48 
        return QString("kh%1.google.com").arg(random);
 - 
 
49 
    }
 - 
 
50 
 
 - 
 
51 
    QPoint GoogleSatMapAdapter::coordinateToDisplay(const QPointF& coordinate) const
 - 
 
52 
    {
 - 
 
53 
        //double x = ((coordinate.x()+180)*(tilesize*numberOfTiles/360));
 - 
 
54 
        //double y = (((coordinate.y()*-1)+90)*(tilesize*numberOfTiles/180));
 - 
 
55 
 
 - 
 
56 
        qreal x = (coordinate.x()+180.) * (numberOfTiles*mytilesize)/360.;              // coord to pixel!
 - 
 
57 
        //double y = -1*(coordinate.y()-90) * (numberOfTiles*tilesize)/180.;    // coord to pixel!
 - 
 
58 
        qreal y = (getMercatorYCoord(coordinate.y())-M_PI) * -1 * (numberOfTiles*mytilesize)/(2*M_PI);  // coord to pixel!
 - 
 
59 
        return QPoint(int(x), int(y));
 - 
 
60 
    }
 - 
 
61 
 
 - 
 
62 
    QPointF GoogleSatMapAdapter::displayToCoordinate(const QPoint& point) const
 - 
 
63 
    {
 - 
 
64 
        //double lon = ((point.x()/tilesize*numberOfTiles)*360)-180;
 - 
 
65 
        //double lat = (((point.y()/tilesize*numberOfTiles)*180)-90)*-1;
 - 
 
66 
 
 - 
 
67 
        qreal lon = (point.x()*(360./(numberOfTiles*mytilesize)))-180.;
 - 
 
68 
        //double lat = -(point.y()*(180./(numberOfTiles*tilesize)))+90;
 - 
 
69 
        //qreal lat = getMercatorLatitude(point.y()*-1*(2*M_PI/(numberOfTiles*tilesize)) + M_PI);
 - 
 
70 
        qreal lat = lat *180./M_PI;
 - 
 
71 
        return QPointF(lon, lat);
 - 
 
72 
    }
 - 
 
73 
 
 - 
 
74 
    qreal GoogleSatMapAdapter::getMercatorLatitude(qreal YCoord) const
 - 
 
75 
    {
 - 
 
76 
        //http://welcome.warnercnr.colostate.edu/class_info/nr502/lg4/projection_mathematics/converting.html
 - 
 
77 
        if (YCoord > M_PI) return 9999.;
 - 
 
78 
        if (YCoord < -M_PI) return -9999.;
 - 
 
79 
 
 - 
 
80 
        qreal t = atan(exp(YCoord));
 - 
 
81 
        qreal res = (2.*(t))-(M_PI/2.);
 - 
 
82 
        return res;
 - 
 
83 
    }
 - 
 
84 
 
 - 
 
85 
    qreal GoogleSatMapAdapter::getMercatorYCoord(qreal lati) const
 - 
 
86 
    {
 - 
 
87 
        qreal lat = lati;
 - 
 
88 
 
 - 
 
89 
        // conversion degre=>radians
 - 
 
90 
        qreal phi = M_PI * lat / 180;
 - 
 
91 
 
 - 
 
92 
        qreal res;
 - 
 
93 
        //double temp = Math.Tan(Math.PI / 4 - phi / 2);
 - 
 
94 
        //res = Math.Log(temp);
 - 
 
95 
        res = 0.5 * log((1 + sin(phi)) / (1 - sin(phi)));
 - 
 
96 
 
 - 
 
97 
        return res;
 - 
 
98 
    }
 - 
 
99 
 
 - 
 
100 
    void GoogleSatMapAdapter::zoom_in()
 - 
 
101 
    {
 - 
 
102 
        current_zoom+=1;
 - 
 
103 
        numberOfTiles = pow(2, current_zoom+0.0);
 - 
 
104 
        coord_per_x_tile = 360. / numberOfTiles;
 - 
 
105 
        coord_per_y_tile = 180. / numberOfTiles;
 - 
 
106 
    }
 - 
 
107 
 
 - 
 
108 
    void GoogleSatMapAdapter::zoom_out()
 - 
 
109 
    {
 - 
 
110 
        current_zoom-=1;
 - 
 
111 
        numberOfTiles = pow(2, current_zoom+0.0);
 - 
 
112 
        coord_per_x_tile = 360. / numberOfTiles;
 - 
 
113 
        coord_per_y_tile = 180. / numberOfTiles;
 - 
 
114 
    }
 - 
 
115 
 
 - 
 
116 
    bool GoogleSatMapAdapter::isValid(int x, int y, int z) const
 - 
 
117 
    {
 - 
 
118 
        if ((x>=0 && x < numberOfTiles) && (y>=0 && y < numberOfTiles) && z>=0)
 - 
 
119 
        {
 - 
 
120 
            return true;
 - 
 
121 
        }
 - 
 
122 
        return false;
 - 
 
123 
    }
 - 
 
124 
    QString GoogleSatMapAdapter::query(int i, int j, int z) const
 - 
 
125 
    {
 - 
 
126 
        return getQ(-180+i*coord_per_x_tile,
 - 
 
127 
                    90-(j+1)*coord_per_y_tile, z);
 - 
 
128 
    }
 - 
 
129 
 
 - 
 
130 
    QString GoogleSatMapAdapter::getQ(qreal longitude, qreal latitude, int zoom) const
 - 
 
131 
    {
 - 
 
132 
        qreal xmin=-180;
 - 
 
133 
        qreal xmax=180;
 - 
 
134 
        qreal ymin=-90;
 - 
 
135 
        qreal ymax=90;
 - 
 
136 
 
 - 
 
137 
        qreal xmoy=0;
 - 
 
138 
        qreal ymoy=0;
 - 
 
139 
        QString location="t";
 - 
 
140 
 
 - 
 
141 
        //Google uses a latitude divided by 2;
 - 
 
142 
        qreal halflat = latitude;
 - 
 
143 
 
 - 
 
144 
        for (int i = 0; i < zoom; i++)
 - 
 
145 
        {
 - 
 
146 
            xmoy = (xmax + xmin) / 2;
 - 
 
147 
            ymoy = (ymax + ymin) / 2;
 - 
 
148 
            if (halflat >= ymoy) //upper part (q or r)
 - 
 
149 
            {
 - 
 
150 
                ymin = ymoy;
 - 
 
151 
                if (longitude < xmoy)
 - 
 
152 
                { /*q*/
 - 
 
153 
                    location+= "q";
 - 
 
154 
                    xmax = xmoy;
 - 
 
155 
                }
 - 
 
156 
                else
 - 
 
157 
                {/*r*/
 - 
 
158 
                    location+= "r";
 - 
 
159 
                    xmin = xmoy;
 - 
 
160 
                }
 - 
 
161 
            }
 - 
 
162 
            else //lower part (t or s)
 - 
 
163 
            {
 - 
 
164 
                ymax = ymoy;
 - 
 
165 
                if (longitude < xmoy)
 - 
 
166 
                { /*t*/
 - 
 
167 
 
 - 
 
168 
                    location+= "t";
 - 
 
169 
                    xmax = xmoy;
 - 
 
170 
                }
 - 
 
171 
                else
 - 
 
172 
                {/*s*/
 - 
 
173 
                    location+= "s";
 - 
 
174 
                    xmin = xmoy;
 - 
 
175 
                }
 - 
 
176 
            }
 - 
 
177 
        }
 - 
 
178 
        return QString("/kh?n=404&v=24&t=%1").arg(location);
 35 
    {
179 
    }
 -