org.jgrasstools.gears.modules.r.carver.OmsCarver 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 riverGeoms = FeatureUtilities.featureCollectionToGeometriesList(inCarverLines,
// true, null);
// List bufferedRiverGeoms = new ArrayList();
// for( Geometry riverGeom : riverGeoms ) {
// Geometry bufferedGeom = riverGeom.buffer(pBuffer);
// bufferedRiverGeoms.add(bufferedGeom);
// }
//
// pm.beginTask("Union...", -1);
// GeometryCollection gc = new
// GeometryCollection(bufferedRiverGeoms.toArray(GeometryUtilities.TYPE_GEOMETRY),
// GeometryUtilities.gf());
// Geometry bufferUnion = gc.union();
// pm.done();
//
// SimpleFeatureCollection outPoints = new DefaultFeatureCollection();
// SimpleFeatureCollection outPolygons = new DefaultFeatureCollection();
// SimpleFeatureTypeBuilder polygonTypeBuilder = new SimpleFeatureTypeBuilder();
// polygonTypeBuilder.setName("cpoints");
// polygonTypeBuilder.setCRS(crs);
// polygonTypeBuilder.add("the_geom", Polygon.class);
// SimpleFeatureType polygonType = polygonTypeBuilder.buildFeatureType();
// SimpleFeatureBuilder polygonBuilder = new SimpleFeatureBuilder(polygonType);
//
// pm.beginTask("Bob it...", bufferUnion.getNumGeometries());
// for( int i = 0; i < bufferUnion.getNumGeometries(); i++ ) {
// Geometry bufferArea = bufferUnion.getGeometryN(i);
// PreparedGeometry preparedArea = PreparedGeometryFactory.prepare(bufferArea);
// List involvedRiversCoordinate = new ArrayList();
// // pm.beginTask("Find involved rivers...", riverGeoms.size());
// for( Geometry riverGeom : riverGeoms ) {
// if (preparedArea.intersects(riverGeom)) {
// Coordinate[] coordinates = riverGeom.getCoordinates();
// for( Coordinate coordinate : coordinates ) {
// involvedRiversCoordinate.add(coordinate);
// }
// }
// // pm.worked(1);
// }
// // pm.done();
//
// if (involvedRiversCoordinate.size() < 4) {
// continue;
// }
//
// // use buffered area and involved rivers to carve - one area at a time
//
// SimpleFeatureCollection controlPoints = createFC4ControlPoints(crs, involvedRiversCoordinate,
// inRaster);
// SimpleFeatureCollection areaFC = FeatureUtilities.featureCollectionFromGeometry(crs,
// bufferArea);
// outPoints.addAll(controlPoints);
//
// OmsBobTheBuilder bob = new OmsBobTheBuilder();
// bob.pm = new DummyProgressMonitor();// pm;
// bob.inRaster = inRaster;
// bob.inArea = areaFC;
// bob.inElevations = controlPoints;
// bob.pMaxbuffer = pBuffer + 10;
// bob.doErode = true;
// bob.doPolygonborder = true;
// bob.doUseOnlyInternal = true;
// bob.fElevation = "elev";
// bob.process();
// inRaster = bob.outRaster;
//
// polygonBuilder.addAll(new Object[]{bufferArea});
// SimpleFeature polygonFeature = polygonBuilder.buildFeature(null);
// outPolygons.add(polygonFeature);
//
// pm.worked(1);
// }
// pm.done();
//
// outRaster = inRaster;
// }
// private SimpleFeatureCollection createFC4ControlPoints( CoordinateReferenceSystem crs,
// List riverCoords,
// GridCoverage2D dtm ) throws Exception {
// SimpleFeatureTypeBuilder b = new SimpleFeatureTypeBuilder();
// b.setName("cpoints");
// b.setCRS(crs);
// b.add("the_geom", Point.class);
// b.add("elev", Double.class);
// SimpleFeatureType type = b.buildFeatureType();
// SimpleFeatureBuilder builder = new SimpleFeatureBuilder(type);
//
// SimpleFeatureCollection fc = new DefaultFeatureCollection();
//
// for( Coordinate coordinate : riverCoords ) {
// double value = CoverageUtilities.getValue(dtm, coordinate);
// Point point = GeometryUtilities.gf().createPoint(coordinate);
// value = value - pDepth;
// Object[] values = new Object[]{point, value};
// builder.addAll(values);
// SimpleFeature feature = builder.buildFeature(null);
// fc.add(feature);
// }
//
// return fc;
// }
}