/* (c) 2014 Open Source Geospatial Foundation - all rights reserved
* (c) 2001 - 2013 OpenPlans
* This code is licensed under the GPL 2.0 license, available at the root
* application directory.
*/
package org.geoserver.test;
import static org.junit.Assert.*;
import org.junit.Test;
import org.geoserver.data.test.SystemTestData;
import org.geotools.factory.Hints;
import org.geotools.geometry.jts.JTS;
import org.geotools.referencing.CRS;
import org.geotools.referencing.CRS.AxisOrder;
import org.opengis.geometry.MismatchedDimensionException;
import org.opengis.referencing.FactoryException;
import org.opengis.referencing.NoSuchAuthorityCodeException;
import org.opengis.referencing.crs.CoordinateReferenceSystem;
import org.opengis.referencing.operation.MathTransform;
import org.opengis.referencing.operation.TransformException;
import org.w3c.dom.Document;
import com.vividsolutions.jts.geom.Coordinate;
import com.vividsolutions.jts.geom.GeometryFactory;
import com.vividsolutions.jts.geom.Point;
import com.vividsolutions.jts.geom.Polygon;
/**
* This is to test encoding of SRS information and reprojection values in app-schema features.
*
* @author Rini Angreani, Curtin University of Technology
*/
public class SRSWfsTest extends AbstractAppSchemaTestSupport {
final String EPSG_4326 = "urn:x-ogc:def:crs:EPSG:4326";
final String EPSG_4283 = "urn:x-ogc:def:crs:EPSG:4283";
final String DIMENSION = "2";
@Override
protected SRSMockData createTestData() {
return new SRSMockData();
}
@Override
protected void setUpTestData(SystemTestData testData) throws Exception {
// use the OGC standard for axis order
//
// must be done *before* super.oneTimeSetUp() to ensure CRS factories
// configured before data is loaded
//
// if this property is null, GeoServerAbstractTestSupport.oneTimeSetUp()
// will blow away our changes
System.setProperty("org.geotools.referencing.forceXY", "false");
// yes, we need this too
Hints.putSystemDefault(Hints.FORCE_LONGITUDE_FIRST_AXIS_ORDER, false);
// if this is set to anything but "http", GeoServerAbstractTestSupport.oneTimeSetUp()
// will blow away our changes
Hints.putSystemDefault(Hints.FORCE_AXIS_ORDER_HONORING, "http");
// apply changes
CRS.reset("all");
super.setUpTestData(testData);
}
@Override
protected void onTearDown(SystemTestData testData) throws Exception {
super.onTearDown(testData);
// undo the changes made for this suite and reset
System.clearProperty("org.geotools.referencing.forceXY");
Hints.removeSystemDefault(Hints.FORCE_LONGITUDE_FIRST_AXIS_ORDER);
Hints.removeSystemDefault(Hints.FORCE_AXIS_ORDER_HONORING);
CRS.reset("all");
}
/**
* Test content of GetFeature response.
*/
@Test
public void testGetFeatureContent() throws NoSuchAuthorityCodeException, FactoryException {
String id = "1";
// make sure we are really working in lat/lon order
CoordinateReferenceSystem crs = CRS.decode("EPSG:4326");
assertEquals(AxisOrder.NORTH_EAST, CRS.getAxisOrder(crs));
Document doc = getAsDOM("wfs?request=GetFeature&typename=ex:geomContainer&version=1.1.0");
LOGGER.info("WFS GetFeature&typename=ex:geomContainer response:\n" + prettyString(doc));
assertXpathEvaluatesTo("2", "/wfs:FeatureCollection/@numberOfFeatures", doc);
assertXpathCount(2, "//ex:geomContainer", doc);
// 1st feature
assertXpathEvaluatesTo(id, "(//ex:geomContainer)[1]/@gml:id", doc);
assertXpathEvaluatesTo("1.0", "(//ex:geomContainer)[1]/ex:length", doc);
assertXpathEvaluatesTo("m", "(//ex:geomContainer)[1]/ex:length/@uom", doc);
// check gml:id for geometries
assertXpathEvaluatesTo("geom_" + id, "(//ex:geomContainer)[1]/ex:geom/gml:Polygon/@gml:id",
doc);
// check srs properties
assertXpathEvaluatesTo(EPSG_4283, "(//ex:geomContainer)[1]/ex:geom/gml:Polygon/@srsName",
doc);
assertXpathEvaluatesTo(DIMENSION,
"(//ex:geomContainer)[1]/ex:geom/gml:Polygon/@srsDimension", doc);
// test geometry values
assertXpathEvaluatesTo(
"-1.2 52.5 -1.2 52.6 -1.1 52.6 -1.1 52.5 -1.2 52.5",
"(//ex:geomContainer)[1]/ex:geom/gml:Polygon/gml:exterior/gml:LinearRing/gml:posList",
doc);
// test second geometry (multi geometry feature) is encoded correctly, with srs info
assertXpathEvaluatesTo("shape_" + id,
"(//ex:geomContainer)[1]/ex:shape/gml:LineString/@gml:id", doc);
assertXpathEvaluatesTo(EPSG_4283,
"(//ex:geomContainer)[1]/ex:shape/gml:LineString/@srsName", doc);
assertXpathEvaluatesTo(DIMENSION,
"(//ex:geomContainer)[1]/ex:shape/gml:LineString/@srsDimension", doc);
assertXpathEvaluatesTo("-1.2 52.5 -1.2 52.6 -1.1 52.6 -1.1 52.5 -1.2 52.5",
"(//ex:geomContainer)[1]/ex:shape/gml:LineString/gml:posList", doc);
// test nested geometry via feature chaining
assertXpathEvaluatesTo("nested_one_2",
"(//ex:geomContainer)[1]/ex:nestedFeature[1]/ex:nestedGeom/@gml:id", doc);
assertXpathEvaluatesTo(
"nested_one_geom_2",
"(//ex:geomContainer)[1]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Point/@gml:id",
doc);
assertXpathEvaluatesTo(
EPSG_4283,
"(//ex:geomContainer)[1]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Point/@srsName",
doc);
assertXpathEvaluatesTo(
DIMENSION,
"(//ex:geomContainer)[1]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Point/@srsDimension",
doc);
assertXpathEvaluatesTo(
"42.58 31.29",
"(//ex:geomContainer)[1]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Point/gml:pos",
doc);
// test nested inline geometry
assertXpathEvaluatesTo("nested_two_" + id,
"(//ex:geomContainer)[1]/ex:nestedFeature[2]/ex:nestedGeom/@gml:id", doc);
assertXpathEvaluatesTo(
"nested_two_geom_" + id,
"(//ex:geomContainer)[1]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Polygon/@gml:id",
doc);
assertXpathEvaluatesTo(
EPSG_4283,
"(//ex:geomContainer)[1]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Polygon/@srsName",
doc);
assertXpathEvaluatesTo(
DIMENSION,
"(//ex:geomContainer)[1]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Polygon/@srsDimension",
doc);
assertXpathEvaluatesTo(
"-1.2 52.5 -1.2 52.6 -1.1 52.6 -1.1 52.5 -1.2 52.5",
"(//ex:geomContainer)[1]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Polygon/gml:exterior/gml:LinearRing/gml:posList",
doc);
// second feature
id = "2";
assertXpathEvaluatesTo(id, "(//ex:geomContainer)[2]/@gml:id", doc);
assertXpathEvaluatesTo("1.0", "(//ex:geomContainer)[2]/ex:length", doc);
assertXpathEvaluatesTo("m", "(//ex:geomContainer)[2]/ex:length/@uom", doc);
// check gml:id for geometries
assertXpathEvaluatesTo("geom_" + id, "(//ex:geomContainer)[2]/ex:geom/gml:Point/@gml:id",
doc);
// check srs properties
assertXpathEvaluatesTo(EPSG_4283, "(//ex:geomContainer)[2]/ex:geom/gml:Point/@srsName", doc);
assertXpathEvaluatesTo(DIMENSION,
"(//ex:geomContainer)[2]/ex:geom/gml:Point/@srsDimension", doc);
// test geometry values
assertXpathEvaluatesTo("42.58 31.29", "(//ex:geomContainer)[2]/ex:geom/gml:Point/gml:pos",
doc);
// test second geometry (multi geometry feature) is encoded correctly, with srs info
assertXpathEvaluatesTo("shape_" + id,
"(//ex:geomContainer)[2]/ex:shape/gml:LineString/@gml:id", doc);
assertXpathEvaluatesTo(EPSG_4283,
"(//ex:geomContainer)[2]/ex:shape/gml:LineString/@srsName", doc);
assertXpathEvaluatesTo(DIMENSION,
"(//ex:geomContainer)[2]/ex:shape/gml:LineString/@srsDimension", doc);
assertXpathEvaluatesTo("42.58 31.29 42.58 31.29",
"(//ex:geomContainer)[2]/ex:shape/gml:LineString/gml:posList", doc);
// test nested geometry via feature chaining
assertXpathEvaluatesTo("nested_one_1",
"(//ex:geomContainer)[2]/ex:nestedFeature[1]/ex:nestedGeom/@gml:id", doc);
assertXpathEvaluatesTo(
"nested_one_geom_1",
"(//ex:geomContainer)[2]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Polygon/@gml:id",
doc);
assertXpathEvaluatesTo(
EPSG_4283,
"(//ex:geomContainer)[2]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Polygon/@srsName",
doc);
assertXpathEvaluatesTo(
DIMENSION,
"(//ex:geomContainer)[2]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Polygon/@srsDimension",
doc);
assertXpathEvaluatesTo(
"-1.2 52.5 -1.2 52.6 -1.1 52.6 -1.1 52.5 -1.2 52.5",
"(//ex:geomContainer)[2]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Polygon/gml:exterior/gml:LinearRing/gml:posList",
doc);
// test nested inline geometry
assertXpathEvaluatesTo("nested_two_" + id,
"(//ex:geomContainer)[2]/ex:nestedFeature[2]/ex:nestedGeom/@gml:id", doc);
assertXpathEvaluatesTo(
"nested_two_geom_" + id,
"(//ex:geomContainer)[2]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Point/@gml:id",
doc);
assertXpathEvaluatesTo(
EPSG_4283,
"(//ex:geomContainer)[2]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Point/@srsName",
doc);
assertXpathEvaluatesTo(
DIMENSION,
"(//ex:geomContainer)[2]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Point/@srsDimension",
doc);
assertXpathEvaluatesTo(
"42.58 31.29",
"(//ex:geomContainer)[2]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Point/gml:pos",
doc);
}
/**
* Test SRS reprojection. Ensure both geometry values and SRS information are encoded correctly.
*/
@Test
public void testReproject() throws NoSuchAuthorityCodeException, FactoryException,
MismatchedDimensionException, TransformException {
// reprojected geometries
CoordinateReferenceSystem sourceCRS = (CoordinateReferenceSystem) CRS.decode(EPSG_4283);
CoordinateReferenceSystem targetCRS = (CoordinateReferenceSystem) CRS.decode(EPSG_4326);
MathTransform transform = CRS.findMathTransform(sourceCRS, targetCRS);
GeometryFactory factory = new GeometryFactory();
Polygon srcPolygon = factory.createPolygon(factory.createLinearRing(factory
.getCoordinateSequenceFactory().create(
new Coordinate[] { new Coordinate(-1.2, 52.5), new Coordinate(-1.2, 52.6),
new Coordinate(-1.1, 52.6), new Coordinate(-1.1, 52.5),
new Coordinate(-1.2, 52.5) })), null);
Polygon targetPolygon = (Polygon) JTS.transform(srcPolygon, transform);
StringBuffer polygonBuffer = new StringBuffer();
for (Coordinate coord : targetPolygon.getCoordinates()) {
polygonBuffer.append(coord.x).append(" ");
polygonBuffer.append(coord.y).append(" ");
}
String targetPolygonCoords = polygonBuffer.toString().trim();
Point targetPoint = (Point) JTS.transform(
factory.createPoint(new Coordinate(42.58, 31.29)), transform);
String targetPointCoord = targetPoint.getCoordinate().x + " "
+ targetPoint.getCoordinate().y;
Document doc = getAsDOM("wfs?request=GetFeature&version=1.1.0&typename=ex:geomContainer" +
"&srsname=urn:x-ogc:def:crs:EPSG::4326");
LOGGER
.info("WFS GetFeature&typename=ex:geomContainer&srsname=urn:x-ogc:def:crs:EPSG::4326 response:\n"
+ prettyString(doc));
assertXpathEvaluatesTo("2", "/wfs:FeatureCollection/@numberOfFeatures", doc);
assertXpathCount(2, "//ex:geomContainer", doc);
// 1st feature
String id = "1";
assertXpathEvaluatesTo(id, "(//ex:geomContainer)[1]/@gml:id", doc);
assertXpathEvaluatesTo("1.0", "(//ex:geomContainer)[1]/ex:length", doc);
assertXpathEvaluatesTo("m", "(//ex:geomContainer)[1]/ex:length/@uom", doc);
// test gml:id for geometries
assertXpathEvaluatesTo("geom_" + id, "(//ex:geomContainer)[1]/ex:geom/gml:Polygon/@gml:id",
doc);
// check srs properties
assertXpathEvaluatesTo(EPSG_4326, "//ex:geomContainer[1]/ex:geom/gml:Polygon/@srsName", doc);
assertXpathEvaluatesTo(DIMENSION,
"(//ex:geomContainer)[1]/ex:geom/gml:Polygon/@srsDimension", doc);
// test values
assertXpathEvaluatesTo(
targetPolygonCoords,
"(//ex:geomContainer)[1]/ex:geom/gml:Polygon/gml:exterior/gml:LinearRing/gml:posList",
doc);
// test second geometry (multi geometry feature) is encoded correctly, with srs info
assertXpathEvaluatesTo("shape_" + id,
"(//ex:geomContainer)[1]/ex:shape/gml:LineString/@gml:id", doc);
assertXpathEvaluatesTo(EPSG_4326,
"(//ex:geomContainer)[1]/ex:shape/gml:LineString/@srsName", doc);
assertXpathEvaluatesTo(DIMENSION,
"(//ex:geomContainer)[1]/ex:shape/gml:LineString/@srsDimension", doc);
assertXpathEvaluatesTo(targetPolygonCoords,
"(//ex:geomContainer)[1]/ex:shape/gml:LineString/gml:posList", doc);
// test nested geometry via feature chaining
assertXpathEvaluatesTo("nested_one_2",
"(//ex:geomContainer)[1]/ex:nestedFeature[1]/ex:nestedGeom/@gml:id", doc);
assertXpathEvaluatesTo(
"nested_one_geom_2",
"(//ex:geomContainer)[1]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Point/@gml:id",
doc);
assertXpathEvaluatesTo(
EPSG_4326,
"(//ex:geomContainer)[1]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Point/@srsName",
doc);
assertXpathEvaluatesTo(
DIMENSION,
"(//ex:geomContainer)[1]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Point/@srsDimension",
doc);
assertXpathEvaluatesTo(
targetPointCoord,
"(//ex:geomContainer)[1]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Point/gml:pos",
doc);
// test nested inline geometry
assertXpathEvaluatesTo("nested_two_" + id,
"(//ex:geomContainer)[1]/ex:nestedFeature[2]/ex:nestedGeom/@gml:id", doc);
assertXpathEvaluatesTo(
"nested_two_geom_" + id,
"(//ex:geomContainer)[1]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Polygon/@gml:id",
doc);
assertXpathEvaluatesTo(
EPSG_4326,
"(//ex:geomContainer)[1]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Polygon/@srsName",
doc);
assertXpathEvaluatesTo(
DIMENSION,
"(//ex:geomContainer)[1]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Polygon/@srsDimension",
doc);
assertXpathEvaluatesTo(
targetPolygonCoords,
"(//ex:geomContainer)[1]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Polygon/gml:exterior/gml:LinearRing/gml:posList",
doc);
// second feature
id = "2";
assertXpathEvaluatesTo(id, "(//ex:geomContainer)[2]/@gml:id", doc);
assertXpathEvaluatesTo("1.0", "(//ex:geomContainer)[2]/ex:length", doc);
assertXpathEvaluatesTo("m", "(//ex:geomContainer)[2]/ex:length/@uom", doc);
// test gml:id for geometries
assertXpathEvaluatesTo("geom_" + id, "(//ex:geomContainer)[2]/ex:geom/gml:Point/@gml:id",
doc);
// check srs properties
assertXpathEvaluatesTo(EPSG_4326, "(//ex:geomContainer)[2]/ex:geom/gml:Point/@srsName", doc);
assertXpathEvaluatesTo(DIMENSION,
"(//ex:geomContainer)[2]/ex:geom/gml:Point/@srsDimension", doc);
// test values
assertXpathEvaluatesTo(targetPointCoord,
"(//ex:geomContainer)[2]/ex:geom/gml:Point/gml:pos", doc);
// test second geometry (multi geometry feature) is encoded correctly, with srs info
assertXpathEvaluatesTo("shape_" + id,
"(//ex:geomContainer)[2]/ex:shape/gml:LineString/@gml:id", doc);
assertXpathEvaluatesTo(EPSG_4326,
"(//ex:geomContainer)[2]/ex:shape/gml:LineString/@srsName", doc);
assertXpathEvaluatesTo(DIMENSION,
"(//ex:geomContainer)[2]/ex:shape/gml:LineString/@srsDimension", doc);
assertXpathEvaluatesTo(targetPointCoord + " " + targetPointCoord,
"(//ex:geomContainer)[2]/ex:shape/gml:LineString/gml:posList", doc);
// test nested geometry via feature chaining
assertXpathEvaluatesTo("nested_one_1",
"(//ex:geomContainer)[2]/ex:nestedFeature[1]/ex:nestedGeom/@gml:id", doc);
assertXpathEvaluatesTo(
"nested_one_geom_1",
"(//ex:geomContainer)[2]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Polygon/@gml:id",
doc);
assertXpathEvaluatesTo(
EPSG_4326,
"(//ex:geomContainer)[2]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Polygon/@srsName",
doc);
assertXpathEvaluatesTo(
DIMENSION,
"(//ex:geomContainer)[2]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Polygon/@srsDimension",
doc);
assertXpathEvaluatesTo(
targetPolygonCoords,
"(//ex:geomContainer)[2]/ex:nestedFeature[1]/ex:nestedGeom/ex:geom/gml:Polygon/gml:exterior/gml:LinearRing/gml:posList",
doc);
// test nested inline geometry
assertXpathEvaluatesTo("nested_two_" + id,
"(//ex:geomContainer)[2]/ex:nestedFeature[2]/ex:nestedGeom/@gml:id", doc);
assertXpathEvaluatesTo(
"nested_two_geom_" + id,
"(//ex:geomContainer)[2]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Point/@gml:id",
doc);
assertXpathEvaluatesTo(
EPSG_4326,
"(//ex:geomContainer)[2]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Point/@srsName",
doc);
assertXpathEvaluatesTo(
DIMENSION,
"(//ex:geomContainer)[2]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Point/@srsDimension",
doc);
assertXpathEvaluatesTo(
targetPointCoord,
"(//ex:geomContainer)[2]/ex:nestedFeature[2]/ex:nestedGeom/ex:geom/gml:Point/gml:pos",
doc);
}
/**
* Ensure filters are still working.
*/
@Test
public void testFilters() {
Document doc = getAsDOM("wfs?request=GetFeature&version=1.1.0&typename=ex:geomContainer" +
"&srsname=urn:x-ogc:def:crs:EPSG::4326&featureid=1");
LOGGER
.info("WFS GetFeature&typename=ex:geomContainer&srsname=urn:x-ogc:def:crs:EPSG::4326&featureid=1"
+ "response:\n" + prettyString(doc));
assertXpathEvaluatesTo("1", "/wfs:FeatureCollection/@numberOfFeatures", doc);
assertXpathCount(1, "//ex:geomContainer", doc);
}
}