Search in sources :

Example 1 with OrientableCurve

use of org.opengis.geometry.primitive.OrientableCurve in project geotoolkit by Geomatys.

the class JTSMultiCurve method applyCRSonChild.

public void applyCRSonChild() {
    for (OrientableCurve child : getElements()) {
        if (child.getCoordinateReferenceSystem() == null && child instanceof JTSCurve) {
            JTSCurve childCurve = (JTSCurve) child;
            childCurve.setCoordinateReferenceSystem(getCoordinateReferenceSystem());
            for (CurveSegment cv : childCurve.getSegments()) {
                if (cv instanceof JTSLineString) {
                    ((JTSLineString) cv).setCoordinateReferenceSystem(getCoordinateReferenceSystem());
                    ((JTSLineString) cv).applyCRSOnChild();
                }
            }
        }
    }
}
Also used : CurveSegment(org.opengis.geometry.primitive.CurveSegment) OrientableCurve(org.opengis.geometry.primitive.OrientableCurve) JTSCurve(org.geotoolkit.geometry.isoonjts.spatialschema.geometry.primitive.JTSCurve) JTSLineString(org.geotoolkit.geometry.isoonjts.spatialschema.geometry.geometry.JTSLineString)

Example 2 with OrientableCurve

use of org.opengis.geometry.primitive.OrientableCurve in project geotoolkit by Geomatys.

the class SurfaceImplTest method testLargeSurface.

/**
 * We need to create a large surface with 7000 points
 */
@Test
public void testLargeSurface() {
    int NUMBER = 100000;
    double delta = 360.0 / (double) NUMBER;
    PointArray points = postitionFactory.createPointArray();
    for (double angle = 0.0; angle < 360.0; angle += delta) {
        double[] coordinates = new double[] { Math.sin(Math.toRadians(angle)), Math.cos(Math.toRadians(angle)) };
        DirectPosition point = postitionFactory.createDirectPosition(coordinates);
        points.add(point);
    }
    List<OrientableCurve> curves = new ArrayList<OrientableCurve>();
    // A curve will be created
    // - The curve will be set as parent curves for the Curve segments
    // - Start and end params for the CurveSegments will be set
    List<CurveSegment> segmentList = new ArrayList<CurveSegment>();
    for (int i = 0; i < points.size(); i++) {
        int start = i;
        int end = (i + 1) % points.size();
        DirectPosition point1 = points.getDirectPosition(start, null);
        DirectPosition point2 = points.getDirectPosition(end, null);
        LineSegment segment = geometryFactory.createLineSegment(point1, point2);
        segmentList.add(segment);
    }
    Curve curve = primitiveFactory.createCurve(segmentList);
    curves.add(curve);
    Ring ring = primitiveFactory.createRing(curves);
    SurfaceBoundary boundary = primitiveFactory.createSurfaceBoundary(ring, new ArrayList());
    JTSSurface surface = (JTSSurface) primitiveFactory.createSurface(boundary);
    Geometry peer = surface.computeJTSPeer();
}
Also used : DirectPosition(org.opengis.geometry.DirectPosition) SurfaceBoundary(org.opengis.geometry.primitive.SurfaceBoundary) ArrayList(java.util.ArrayList) Curve(org.opengis.geometry.primitive.Curve) OrientableCurve(org.opengis.geometry.primitive.OrientableCurve) OrientableCurve(org.opengis.geometry.primitive.OrientableCurve) PointArray(org.opengis.geometry.coordinate.PointArray) Geometry(org.locationtech.jts.geom.Geometry) CurveSegment(org.opengis.geometry.primitive.CurveSegment) Ring(org.opengis.geometry.primitive.Ring) LineSegment(org.opengis.geometry.coordinate.LineSegment) Test(org.junit.Test)

Example 3 with OrientableCurve

use of org.opengis.geometry.primitive.OrientableCurve in project geotoolkit by Geomatys.

