All Downloads are FREE. Search and download functionalities are using the official Maven repository.
Many resources are needed to download a project. Please understand that we have to compensate our server costs. Thank you in advance.
Project price only 1 $
You can buy this project and download/modify it how often you want.
org.hortonmachine.gears.modules.v.marchingsquares.OmsMarchingSquaresVectorializer Maven / Gradle / Ivy
/*
* This file is part of HortonMachine (http://www.hortonmachine.org)
* (C) HydroloGIS - www.hydrologis.com
*
* The HortonMachine 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 awtGeometriesList;
private RandomIter iter = null;
private double xRes;
private double yRes;
private GridGeometry2D gridGeometry;
private int height;
private int width;
private BitSet bitSet = null;
private CoordinateReferenceSystem crs;
@Execute
public void process() throws Exception {
if (!concatOr(outGeodata == null, doReset)) {
return;
}
if (iter == null) {
RenderedImage inputRI = inGeodata.getRenderedImage();
iter = RandomIterFactory.create(inputRI, null);
HashMap regionMap = getRegionParamsFromGridCoverage(inGeodata);
height = regionMap.get(ROWS).intValue();
width = regionMap.get(COLS).intValue();
xRes = regionMap.get(XRES);
yRes = regionMap.get(YRES);
crs = inGeodata.getCoordinateReferenceSystem();
bitSet = new BitSet(width * height);
gridGeometry = inGeodata.getGridGeometry();
}
List geometriesList = new ArrayList();
awtGeometriesList = new ArrayList();
pm.beginTask("Extracting vectors...", height);
/*
* a List where to store the different value of the raster if pValue is
* null, otherwise only pValue is putted into.
*/
ArrayList valueRaster = new ArrayList();
/*
* if pValue is a number then extract the polygon from the raster (if
* the pixel value is equals to pValue).
*/
if (pValue != null) {
valueRaster.add(0, pValue);
for( int row = 0; row < height; row++ ) {
for( int col = 0; col < width; col++ ) {
double value = iter.getSampleDouble(col, row, 0);
if (!isNovalue(value) && !bitSet.get(row * width + col)) {
if (abs(value - pValue) < .0000001) {
java.awt.Polygon awtPolygon = new java.awt.Polygon();
Polygon polygon = identifyPerimeter(col, row, awtPolygon);
if (polygon != null) {
geometriesList.add(polygon);
awtGeometriesList.add(awtPolygon);
}
}
}
}
pm.worked(1);
}
} else {
/*
*
*/
pValue = doubleNovalue;
for( int row = 0; row < height; row++ ) {
for( int col = 0; col < width; col++ ) {
double value = iter.getSampleDouble(col, row, 0);
if (value != pValue) {
pValue = value;
if (!isNovalue(value) && !bitSet.get(row * width + col)) {
if (value == pValue) {
java.awt.Polygon awtPolygon = new java.awt.Polygon();
Polygon polygon = identifyPerimeter(col, row, awtPolygon);
if (polygon != null) {
valueRaster.add(pValue);
geometriesList.add(polygon);
awtGeometriesList.add(awtPolygon);
}
}
}
}
}
pm.worked(1);
}
}
pm.done();
SimpleFeatureTypeBuilder b = new SimpleFeatureTypeBuilder();
b.setName("raster2vector");
b.setCRS(crs);
b.add("the_geom", Polygon.class);
b.add("cat", Integer.class);
b.add(defaultFeatureField, Double.class);
SimpleFeatureType type = b.buildFeatureType();
outGeodata = new DefaultFeatureCollection();
int index = 0;
for( Polygon polygon : geometriesList ) {
Double tmpValue = valueRaster.get(index);
SimpleFeatureBuilder builder = new SimpleFeatureBuilder(type);
Object[] values = new Object[]{polygon, index, tmpValue};
builder.addAll(values);
SimpleFeature feature = builder.buildFeature(type.getTypeName() + "." + index);
index++;
((DefaultFeatureCollection) outGeodata).add(feature);
}
}
private Polygon identifyPerimeter( int initialX, int initialY, java.awt.Polygon awtPolygon ) throws TransformException {
if (initialX < 0 || initialX > width - 1 || initialY < 0 || initialY > height - 1)
throw new IllegalArgumentException("Coordinate outside the bounds.");
int initialValue = value(initialX, initialY);
if (initialValue == 0) {
throw new IllegalArgumentException(String.format("Supplied initial coordinates (%d, %d) do not lie on a perimeter.",
initialX, initialY));
}
final Point2D worldPosition = new Point2D.Double(initialX, initialY);
gridGeometry.getGridToCRS2D().transform(worldPosition, worldPosition);
if (initialValue == 15) {
// not a border pixel
return null;
}
Coordinate startCoordinate = new Coordinate(worldPosition.getX() + xRes / 2.0, worldPosition.getY() - yRes / 2.0);
List coordinateList = new ArrayList(200);
coordinateList.add(startCoordinate);
double currentX = startCoordinate.x;
double currentY = startCoordinate.y;
int x = initialX;
int y = initialY;
awtPolygon.addPoint(x, y);
boolean previousWentNorth = false;
boolean previousWentEast = false;
do {
Coordinate direction = null;
int dx = 0;
int dy = 0;
int v = value(x, y);
switch( v ) {
case 1:
dy = -1; // N
currentY = currentY + yRes;
previousWentNorth = true;
break;
case 2:
dx = +1; // E
currentX = currentX + xRes;
previousWentEast = true;
break;
case 3:
dx = +1; // E
currentX = currentX + xRes;
previousWentEast = true;
break;
case 4:
dx = -1; // W
currentX = currentX - xRes;
previousWentEast = false;
break;
case 5:
dy = -1; // N
currentY = currentY + yRes;
previousWentNorth = true;
break;
case 6:
if (!previousWentNorth) {// W
dx = -1;
currentX = currentX - xRes;
previousWentEast = false;
} else {
dx = +1;// E
currentX = currentX + xRes;
previousWentEast = true;
}
break;
case 7:
dx = +1; // E
currentX = currentX + xRes;
previousWentEast = true;
break;
case 8:
dy = +1; // S
currentY = currentY - yRes;
previousWentNorth = false;
break;
case 9:
if (previousWentEast) {
dy = -1; // N
currentY = currentY + yRes;
previousWentNorth = true;
} else {
dy = +1; // S
currentY = currentY - yRes;
previousWentNorth = false;
}
break;
case 10:
dy = +1; // S
currentY = currentY - yRes;
previousWentNorth = false;
break;
case 11:
dy = +1; // S
currentY = currentY - yRes;
previousWentNorth = false;
break;
case 12:
dx = -1; // W
currentX = currentX - xRes;
previousWentEast = false;
break;
case 13:
dy = -1; // N
currentY = currentY + yRes;
previousWentNorth = true;
break;
case 14:
dx = -1; // W
currentX = currentX - xRes;
previousWentEast = false;
break;
default:
throw new IllegalStateException("Illegat state: " + v);
}
direction = new Coordinate(currentX, currentY);
coordinateList.add(direction);
x = x + dx;
y = y + dy;
awtPolygon.addPoint(x, y);
} while( x != initialX || y != initialY );
double polygonArea = GeometryUtilities.getPolygonArea(awtPolygon.xpoints, awtPolygon.ypoints, coordinateList.size() - 1);
if (polygonArea < pThres) {
return null;
}
GeometryFactory gf = GeometryUtilities.gf();
coordinateList.add(startCoordinate);
Coordinate[] coordinateArray = (Coordinate[]) coordinateList.toArray(new Coordinate[coordinateList.size()]);
LinearRing linearRing = gf.createLinearRing(coordinateArray);
Polygon polygon = gf.createPolygon(linearRing, null);
return polygon;
}
private int value( int x, int y ) {
int sum = 0;
if (isSet(x, y)) // UL
sum |= 1;
if (isSet(x + 1, y)) // UR
sum |= 2;
if (isSet(x, y + 1)) // LL
sum |= 4;
if (isSet(x + 1, y + 1)) // LR
sum |= 8;
if (sum == 0) {
System.out.println(x + "/" + y);
}
// mark the used position
bitSet.set(y * width + x);
return sum;
}
private boolean isSet( int x, int y ) {
boolean isOutsideGrid = x < 0 || x >= width || y < 0 || y >= height;
if (isOutsideGrid) {
return false;
}
double value = iter.getSampleDouble(x, y, 0);
if (isNovalue(value)) {
return false;
}
if (value == pValue) {
return true;
}
return false;
}
}
