Index: src/org/openstreetmap/josm/data/projection/Projections.java IDEA additional info: Subsystem: com.intellij.openapi.diff.impl.patch.CharsetEP <+>UTF-8 =================================================================== --- src/org/openstreetmap/josm/data/projection/Projections.java (date 1456352671000) +++ src/org/openstreetmap/josm/data/projection/Projections.java (revision ) @@ -40,6 +40,7 @@ import org.openstreetmap.josm.data.projection.proj.PolarStereographic; import org.openstreetmap.josm.data.projection.proj.Proj; import org.openstreetmap.josm.data.projection.proj.ProjFactory; +import org.openstreetmap.josm.data.projection.proj.Sinusoidal; import org.openstreetmap.josm.data.projection.proj.SwissObliqueMercator; import org.openstreetmap.josm.data.projection.proj.TransverseMercator; import org.openstreetmap.josm.gui.preferences.projection.ProjectionChoice; @@ -94,6 +95,7 @@ registerBaseProjection("merc", Mercator.class, "core"); registerBaseProjection("omerc", ObliqueMercator.class, "core"); registerBaseProjection("somerc", SwissObliqueMercator.class, "core"); + registerBaseProjection("sinu", Sinusoidal.class, "core"); registerBaseProjection("stere", PolarStereographic.class, "core"); registerBaseProjection("sterea", DoubleStereographic.class, "core"); registerBaseProjection("tmerc", TransverseMercator.class, "core"); Index: data_nodist/projection/josm-epsg IDEA additional info: Subsystem: com.intellij.openapi.diff.impl.patch.CharsetEP <+>UTF-8 =================================================================== --- data_nodist/projection/josm-epsg (date 1456352671000) +++ data_nodist/projection/josm-epsg (revision ) @@ -417,3 +417,8 @@ <6811> +proj=omerc +lat_0=45.91666666666666 +lonc=-123 +alpha=295 +k=1 +x_0=7000000.00000248 +y_0=-2999999.999988 +no_uoff +gamma=295 +ellps=GRS80 +units=ft +no_defs +bounds=-125,45,-121,47 <> # NAD83(2011) / Oregon Columbia River West zone (m) <6810> +proj=omerc +lat_0=45.91666666666666 +lonc=-123 +alpha=295 +k=1 +x_0=7000000 +y_0=-3000000 +no_uoff +gamma=295 +ellps=GRS80 +units=m +no_defs +bounds=-125,45,-121,47 <> +## +## Following entries use Sinusoidal projection +## +# ESRI:54008 / World Sinusoidal +<54008> +proj=sinu +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs <> \ No newline at end of file Index: src/org/openstreetmap/josm/tools/Geometry.java IDEA additional info: Subsystem: com.intellij.openapi.diff.impl.patch.CharsetEP <+>UTF-8 =================================================================== --- src/org/openstreetmap/josm/tools/Geometry.java (date 1456352671000) +++ src/org/openstreetmap/josm/tools/Geometry.java (revision ) @@ -30,6 +30,8 @@ import org.openstreetmap.josm.data.osm.Relation; import org.openstreetmap.josm.data.osm.RelationMember; import org.openstreetmap.josm.data.osm.Way; +import org.openstreetmap.josm.data.projection.Projection; +import org.openstreetmap.josm.data.projection.Projections; /** * Some tools for geometry related tasks. @@ -629,64 +631,14 @@ /** * Returns area of a closed way in square meters. - * (approximate(?), but should be OK for small areas) * - * Relies on the current projection: Works correctly, when - * one unit in projected coordinates corresponds to one meter. - * This is true for most projections, but not for WGS84 and - * Mercator (EPSG:3857). - * * @param way Way to measure, should be closed (first node is the same as last node) * @return area of the closed way. */ public static double closedWayArea(Way way) { - - //http://local.wasp.uwa.edu.au/~pbourke/geometry/polyarea/ - double area = 0; - Node lastN = null; - for (Node n : way.getNodes()) { - if (lastN != null) { - area += (calcX(n) * calcY(lastN)) - (calcY(n) * calcX(lastN)); + return getAreaAndPerimeter(way.getNodes(), Projections.getProjectionByCode("EPSG:54008")).getArea(); - } + } - lastN = n; - } - return Math.abs(area/2); - } - protected static double calcX(Node p1) { - double lat1, lon1, lat2, lon2; - double dlon, dlat; - - lat1 = p1.getCoor().lat() * Math.PI / 180.0; - lon1 = p1.getCoor().lon() * Math.PI / 180.0; - lat2 = lat1; - lon2 = 0; - - dlon = lon2 - lon1; - dlat = lat2 - lat1; - - double a = Math.pow(Math.sin(dlat/2), 2) + Math.cos(lat1) * Math.cos(lat2) * Math.pow(Math.sin(dlon/2), 2); - double c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a)); - return 6367000 * c; - } - - protected static double calcY(Node p1) { - double lat1, lon1, lat2, lon2; - double dlon, dlat; - - lat1 = p1.getCoor().lat() * Math.PI / 180.0; - lon1 = p1.getCoor().lon() * Math.PI / 180.0; - lat2 = 0; - lon2 = lon1; - - dlon = lon2 - lon1; - dlat = lat2 - lat1; - - double a = Math.pow(Math.sin(dlat/2), 2) + Math.cos(lat1) * Math.cos(lat2) * Math.pow(Math.sin(dlon/2), 2); - double c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a)); - return 6367000 * c; - } - /** * Determines whether a way is oriented clockwise. * @@ -985,14 +937,27 @@ * @return area and perimeter */ public static AreaAndPerimeter getAreaAndPerimeter(List nodes) { + return getAreaAndPerimeter(nodes, Main.getProjection()); + } + + /** + * Calculate area and perimeter length of a polygon in the given projection. + * + * @param nodes the list of nodes representing the polygon + * @param projection the projection to use for the calculation + * @return area and perimeter + */ + public static AreaAndPerimeter getAreaAndPerimeter(List nodes, Projection projection) { + CheckParameterUtil.ensureParameterNotNull(nodes, "nodes"); + CheckParameterUtil.ensureParameterNotNull(projection, "projection"); double area = 0; double perimeter = 0; if (!nodes.isEmpty()) { boolean closed = nodes.get(0) == nodes.get(nodes.size() - 1); int numSegments = closed ? nodes.size() - 1 : nodes.size(); - EastNorth p1 = nodes.get(0).getEastNorth(); + EastNorth p1 = projection.latlon2eastNorth(nodes.get(0).getCoor()); for (int i = 1; i <= numSegments; i++) { - EastNorth p2 = nodes.get(i == numSegments ? 0 : i).getEastNorth(); + final EastNorth p2 = projection.latlon2eastNorth(nodes.get(i == numSegments ? 0 : i).getCoor()); area += p1.east() * p2.north() - p2.east() * p1.north(); perimeter += p1.distance(p2); p1 = p2; Index: src/org/openstreetmap/josm/data/projection/proj/Sinusoidal.java IDEA additional info: Subsystem: com.intellij.openapi.diff.impl.patch.CharsetEP <+>UTF-8 =================================================================== --- src/org/openstreetmap/josm/data/projection/proj/Sinusoidal.java (revision ) +++ src/org/openstreetmap/josm/data/projection/proj/Sinusoidal.java (revision ) @@ -0,0 +1,37 @@ +// License: GPL. For details, see LICENSE file. +package org.openstreetmap.josm.data.projection.proj; + +import org.openstreetmap.josm.data.Bounds; + +import static org.openstreetmap.josm.tools.I18n.tr; + +/** + * Sinusoidal projection (aka. Sanson–Flamsteed, Mercator equal-area projection) + */ +public class Sinusoidal extends AbstractProj { + + @Override + public String getName() { + return tr("Sinusoidal"); + } + + @Override + public String getProj4Id() { + return "sinu"; + } + + @Override + public double[] project(double lat_rad, double lon_rad) { + return new double[]{lon_rad * Math.cos(lat_rad), lat_rad}; + } + + @Override + public double[] invproject(double east, double north) { + return new double[]{east / Math.cos(north), north}; + } + + @Override + public Bounds getAlgorithmBounds() { + return new Bounds(-90, -180, 90, 180, false); + } +}