use of org.geotoolkit.gml.xml.WithCoordinates in project geotoolkit by Geomatys.
the class GeometryTransformer method getCoordinates.
/**
* Check that given geometry is a primitive geometry (implements {@link WithCoordinates}), and
* get its points.
* @param source The geometry to extract points from.
* @return Found points, never null, but can be empty.
* @throws UnconvertibleObjectException If the given geometry does not implement {@link WithCoordinates}.
*/
private Spliterator<Coordinate> getCoordinates(final Object source) throws UnconvertibleObjectException {
List<Double> values = null;
if (source instanceof WithCoordinates) {
final Coordinates coords = ((WithCoordinates) source).getCoordinates();
if (coords != null) {
/* HACK : In GML 3, coordinates are just a list of decimal values.
* The grouping by coordinate is done using "srsDimension" attribute
* on parent geometry type. However, with GML 2, there's another
* possibility : Coordinates use two distinct separators : one for
* decimal value separation, and another for coordinate separation.
* To manage both ways, we first check if coordinates object has
* succeeded in splitting underlying decimals in coordinates. If
* it does (return arrays with more than one element), we use
* that approach. Otherwise, we fallback on standard way which
* will try to determine manually the number of dimensions.
*/
Iterator<double[]> it = coords.points().iterator();
if (it.hasNext() && it.next().length > 1) {
return coords.points().map(GeometryTransformer::toCoordinate).spliterator();
} else {
values = coords.getValues();
}
} else if (source instanceof Point) {
DirectPosition dp = ((Point) source).getPos();
if (dp != null) {
return Stream.of(convertDirectPosition(dp)).spliterator();
} else
// recognized object, but no value. Empty geometry
return Spliterators.emptySpliterator();
}
}
// TODO : below conditions should be removed when proper abstraction is setup on GML geometry definition.
if (values == null) {
final boolean isLineString = source instanceof org.geotoolkit.gml.xml.LineString;
final boolean isLineStringSegment = source instanceof org.geotoolkit.gml.xml.LineStringSegment;
if (isLineString || isLineStringSegment) {
final DirectPositionList posList;
if (isLineString) {
posList = ((org.geotoolkit.gml.xml.LineString) source).getPosList();
} else {
posList = ((org.geotoolkit.gml.xml.LineStringSegment) source).getPosList();
}
if (posList != null) {
values = posList.getValue();
} else {
final List<? extends DirectPosition> pList;
if (isLineString) {
pList = ((org.geotoolkit.gml.xml.LineString) source).getPos();
} else {
pList = ((org.geotoolkit.gml.xml.LineStringSegment) source).getPos();
}
if (pList != null) {
return pList.stream().map(GeometryTransformer::convertDirectPosition).filter(Objects::nonNull).spliterator();
} else {
// We've identified a line, but there's no data in it
values = Collections.EMPTY_LIST;
}
}
} else if (source instanceof org.geotoolkit.gml.xml.LinearRing) {
// Note : do not check "getCoordinates", because it should have been done above.
values = asDoubles(() -> ((org.geotoolkit.gml.xml.LinearRing) source).getPosList());
} else if (source instanceof org.geotoolkit.gml.xml.v311.GeodesicStringType) {
values = asDoubles(() -> ((org.geotoolkit.gml.xml.v311.GeodesicStringType) source).getPosList());
} else if (source instanceof org.geotoolkit.gml.xml.v321.GeodesicStringType) {
values = asDoubles(() -> ((org.geotoolkit.gml.xml.v321.GeodesicStringType) source).getPosList());
} else if (source instanceof Curve) {
CurveSegmentArrayProperty segments = ((Curve) source).getSegments();
if (segments != null) {
final List<? extends AbstractCurveSegment> curveSegments = segments.getAbstractCurveSegment();
if (curveSegments != null) {
return curveSegments.stream().flatMap(seg -> StreamSupport.stream(getCoordinates(seg), false)).spliterator();
}
}
// If we arrive here, we've got an empty curve.
values = Collections.EMPTY_LIST;
} else if (source instanceof ArcByCenterPointType) {
final ArcByCenterPointType arc = (ArcByCenterPointType) source;
org.opengis.geometry.DirectPosition dp = arc.getPos();
if (dp == null) {
PointPropertyType pp = arc.getPointProperty();
if (pp == null) {
pp = arc.getPointRep();
}
if (pp == null) {
throw new UnconvertibleObjectException("Not enough information to build an arc.");
}
final Geometry point = new GeometryTransformer(pp.getPoint(), this).get();
dp = JTS.toEnvelope(point).getLowerCorner();
}
CoordinateReferenceSystem crs = dp.getCoordinateReferenceSystem();
if (crs == null) {
crs = getSrsName().map(this::findCRS).orElseThrow(() -> new UnconvertibleObjectException("Cannot create an arc without its coordinate reference system"));
final GeneralDirectPosition gdp = new GeneralDirectPosition(dp);
gdp.setCoordinateReferenceSystem(crs);
dp = gdp;
}
try {
final Measure startAngle, endAngle;
// If we miss a start or end angle, it means we're in presence of a circle.
if (arc.getStartAngle() == null || arc.getEndAngle() == null) {
startAngle = new Measure(0, Units.DEGREE);
endAngle = new Measure(360, Units.DEGREE);
} else {
startAngle = asMeasure(arc.getStartAngle());
endAngle = asMeasure(arc.getEndAngle());
}
final Coordinate[] coordinates = drawArc(dp, asMeasure(arc.getRadius()), startAngle, endAngle, ARC_PRECISION);
return Spliterators.spliterator(coordinates, Spliterator.ORDERED);
} catch (TransformException ex) {
throw new UnconvertibleObjectException("Cannot draw an arc.", ex);
}
} else if (source instanceof org.geotoolkit.gml.xml.v311.ArcByCenterPointType) {
// TODO : factorize with above case
final org.geotoolkit.gml.xml.v311.ArcByCenterPointType arc = (org.geotoolkit.gml.xml.v311.ArcByCenterPointType) source;
org.opengis.geometry.DirectPosition dp = arc.getPos();
if (dp == null) {
org.geotoolkit.gml.xml.v311.PointPropertyType pp = arc.getPointProperty();
if (pp == null) {
pp = arc.getPointRep();
}
if (pp == null) {
throw new UnconvertibleObjectException("Not enough information to build an arc.");
}
final Geometry point = new GeometryTransformer(pp.getPoint(), this).get();
dp = JTS.toEnvelope(point).getLowerCorner();
}
CoordinateReferenceSystem crs = dp.getCoordinateReferenceSystem();
if (crs == null) {
crs = getSrsName().map(this::findCRS).orElseThrow(() -> new UnconvertibleObjectException("Cannot create an arc without its coordinate reference system"));
final GeneralDirectPosition gdp = new GeneralDirectPosition(dp);
gdp.setCoordinateReferenceSystem(crs);
dp = gdp;
}
try {
final Measure startAngle, endAngle;
// If we miss a start or end angle, it means we're in presence of a circle.
if (arc.getStartAngle() == null || arc.getEndAngle() == null) {
startAngle = new Measure(0, Units.DEGREE);
endAngle = new Measure(360, Units.DEGREE);
} else {
startAngle = asMeasure(arc.getStartAngle());
endAngle = asMeasure(arc.getEndAngle());
}
final Coordinate[] coordinates = drawArc(dp, asMeasure(arc.getRadius()), startAngle, endAngle, ARC_PRECISION);
return Spliterators.spliterator(coordinates, Spliterator.ORDERED);
} catch (TransformException ex) {
throw new UnconvertibleObjectException("Cannot draw an arc.", ex);
}
}
}
if (values != null) {
if (values.isEmpty()) {
return Spliterators.emptySpliterator();
}
return new CoordinateSpliterator(values, getCoordinateDimension());
}
throw new UnconvertibleObjectException("Cannot extract coordinates from source geometry.");
}
Aggregations