/* Spatial Operations & Editing Tools for uDig
*
* Axios Engineering under a funding contract with:
* Diputación Foral de Gipuzkoa, Ordenación Territorial
*
* http://b5m.gipuzkoa.net
* http://www.axios.es
*
* (C) 2006, Diputación Foral de Gipuzkoa, Ordenación Territorial (DFG-OT).
* DFG-OT agrees to license under Lesser General Public License (LGPL).
*
* You can redistribute it and/or modify it under the terms of the
* GNU Lesser General Public License as published by the Free Software
* Foundation; version 2.1 of the License.
*
* This library is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* Lesser General Public License for more details.
*/
package es.axios.geotools.util;
import java.util.HashMap;
import java.util.Map;
import javax.measure.unit.Unit;
import org.geotools.geometry.jts.CoordinateSequenceTransformer;
import org.geotools.geometry.jts.GeometryCoordinateSequenceTransformer;
import org.geotools.referencing.CRS;
import org.geotools.resources.CRSUtilities;
import org.opengis.feature.IllegalAttributeException;
import org.opengis.feature.simple.SimpleFeature;
import org.opengis.feature.simple.SimpleFeatureType;
import org.opengis.feature.type.AttributeDescriptor;
import org.opengis.feature.type.GeometryDescriptor;
import org.opengis.referencing.FactoryException;
import org.opengis.referencing.crs.CoordinateReferenceSystem;
import org.opengis.referencing.cs.CoordinateSystem;
import org.opengis.referencing.cs.CoordinateSystemAxis;
import org.opengis.referencing.operation.MathTransform;
import org.opengis.referencing.operation.OperationNotFoundException;
import org.opengis.referencing.operation.TransformException;
import com.vividsolutions.jts.geom.Coordinate;
import com.vividsolutions.jts.geom.CoordinateSequenceFactory;
import com.vividsolutions.jts.geom.Geometry;
import com.vividsolutions.jts.geom.GeometryFactory;
import com.vividsolutions.jts.geom.LineSegment;
import es.axios.lib.geometry.util.GeometryUtil;
/**
* GeoTools Utils
* <p>
* This class has util and convenient methods to work with GeoTools
* </p>
*
*
* @author Mauricio Pazos (www.axios.es)
* @author Aritz Davila (www.axios.es)
* @since 1.1.0
*/
public class GeoToolsUtils {
private GeoToolsUtils() {
// util class
}
/**
* Adds the matching attributes from <code>source</code> to
* <code>target</code>
* <p>
* Two attributes match if they have the same name and type.
* </p>
*
* @param source
* @param target
* @throws IllegalAttributeException
*/
public static void copyAttributes(SimpleFeature source, SimpleFeature target) throws IllegalAttributeException {
Map<String, Class<?>> sourceTypes = new HashMap<String, Class<?>>();
for (AttributeDescriptor att : source.getFeatureType().getAttributeDescriptors()) {
sourceTypes.put(att.getLocalName(), att.getType().getBinding());
}
for (AttributeDescriptor att : target.getFeatureType().getAttributeDescriptors()) {
String name = att.getLocalName();
Class<?> sourceBinding = sourceTypes.get(name);
if (sourceBinding != null && sourceBinding == att.getType().getBinding()) {
Object attribute = source.getAttribute(name);
target.setAttribute(name, attribute);
}
}
}
/**
* Returns a representative (first axis with a length) unit of the given
* crs.
*
* @param crs
* @return a representative unit of the given crs.
*/
public static Unit<?> getDefaultCRSUnit(CoordinateReferenceSystem crs) {
assert crs != null;
CoordinateSystem coordinateSystem = crs.getCoordinateSystem();
Unit<?> unit = CRSUtilities.getUnit(coordinateSystem);
if (unit == null) {
CoordinateSystemAxis axis = coordinateSystem.getAxis(0);
unit = axis.getUnit();
}
return unit;
}
/**
* @param gFactory
* @param geomCrs
* @param reprojectCrs
* @return a GeometryCoordinateSequenceTransformer configured to transform
* geometries from <code>geomCrs</code> to <code>reprojectCrs</code>
* @throws OperationNotFoundException
*/
public static GeometryCoordinateSequenceTransformer getTransformer( final GeometryFactory gFactory,
final CoordinateReferenceSystem geomCrs,
final CoordinateReferenceSystem reprojectCrs)
throws OperationNotFoundException {
assert geomCrs != null;
assert reprojectCrs != null;
CoordinateSequenceFactory csFactory;
CoordinateSequenceTransformer csTransformer;
GeometryCoordinateSequenceTransformer transformer;
csFactory = gFactory.getCoordinateSequenceFactory();
csTransformer = new CoordSeqFactoryPreservingCoordinateSequenceTransformer(csFactory);
transformer = new GeometryCoordinateSequenceTransformer(csTransformer);
MathTransform mathTransform;
try {
mathTransform = CRS.findMathTransform(geomCrs, reprojectCrs, true);
} catch (FactoryException e) {
throw new OperationNotFoundException(e.getMessage());
}
transformer.setMathTransform(mathTransform);
return transformer;
}
/**
* Reprojects <code>geom</code> from <code>geomCrs</code> to
* <code>reprojectCrs</code>
*
* @param geom
* the geometry to reproject
* @param geomCrs
* the original CRS of <code>geom</code>
* @param reprojectCrs
* the CRS to reproject <code>geom</code> onto
* @return <code>geom</code> reprojected from <code>geomCrs</code> to
* <code>reprojectCrs</code>
* @throws OperationNotFoundException
* if there isn't a transformation in GeoTools to convert from
* <code>geomCrs</code> to <code>reprojectCrs</code>
* @throws TransformException
*/
public static Geometry reproject( final Geometry geom,
final CoordinateReferenceSystem geomCrs,
final CoordinateReferenceSystem reprojectCrs)
throws OperationNotFoundException, TransformException {
assert geom != null : "geom cannot be null";
assert geomCrs != null : "geomCrs cannot be null";
assert reprojectCrs != null : "reprojectCrs cannot be null";
if (geomCrs.equals(reprojectCrs)) {
return geom;
}
if (CRS.equalsIgnoreMetadata(geomCrs, reprojectCrs)) {
return geom;
}
GeometryFactory gFactory = geom.getFactory();
GeometryCoordinateSequenceTransformer transformer;
transformer = getTransformer(gFactory, geomCrs, reprojectCrs);
Geometry geometry;
try {
geometry = transformer.transform(geom);
} catch (TransformException e) {
throw e;
}
return geometry;
}
/**
* @param featureType
* @return the dimension of default geometry
*/
public static int getDimensionOf(final SimpleFeatureType featureType) {
GeometryDescriptor geomAttr = featureType.getGeometryDescriptor();
Class<?> geomClass = geomAttr.getType().getBinding();
int dim = GeometryUtil.getDimension(geomClass);
return dim;
}
/**
* Utility method to easily reproject a line segment as LineSegment is not a
* Geometry
*
* @param segment
* the segment to reproject
* @param segmentCrs
* the CRS the segment coordinates are in
* @param reprojectCrs
* the CRS to reproject the segment to
* @return a new line segment built from the reprojected coordinates of
* <code>segment</code> from <code>segmentCrs</code> to
* <code>reprojectCrs</code>
*/
public static LineSegment reproject(LineSegment segment,
CoordinateReferenceSystem segmentCrs,
CoordinateReferenceSystem reprojectCrs) {
assert segment != null;
assert segmentCrs != null;
assert reprojectCrs != null;
if (segmentCrs.equals(reprojectCrs)) {
return segment;
}
if (CRS.equalsIgnoreMetadata(segmentCrs, reprojectCrs)) {
return segment;
}
MathTransform mathTransform;
try {
mathTransform = CRS.findMathTransform(segmentCrs, reprojectCrs, true);
} catch (FactoryException e) {
throw new RuntimeException(e.getMessage());
}
double[] src = { segment.p0.x, segment.p0.y, segment.p1.x, segment.p1.y };
double[] dst = new double[4];
try {
mathTransform.transform(src, 0, dst, 0, 2);
} catch (TransformException e) {
throw new RuntimeException(e.getMessage());
}
Coordinate p0 = new Coordinate(dst[0], dst[1]);
Coordinate p1 = new Coordinate(dst[2], dst[3]);
LineSegment lineSegment = new LineSegment(p0, p1);
return lineSegment;
}
}