/* * Licensed to the Apache Software Foundation (ASF) under one or more * contributor license agreements. See the NOTICE file distributed with * this work for additional information regarding copyright ownership. * The ASF licenses this file to You under the Apache License, Version 2.0 * (the "License"); you may not use this file except in compliance with * the License. You may obtain a copy of the License at * * http://www.apache.org/licenses/LICENSE-2.0 * * Unless required by applicable law or agreed to in writing, software * distributed under the License is distributed on an "AS IS" BASIS, * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. * See the License for the specific language governing permissions and * limitations under the License. */ package org.apache.lucene.spatial3d.geom; /** * Degenerate longitude slice. * * @lucene.internal */ class GeoDegenerateLongitudeSlice extends GeoBaseBBox { /** The longitude of the slice */ protected final double longitude; /** The bounding plane for the slice (through both poles, perpendicular to the slice) */ protected final SidedPlane boundingPlane; /** The plane of the slice */ protected final Plane plane; /** A point on the slice */ protected final GeoPoint interiorPoint; /** An array consisting of the one point chosen on the slice */ protected final GeoPoint[] edgePoints; /** Notable points for the slice (north and south poles) */ protected final GeoPoint[] planePoints; /** * Accepts only values in the following ranges: lon: {@code -PI -> PI} */ public GeoDegenerateLongitudeSlice(final PlanetModel planetModel, final double longitude) { super(planetModel); // Argument checking if (longitude < -Math.PI || longitude > Math.PI) throw new IllegalArgumentException("Longitude out of range"); this.longitude = longitude; final double sinLongitude = Math.sin(longitude); final double cosLongitude = Math.cos(longitude); this.plane = new Plane(cosLongitude, sinLongitude); // We need a bounding plane too, which is perpendicular to the longitude plane and sided so that the point (0.0, longitude) is inside. this.interiorPoint = new GeoPoint(planetModel, 0.0, sinLongitude, 1.0, cosLongitude); this.boundingPlane = new SidedPlane(interiorPoint, -sinLongitude, cosLongitude); this.edgePoints = new GeoPoint[]{interiorPoint}; this.planePoints = new GeoPoint[]{planetModel.NORTH_POLE, planetModel.SOUTH_POLE}; } @Override public GeoBBox expand(final double angle) { // Figuring out when we escalate to a special case requires some prefiguring double newLeftLon = longitude - angle; double newRightLon = longitude + angle; double currentLonSpan = 2.0 * angle; if (currentLonSpan + 2.0 * angle >= Math.PI * 2.0) { newLeftLon = -Math.PI; newRightLon = Math.PI; } return GeoBBoxFactory.makeGeoBBox(planetModel, Math.PI * 0.5, -Math.PI * 0.5, newLeftLon, newRightLon); } @Override public boolean isWithin(final double x, final double y, final double z) { return plane.evaluateIsZero(x, y, z) && boundingPlane.isWithin(x, y, z); } @Override public double getRadius() { return Math.PI * 0.5; } @Override public GeoPoint getCenter() { return interiorPoint; } @Override public GeoPoint[] getEdgePoints() { return edgePoints; } @Override public boolean intersects(final Plane p, final GeoPoint[] notablePoints, final Membership... bounds) { return p.intersects(planetModel, plane, notablePoints, planePoints, bounds, boundingPlane); } @Override public void getBounds(Bounds bounds) { super.getBounds(bounds); bounds .addVerticalPlane(planetModel, longitude, plane, boundingPlane) .addPoint(planetModel.NORTH_POLE).addPoint(planetModel.SOUTH_POLE); } @Override public int getRelationship(final GeoShape path) { // Look for intersections. if (path.intersects(plane, planePoints, boundingPlane)) return OVERLAPS; if (path.isWithin(interiorPoint)) return CONTAINS; return DISJOINT; } @Override protected double outsideDistance(final DistanceStyle distanceStyle, final double x, final double y, final double z) { final double distance = distanceStyle.computeDistance(planetModel, plane, x,y,z, boundingPlane); final double northDistance = distanceStyle.computeDistance(planetModel.NORTH_POLE, x,y,z); final double southDistance = distanceStyle.computeDistance(planetModel.SOUTH_POLE, x,y,z); return Math.min( distance, Math.min(northDistance, southDistance)); } @Override public boolean equals(Object o) { if (!(o instanceof GeoDegenerateLongitudeSlice)) return false; GeoDegenerateLongitudeSlice other = (GeoDegenerateLongitudeSlice) o; return super.equals(other) && other.longitude == longitude; } @Override public int hashCode() { int result = super.hashCode(); long temp = Double.doubleToLongBits(longitude); result = result * 31 + (int) (temp ^ (temp >>> 32)); return result; } @Override public String toString() { return "GeoDegenerateLongitudeSlice: {planetmodel="+planetModel+", longitude=" + longitude + "(" + longitude * 180.0 / Math.PI + ")}"; } }