0

次のコードでは、期待した結果が得られません。ポリゴンのノードとしてすべてのポイントがあり、ポリゴンにはすべてのポイントが含まれていると期待しています。API を正しく使用するには? また、thresold を 1 未満に設定すると、プログラムはある種の無限ループに陥ります。

import com.vividsolutions.jts.geom.Coordinate;
import com.vividsolutions.jts.geom.Geometry;
import com.vividsolutions.jts.geom.GeometryCollection;
import com.vividsolutions.jts.geom.GeometryFactory;
import com.vividsolutions.jts.geom.Point;
import java.util.Set;
import org.geotools.geometry.jts.JTSFactoryFinder;
import org.opensphere.geometry.algorithm.ConcaveHull;


public class ConcaveHullTrial {

    public static void main(String args[]) {

        ConcaveHull cch = null;
        com.vividsolutions.jts.geom.Point point[] = new Point[5];
        Set<GeometryFactory> testing = JTSFactoryFinder.getGeometryFactories();
        for (GeometryFactory f : testing) {
            System.out.println(f);
        }
        GeometryFactory gf = JTSFactoryFinder.getGeometryFactory();
        point[0] = gf.createPoint(new Coordinate(0.0, 0.0));
        point[1] = gf.createPoint(new Coordinate(0.0, 2.0));
        point[2] = gf.createPoint(new Coordinate(2.0, 0.0));
        point[3] = gf.createPoint(new Coordinate(2.0, 2.0));
        point[4] = gf.createPoint(new Coordinate(1.8, 1.5));

        GeometryCollection gc = new GeometryCollection(point, gf);
        double threshold = 4;
        cch = new ConcaveHull(gc, threshold);
        Geometry hull = cch.getConcaveHull();
        for (int i = 0; i < 5; i++) {
            System.out.println(hull.covers(point[i]));
        }
        for (Coordinate c: hull.getCoordinates()) {
            System.out.println(c);
        }
    }
}

結果:

false
true
false
true
false
(0.0, 2.0, NaN)
(2.0, 2.0, NaN)
4

1 に答える 1