the class GMLUtilities method getGMLFromISO.

/**
 * @param geometry The ISO geometry to convert.
 * @return A GML 3.1.1 geometry matching given geometry definition.
 * @deprecated This method should not be used for multiple reasons:
 * <ol>
 * <li>OpenGIS definition is based on an obsolete ISO-19107 draft</li>
 * <li>The GML version is fixed and obsolete</li>
 * <li>Only partial management of geometries</li>
 * </ol>
 */
public static AbstractGeometryType getGMLFromISO(final org.opengis.geometry.Geometry geometry) {
    if (geometry instanceof Point) {
        Point point = (Point) geometry;
        PointType gmlPoint = new PointType(null, point.getDirectPosition());
        return gmlPoint;
    } else if (geometry instanceof OrientableSurface) {
        OrientableSurface surface = (OrientableSurface) geometry;
        SurfaceBoundary boundary = surface.getBoundary();
        Ring exterior = boundary.getExterior();
        List<CurvePropertyType> curves = new ArrayList<CurvePropertyType>();
        for (Primitive p : exterior.getElements()) {
            curves.add(new CurvePropertyType((CurveType) getGMLFromISO(p)));
        }
        RingType gmlExterior = new RingType();
        gmlExterior.getCurveMember().addAll(curves);
        List<Ring> interiors = boundary.getInteriors();
        List<RingType> gmlInteriors = new ArrayList<RingType>();
        for (Ring interior : interiors) {
            List<CurvePropertyType> intcurves = new ArrayList<CurvePropertyType>();
            for (Primitive p : interior.getElements()) {
                intcurves.add(new CurvePropertyType((CurveType) getGMLFromISO(p)));
            }
            RingType gmlinterior = new RingType();
            gmlinterior.getCurveMember().addAll(intcurves);
            gmlInteriors.add(gmlinterior);
        }
        PolygonType poly = new PolygonType(gmlExterior, gmlInteriors);
        return poly;
    } else if (geometry instanceof MultiSurface) {
        MultiSurface multiPrim = (MultiSurface) geometry;
        List<PolygonPropertyType> geometries = new ArrayList<PolygonPropertyType>();
        for (Geometry prim : multiPrim.getElements()) {
            PolygonType element = (PolygonType) getGMLFromISO(prim);
            PolygonPropertyType gp = new PolygonPropertyType((PolygonType) element);
            geometries.add(gp);
        }
        MultiPolygonType gmlMulti = new MultiPolygonType(null, geometries);
        return gmlMulti;
    } else if (geometry instanceof MultiCurve) {
        MultiCurve multiPrim = (MultiCurve) geometry;
        List<CurvePropertyType> geometries = new ArrayList<CurvePropertyType>();
        for (OrientableCurve prim : multiPrim.getElements()) {
            AbstractCurveType element = (AbstractCurveType) getGMLFromISO(prim);
            CurvePropertyType gp = new CurvePropertyType((AbstractCurveType) element);
            geometries.add(gp);
        }
        MultiCurveType gmlMulti = new MultiCurveType(geometries);
        return gmlMulti;
    } else if (geometry instanceof MultiPoint) {
        MultiPoint multiPrim = (MultiPoint) geometry;
        List<PointPropertyType> geometries = new ArrayList<PointPropertyType>();
        for (Point prim : multiPrim.getElements()) {
            PointType element = (PointType) getGMLFromISO(prim);
            PointPropertyType gp = new PointPropertyType((PointType) element);
            geometries.add(gp);
        }
        MultiPointType gmlMulti = new MultiPointType(null, geometries);
        return gmlMulti;
    } else if (geometry instanceof MultiPrimitive) {
        MultiPrimitive multiPrim = (MultiPrimitive) geometry;
        List<GeometryPropertyType> geometries = new ArrayList<GeometryPropertyType>();
        for (Primitive prim : multiPrim.getElements()) {
            AbstractGMLType element = getGMLFromISO(prim);
            GeometryPropertyType gp = new GeometryPropertyType((AbstractGeometryType) element);
            geometries.add(gp);
        }
        MultiGeometryType gmlMulti = new MultiGeometryType(geometries);
        return gmlMulti;
    } else if (geometry instanceof Curve) {
        Curve curve = (Curve) geometry;
        List<? extends CurveSegment> segments = curve.getSegments();
        List<LineStringSegmentType> gmlSegments = new ArrayList<LineStringSegmentType>();
        for (CurveSegment segment : segments) {
            CurveInterpolationType interpolation = CurveInterpolationType.fromValue(segment.getInterpolation().identifier());
            PointArray array = GeometricUtilities.getSamplePoints(segment);
            List<DirectPosition> positions = new ArrayList<DirectPosition>();
            for (int i = 0; i < array.size(); i++) {
                positions.add(array.getDirectPosition(i, null));
            }
            LineStringSegmentType gmlSegment = new LineStringSegmentType(segment.getNumDerivativesAtStart(), segment.getNumDerivativesAtEnd(), segment.getNumDerivativesInterior(), interpolation, positions);
            gmlSegments.add(gmlSegment);
        }
        CurveType gmlCurve = new CurveType(gmlSegments);
        return gmlCurve;
    } else if (geometry instanceof LineString) {
        LineString line = (LineString) geometry;
        PointArray array = GeometricUtilities.getSamplePoints(line);
        List<DirectPosition> positions = new ArrayList<DirectPosition>();
        for (int i = 0; i < array.size(); i++) {
            positions.add(array.getDirectPosition(i, null));
        }
        LineStringType gmlLine = new LineStringType(positions);
        return gmlLine;
    } else if (geometry instanceof Polygon) {
        Polygon polygon = (Polygon) geometry;
        SurfaceBoundary boundary = polygon.getBoundary();
        Ring exterior = boundary.getExterior();
        List<CurvePropertyType> curves = new ArrayList<CurvePropertyType>();
        for (Primitive p : exterior.getElements()) {
            curves.add(new CurvePropertyType((CurveType) getGMLFromISO(p)));
        }
        RingType gmlExterior = new RingType();
        gmlExterior.getCurveMember().addAll(curves);
        List<Ring> interiors = boundary.getInteriors();
        List<RingType> gmlInteriors = new ArrayList<RingType>();
        for (Ring interior : interiors) {
            List<CurvePropertyType> intcurves = new ArrayList<CurvePropertyType>();
            for (Primitive p : interior.getElements()) {
                intcurves.add(new CurvePropertyType((CurveType) getGMLFromISO(p)));
            }
            RingType gmlinterior = new RingType();
            gmlinterior.getCurveMember().addAll(intcurves);
            gmlInteriors.add(gmlinterior);
        }
        PolygonType gmlPolygon = new PolygonType(gmlExterior, gmlInteriors);
        return gmlPolygon;
    } else if (geometry instanceof PolyhedralSurface) {
        PolyhedralSurface polySurface = (PolyhedralSurface) geometry;
        List<PolygonPatchType> gmlPatches = new ArrayList<PolygonPatchType>();
        List<? extends Polygon> patches = polySurface.getPatches();
        for (Polygon polygon : patches) {
            SurfaceInterpolationType interpolation = SurfaceInterpolationType.fromValue(polygon.getInterpolation().identifier());
            SurfaceBoundary boundary = polygon.getBoundary();
            Ring exterior = boundary.getExterior();
            List<CurvePropertyType> curves = new ArrayList<CurvePropertyType>();
            for (Primitive p : exterior.getElements()) {
                curves.add(new CurvePropertyType((CurveType) getGMLFromISO(p)));
            }
            RingType gmlExterior = new RingType();
            gmlExterior.getCurveMember().addAll(curves);
            List<Ring> interiors = boundary.getInteriors();
            List<RingType> gmlInteriors = new ArrayList<RingType>();
            for (Ring interior : interiors) {
                List<CurvePropertyType> intcurves = new ArrayList<CurvePropertyType>();
                for (Primitive p : interior.getElements()) {
                    intcurves.add(new CurvePropertyType((CurveType) getGMLFromISO(p)));
                }
                RingType gmlinterior = new RingType();
                gmlinterior.getCurveMember().addAll(intcurves);
                gmlInteriors.add(gmlinterior);
            }
            PolygonPatchType patche = new PolygonPatchType(interpolation, gmlExterior, gmlInteriors);
            gmlPatches.add(patche);
        }
        PolygonPatchArrayPropertyType pathArray = new PolygonPatchArrayPropertyType(gmlPatches);
        PolyhedralSurfaceType gmlPolySurface = new PolyhedralSurfaceType(pathArray);
        return gmlPolySurface;
    } else {
        System.out.println("unexpected iso geometry type:" + geometry.getClass().getName());
    }
    return null;
}
Also used : MultiPoint(org.opengis.geometry.aggregate.MultiPoint) DirectPosition(org.opengis.geometry.DirectPosition) MultiCurveType(org.geotoolkit.gml.xml.v311.MultiCurveType) OrientableSurface(org.opengis.geometry.primitive.OrientableSurface) ArrayList(java.util.ArrayList) MultiGeometryType(org.geotoolkit.gml.xml.v311.MultiGeometryType) MultiPrimitive(org.opengis.geometry.aggregate.MultiPrimitive) CurveSegment(org.opengis.geometry.primitive.CurveSegment) MultiCurve(org.opengis.geometry.aggregate.MultiCurve) ArrayList(java.util.ArrayList) List(java.util.List) CurvePropertyType(org.geotoolkit.gml.xml.v311.CurvePropertyType) Polygon(org.opengis.geometry.coordinate.Polygon) SurfaceBoundary(org.opengis.geometry.primitive.SurfaceBoundary) AbstractGeometryType(org.geotoolkit.gml.xml.v311.AbstractGeometryType) Curve(org.opengis.geometry.primitive.Curve) MultiCurve(org.opengis.geometry.aggregate.MultiCurve) OrientableCurve(org.opengis.geometry.primitive.OrientableCurve) PolygonType(org.geotoolkit.gml.xml.v311.PolygonType) MultiPolygonType(org.geotoolkit.gml.xml.v311.MultiPolygonType) OrientableCurve(org.opengis.geometry.primitive.OrientableCurve) LineStringType(org.geotoolkit.gml.xml.v311.LineStringType) MultiPointType(org.geotoolkit.gml.xml.v311.MultiPointType) CurveInterpolationType(org.geotoolkit.gml.xml.v311.CurveInterpolationType) RingType(org.geotoolkit.gml.xml.v311.RingType) PolyhedralSurfaceType(org.geotoolkit.gml.xml.v311.PolyhedralSurfaceType) LineStringSegmentType(org.geotoolkit.gml.xml.v311.LineStringSegmentType) MultiSurface(org.opengis.geometry.aggregate.MultiSurface) SurfaceInterpolationType(org.geotoolkit.gml.xml.v311.SurfaceInterpolationType) AbstractCurveType(org.geotoolkit.gml.xml.v311.AbstractCurveType) MultiCurveType(org.geotoolkit.gml.xml.v311.MultiCurveType) CurveType(org.geotoolkit.gml.xml.v311.CurveType) AbstractGMLType(org.geotoolkit.gml.xml.v311.AbstractGMLType) PolygonPatchArrayPropertyType(org.geotoolkit.gml.xml.v311.PolygonPatchArrayPropertyType) MultiPrimitive(org.opengis.geometry.aggregate.MultiPrimitive) Primitive(org.opengis.geometry.primitive.Primitive) MultiPolygonType(org.geotoolkit.gml.xml.v311.MultiPolygonType) PolygonPatchType(org.geotoolkit.gml.xml.v311.PolygonPatchType) Point(org.opengis.geometry.primitive.Point) MultiPoint(org.opengis.geometry.aggregate.MultiPoint) PointArray(org.opengis.geometry.coordinate.PointArray) Point(org.opengis.geometry.primitive.Point) MultiPoint(org.opengis.geometry.aggregate.MultiPoint) PolyhedralSurface(org.opengis.geometry.coordinate.PolyhedralSurface) Geometry(org.opengis.geometry.Geometry) AbstractCurveType(org.geotoolkit.gml.xml.v311.AbstractCurveType) LineString(org.opengis.geometry.coordinate.LineString) Ring(org.opengis.geometry.primitive.Ring) PolygonPropertyType(org.geotoolkit.gml.xml.v311.PolygonPropertyType) PointType(org.geotoolkit.gml.xml.v311.PointType) MultiPointType(org.geotoolkit.gml.xml.v311.MultiPointType) PointPropertyType(org.geotoolkit.gml.xml.v311.PointPropertyType) GeometryPropertyType(org.geotoolkit.gml.xml.v311.GeometryPropertyType)

