content: /**
   * 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.spatial.geometry.shape;
  
  import org.apache.lucene.spatial.geometry.FloatLatLng;
  import org.apache.lucene.spatial.geometry.LatLng;
  
  
  
  /**
   * Lat-long rect.  Instances are mutable.
   *
   * <p><font color="red"><b>NOTE:</b> This API is still in
   * flux and might change in incompatible ways in the next
   * release.</font>
   */
  public class LLRect {
    private LatLng ll, ur;
    
    public LLRect(LatLng ll, LatLng ur) {
      this.ll=ll;
      this.ur=ur;
    }
    
    public LLRect(LLRect other) {
      this.ll=other.ll;
      this.ur=other.ur;
    }
    
    /**
     * Return the area in units of lat-lng squared.  This is a contrived unit
     * that only has value when comparing to something else.
     */
    public double area() {
      return Math.abs((ll.getLat()-ur.getLat()) * (ll.getLng()-ur.getLng()));
    }
  
    public LatLng getLowerLeft() {
      return ll;
    }
    
    public LatLng getUpperRight() {
      return ur;
    }
    
    @Override
    public String toString() {
      return "{" + ll + ", " + ur + "}";
    }
  
    public LatLng getMidpoint() {
      return ll.calculateMidpoint(ur);
    }
  
    /**
     * Approximates a box centered at the given point with the given width and height in miles.
     * @param center
     * @param widthMi
     * @param heightMi
     */
    public static LLRect createBox(LatLng center, double widthMi, double heightMi) {
      double d = widthMi;
      LatLng ur = boxCorners(center, d, 45.0); // assume right angles
      LatLng ll = boxCorners(center, d, 225.0);
  
      //System.err.println("boxCorners: ur " + ur.getLat() + ',' + ur.getLng());
      //System.err.println("boxCorners: cnt " + center.getLat() + ',' + center.getLng());
      //System.err.println("boxCorners: ll " + ll.getLat() + ',' + ll.getLng());
      return new LLRect(ll, ur);
    }
    
    /**
     * Returns a rectangle shape for the bounding box
     */
    public Rectangle toRectangle() {
      return new Rectangle(ll.getLng(), ll.getLat(), ur.getLng(), ur.getLat());
    }
  
    private static LatLng boxCorners(LatLng center, double d, double brngdeg) {
      double a = center.getLat();
      double b = center.getLng();
      double R = 3963.0; // radius of earth in miles
      double brng = (Math.PI*brngdeg/180);
      double lat1 = (Math.PI*a/180);
      double lon1 = (Math.PI*b/180);
  
      // Haversine formula
      double lat2 = Math.asin( Math.sin(lat1)*Math.cos(d/R) +
                               Math.cos(lat1)*Math.sin(d/R)*Math.cos(brng) );
      double lon2 = lon1 + Math.atan2(Math.sin(brng)*Math.sin(d/R)*Math.cos(lat1),
                                      Math.cos(d/R)-Math.sin(lat1)*Math.sin(lat2));
  
      lat2 = (lat2*180)/Math.PI;
      lon2 = (lon2*180)/Math.PI;
  
      // normalize long first
      LatLng ll = normLng(lat2,lon2);
  
      // normalize lat - could flip poles
      ll = normLat(ll.getLat(),ll.getLng());
  
      return ll;
  }
  
    /**
     * Returns a normalized Lat rectangle shape for the bounding box
     * If you go over the poles, you need to flip the lng value too
     */
    private static LatLng normLat(double lat, double lng) {
      if (lat > 90.0) {
          lat = 90.0 - (lat - 90.0);
          if (lng < 0) {
                  lng = lng+180;
          } else {
                  lng=lng-180;
          }
      }
      else if (lat < -90.0) {
          lat = -90.0 - (lat + 90.0);
          if (lng < 0) {
                  lng = lng+180;
          } else {
                  lng=lng-180;
          }
      }
      LatLng ll=new FloatLatLng(lat, lng);
      return ll;
    }
  
    /**
     * Returns a normalized Lng rectangle shape for the bounding box
     */
    private static LatLng normLng(double lat,double lng) {
      if (lng > 180.0) {
          lng = -1.0*(180.0 - (lng - 180.0));
      }
      else if (lng < -180.0) {
          lng = (lng + 180.0)+180.0;
      }
      LatLng ll=new FloatLatLng(lat, lng);
      return ll;
    }
  
    @Override
    public int hashCode() {
      final int prime = 31;
      int result = 1;
      result = prime * result + ((ll == null) ? 0 : ll.hashCode());
      result = prime * result + ((ur == null) ? 0 : ur.hashCode());
      return result;
    }
  
    @Override
    public boolean equals(Object obj) {
      if (this == obj)
        return true;
      if (obj == null)
        return false;
      if (getClass() != obj.getClass())
        return false;
      LLRect other = (LLRect) obj;
      if (ll == null) {
        if (other.ll != null)
          return false;
      } else if (!ll.equals(other.ll))
        return false;
      if (ur == null) {
        if (other.ur != null)
          return false;
      } else if (!ur.equals(other.ur))
        return false;
      return true;
    }
    
    
  }
  