Subversion Repositories splitter

Rev

Rev 529 | Blame | Compare with Previous | Last modification | View Log | RSS feed

/*
 * Copyright (c) 2009, Chris Miller
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License version 3 as
 * published by the Free Software Foundation.
 *
 * This program 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.
 */


package uk.me.parabola.splitter;

/**
 * Utility methods for rounding numbers and areas
 *
 * @author Chris Miller
 */

public class RoundingUtils {
        /**
         * Rounds an integer down to the nearest multiple of {@code 2^shift}.
         * Works with both positive and negative integers.
         * @param val the integer to round down.
         * @param shift the power of two to round down to.
         * @return the rounded integer.
         */

        public static int roundDown(int val, int shift) {
                return val >>> shift << shift;
        }

        /**
         * Rounds an integer up to the nearest multiple of {@code 2^shift}.
         * Works with both positive and negative integers.
         * @param val the integer to round up.
         * @param shift the power of two to round up to.
         * @return the rounded integer.
         */

        public static int roundUp(int val, int shift) {
                return (val + (1 << shift) - 1) >>> shift << shift;
        }

        /**
         * Rounds an integer up or down to the nearest multiple of {@code 2^shift}.
         * Works with both positive and negative integers.
         * @param val the integer to round.
         * @param shift the power of two to round to.
         * @return the rounded integer.
         */

        public static int round(int val, int shift) {
                return (val + (1 << (shift - 1))) >>> shift << shift;
        }

        /**
         * Rounds an area's borders to suit the supplied resolution. This
         * means edges are aligned at 2 ^ (24 - resolution) boundaries
         *
         * @param b the area to round
         * @param resolution the map resolution to align the borders at
         * @return the rounded area
         */

        public static Area round(Area b, int resolution) {
                int shift = 24 - resolution;
                int alignment = 1 << shift;

                // Avoid pathological behaviour near the poles by discarding anything
                // greater than +/-85 degrees latitude.
                int minLat = Math.max(b.getMinLat(), Utils.toMapUnit(-85.0d));
                int maxLat = Math.min(b.getMaxLat(), Utils.toMapUnit(85.0d));

                int roundedMinLat = roundDown(minLat, shift);
                int roundedMaxLat = roundUp(maxLat, shift);
                assert roundedMinLat % alignment == 0 : "The area's min latitude is not aligned to a multiple of " + alignment;
                assert roundedMaxLat % alignment == 0 : "The area's max latitude is not aligned to a multiple of " + alignment;

                int roundedMinLon = roundDown(b.getMinLong(), shift);
                int roundedMaxLon = roundUp(b.getMaxLong(), shift);
                // don't produce illegal values
                if (roundedMinLon < -0x800000)
                        roundedMinLon = -0x800000;
                if (roundedMaxLon > 0x800000)
                        roundedMaxLon = 0x800000;
                assert roundedMinLon % alignment == 0 : "The area's min longitude is not aligned to a multiple of " + alignment;
                assert roundedMaxLon % alignment == 0 : "The area's max longitude is not aligned to a multiple of " + alignment;

                return new Area(roundedMinLat, roundedMinLon, roundedMaxLat, roundedMaxLon);
        }
}