Question

我从WorldWind的Sphere - 线交叉逻辑中看到了看似矛盾的行为。我创建了一个Sphere和Line并且它们相交但是然后交集返回null（扫描代码用于注释：// ***这是它变得糟糕的地方）。

这是视觉上发生的事情（线条是灰色的，但很难看到）：

public class WWTest extends ApplicationTemplate {

    public static class VisualizationFrame extends ApplicationTemplate.AppFrame {

        public VisualizationFrame() {
            super(new Dimension(1200, 1024));
            final Globe globe = getWwd().getModel().getGlobe();

            //Create a sphere at 0,0 on the surface of the Earth wtih a 60 NMi radius
            final Vec4 sphereCenter = globe.computePointFromLocation(LatLon.ZERO);
            final Sphere sphere = new Sphere(sphereCenter, 111120);
            // Draw the sphere
            final RenderableLayer sphereLayer = new RenderableLayer();
            sphereLayer.addRenderable(sphere);

            final RenderableLayer pathLayer = new RenderableLayer();
            // Create a line at 10k feet (3048 meters) that starts outside the sphere at (2,-2) and proceeds into the sphere at (0.5, 0.5)
            final Position lineStart = Position.fromDegrees(2, -2, 3048);
            final Position lineEnds = Position.fromDegrees(0.5, 0.5, 3048);
            final Path asPath = new Path(lineStart, lineEnds);
            pathLayer.addRenderable(asPath);

            // Now that we've visualized the line, let's do some intersection math
            final Vec4 lineStartsAsVec = globe.computePointFromPosition(lineStart);
            final Vec4 lineEndsAsVec = globe.computePointFromPosition(lineEnds);
            final Line asLine = Line.fromSegment(lineStartsAsVec, lineEndsAsVec);

            // *** This is where it gets whacky - true, but no intersection?
            final boolean doesIntersect = sphere.intersects(asLine);
            final Intersection[] intersection = sphere.intersect(asLine);
            //outputs: Intersection found: null
            System.out.println(doesIntersect ? "Intersection found: " + Arrays.toString(intersection) : "No intersection, Why Not!?!?"); 

            insertBeforeCompass(getWwd(), sphereLayer);
            insertBeforeCompass(getWwd(), pathLayer);
            getWwd().getView().setEyePosition(Position.fromDegrees(0, 0, 500_000));
            getLayerPanel().update(getWwd());
        }
    }

    public static void main(String[] args) {
        ApplicationTemplate.start("World Wind Sphere-Line Intersection", VisualizationFrame.class);

    }
}

以下是我宣布将WorldWind纳入我的maven项目的依赖项（我也尝试过版本＆＃39; 2.0.0-986＆＃39;，但这似乎没有帮助）：

<dependency>
    <groupId>gov.nasa</groupId>
    <artifactId>worldwind</artifactId>
    <version>2.0.0</version>
</dependency>
<dependency>
    <groupId>gov.nasa</groupId>
    <artifactId>worldwindx</artifactId>
    <version>2.0.0</version>
</dependency>
<dependency>
    <groupId>org.jogamp.gluegen</groupId>
    <artifactId>gluegen-rt-main</artifactId>
    <version>2.2.4</version>
</dependency>
<dependency>
    <groupId>org.jogamp.jogl</groupId>
    <artifactId>jogl-all-main</artifactId>
    <version>2.2.4</version>
</dependency>

要完全彻底，以下是代码导入：

import gov.nasa.worldwind.geom.Intersection;
import gov.nasa.worldwind.geom.LatLon;
import gov.nasa.worldwind.geom.Line;
import gov.nasa.worldwind.geom.Position;
import gov.nasa.worldwind.geom.Sphere;
import gov.nasa.worldwind.geom.Vec4;
import gov.nasa.worldwind.globes.Globe;
import gov.nasa.worldwind.layers.RenderableLayer;
import gov.nasa.worldwind.render.Path;
import gov.nasa.worldwindx.examples.ApplicationTemplate;
import static gov.nasa.worldwindx.examples.ApplicationTemplate.insertBeforeCompass;
import java.awt.Dimension;
import java.util.Arrays;

Answer 1

如果你看一下Sphere#intersect()的实现，它期望坐标中的线以球体的起源（不是地球的）为中心，这几乎肯定是一个错误。你应该能够做到：

