/*
 * 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;

import java.io.InputStream;
import java.io.OutputStream;
import java.io.IOException;

Degenerate bounding box wider than PI and limited on two sides (left lon, right lon).
@lucene.internal
/** * Degenerate bounding box wider than PI and limited on two sides (left lon, right lon). * * @lucene.internal */
class GeoWideDegenerateHorizontalLine extends GeoBaseBBox {
The latitude of the line
/** The latitude of the line */
protected final double latitude;
The left longitude cutoff of the line
/** The left longitude cutoff of the line */
protected final double leftLon;
The right longitude cutoff of the line
/** The right longitude cutoff of the line */
protected final double rightLon;
The left end of the line
/** The left end of the line */
protected final GeoPoint LHC;
The right end of the line
/** The right end of the line */
protected final GeoPoint RHC;
The plane the line is in
/** The plane the line is in */
protected final Plane plane;
The left cutoff plane
/** The left cutoff plane */
protected final SidedPlane leftPlane;
The right cutoff plane
/** The right cutoff plane */
protected final SidedPlane rightPlane;
Notable points for the line
/** Notable points for the line */
protected final GeoPoint[] planePoints;
Center point for the line
/** Center point for the line */
protected final GeoPoint centerPoint;
Left/right combination bound
/** Left/right combination bound */
protected final EitherBound eitherBound;
A point on the line
/** A point on the line */
protected final GeoPoint[] edgePoints;
Accepts only values in the following ranges: lat: -PI/2 -> PI/2, lon: -PI -> PI. Horizontal angle must be greater than or equal to PI.
Params:
  • planetModel – is the planet model.
  • latitude – is the line latitude.
  • leftLon – is the left cutoff longitude.
  • rightLon – is the right cutoff longitude.
/** * Accepts only values in the following ranges: lat: {@code -PI/2 -> PI/2}, lon: {@code -PI -> PI}. * Horizontal angle must be greater than or equal to PI. *@param planetModel is the planet model. *@param latitude is the line latitude. *@param leftLon is the left cutoff longitude. *@param rightLon is the right cutoff longitude. */
public GeoWideDegenerateHorizontalLine(final PlanetModel planetModel, final double latitude, final double leftLon, double rightLon) { super(planetModel); // Argument checking if (latitude > Math.PI * 0.5 || latitude < -Math.PI * 0.5) throw new IllegalArgumentException("Latitude out of range"); if (leftLon < -Math.PI || leftLon > Math.PI) throw new IllegalArgumentException("Left longitude out of range"); if (rightLon < -Math.PI || rightLon > Math.PI) throw new IllegalArgumentException("Right longitude out of range"); double extent = rightLon - leftLon; if (extent < 0.0) { extent += 2.0 * Math.PI; } if (extent < Math.PI) throw new IllegalArgumentException("Width of rectangle too small"); this.latitude = latitude; this.leftLon = leftLon; this.rightLon = rightLon; final double sinLatitude = Math.sin(latitude); final double cosLatitude = Math.cos(latitude); final double sinLeftLon = Math.sin(leftLon); final double cosLeftLon = Math.cos(leftLon); final double sinRightLon = Math.sin(rightLon); final double cosRightLon = Math.cos(rightLon); // Now build the two points this.LHC = new GeoPoint(planetModel, sinLatitude, sinLeftLon, cosLatitude, cosLeftLon, latitude, leftLon); this.RHC = new GeoPoint(planetModel, sinLatitude, sinRightLon, cosLatitude, cosRightLon, latitude, rightLon); this.plane = new Plane(planetModel, sinLatitude); // Normalize while (leftLon > rightLon) { rightLon += Math.PI * 2.0; } double middleLon = (leftLon + rightLon) * 0.5; double sinMiddleLon = Math.sin(middleLon); double cosMiddleLon = Math.cos(middleLon); this.centerPoint = new GeoPoint(planetModel, sinLatitude, sinMiddleLon, cosLatitude, cosMiddleLon); this.leftPlane = new SidedPlane(centerPoint, cosLeftLon, sinLeftLon); this.rightPlane = new SidedPlane(centerPoint, cosRightLon, sinRightLon); this.planePoints = new GeoPoint[]{LHC, RHC}; this.eitherBound = new EitherBound(); this.edgePoints = new GeoPoint[]{centerPoint}; }
Constructor for deserialization.
Params:
  • planetModel – is the planet model.
  • inputStream – is the input stream.
/** * Constructor for deserialization. * @param planetModel is the planet model. * @param inputStream is the input stream. */
public GeoWideDegenerateHorizontalLine(final PlanetModel planetModel, final InputStream inputStream) throws IOException { this(planetModel, SerializableObject.readDouble(inputStream), SerializableObject.readDouble(inputStream), SerializableObject.readDouble(inputStream)); } @Override public void write(final OutputStream outputStream) throws IOException { SerializableObject.writeDouble(outputStream, latitude); SerializableObject.writeDouble(outputStream, leftLon); SerializableObject.writeDouble(outputStream, rightLon); } @Override public GeoBBox expand(final double angle) { final double newTopLat = latitude + angle; final double newBottomLat = latitude - angle; // Figuring out when we escalate to a special case requires some prefiguring double currentLonSpan = rightLon - leftLon; if (currentLonSpan < 0.0) currentLonSpan += Math.PI * 2.0; double newLeftLon = leftLon - angle; double newRightLon = rightLon + angle; if (currentLonSpan + 2.0 * angle >= Math.PI * 2.0) { newLeftLon = -Math.PI; newRightLon = Math.PI; } return GeoBBoxFactory.makeGeoBBox(planetModel, newTopLat, newBottomLat, newLeftLon, newRightLon); } @Override public boolean isWithin(final double x, final double y, final double z) { return plane.evaluateIsZero(x, y, z) && (leftPlane.isWithin(x, y, z) || rightPlane.isWithin(x, y, z)); } @Override public double getRadius() { // Here we compute the distance from the middle point to one of the corners. However, we need to be careful // to use the longest of three distances: the distance to a corner on the top; the distnace to a corner on the bottom, and // the distance to the right or left edge from the center. final double topAngle = centerPoint.arcDistance(RHC); final double bottomAngle = centerPoint.arcDistance(LHC); return Math.max(topAngle, bottomAngle); } @Override public GeoPoint getCenter() { return centerPoint; } @Override public GeoPoint[] getEdgePoints() { return edgePoints; } @Override public boolean intersects(final Plane p, final GeoPoint[] notablePoints, final Membership... bounds) { // Right and left bounds are essentially independent hemispheres; crossing into the wrong part of one // requires crossing into the right part of the other. So intersection can ignore the left/right bounds. return p.intersects(planetModel, plane, notablePoints, planePoints, bounds, eitherBound); } @Override public boolean intersects(final GeoShape geoShape) { // Right and left bounds are essentially independent hemispheres; crossing into the wrong part of one // requires crossing into the right part of the other. So intersection can ignore the left/right bounds. return geoShape.intersects(plane, planePoints, eitherBound); } @Override public void getBounds(Bounds bounds) { super.getBounds(bounds); bounds.isWide() .addHorizontalPlane(planetModel, latitude, plane, eitherBound) .addPoint(LHC) .addPoint(RHC); } @Override public int getRelationship(final GeoShape path) { if (intersects(path)) { return OVERLAPS; } if (path.isWithin(centerPoint)) { 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, eitherBound); final double LHCDistance = distanceStyle.computeDistance(LHC, x,y,z); final double RHCDistance = distanceStyle.computeDistance(RHC, x,y,z); return Math.min( distance, Math.min(LHCDistance, RHCDistance)); } @Override public boolean equals(Object o) { if (!(o instanceof GeoWideDegenerateHorizontalLine)) return false; GeoWideDegenerateHorizontalLine other = (GeoWideDegenerateHorizontalLine) o; return super.equals(other) && other.LHC.equals(LHC) && other.RHC.equals(RHC); } @Override public int hashCode() { int result = super.hashCode(); result = 31 * result + LHC.hashCode(); result = 31 * result + RHC.hashCode(); return result; } @Override public String toString() { return "GeoWideDegenerateHorizontalLine: {planetmodel="+planetModel+", latitude=" + latitude + "(" + latitude * 180.0 / Math.PI + "), leftlon=" + leftLon + "(" + leftLon * 180.0 / Math.PI + "), rightLon=" + rightLon + "(" + rightLon * 180.0 / Math.PI + ")}"; }
Membership implementation representing a wide cutoff (more than 180 degrees).
/** Membership implementation representing a wide cutoff (more than 180 degrees). */
protected class EitherBound implements Membership {
Constructor.
/** Constructor. */
public EitherBound() { } @Override public boolean isWithin(final double x, final double y, final double z) { return leftPlane.isWithin(x, y, z) || rightPlane.isWithin(x, y, z); } } }