Example 4 with OrientableCurve

use of org.opengis.geometry.primitive.OrientableCurve in project geotoolkit by Geomatys.

the class ISOGeometryJ2D method createRectangle.

/**
 * Creates a jts Geometry object representing a rectangle with the given
 * parameters
 *
 * @param x left coordinate
 * @param y bottom coordinate
 * @param w width
 * @param h height
 * @return a rectangle with the specified position and size
 */
private Geometry createRectangle(final double x, final double y, final double w, final double h) {
    final GeometryFactory gf = new JTSGeometryFactory(geometry.getCoordinateReferenceSystem());
    final PrimitiveFactory pf = new JTSPrimitiveFactory(geometry.getCoordinateReferenceSystem());
    final List<Position> points = new ArrayList<Position>();
    points.add(new DirectPosition2D(x, y));
    points.add(new DirectPosition2D(x + w, y));
    points.add(new DirectPosition2D(x + w, y + h));
    points.add(new DirectPosition2D(x, y + h));
    points.add(new DirectPosition2D(x, y));
    CurveSegment segment = gf.createLineString(points);
    OrientableCurve curve = pf.createCurve(Collections.singletonList(segment));
    Ring exterior = pf.createRing(Collections.singletonList(curve));
    SurfaceBoundary boundary = pf.createSurfaceBoundary(exterior, null);
    return boundary;
}
Also used : JTSGeometryFactory(org.geotoolkit.geometry.isoonjts.spatialschema.geometry.geometry.JTSGeometryFactory) GeometryFactory(org.opengis.geometry.coordinate.GeometryFactory) JTSGeometryFactory(org.geotoolkit.geometry.isoonjts.spatialschema.geometry.geometry.JTSGeometryFactory) CurveSegment(org.opengis.geometry.primitive.CurveSegment) SurfaceBoundary(org.opengis.geometry.primitive.SurfaceBoundary) Position(org.opengis.geometry.coordinate.Position) Ring(org.opengis.geometry.primitive.Ring) ArrayList(java.util.ArrayList) JTSPrimitiveFactory(org.geotoolkit.geometry.isoonjts.spatialschema.geometry.primitive.JTSPrimitiveFactory) OrientableCurve(org.opengis.geometry.primitive.OrientableCurve) JTSPrimitiveFactory(org.geotoolkit.geometry.isoonjts.spatialschema.geometry.primitive.JTSPrimitiveFactory) PrimitiveFactory(org.opengis.geometry.primitive.PrimitiveFactory) DirectPosition2D(org.apache.sis.geometry.DirectPosition2D)

Example 5 with OrientableCurve

use of org.opengis.geometry.primitive.OrientableCurve in project geotoolkit by Geomatys.

the class GeometryParser method readPolygonText.

/**
 * Creates a <code>SurfaceBoundary</code> using the next token in the stream.
 *
 * @param tokenizer tokenizer over a stream of text in Well-known Text
 *                  format. The next tokens must form a &lt;Polygon Text&gt;.
 * @return a <code>Surface</code> specified by the vertices in the stream
 * @throws ParseException if the coordinates used to create the <code>Polygon</code>
 *                        shell and holes do not form closed linestrings, or if an unexpected
 *                        token was encountered.
 * @throws IOException    if an I/O error occurs
 */
private Surface readPolygonText(final StreamTokenizer tokenizer) throws IOException, ParseException {
    String nextToken = getNextEmptyOrOpener(tokenizer);
    if (nextToken.equals(EMPTY)) {
        return null;
    }
    Curve curve = readLinearRingText(tokenizer);
    List<OrientableCurve> curveList = Collections.singletonList((OrientableCurve) curve);
    Ring shell = primitiveFactory.createRing(curveList);
    // Ring shell = readLinearRingText(tokenizer);
    List<Ring> holes = new ArrayList<Ring>();
    nextToken = getNextCloserOrComma(tokenizer);
    while (nextToken.equals(COMMA)) {
        Curve holecurve = readLinearRingText(tokenizer);
        List<OrientableCurve> holeList = Collections.singletonList((OrientableCurve) holecurve);
        Ring hole = primitiveFactory.createRing(holeList);
        // Ring hole = readLinearRingText(tokenizer);
        holes.add(hole);
        nextToken = getNextCloserOrComma(tokenizer);
    }
    SurfaceBoundary sb = primitiveFactory.createSurfaceBoundary(shell, holes);
    return primitiveFactory.createSurface(sb);
}
Also used : SurfaceBoundary(org.opengis.geometry.primitive.SurfaceBoundary) Ring(org.opengis.geometry.primitive.Ring) Curve(org.opengis.geometry.primitive.Curve) OrientableCurve(org.opengis.geometry.primitive.OrientableCurve) ArrayList(java.util.ArrayList) LineString(org.opengis.geometry.coordinate.LineString) OrientableCurve(org.opengis.geometry.primitive.OrientableCurve)

Aggregations

OrientableCurve (org.opengis.geometry.primitive.OrientableCurve)5 ArrayList (java.util.ArrayList)4 CurveSegment (org.opengis.geometry.primitive.CurveSegment)4 Ring (org.opengis.geometry.primitive.Ring)4 SurfaceBoundary (org.opengis.geometry.primitive.SurfaceBoundary)4 Curve (org.opengis.geometry.primitive.Curve)3 DirectPosition (org.opengis.geometry.DirectPosition)2 PointArray (org.opengis.geometry.coordinate.PointArray)2 List (java.util.List)1 DirectPosition2D (org.apache.sis.geometry.DirectPosition2D)1 JTSGeometryFactory (org.geotoolkit.geometry.isoonjts.spatialschema.geometry.geometry.JTSGeometryFactory)1 JTSLineString (org.geotoolkit.geometry.isoonjts.spatialschema.geometry.geometry.JTSLineString)1 JTSCurve (org.geotoolkit.geometry.isoonjts.spatialschema.geometry.primitive.JTSCurve)1 JTSPrimitiveFactory (org.geotoolkit.geometry.isoonjts.spatialschema.geometry.primitive.JTSPrimitiveFactory)1 AbstractCurveType (org.geotoolkit.gml.xml.v311.AbstractCurveType)1 AbstractGMLType (org.geotoolkit.gml.xml.v311.AbstractGMLType)1 AbstractGeometryType (org.geotoolkit.gml.xml.v311.AbstractGeometryType)1 CurveInterpolationType (org.geotoolkit.gml.xml.v311.CurveInterpolationType)1 CurvePropertyType (org.geotoolkit.gml.xml.v311.CurvePropertyType)1 CurveType (org.geotoolkit.gml.xml.v311.CurveType)1