final Vec4 pa = lineStartsAsVec.subtract3(sphereCenter);
final Vec4 pb = lineEndsAsVec.subtract3(sphereCenter);
final Line asLine2 = Line.fromSegment(pa, pb);
final Intersection[] intersection = sphere.intersect(asLine2);

请记住，返回的交叉点仍然是以球体原点为中心的笛卡尔坐标，所以要将它们转换回你需要做的World Wind Cartesian：

final Vec4 intersectionPos = intersection[0].getIntersectionPoint().add3(sphereCenter);

该实现还认为该行无限长，因此它将返回两个点，而不是一个。

实现您自己的intersect()版本并在正常坐标中工作并考虑行的长度非常简单，请参阅here。

Answer 2

如果查看Worldwind source code，特别是Sphere类的intersects（）和intersect（）方法，并使用代码作为输入逐步执行它们，您会看到以下内容：

方法：

public boolean intersects(Line line)

返回true，因为从球体中心到线的距离小于球体的半径，正如预期的那样。

方法：

public final Intersection[] intersect(Line line)

事实证明，二次方的判别式小于零（即二次方程没有真正的根 - 两个不同的复根）。

WorldWind API参考是here。

涉及的具体方法是：

    /**
 * Tests for intersection with a <code>Line</code>.
 *
 * @param line the <code>Line</code> with which to test for intersection
 *
 * @return true if <code>line</code> intersects or makes a tangent with the surface of this <code>Sphere</code>
 *
 * @throws IllegalArgumentException if <code>line</code> is null
 */
public boolean intersects(Line line)
{
    if (line == null)
    {
        String msg = Logging.getMessage("nullValue.LineIsNull");
        Logging.logger().severe(msg);
        throw new IllegalArgumentException(msg);
    }
    return line.distanceTo(this.center) <= this.radius;
}

和

    /**
 * Obtains the intersections of this sphere with a line. The returned array may be either null or of zero length if
 * no intersections are discovered. It does not contain null elements and will have a size of 2 at most. Tangential
 * intersections are marked as such. <code>line</code> is considered to have infinite length in both directions.
 *
 * @param line the <code>Line</code> with which to intersect this <code>Sphere</code>
 *
 * @return an array containing all the intersections of this <code>Sphere</code> and <code>line</code>
 *
 * @throws IllegalArgumentException if <code>line</code> is null
 */
public final Intersection[] intersect(Line line)
{
    if (line == null)
    {
        String message = Logging.getMessage("nullValue.LineIsNull");
        Logging.logger().severe(message);
        throw new IllegalArgumentException(message);
    }

    double a = line.getDirection().getLengthSquared3();
    double b = 2 * line.selfDot();
    double c = line.getOrigin().getLengthSquared3() - this.radius * this.radius;

    double discriminant = Sphere.discriminant(a, b, c);
    if (discriminant < 0)
        return null;

    double discriminantRoot = Math.sqrt(discriminant);
    if (discriminant == 0)
    {
        Vec4 p = line.getPointAt((-b - discriminantRoot) / (2 * a));
        return new Intersection[] {new Intersection(p, true)};
    }
    else // (discriminant > 0)
    {
        Vec4 near = line.getPointAt((-b - discriminantRoot) / (2 * a));
        Vec4 far = line.getPointAt((-b + discriminantRoot) / (2 * a));
        return new Intersection[] {new Intersection(near, false), new Intersection(far, false)};
    }
}

使用：

    /**
 * Calculates a discriminant. A discriminant is useful to determine the number of roots to a quadratic equation. If
 * the discriminant is less than zero, there are no roots. If it equals zero, there is one root. If it is greater
 * than zero, there are two roots.
 *
 * @param a the coefficient of the second order pronumeral
 * @param b the coefficient of the first order pronumeral
 * @param c the constant parameter in the quadratic equation
 *
 * @return the discriminant "b squared minus 4ac"
 */
private static double discriminant(double a, double b, double c)
{
    return b * b - 4 * a * c;
}

在这种情况下，您的代码将失败：

if (discriminant < 0)

测试。

看起来我回答这个问题的速度有点慢，与此同时Chris K已经指出这是因为期望线坐标居中的intersect（）方法球体的起源，而不是地球的起源。

正如Chris K所说，这似乎是一个错误，应该与这个源代码的维护者一起记录。

WorldWind Sphere Line Intersection Bug？

2 个答案: