org.jgrasstools.gears.modules.v.smoothing.OmsLineSmootherMcMaster Maven / Gradle / Ivy
The newest version!
/*
* This file is part of JGrasstools (http://www.jgrasstools.org)
* (C) HydroloGIS - www.hydrologis.com
*
* JGrasstools is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* This program 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 General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program. If not, see linesList;
@Execute
public void process() throws Exception {
if (!concatOr(outVector == null, doReset)) {
return;
}
if (pDensify != null) {
densify = pDensify;
}
if (pSimplify != null) {
simplify = pSimplify;
}
outVector = new DefaultFeatureCollection();
pm.message("Collecting geometries...");
linesList = FeatureUtilities.featureCollectionToList(inVector);
int size = inVector.size();
FeatureGeometrySubstitutor fGS = new FeatureGeometrySubstitutor(inVector.getSchema());
pm.beginTask("Smoothing features...", size);
for( SimpleFeature line : linesList ) {
Geometry geometry = (Geometry) line.getDefaultGeometry();
List lsList = smoothGeometries(geometry);
if (lsList.size() != 0) {
LineString[] lsArray = (LineString[]) lsList.toArray(new LineString[lsList.size()]);
MultiLineString multiLineString = gF.createMultiLineString(lsArray);
SimpleFeature newFeature = fGS.substituteGeometry(line, multiLineString);
((DefaultFeatureCollection) outVector).add(newFeature);
}
pm.worked(1);
}
pm.done();
}
private List smoothGeometries( Geometry geometry ) {
List lsList = new ArrayList();
int numGeometries = geometry.getNumGeometries();
for( int i = 0; i < numGeometries; i++ ) {
Geometry geometryN = geometry.getGeometryN(i);
double length = geometryN.getLength();
Coordinate[] smoothedArray = geometryN.getCoordinates();
Coordinate first = smoothedArray[0];
Coordinate last = smoothedArray[smoothedArray.length - 1];
if (length <= pLimit) {
// if it is circle remove it, else just do not smooth it
if (first.distance(last) < SAMEPOINTTHRESHOLD) {
continue;
}
// check if the line is an error lying around somewhere
if (isAlone(geometryN)) {
continue;
}
} else {
if (densify != -1) {
geometryN = Densifier.densify(geometryN, pDensify);
}
List smoothedCoords = Collections.emptyList();;
FeatureSlidingAverage fSA = new FeatureSlidingAverage(geometryN);
smoothedCoords = fSA.smooth(pLookahead, false, pSlide);
if (smoothedCoords != null) {
smoothedArray = (Coordinate[]) smoothedCoords.toArray(new Coordinate[smoothedCoords.size()]);
} else {
smoothedArray = geometryN.getCoordinates();
}
}
LineString lineString = gF.createLineString(smoothedArray);
if (simplify != -1) {
TopologyPreservingSimplifier tpSimplifier = new TopologyPreservingSimplifier(lineString);
tpSimplifier.setDistanceTolerance(pSimplify);
lineString = (LineString) tpSimplifier.getResultGeometry();
}
lsList.add(lineString);
}
return lsList;
}
/**
* Checks if the given geometry is connected to any other line.
*
* @param geometryN the geometry to test.
* @return true if the geometry is alone in the space, i.e. not connected at
* one of the ends to any other geometry.
*/
private boolean isAlone( Geometry geometryN ) {
Coordinate[] coordinates = geometryN.getCoordinates();
if (coordinates.length > 1) {
Coordinate first = coordinates[0];
Coordinate last = coordinates[coordinates.length - 1];
for( SimpleFeature line : linesList ) {
Geometry lineGeom = (Geometry) line.getDefaultGeometry();
int numGeometries = lineGeom.getNumGeometries();
for( int i = 0; i < numGeometries; i++ ) {
Geometry subGeom = lineGeom.getGeometryN(i);
Coordinate[] lineCoordinates = subGeom.getCoordinates();
if (lineCoordinates.length < 2) {
continue;
} else {
Coordinate tmpFirst = lineCoordinates[0];
Coordinate tmpLast = lineCoordinates[lineCoordinates.length - 1];
if (tmpFirst.distance(first) < SAMEPOINTTHRESHOLD || tmpFirst.distance(last) < SAMEPOINTTHRESHOLD
|| tmpLast.distance(first) < SAMEPOINTTHRESHOLD || tmpLast.distance(last) < SAMEPOINTTHRESHOLD) {
return false;
}
}
}
}
}
// 1 point line or no connection, mark it as alone for removal
return true;
}
/**
* An utility method to use the module with default values and shapefiles.
*
*
* This will use the windowed average and a density of 0.2, simplification threshold of 0.1
* and a lookahead of 13, as well as a length filter of 10.
*
*
* @param shapePath the input file.
* @param outPath the output smoothed path.
* @throws Exception
*/
public static void defaultSmoothShapefile( String shapePath, String outPath ) throws Exception {
PrintStreamProgressMonitor pm = new PrintStreamProgressMonitor(System.out, System.err);
SimpleFeatureCollection initialFC = OmsShapefileFeatureReader.readShapefile(shapePath);
OmsLineSmootherMcMaster smoother = new OmsLineSmootherMcMaster();
smoother.pm = pm;
smoother.pLimit = 10;
smoother.inVector = initialFC;
smoother.pLookahead = 13;
// smoother.pSlide = 1;
smoother.pDensify = 0.2;
smoother.pSimplify = 0.01;
smoother.process();
SimpleFeatureCollection smoothedFeatures = smoother.outVector;
OmsShapefileFeatureWriter.writeShapefile(outPath, smoothedFeatures, pm);
}
}