Method _computeBoundingBox

CesiumGeospatial/src/BoundingRegion.cpp:274–474

static*/

} // namespace

/*static*/ OrientedBoundingBox BoundingRegion::_computeBoundingBox(
    const GlobeRectangle& rectangle,
    double minimumHeight,
    double maximumHeight,
    const Ellipsoid& ellipsoid) {
  //>>includeStart('debug', pragmas.debug);
  if (!Math::equalsEpsilon(
          ellipsoid.getRadii().x,
          ellipsoid.getRadii().y,
          Math::Epsilon15)) {
    throw std::runtime_error(
        "Ellipsoid must be an ellipsoid of revolution (radii.x == radii.y)");
  }
  //>>includeEnd('debug');

  double minX, maxX, minY, maxY, minZ, maxZ;
  Plane plane(glm::dvec3(0.0, 0.0, 1.0), 0.0);

  if (rectangle.computeWidth() <= Math::OnePi) {
    // The bounding box will be aligned with the tangent plane at the center of
    // the rectangle.
    const Cartographic tangentPointCartographic = rectangle.computeCenter();
    const glm::dvec3 tangentPoint =
        ellipsoid.cartographicToCartesian(tangentPointCartographic);
    const EllipsoidTangentPlane tangentPlane(tangentPoint, ellipsoid);
    plane = tangentPlane.getPlane();

    // If the rectangle spans the equator, CW is instead aligned with the
    // equator (because it sticks out the farthest at the equator).
    const double lonCenter = tangentPointCartographic.longitude;
    const double latCenter =
        rectangle.getSouth() < 0.0 && rectangle.getNorth() > 0.0
            ? 0.0
            : tangentPointCartographic.latitude;

    // Compute XY extents using the rectangle at maximum height
    const Cartographic perimeterCartographicNC(
        lonCenter,
        rectangle.getNorth(),
        maximumHeight);
    Cartographic perimeterCartographicNW(
        rectangle.getWest(),
        rectangle.getNorth(),
        maximumHeight);
    const Cartographic perimeterCartographicCW(
        rectangle.getWest(),
        latCenter,
        maximumHeight);
    Cartographic perimeterCartographicSW(
        rectangle.getWest(),
        rectangle.getSouth(),
        maximumHeight);
    const Cartographic perimeterCartographicSC(
        lonCenter,
        rectangle.getSouth(),
        maximumHeight);

    const glm::dvec3 perimeterCartesianNC =