org.jgrasstools.lesto.modules.utilities.LasSlicer Maven / Gradle / Ivy
/*
* 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 > recordsMap = new LinkedHashMap<>();
pm.beginTask("Producing slices...", (xRange.length - 1));
for( int x = 0; x < xRange.length - 1; x++ ) {
for( int y = 0; y < yRange.length - 1; y++ ) {
tilesCount++;
Envelope env = new Envelope(xRange[x], xRange[x + 1], yRange[y], yRange[y + 1]);
Polygon polygon = GeometryUtilities.createPolygonFromEnvelope(env);
List pointsInGeometry = dataManager.getPointsInGeometry(polygon, true);
pm.message("Points in tile " + x + "/" + y + "(" + tilesCount + " of " + tilesNum + "): "
+ pointsInGeometry.size());
if (pointsInGeometry.size() == 0) {
continue;
}
for( double z = minZ + pInterval; z < maxZ; z = z + pInterval ) {
String key = String.valueOf(z);
List pointsInSlice = recordsMap.get(key);
if (pointsInSlice == null) {
pointsInSlice = new ArrayList();
recordsMap.put(key, pointsInSlice);
}
double height = z - minZ;
double low = height - pThickness / 2.0;
double high = height + pThickness / 2.0;
List pointsInHeightRange = ALasDataManager.getPointsInHeightRange(pointsInGeometry, low, high);
if (pointsInHeightRange.size() > 0) {
pointsInSlice.addAll(pointsInHeightRange);
// pm.message("Added points: " + pointsInHeightRange.size());
if (!doRaster)
chartPoints(outputFolder, height, pointsInSlice, chartWidth, chartHeigth, minX, maxX, minY, maxY);
}
}
}
pm.worked(1);
}
pm.done();
if (doRaster) {
Set>> entrySet = recordsMap.entrySet();
int size = entrySet.size();
pm.beginTask("Generating rasters...", size);
for( Entry> entry : entrySet ) {
double z = Double.parseDouble(entry.getKey());
double height = z - minZ;
List points = entry.getValue();
dumpRaster(outputFolder, height, points, pResolution, minX, maxX, minY, maxY, crs);
pm.worked(1);
}
pm.done();
}
}
}
private void dumpRaster( File outputFolder, double height, List points, double resolution, double minX,
double maxX, double minY, double maxY, CoordinateReferenceSystem crs ) throws Exception {
File rasterFile = new File(outputFolder, "slice_" + height + ".asc");
int rows = (int) round((maxY - minY) / resolution);
int cols = (int) round((maxX - minX) / resolution);
GridGeometry2D gridGeometry = CoverageUtilities.gridGeometryFromRegionValues(maxY, minY, maxX, minX, cols, rows, crs);
RegionMap regionMap = CoverageUtilities.gridGeometry2RegionParamsMap(gridGeometry);
WritableRaster outWR = CoverageUtilities.createDoubleWritableRaster(cols, rows, null, null, JGTConstants.doubleNovalue);
WritableRandomIter outIter = CoverageUtilities.getWritableRandomIterator(outWR);
Point point = new Point();
for( LasRecord dot : points ) {
CoverageUtilities.colRowFromCoordinate(new Coordinate(dot.x, dot.y), gridGeometry, point);
outIter.setSample(point.x, point.y, 0, height);
}
outIter.done();
GridCoverage2D outRaster = CoverageUtilities.buildCoverage(rasterFile.getName(), outWR, regionMap, crs);
dumpRaster(outRaster, rasterFile.getAbsolutePath());
}
private void chartPoints( File chartFolder, double z, List pointsInSlice, int width, int height, double minX,
double maxX, double minY, double maxY ) throws IOException {
File chartFile = new File(chartFolder, "slice_" + z + ".png");
int size = pointsInSlice.size();
double[] xPlanim = new double[size];
double[] yPlanim = new double[size];
for( int i = 0; i < size; i++ ) {
LasRecord dot = pointsInSlice.get(i);
xPlanim[i] = dot.x;
yPlanim[i] = dot.y;
}
Scatter scatterPlanim = new Scatter("Slice " + z);
scatterPlanim.addSeries("planimetry", xPlanim, yPlanim);
scatterPlanim.setShowLines(false);
scatterPlanim.setXLabel("longitude");
scatterPlanim.setYLabel("latitude");
scatterPlanim.setXRange(minX, maxX);
scatterPlanim.setYRange(minY, maxY);
BufferedImage imagePlanim = scatterPlanim.getImage(width, height);
ImageIO.write(imagePlanim, "png", chartFile);
}
}
© 2015 - 2025 Weber Informatics LLC | Privacy Policy