dive4elements/river: flys-artifacts/src/main/java/de/intevation/flys/utils/GeometryUtils.java comparison

comparison flys-artifacts/src/main/java/de/intevation/flys/utils/GeometryUtils.java @ 3301:c8f670ae96e8

Create extent of a map layer in the correct coordinate reference system. flys-artifacts/trunk@4986 c6561f87-3c4e-4783-a992-168aeb5c3f6f

author	Ingo Weinzierl <ingo.weinzierl@intevation.de>
date	Fri, 13 Jul 2012 10:25:24 +0000
parents	f062b5a90e26
children	96daa8a4c48f

comparison

equal deleted inserted replaced

-:e1cf76b3ecb4
+:c8f670ae96e8
 import org.geotools.feature.FeatureCollection;
 import org.geotools.feature.FeatureIterator;
 import org.geotools.feature.simple.SimpleFeatureTypeBuilder;
 import org.geotools.geojson.feature.FeatureJSON;
 import org.geotools.geometry.jts.JTS;
+import org.geotools.geometry.jts.ReferencedEnvelope;
 import org.geotools.referencing.CRS;
 import org.opengis.feature.simple.SimpleFeature;
 import org.opengis.feature.simple.SimpleFeatureType;
 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;
 public class GeometryUtils {
 private static final Logger logger = Logger.getLogger(GeometryUtils.class);
+public static final String PREFIX_EPSG  = "EPSG:";
 public static final String DEFAULT_EPSG = "EPSG:31467";
 private GeometryUtils() {
 }
 public static Envelope getRiverBoundary(String rivername) {
 List<RiverAxis> axes = RiverAxis.getRiverAxis(rivername);
 if (axes != null && axes.size() > 0) {
 Envelope max = null;
 }
 return null;
 }
 public static String getRiverBounds(String rivername) {
 Envelope env = getRiverBoundary(rivername);
 if (env == null) {
 return jtsBoundsToOLBounds(env);
 }
 return null;
 }
 /**
-* Returns the boundary of Envelope <i>env</i> in OpenLayers
+* Returns the boundary of Envelope <i>env</i> in OpenLayers representation.
-* representation.
 *
 * @param env The envelope of a geometry.
 *
 * @return the OpenLayers boundary of <i>env</i>.
 */
 env.getMinY() + " " +
 env.getMaxX() + " " +
 env.getMaxY();
 }
 public static String createOLBounds(Geometry a, Geometry b) {
 Coordinate[] ca = a.getCoordinates();
 Coordinate[] cb = b.getCoordinates();
 double lowerX = Double.MAX_VALUE;
 }
 return "" + lowerX + " " + lowerY + " " + upperX + " " + upperY;
 }
 public static SimpleFeatureType buildFeatureType(
 String name, String srs, Class<?> geometryType)
 {
 return buildFeatureType(name, srs, geometryType, null);
 }
 /**
 * Creates a new SimpleFeatureType using a SimpleFeatureTypeBuilder.
 *
 * @param name The name of the FeatureType.
 * @param srs The SRS (e.g. "EPSG:31466").
 * @param geometryType The geometry type's class (e.g. Polygon.class).
-* @param attrs Optional. An object with attribute-name/attribute-class
+* @param attrs Optional. An object with attribute-name/attribute-class pairs
-* pairs where index 0 specifies the name as string and index 1 the type
+* where index 0 specifies the name as string and index 1 the
-* as class.
+* ype as class.
 *
 * @return a new SimpleFeatureType.
 */
-public static SimpleFeatureType buildFeatureType(
+public static SimpleFeatureType buildFeatureType(String name, String srs,
-String name, String srs, Class<?> geometryType, Object[][] attrs)
+Class<?> geometryType, Object[][] attrs) {
-{
 try {
 SimpleFeatureTypeBuilder builder = new SimpleFeatureTypeBuilder();
 CoordinateReferenceSystem crs    = CRS.decode(srs);
 builder.setName("flys");
 }
 return null;
 }
 public static List<SimpleFeature> parseGeoJSON(
 String geojson, SimpleFeatureType ft
 ) {
 List<SimpleFeature> collection = new ArrayList<SimpleFeature>();
 return collection;
 }
-public static boolean writeShapefile(
+/**
-File              shape,
+* This method returns the {@link CoordinateReferenceSystem} by the
-SimpleFeatureType featureType,
+* {@link String} <i>epsg</i>.
-FeatureCollection collection
+*
+* @param epsg An EPSG code like <b>EPSG:31466</b>
+*
+* @return the {@link CoordinateReferenceSystem} specified by <i>epsg</i>.
+*/
+public static CoordinateReferenceSystem getCoordinateReferenceSystem(
+String epsg
 ) {
-return writeShapefile(
+if (epsg == null) {
-shape, featureType, collection, featureType.getCoordinateReferenceSystem());
+logger.warn("cannot create CoordinateReferenceSystem with null");
-}
+return null;
+}
-public static boolean writeShapefile(
+if (!epsg.startsWith(PREFIX_EPSG)) {
-File                      shape,
+epsg = PREFIX_EPSG + epsg;
-SimpleFeatureType         featureType,
+}
-FeatureCollection         collection,
-CoordinateReferenceSystem crs
+try {
+return CRS.decode(epsg);
+}
+catch (FactoryException fe) {
+logger.error(
+"unable to get CoordinateReferenceSystem for: " + epsg,
+fe);
+}
+return null;
+}
+public static Envelope transform(Envelope orig, String targetSrs) {
+return transform(orig, targetSrs, DEFAULT_EPSG);
+}
+public static Envelope transform(
+Envelope orig,
+String   targetSrs,
+String   origSrs
 ) {
+if (targetSrs == null || orig == null || origSrs == null) {
+logger.warn("unable to transform envelope: empty parameters");
+return orig;
+}
+logger.debug("Transform envlope to '" + targetSrs + "'");
+try {
+CoordinateReferenceSystem sourceCRS =
+getCoordinateReferenceSystem(origSrs);
+CoordinateReferenceSystem targetCRS =
+getCoordinateReferenceSystem(targetSrs);
+if (sourceCRS != null && targetCRS != null) {
+ReferencedEnvelope tmpEnv =
+new ReferencedEnvelope(orig, CRS.decode(origSrs));
+Envelope target = tmpEnv.transform(targetCRS, false);
+if (logger.isDebugEnabled()) {
+logger.debug("   orig envelope       : " + orig);
+logger.debug("   transformed envelope: " + target);
+}
+return target;
+}
+}
+catch (NoSuchAuthorityCodeException nsae) {
+logger.error("Cannot get CoordinateReferenceSystem!", nsae);
+}
+catch (FactoryException fe) {
+logger.error("Cannot get CoordinateReferenceSystem!", fe);
+}
+catch (TransformException te) {
+logger.error("Cannot transform envelope from source "
++ origSrs + " to target srs " + targetSrs);
+}
+return null;
+}
+public static boolean writeShapefile(File shape,
+SimpleFeatureType featureType, FeatureCollection collection) {
+return writeShapefile(shape, featureType, collection,
+featureType.getCoordinateReferenceSystem());
+}
+public static boolean writeShapefile(File shape,
+SimpleFeatureType featureType, FeatureCollection collection,
+CoordinateReferenceSystem crs) {
 if (collection.isEmpty()) {
 logger.warn("Shapefile is not written - no features given!");
 return false;
 }

Mercurial > dive4elements > river

comparison flys-artifacts/src/main/java/de/intevation/flys/utils/GeometryUtils.java @ 3301:c8f670ae96e8