2 * Copyright © 2014 Keith Packard <keithp@keithp.com>
4 * This program is free software; you can redistribute it and/or modify
5 * it under the terms of the GNU General Public License as published by
6 * the Free Software Foundation; either version 2 of the License, or
7 * (at your option) any later version.
9 * This program is distributed in the hope that it will be useful, but
10 * WITHOUT ANY WARRANTY; without even the implied warranty of
11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 * General Public License for more details.
14 * You should have received a copy of the GNU General Public License along
15 * with this program; if not, write to the Free Software Foundation, Inc.,
16 * 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA.
19 package org.altusmetrum.altoslib_14;
22 import java.lang.Math;
24 import java.util.concurrent.*;
26 public class AltosMapTransform {
28 double scale_x, scale_y;
30 double offset_x, offset_y;
32 public AltosLatLon lat_lon (AltosPointDouble point) {
36 lon = point.x/scale_x;
37 rads = 2 * Math.atan(Math.exp(-point.y/scale_y));
38 lat = Math.toDegrees(rads - Math.PI/2);
40 return new AltosLatLon(lat,lon);
43 public AltosLatLon lat_lon (AltosPointInt point) {
44 return lat_lon(new AltosPointDouble(point.x, point.y));
47 public AltosPointDouble screen_point(AltosPointInt screen) {
48 return new AltosPointDouble(screen.x + offset_x, screen.y + offset_y);
51 public AltosPointDouble screen_point(AltosPointDouble screen) {
52 return new AltosPointDouble(screen.x + offset_x, screen.y + offset_y);
55 public double hypot(AltosLatLon a, AltosLatLon b) {
56 AltosPointDouble a_pt = point(a);
57 AltosPointDouble b_pt = point(b);
59 return Math.hypot(a_pt.x - b_pt.x, a_pt.y - b_pt.y);
62 public AltosLatLon screen_lat_lon(AltosPointInt screen) {
63 return lat_lon(screen_point(screen));
66 public AltosLatLon screen_lat_lon(AltosPointDouble screen) {
67 return lat_lon(screen_point(screen));
70 public AltosPointDouble point(double lat, double lon) {
76 e = Math.sin(Math.toRadians(lat));
77 e = Math.max(e,-(1-1.0E-15));
78 e = Math.min(e, 1-1.0E-15 );
80 y = 0.5*Math.log((1+e)/(1-e))*-scale_y;
82 return new AltosPointDouble(x, y);
85 public AltosPointDouble point(AltosLatLon lat_lon) {
86 return point(lat_lon.lat, lat_lon.lon);
89 public AltosPointDouble screen(AltosPointDouble point) {
90 return new AltosPointDouble(point.x - offset_x, point.y - offset_y);
93 public AltosPointInt screen(AltosPointInt point) {
94 return new AltosPointInt((int) (point.x - offset_x + 0.5),
95 (int) (point.y - offset_y + 0.5));
98 public AltosRectangle screen(AltosMapRectangle map_rect) {
99 AltosPointDouble ul = screen(map_rect.ul);
100 AltosPointDouble lr = screen(map_rect.lr);
102 return new AltosRectangle((int) ul.x, (int) ul.y, (int) (lr.x - ul.x), (int) (lr.y - ul.y));
105 public AltosPointDouble screen(AltosLatLon lat_lon) {
106 return screen(point(lat_lon));
109 public AltosPointDouble screen(double lat, double lon) {
110 return screen(point(lat, lon));
113 private boolean has_location;
115 public boolean has_location() {
119 public AltosMapTransform(int width, int height, int zoom, AltosLatLon centre_lat_lon) {
120 scale_x = 256/360.0 * Math.pow(2, zoom);
121 scale_y = 256/(2.0*Math.PI) * Math.pow(2, zoom);
123 AltosPointDouble centre_pt = point(centre_lat_lon);
125 has_location = (centre_lat_lon.lat != 0 || centre_lat_lon.lon != 0);
126 offset_x = centre_pt.x - width / 2.0;
127 offset_y = centre_pt.y - height / 2.0;
130 public static double lon_from_distance(double lat, double distance) {
131 double c = AltosGreatCircle.earth_radius * Math.cos(lat * Math.PI / 180) * 2 * Math.PI;
135 return distance/c * 360.0;