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
|
/* reprojection.h
*
* Convert OSM lattitude / longitude from degrees to mercator
* so that Mapnik does not have to project the data again
*
*/
#ifndef REPROJECTION_H
#define REPROJECTION_H
#include <boost/noncopyable.hpp>
#include <osmium/geom/projection.hpp>
#include <osmium/osm/location.hpp>
enum Projection { PROJ_LATLONG = 4326, PROJ_SPHERE_MERC = 3857 };
class reprojection : public boost::noncopyable
{
public:
virtual ~reprojection() = default;
/**
* Reproject from the source projection lat/lon (EPSG:4326)
* to target projection.
*/
virtual osmium::geom::Coordinates reproject(osmium::Location loc) const = 0;
/**
* Converts coordinates from target projection to tile projection (EPSG:3857)
*
* Do not confuse with coords_to_tile which does *not* calculate coordinates in the
* tile projection, but tile coordinates.
*/
virtual void target_to_tile(double *lat, double *lon) const = 0;
/**
* Converts from target coordinates to tile coordinates.
*
* The zoom level for the coordinates is explicitly given in the
* variable map_width.
*/
void coords_to_tile(double *tilex, double *tiley,
double lon, double lat, int map_width);
virtual int target_srs() const = 0;
virtual const char *target_desc() const = 0;
bool target_latlon() const
{
return target_srs() == PROJ_LATLONG;
}
/**
* Create a reprojection object with target srs `srs`.
*
* The target projection (used in the PostGIS tables).
* Controlled by the -l/-m/-E options.
*/
static reprojection *create_projection(int srs);
};
#endif
|
