org.jgrasstools.lesto.modules.vegetation.rastermaxima.OmsRasterMaximaFinder 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 surroundingNodes = node.getSurroundingNodes();
int topBufferThresViolationCount = 0;
for( GridNode gridNode : surroundingNodes ) {
if (gridNode != null) {
double elev = gridNode.elevation;
double elevDiff = elevation - elev;
if (elevDiff > pTopBufferThres) {
topBufferThresViolationCount++;
}
} else {
nonValidCells++;
}
if (nonValidCells > 1) {
note = "exclude: found invalid neighbor cells = " + nonValidCells;
isMax = false;
break;
}
}
if (isMax && pTopBufferThresCellCount > 0 && topBufferThresViolationCount >= pTopBufferThresCellCount) {
isMax = false;
note = "exclude: elevation diff of neighbors from top violates thres (" + pTopBufferThres + "/"
+ topBufferThresViolationCount + ")";
}
// create circle
if (isMax) {
if (elevation >= pThreshold) {
Coordinate coordinate = CoverageUtilities.coordinateFromColRow(c, r, gridGeometry);
Point point = gf.createPoint(coordinate);
String elevStr = formatter.format(elevation);
elevStr = elevStr.replace(',', '.');
elevation = Double.parseDouble(elevStr);
Object[] values = new Object[]{point, id, elevation, note};
maximaBuilder.addAll(values);
SimpleFeature feature = maximaBuilder.buildFeature(null);
((DefaultFeatureCollection) outMaxima).add(feature);
double radius = (size * xRes) / 2.0;
Geometry buffer = point.buffer(radius);
values = new Object[]{buffer, id, buffer.getArea()};
circleBuilder.addAll(values);
feature = circleBuilder.buildFeature(null);
((DefaultFeatureCollection) outCircles).add(feature);
id++;
}
}
}
}
pm.worked(1);
}
pm.done();
elevIter.done();
if (pBorderDistanceThres > 0) {
OmsVectorizer vectorizer = new OmsVectorizer();
vectorizer.pm = pm;
vectorizer.inRaster = inDsmDtmDiff;
vectorizer.pValue = null;
vectorizer.pThres = 0;
vectorizer.doMask = true;
vectorizer.pMaskThreshold = pThreshold;
vectorizer.fDefault = "rast";
vectorizer.process();
SimpleFeatureCollection diffPolygons = vectorizer.outVector;
List diffGeoms = FeatureUtilities.featureCollectionToGeometriesList(diffPolygons, true, null);
List bordersGeoms = new ArrayList();
for( Geometry geometry : diffGeoms ) {
if (geometry instanceof Polygon) {
Polygon polygon = (Polygon) geometry;
LineString exteriorRing = polygon.getExteriorRing();
bordersGeoms.add(exteriorRing);
int numInteriorRing = polygon.getNumInteriorRing();
for( int i = 0; i < numInteriorRing; i++ ) {
LineString interiorRingN = polygon.getInteriorRingN(i);
bordersGeoms.add(interiorRingN);
}
}
}
MultiLineString allBorders = gf.createMultiLineString(bordersGeoms.toArray(GeometryUtilities.TYPE_LINESTRING));
outBorders = FeatureUtilities.featureCollectionFromGeometry(inDsmDtmDiff.getCoordinateReferenceSystem(), allBorders);
SimpleFeatureIterator maximaIter = outMaxima.features();
while( maximaIter.hasNext() ) {
SimpleFeature maxima = maximaIter.next();
Geometry maximaGeometry = (Geometry) maxima.getDefaultGeometry();
double distance = DistanceOp.distance(allBorders, maximaGeometry);
if (distance < pBorderDistanceThres) {
maxima.setAttribute(NOTE, "exclude: near border: " + distance);
}
}
maximaIter.close();
}
}
}
© 2015 - 2025 Weber Informatics LLC | Privacy Policy