org.jgrasstools.lesto.modules.buildings.LasOnDtmBuildingsExtractor 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 pointsInGeometry = lasHandler.getPointsInGeometry(regionPolygon, true);
pm.done();
pm.beginTask("Buildings filtering...", (int) pointsInGeometry.size());
for( LasRecord dot : pointsInGeometry ) {
Coordinate c = new Coordinate(dot.x, dot.y);
if (!envelope2d.contains(dot.x, dot.y)) {
continue;
}
double dtmValue = CoverageUtilities.getValue(inDtmGC, c);
double height = abs(dot.z - dtmValue);
CoverageUtilities.colRowFromCoordinate(c, gridGeometry, p);
if (height < pDtmThres) {
buildingsHolder[0].setSample(p.x, p.y, 0, doubleNovalue);
}
pm.worked(1);
}
pm.done();
OmsMorpher morpher = new OmsMorpher();
morpher.pm = pm;
morpher.inMap = newBuildingsRaster;
morpher.pMode = Variables.OPEN;
morpher.process();
newBuildingsRaster = morpher.outMap;
OmsVectorizer buildingsVectorizer = new OmsVectorizer();
buildingsVectorizer.pm = pm;
buildingsVectorizer.inRaster = newBuildingsRaster;
buildingsVectorizer.doRemoveHoles = true;
double minAreaInCells = pMinArea / pRasterResolution / pRasterResolution;
buildingsVectorizer.pThres = minAreaInCells;
buildingsVectorizer.fDefault = "rast";
buildingsVectorizer.process();
buildingsFC = buildingsVectorizer.outVector;
dumpVector(buildingsFC, outBuildings);
if (doSmoothing && outCleanBuildings != null) {
// smooth
SimpleFeatureCollection smoothedFC = smoothBuildings(pDensifyResolution, buildingsFC);
// remove trees
DefaultFeatureCollection removedAndSmoothedFC = removeNonBuildings(lasHandler, smoothedFC, inDtmGC,
pBuildingsBuffer);
dumpVector(removedAndSmoothedFC, outCleanBuildings);
}
}
}
private DefaultFeatureCollection removeNonBuildings( ALasDataManager lasHandler, SimpleFeatureCollection buildingsFC,
GridCoverage2D dem, double buildingsBuffer ) throws Exception {
final List buildingsList = FeatureUtilities.featureCollectionToList(buildingsFC);
final List checkedBuildings = new ArrayList();
pm.beginTask("Removing buildings...", buildingsList.size());
for( int i = 0; i < buildingsList.size(); i++ ) {
SimpleFeature building = buildingsList.get(i);
Geometry buildingGeom = (Geometry) building.getDefaultGeometry();
Geometry bufferedGeom = buildingGeom.buffer(buildingsBuffer);
List points = lasHandler.getPointsInGeometry(bufferedGeom, false);
int percOfOnes = checkReturnNum(points, bufferedGeom);
if (percOfOnes >= 96) {
checkedBuildings.add(building);
}
pm.worked(1);
}
pm.done();
DefaultFeatureCollection fc = new DefaultFeatureCollection();
fc.addAll(checkedBuildings);
return fc;
}
private int checkReturnNum( List points, Geometry bufferedGeom ) {
int count = 0;
for( LasRecord l : points ) {
// System.out.print(l.returnNumber);
if (l.returnNumber == 1) {
count++;
}
}
double perc = count * 100.0 / points.size();
return (int) round(perc);
}
private SimpleFeatureCollection smoothBuildings( double densifyResolution, SimpleFeatureCollection buildingsFC )
throws Exception {
final List buildingsList = FeatureUtilities.featureCollectionToList(buildingsFC);
List newPolygons = new ArrayList();
pm.beginTask("Smoothing buildings...", buildingsList.size());
for( final SimpleFeature building : buildingsList ) {
Geometry buildingGeom = (Geometry) building.getDefaultGeometry();
buildingGeom = Densifier.densify(buildingGeom, densifyResolution);
Coordinate[] coordinates = buildingGeom.getCoordinates();
List newCoords = new ArrayList();
for( int i = 0; i < coordinates.length - 1; i++ ) {
Coordinate from = coordinates[i];
Coordinate to = coordinates[i + 1];
LineSegment l = new LineSegment(from, to);
Coordinate newCoord = l.pointAlong(0.5);
newCoords.add(newCoord);
}
newCoords.add(newCoords.get(0));
Polygon newPolygon = gf.createPolygon(newCoords.toArray(new Coordinate[0]));
// newPolygon = (Polygon) DouglasPeuckerSimplifier.simplify(newPolygon, 0.1);
newPolygons.add(newPolygon);
pm.worked(1);
}
pm.done();
SimpleFeatureCollection fc = FeatureUtilities.featureCollectionFromGeometry(buildingsFC.getBounds()
.getCoordinateReferenceSystem(), newPolygons.toArray(GeometryUtilities.TYPE_POLYGON));
OmsLineSmootherMcMaster smoother = new OmsLineSmootherMcMaster();
smoother.pm = pm;
smoother.pLimit = 10;
smoother.inVector = fc;
smoother.pLookahead = 5;
// smoother.pSlide = 1;
// smoother.pDensify = 0.2;
// smoother.pSimplify = 0.1;
smoother.process();
fc = smoother.outVector;
List geomsList = FeatureUtilities.featureCollectionToGeometriesList(fc, true, null);
List newGeomsList = new ArrayList();
for( Geometry geometry : geomsList ) {
Polygon polygon = gf.createPolygon(geometry.getCoordinates());
polygon = (Polygon) TopologyPreservingSimplifier.simplify(polygon, 0.4);
newGeomsList.add(polygon);
}
fc = FeatureUtilities.featureCollectionFromGeometry(buildingsFC.getBounds().getCoordinateReferenceSystem(),
newGeomsList.toArray(GeometryUtilities.TYPE_POLYGON));
return fc;
}
}
© 2015 - 2025 Weber Informatics LLC | Privacy Policy