org.jgrasstools.gears.modules.r.scanline.OmsScanLineRasterizer 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 featureIterator = inVector.features();
ExecutorService fixedThreadPool = Executors.newFixedThreadPool(getDefaultThreadsNum());
while( featureIterator.hasNext() ) {
final SimpleFeature feature = featureIterator.next();
// extract the value to put into the raster.
double tmpValue = -1.0;
if (pValue == null) {
tmpValue = ((Number) feature.getAttribute(fCat)).doubleValue();
} else {
tmpValue = pValue;
}
final double value = tmpValue;
final double delta = xRes / 4.0;
Runnable runner = new Runnable(){
public void run() {
try {
Geometry geometry = (Geometry) feature.getDefaultGeometry();
int numGeometries = geometry.getNumGeometries();
for( int i = 0; i < numGeometries; i++ ) {
final Geometry geometryN = geometry.getGeometryN(i);
// PreparedGeometry preparedGeometryN =
// PreparedGeometryFactory.prepare(geometryN);
for( int r = 0; r < height; r++ ) {
// do scan line to fill the polygon
double[] westPos = gridGeometry.gridToWorld(new GridCoordinates2D(0, r)).getCoordinate();
double[] eastPos = gridGeometry.gridToWorld(new GridCoordinates2D(width - 1, r)).getCoordinate();
Coordinate west = new Coordinate(westPos[0], westPos[1]);
Coordinate east = new Coordinate(eastPos[0], eastPos[1]);
LineString line = gf.createLineString(new Coordinate[]{west, east});
if (geometryN.intersects(line)) {
Geometry internalLines = geometryN.intersection(line);
int lineNums = internalLines.getNumGeometries();
for( int l = 0; l < lineNums; l++ ) {
Coordinate[] coords = internalLines.getGeometryN(l).getCoordinates();
if (coords.length == 2) {
for( int j = 0; j < coords.length; j = j + 2 ) {
Coordinate startC = new Coordinate(coords[j].x + delta, coords[j].y);
Coordinate endC = new Coordinate(coords[j + 1].x - delta, coords[j + 1].y);
DirectPosition2D startDP;
DirectPosition2D endDP;
if (startC.x < endC.x) {
startDP = new DirectPosition2D(startC.x, startC.x);
endDP = new DirectPosition2D(endC.x, endC.x);
} else {
startDP = new DirectPosition2D(endC.x, endC.x);
endDP = new DirectPosition2D(startC.x, startC.x);
}
GridCoordinates2D startGridCoord = gridGeometry.worldToGrid(startDP);
GridCoordinates2D endGridCoord = gridGeometry.worldToGrid(endDP);
/*
* the part in between has to be filled
*/
for( int k = startGridCoord.x; k <= endGridCoord.x; k++ ) {
if (inIter != null && fCat == null) {
double v = inIter.getSampleDouble(k, r, 0);
outWR.setSample(k, r, 0, v);
} else {
outWR.setSample(k, r, 0, value);
}
}
}
} else {
if (coords.length == 1) {
pm.errorMessage(MessageFormat.format("Found a cusp in: {0}/{1}", coords[0].x,
coords[0].y));
} else {
throw new ModelsIOException(MessageFormat.format(
"Found intersection with more than 2 points in: {0}/{1}", coords[0].x,
coords[0].y), this);
}
}
}
}
}
}
pm.worked(1);
} catch (Exception e) {
pm.errorMessage(e.getLocalizedMessage());
e.printStackTrace();
}
}
};
fixedThreadPool.execute(runner);
}
try {
fixedThreadPool.shutdown();
fixedThreadPool.awaitTermination(30, TimeUnit.DAYS);
fixedThreadPool.shutdownNow();
} catch (InterruptedException e) {
e.printStackTrace();
}
pm.done();
featureIterator.close();
}
}