1 | // License: GPL. Copyright 2007 by Immanuel Scholz and others
|
---|
2 | package org.openstreetmap.josm.data.projection;
|
---|
3 |
|
---|
4 | import static org.openstreetmap.josm.tools.I18n.tr;
|
---|
5 |
|
---|
6 | import org.openstreetmap.josm.data.Bounds;
|
---|
7 | import org.openstreetmap.josm.data.coor.EastNorth;
|
---|
8 | import org.openstreetmap.josm.data.coor.LatLon;
|
---|
9 |
|
---|
10 | /**
|
---|
11 | * Implement Mercator Projection code, coded after documentation
|
---|
12 | * from wikipedia.
|
---|
13 | *
|
---|
14 | * The center of the mercator projection is always the 0 grad
|
---|
15 | * coordinate.
|
---|
16 | *
|
---|
17 | * See also USGS Bulletin 1532
|
---|
18 | * (http://egsc.usgs.gov/isb/pubs/factsheets/fs08799.html)
|
---|
19 | *
|
---|
20 | * @author imi
|
---|
21 | */
|
---|
22 | public class Mercator implements Projection {
|
---|
23 |
|
---|
24 | final double radius = 6378137.0;
|
---|
25 |
|
---|
26 | public EastNorth latlon2eastNorth(LatLon p) {
|
---|
27 | return new EastNorth(
|
---|
28 | p.lon()*Math.PI/180*radius,
|
---|
29 | Math.log(Math.tan(Math.PI/4+p.lat()*Math.PI/360))*radius);
|
---|
30 | }
|
---|
31 |
|
---|
32 | public LatLon eastNorth2latlon(EastNorth p) {
|
---|
33 | return new LatLon(
|
---|
34 | Math.atan(Math.sinh(p.north()/radius))*180/Math.PI,
|
---|
35 | p.east()/radius*180/Math.PI);
|
---|
36 | }
|
---|
37 |
|
---|
38 | @Override public String toString() {
|
---|
39 | return tr("Mercator");
|
---|
40 | }
|
---|
41 |
|
---|
42 | public String toCode() {
|
---|
43 | return "EPSG:3857"; /* initially they used 3785 but that has been superseded, see http://www.epsg-registry.org/ */
|
---|
44 | }
|
---|
45 |
|
---|
46 | @Override
|
---|
47 | public int hashCode() {
|
---|
48 | return getClass().getName().hashCode(); // we have no variables
|
---|
49 | }
|
---|
50 |
|
---|
51 | public String getCacheDirectoryName() {
|
---|
52 | return "mercator";
|
---|
53 | }
|
---|
54 |
|
---|
55 | public Bounds getWorldBoundsLatLon()
|
---|
56 | {
|
---|
57 | return new Bounds(
|
---|
58 | new LatLon(-85.05112877980659, -180.0),
|
---|
59 | new LatLon(85.05112877980659, 180.0));
|
---|
60 | }
|
---|
61 |
|
---|
62 | public double getDefaultZoomInPPD() {
|
---|
63 | // This will set the scale bar to about 100 km
|
---|
64 | return 1000.0;/*0.000158*/
|
---|
65 | }
|
---|
66 | }
|
---|