本文整理了Java中org.geotools.geometry.jts.JTS类的一些代码示例，展示了JTS类的具体用法。这些代码示例主要来源于Github/Stackoverflow/Maven等平台，是从一些精选项目中提取出来的代码，具有较强的参考意义，能在一定程度帮忙到你。JTS类的具体详情如下：
包路径：org.geotools.geometry.jts.JTS
类名称：JTS
[英]JTS Geometry utility methods, bringing Geotools to JTS.
Offers geotools based services such as reprojection.
Responsibilities:
代码示例来源：origin: opentripplanner/OpenTripPlanner
CoordinateReferenceSystem crs = CRS.decode(crsCode);
MathTransform toMeters = CRS.findMathTransform(wgs, crs, false);
MathTransform fromMeters = CRS.findMathTransform(crs, wgs, false);
GeometryFactory geomf = JTSFactoryFinder.getGeometryFactory();
if (remainingSeconds > 60) { // avoid degenerate geometries
double remainingMeters = remainingSeconds * request.walkSpeed;
Geometry point = geomf.createPoint(vertexSeconds.getKey().getCoordinate());
point = JTS.transform(point, toMeters);
Geometry buffer = point.buffer(remainingMeters);
bufferLists.put(thresholdSeconds, buffer);
for (Integer threshold : bufferLists.keys()) {
Collection<Geometry> buffers = bufferLists.get(threshold);
Geometry geom = geomf.buildGeometry(buffers); // make geometry collection
geom = geom.union(); // combine all individual buffers in this contour into one
if ( ! resultsProjected) geom = JTS.transform(geom, fromMeters);
contours.put(threshold, geom);代码示例来源：origin: geoserver/geoserver
JTS.toGeometry((Envelope) new ReferencedEnvelope(grid.getEnvelope2D()));
if (coverageBounds.intersects(rasterFilter)) {
final ParameterValueGroup param = cropParams.clone();代码示例来源：origin: geotools/geotools
/**
* Creates a {@link CoordinateSequence} using the provided factory confirming the provided size
* and dimension are respected.
*
* <p>If the requested dimension is larger than the CoordinateSequence implementation can
* provide, then a sequence of maximum possible dimension should be created. An error should not
* be thrown.
*
* <p>This method is functionally identical to calling csFactory.create(size,dim) - it contains
* additional logic to work around a limitation on the commonly used
* CoordinateArraySequenceFactory.
*
* @param size the number of coordinates in the sequence
* @param dimension the dimension of the coordinates in the sequence
*/
public static CoordinateSequence createCS(
CoordinateSequenceFactory csFactory, int size, int dimension) {
// the coordinates don't have measures
return createCS(csFactory, size, dimension, 0);
}代码示例来源：origin: au.gov.amsa.risky/geo-analyzer
public static Point createPoint(double lat, double lon, String srsName) {
GeometryFactory geometryFactory = new GeometryFactory();
Coordinate coordinate = new Coordinate(lon, lat);
Point point = geometryFactory.createPoint(coordinate);
try {
if (!srsName.equals(EPSG_4326)) {
MathTransform transform = CRS.findMathTransform(getCrs(EPSG_4326), getCrs(srsName));
point = (Point) JTS.transform(point, transform);
}
return point;
} catch (NoSuchAuthorityCodeException e) {
throw new RuntimeException(e);
} catch (FactoryException e) {
throw new RuntimeException(e);
} catch (MismatchedDimensionException e) {
throw new RuntimeException(e);
} catch (TransformException e) {
throw new RuntimeException(e);
}
}代码示例来源：origin: org.geotools/gt-process-raster
throw new ProcessException("height and width parameters must be greater than 0");
if (bounds.getCoordinateReferenceSystem() == null) {
throw new ProcessException("Envelope CRS must not be null");
GeometryFactory geomFactory = new GeometryFactory();
try {
Polygon polygon = null;
CoordinateReferenceSystem targetCRS = CRS.parseWKT(targetCRSWKT);
MathTransform transform = CRS.findMathTransform(sourceCRS, targetCRS);
double pX = bounds.getMinX();
double pY = bounds.getMaxY();
double stepX = (bounds.getMaxX() - bounds.getMinX()) / width;
double stepY = (bounds.getMaxY() - bounds.getMinY()) / height;
tempCoordinates[2].x = nX; tempCoordinates[2].y = nY;
tempCoordinates[3].x = pX; tempCoordinates[3].y = nY;
polygon.geometryChanged();
Geometry targetGeometry = JTS.transform(polygon, transform);
matrix[i][j] = (float) targetGeometry.getArea();代码示例来源：origin: geotools/geotools
@Test
public void testUTMDatelineWrapping() throws Exception {
CoordinateReferenceSystem crs = CRS.decode("EPSG:32601", true);
ReferencedEnvelope re = new ReferencedEnvelope(300000, 409800, 5890200, 6000000, crs);
MathTransform mt = CRS.findMathTransform(crs, WGS84);
Geometry geom = JTS.toGeometry(re);
ReferencedEnvelope targetReferenceEnvelope =
new ReferencedEnvelope(-180, 180, -90, 90, WGS84);
ProjectionHandler ph =
ProjectionHandlerFinder.getHandler(targetReferenceEnvelope, crs, true);
Geometry preProcessed = ph.preProcess(geom);
Geometry transformed = JTS.transform(preProcessed, mt);
Geometry postProcessed = ph.postProcess(mt.inverse(), transformed);
// sits across the dateline and it's "small" (used to cover the entire planet)
Envelope ppEnvelope = postProcessed.getGeometryN(0).getEnvelopeInternal();
assertTrue(ppEnvelope.contains(180, 54));
// the original width is 109km, at this latitude one degree of longitude is only 65km
assertEquals(1.7, ppEnvelope.getWidth(), 0.1);
}代码示例来源：origin: org.geotools/gt-render
final CoordinateReferenceSystem tempCRS = CRS.getHorizontalCRS(coordinateReferenceSystem);
if(tempCRS==null)
throw new TransformException(Errors.format(
final Coordinate p1 = new Coordinate(envelope.getMinX(), envelope
.getMinY());
final Coordinate p2 = new Coordinate(envelope.getMaxX(), envelope
.getMaxY());
cs[0] = p1.x;
final MathTransform transform = CRS.findMathTransform(tempCRS,
DefaultGeographicCRS.WGS84, true);
transform.transform(cs, 0, csLatLong, 0, 2);
/ envelope.getHeight() * imageHeight;
double distance_ground = JTS.orthodromicDistance(new Coordinate(
newCsLatLong[0], newCsLatLong[1]), new Coordinate(
newCsLatLong[2], newCsLatLong[3]),
double diagonalGroundDistance = JTS.orthodromicDistance(new Coordinate(
csLatLong[0], csLatLong[1]), new Coordinate(csLatLong[2],
csLatLong[3]), DefaultGeographicCRS.WGS84);代码示例来源：origin: senbox-org/s2tbx
private ReferencedEnvelope transformEnvelope(CoordinateReferenceSystem fromCRS,
CoordinateReferenceSystem toCRS,
Rectangle2D bounds) {
try {
MathTransform transform = CRS.findMathTransform(fromCRS, toCRS);
Envelope sourceEnvelope = new Envelope(bounds.getMinX(), bounds.getMaxX(),
bounds.getMinY(), bounds.getMaxY());
final Envelope envelope = JTS.transform(sourceEnvelope, transform);
return new ReferencedEnvelope(envelope.getMinX(), envelope.getMaxX(),
envelope.getMinY(), envelope.getMaxY(),
toCRS);
} catch (FactoryException | TransformException e) {
e.printStackTrace();
return null;
}
}代码示例来源：origin: org.integratedmodelling/klab-engine
if (dataEnvelope != null) {
try {
geometry = JTS.transform(geometry, CRS
.findMathTransform(dataEnvelope.getCoordinateReferenceSystem(), ((Grid)extent).getCRS()));
} catch (Exception e) {
return;
for (int i = 0; i < geometry.getNumGeometries(); i++) {
Polygon poly = (Polygon) geometry.getGeometryN(i);
LinearRing lr = geoFactory.createLinearRing(poly.getExteriorRing().getCoordinates());
Polygon part = geoFactory.createPolygon(lr, null);
drawGeometry(part, false, swapAxis);
for (int j = 0; j < poly.getNumInteriorRing(); j++) {
lr = geoFactory.createLinearRing(poly.getInteriorRingN(j).getCoordinates());
part = geoFactory.createPolygon(lr, null);
drawGeometry(part, true, swapAxis);代码示例来源：origin: mapplus/spatial_statistics_for_geotools_udig
if (cropShape == null || cropShape.isEmpty()) {
roi = null;
} else {
if (geomBinding.isAssignableFrom(MultiPolygon.class)
|| geomBinding.isAssignableFrom(Polygon.class)) {
if (!cropShape.isValid() || !cropShape.isSimple()) {
cropShape = cropShape.buffer(0);
cropShape.setUserData(cropShape.getUserData());
roi = (GeometryCollection) cropShape;
} else {
roi = cropShape.getFactory().createGeometryCollection(new Geometry[] { cropShape });
if (extent == null || extent.isEmpty() || extent.isNull()) {
if (roi != null) {
bounds = new GeneralEnvelope(new ReferencedEnvelope(roi.getEnvelopeInternal(), crs));
CoordinateReferenceSystem extCrs = extent.getCoordinateReferenceSystem();
if (!CRS.equalsIgnoreMetadata(crs, extCrs)) {
try {
MathTransform transform = CRS.findMathTransform(extCrs, crs, true);
Envelope env = JTS.transform(extent, transform);
bounds = new GeneralEnvelope(new ReferencedEnvelope(env, crs));
} catch (FactoryException e) {代码示例来源：origin: geotools/geotools
@Test
public void testCutGeometryLambertConformal() throws Exception {
// get a lambert conformal conic with
ReferencedEnvelope wgs84South = new ReferencedEnvelope(-180, -90, -40, 0, WGS84);
ReferencedEnvelope laeSouth = wgs84South.transform(CRS.decode("EPSG:2194"), true);
ProjectionHandler handler = ProjectionHandlerFinder.getHandler(laeSouth, WGS84, true);
// a Geometry that crosses the opposite of the central meridian
Polygon geometry = JTS.toGeometry(new Envelope(5, 15, 0, 10));
Geometry preProcessed = handler.preProcess(geometry);
assertTrue(
"Should have sliced the geometry in two parts",
preProcessed instanceof MultiPolygon);
}代码示例来源：origin: org.geoserver/gwc
return JTS.toGeometry(new Envelope(prescribedBounds.getMinX(), prescribedBounds
.getMaxX(), prescribedBounds.getMinY(), prescribedBounds.getMaxY()));
final org.opengis.geometry.Envelope envelope = CRS.getEnvelope(targetCrs);
if (envelope == null) {
return null;
if (envelope.getSpan(0) < tolerance || envelope.getSpan(1) < tolerance) {
GeographicBoundingBox latLonBBox = CRS.getGeographicBoundingBox(targetCrs);
ReferencedEnvelope bbox = new ReferencedEnvelope(new GeneralEnvelope(latLonBBox));
Polygon geometry = JTS.toGeometry(bbox);
double distanceTolerance = Math.max(bbox.getSpan(0), bbox.getSpan(1)) / 2E5;
Geometry densifiedGeom = Densifier.densify(geometry, distanceTolerance);
MathTransform mathTransform;
CoordinateReferenceSystem sourceCRS = bbox.getCoordinateReferenceSystem();
mathTransform = CRS.findMathTransform(sourceCRS, targetCrs);
aovGeom = JTS.transform(densifiedGeom, mathTransform);
} catch (Exception e) {
throw new RuntimeException(e);
aovGeom = JTS.toGeometry(new Envelope(envelope.getMinimum(0), envelope.getMaximum(0),
envelope.getMinimum(1), envelope.getMaximum(1)));
aovGeom.setUserData(targetCrs);
return aovGeom;代码示例来源：origin: geotools/geotools
@Test
public void testWrapGeometryReprojectToLatLonED50() throws Exception {
ReferencedEnvelope world = new ReferencedEnvelope(-80, 80, -180, 180, ED50_LATLON);
// make sure the geometry is not wrapped, but it is preserved
ProjectionHandler handler = ProjectionHandlerFinder.getHandler(world, WGS84, true);
// a geometry that will cross the dateline and sitting in the same area as the
// rendering envelope (with wgs84 lon/latcoordinates)
String wkt = "POLYGON((178 -80, 178 80, 182 80, 182 80, 178 -80))";
Geometry g = new WKTReader().read(wkt);
Geometry original = new WKTReader().read(wkt);
MathTransform mt = CRS.findMathTransform(WGS84, ED50_LATLON);
MathTransform prepared =
handler.getRenderingTransform(CRS.findMathTransform(WGS84, ED50_LATLON));
Geometry reprojected = JTS.transform(original, prepared);
assertTrue(handler.requiresProcessing(g));
Geometry preProcessed = handler.preProcess(g);
// no cutting expected
assertEquals(original, preProcessed);
// post process, this should wrap the geometry and clone it
Geometry postProcessed = handler.postProcess(prepared, reprojected);
assertTrue(postProcessed instanceof MultiPolygon);
}代码示例来源：origin: org.geowebcache/gwc-georss
try {
Envelope coveredLevelEnvelope;
coveredLevelEnvelope = JTS.transform(aggregatedGeomBounds, worldToGrid);
Geometry bufferedEnvelopeInGridCrs;
bufferedEnvelopeInGridCrs = JTS.toGeometry(coveredLevelEnvelope).buffer(
ENVELOPE_BUFFER_RATIO);
coveredLevelEnvelope = bufferedEnvelopeInGridCrs.getEnvelopeInternal();
MathTransform gridToWorld = worldToGrid.inverse();
Envelope bufferedEnvelope = JTS.transform(coveredLevelEnvelope, gridToWorld);
expandedBounds = new BoundingBox(bufferedEnvelope.getMinX(),
bufferedEnvelope.getMinY(), bufferedEnvelope.getMaxX(),
bufferedEnvelope.getMaxY());
} catch (TransformException e) {代码示例来源：origin: org.geotools/gt-widgets-swing-pending
GeometryFactory geofact = new GeometryFactory();
Coordinate coord = toMapCoord(event.getX(), event.getY(), env.getWidth(), env.getHeight(), rec,env);
Point point = geofact.createPoint(coord);
MathTransform transform = CRS.findMathTransform(sourceCRS, targetCRS);
Geometry targetGeometry = JTS.transform( point, transform);
coord = targetGeometry.getCoordinate();
jtf_coord.setText("X= "+ coord.x +" / Y= "+coord.y);
}catch(Exception e){代码示例来源：origin: iandees/shp-to-osm
transform = CRS.findMathTransform(sourceCRS, targetCRS, true);
} catch (MalformedURLException e) {
throw new ShpToOsmException("URL could not be created for input file.", e);
String geometryType = rawGeom.getGeometryType();
geometry = JTS.transform(rawGeom, transform);
} catch (TransformException e) {
throw new ShpToOsmException("Could not transform to spherical mercator.", e);
for (int i = 0; i < geometry.getNumGeometries(); i++) {
LineString geometryN = (LineString) geometry
.getGeometryN(i);
LineString outerLine = geometryN.getExteriorRing();
if (geometryN.getNumInteriorRing() > 0) {
Relation r = new Relation();
r.addTag(new Tag("type", "multipolygon"));
for (int j = 0; j < geometryN.getNumInteriorRing(); j++) {
LineString innerLine = geometryN.getInteriorRingN(j);代码示例来源：origin: org.geotools/gt-main
/**
* Helper method to reproject a geometry.
*/
protected Geometry reproject( Object value, CoordinateReferenceSystem propertyCrs) {
if ( value == null ) {
return null;
}
if (!(value instanceof Geometry))
throw new IllegalArgumentException("Binary geometry filter, but second expression "
+ "is not a geometry literal? (it's a " + value.getClass() + ")");
Geometry geom = (Geometry) value;
// does it make sense to proceed?
if (geom.getUserData() == null
|| !(geom.getUserData() instanceof CoordinateReferenceSystem))
return geom;
try {
// reproject
CoordinateReferenceSystem geomCRS = (CoordinateReferenceSystem) geom.getUserData();
Geometry transformed = JTS.transform(geom, CRS.findMathTransform(geomCRS, propertyCrs, true));
transformed.setUserData(propertyCrs);
return transformed;
} catch(Exception e) {
throw new RuntimeException("Could not reproject geometry " + value, e);
}
}代码示例来源：origin: org.geotools/gt-render
if (CRS.equalsIgnoreMetadata(geomCRS, renderingEnvelope.getCoordinateReferenceSystem())) {
return geometry;
ReferencedEnvelope ge = new ReferencedEnvelope(geometry.getEnvelopeInternal(), geomCRS);
ReferencedEnvelope geWGS84 = ge.transform(WGS84, true);
if (validArea.contains((Envelope) geWGS84)) {
return geometry;
Polygon envelopeGeometry = JTS.toGeometry((Envelope) envInt);
Geometry result;
try {
result = geometry.intersection(envelopeGeometry);
} catch(Exception e1) {
try {代码示例来源：origin: org.locationtech.geogig/geogig-datastore
Polygon geometry = JTS.toGeometry(queryBounds);
filter = ff.intersects(ff.property(pointsType.getGeometryDescriptor().getLocalName()),
ff.literal(geometry));
CoordinateReferenceSystem queryCrs = CRS.decode("EPSG:3857");
transformedQueryBounds = queryBounds.transform(queryCrs, true);
geometry = JTS.toGeometry(transformedQueryBounds);
geometry.setUserData(queryCrs);代码示例来源：origin: org.locationtech.geogig/geogig-core
public Envelope getQuadCenter(Quadrant... quadrants) {
Envelope quadBounds = getMaxBounds();
if (quadrants != null) {
for (Quadrant quad : quadrants) {
quadBounds = quad.slice(quadBounds);
}
}
Coordinate center = quadBounds.centre();
Envelope nodeBounds = Node.makePrecise(new Envelope(center));
if (!quadBounds.contains(nodeBounds)) {
GeometryFactory gf = new GeometryFactory();
Polygon qgeom = JTS.toGeometry(quadBounds, gf);
ArrayList<Geometry> pointGeoms = Lists.newArrayList(gf.createPoint(center),
JTS.toGeometry(nodeBounds, gf));
Geometry point = gf.buildGeometry(pointGeoms);
String msg = qgeom + " does not contain " + point;
Assert.fail(msg);
}
return nodeBounds;
}内容来源于网络，如有侵权，请联系作者删除！