git.gag.com Git - fw/altos/blob - altoslib/AltosMapTransform.java

   1 /*
   2  * Copyright © 2014 Keith Packard <keithp@keithp.com>
   3  *
   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; version 2 of the License.
   7  *
   8  * This program is distributed in the hope that it will be useful, but
   9  * WITHOUT ANY WARRANTY; without even the implied warranty of
  10  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11  * General Public License for more details.
  12  *
  13  * You should have received a copy of the GNU General Public License along
  14  * with this program; if not, write to the Free Software Foundation, Inc.,
  15  * 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA.
  16  */
  17
  18 package org.altusmetrum.altoslib_6;
  19
  20 import java.io.*;
  21 import java.lang.Math;
  22 import java.util.*;
  23 import java.util.concurrent.*;
  24
  25 public class AltosMapTransform {
  26
  27         double  scale_x, scale_y;
  28
  29         double  offset_x, offset_y;
  30
  31         public AltosLatLon lat_lon (AltosPointDouble point) {
  32                 double lat, lon;
  33                 double rads;
  34
  35                 lon = point.x/scale_x;
  36                 rads = 2 * Math.atan(Math.exp(-point.y/scale_y));
  37                 lat = Math.toDegrees(rads - Math.PI/2);
  38
  39                 return new AltosLatLon(lat,lon);
  40         }
  41
  42         public AltosPointDouble screen_point(AltosPointDouble screen) {
  43                 return new AltosPointDouble(screen.x + offset_x, screen.y + offset_y);
  44         }
  45
  46         public AltosLatLon screen_lat_lon(AltosPointDouble screen) {
  47                 return lat_lon(screen_point(screen));
  48         }
  49
  50         public AltosPointDouble point(AltosLatLon lat_lon) {
  51                 double x, y;
  52                 double e;
  53
  54                 x = lat_lon.lon * scale_x;
  55
  56                 e = Math.sin(Math.toRadians(lat_lon.lat));
  57                 e = Math.max(e,-(1-1.0E-15));
  58                 e = Math.min(e,  1-1.0E-15 );
  59
  60                 y = 0.5*Math.log((1+e)/(1-e))*-scale_y;
  61
  62                 return new AltosPointDouble(x, y);
  63         }
  64
  65         public AltosPointDouble screen(AltosPointDouble point) {
  66                 return new AltosPointDouble(point.x - offset_x, point.y - offset_y);
  67         }
  68
  69         public AltosPointInt screen(AltosPointInt point) {
  70                 return new AltosPointInt((int) (point.x - offset_x + 0.5),
  71                                          (int) (point.y - offset_y + 0.5));
  72         }
  73
  74 //      public Rectangle screen(AltosMapRectangle map_rect) {
  75 //              AltosPoint2D    ul = screen(map_rect.ul);
  76 //              AltosPoint2D    lr = screen(map_rect.lr);
  77 //
  78 //              return new Rectangle((int) ul.x, (int) ul.y, (int) (lr.x - ul.x), (int) (lr.y - ul.y));
  79 //      }
  80
  81         public AltosPointDouble screen(AltosLatLon lat_lon) {
  82                 return screen(point(lat_lon));
  83         }
  84
  85         private boolean has_location;
  86
  87         public boolean has_location() {
  88                 return has_location;
  89         }
  90
  91         public AltosMapTransform(int width, int height, int zoom, AltosLatLon centre_lat_lon) {
  92                 scale_x = 256/360.0 * Math.pow(2, zoom);
  93                 scale_y = 256/(2.0*Math.PI) * Math.pow(2, zoom);
  94
  95                 AltosPointDouble centre_pt = point(centre_lat_lon);
  96
  97                 has_location = (centre_lat_lon.lat != 0 || centre_lat_lon.lon != 0);
  98                 offset_x = centre_pt.x - width / 2.0;
  99                 offset_y = centre_pt.y - height / 2.0;
 100         }
 101 }