/*
* Copyright (C) 2010, Emergya (http://www.emergya.es)
*
* @author <a href="mailto:jlrodriguez@emergya.es">Juan Luís Rodríguez</a>
* @author <a href="mailto:marias@emergya.es">María Arias</a>
*
* This file is part of GoFleet
*
* This software is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* This software is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this library; if not, write to the Free Software
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
*
* As a special exception, if you link this library with other files to
* produce an executable, this library does not by itself cause the
* resulting executable to be covered by the GNU General Public License.
* This exception does not however invalidate any other reasons why the
* executable file might be covered by the GNU General Public License.
*/
//License: GPL. For details, see LICENSE file.
package es.emergya.geo.util;
import static org.openstreetmap.josm.tools.I18n.tr;
import org.openstreetmap.josm.data.Bounds;
import org.openstreetmap.josm.data.coor.EastNorth;
import org.openstreetmap.josm.data.coor.LatLon;
/**
* Projection for the SwissGrid, see http://de.wikipedia.org/wiki/Swiss_Grid.
*
* Calculations are based on formula from
* http://www.swisstopo.admin.ch/internet/swisstopo/en/home/topics/survey/sys/refsys/switzerland.parsysrelated1.37696.downloadList.12749.DownloadFile.tmp/ch1903wgs84en.pdf
*
*/
public class SwissGrid implements Projection {
/**
* @param wgs WGS84 lat/lon (ellipsoid GRS80) (in degree)
* @return eastnorth projection in Swiss National Grid (ellipsoid Bessel)
*/
public EastNorth latlon2eastNorth(LatLon wgs) {
double phi = 3600d * wgs.lat();
double lambda = 3600d * wgs.lon();
double phiprime = (phi - 169028.66d) / 10000d;
double lambdaprime = (lambda - 26782.5d) / 10000d;
// precompute squares for lambdaprime and phiprime
//
double lambdaprime_2 = Math.pow(lambdaprime,2);
double phiprime_2 = Math.pow(phiprime,2);
double north =
200147.07d
+ 308807.95d * phiprime
+ 3745.25d * lambdaprime_2
+ 76.63d * phiprime_2
- 194.56d * lambdaprime_2 * phiprime
+ 119.79d * Math.pow(phiprime,3);
double east =
600072.37d
+ 211455.93d * lambdaprime
- 10938.51d * lambdaprime * phiprime
- 0.36d * lambdaprime * phiprime_2
- 44.54d * Math.pow(lambdaprime,3);
return new EastNorth(east, north);
}
/**
* @param xy SwissGrid east/north (in meters)
* @return LatLon WGS84 (in degree)
*/
public LatLon eastNorth2latlon(EastNorth xy) {
double yp = (xy.east() - 600000d) / 1000000d;
double xp = (xy.north() - 200000d) / 1000000d;
// precompute squares of xp and yp
//
double xp_2 = Math.pow(xp,2);
double yp_2 = Math.pow(yp,2);
// lambda = latitude, phi = longitude
double lmbp = 2.6779094d
+ 4.728982d * yp
+ 0.791484d * yp * xp
+ 0.1306d * yp * xp_2
- 0.0436d * Math.pow(yp,3);
double phip = 16.9023892d
+ 3.238272d * xp
- 0.270978d * yp_2
- 0.002528d * xp_2
- 0.0447d * yp_2 * xp
- 0.0140d * Math.pow(xp,3);
double lmb = lmbp * 100.0d / 36.0d;
double phi = phip * 100.0d / 36.0d;
return new LatLon(phi,lmb);
}
@Override public String toString() {
return tr("Swiss Grid (Switzerland)");
}
public String toCode() {
return "EPSG:21781";
}
public String getCacheDirectoryName() {
return "swissgrid";
}
public ProjectionBounds getWorldBounds()
{
Bounds b = getWorldBoundsLatLon();
return new ProjectionBounds(latlon2eastNorth(b.min), latlon2eastNorth(b.max));
}
public Bounds getWorldBoundsLatLon()
{
return new Bounds(
new LatLon(45.7, 5.7),
new LatLon(47.9, 10.6));
}
}