polygonizer.add(g);
}
ArrayList polys = (ArrayList) polygonizer.getPolygons();
final SpatialIndex idx = new STRtree();
FeatureIterator mt = areaMainCollection.features();
try {
while (mt.hasNext()) {
try {
SimpleFeature feat = (SimpleFeature) mt.next();
Geometry point = (Geometry) feat.getAttribute(geomName
+ "_point");
idx.insert(point.getEnvelopeInternal(), feat);
} catch (NullPointerException e) {
// Gibts anscheinend in anderen Kantonen.
logger.warn("Empty point found.");
}
}
} finally {
mt.close();
}
ArrayList<Polygon> holes = new ArrayList();
for (int i = 0; i < polys.size(); i++) {
boolean found = false;
Polygon p = (Polygon) polys.get(i);
PreparedPolygon ppoly = new PreparedPolygon(p);
for (final Object o : idx.query(p.getEnvelopeInternal())) {
SimpleFeature f = (SimpleFeature) o;
// Funktioniert auch intersect? Schneller?
if (ppoly.contains(((Geometry) f.getAttribute(geomName
+ "_point")))) {
f.setAttribute(geomName, p);
features.add(f);
found = true;
break;
}
}
if (!found) {
// Falls wirklich bewusste Löcher im Datensatz sind, führt
// das zu falschen
// Resultaten. Nur zu den "holes" hinzufügen, falls kleiner
// als ein bestimmter Wert.
if (p.getArea() < 1.0) {
holes.add(p);
}
}
}
// Restflächen, die aufgrund des neu Verknotens entstanden sind,
// werden dem Polygon mit der längsten gemeinsamen Kanten
// zugewiesen.
// Mal schauen obs schon hinhaut.
boolean repair = true;
if (repair) {
final SpatialIndex spatialIndex = new STRtree();
Iterator nt = features.iterator();
while (nt.hasNext()) {
SimpleFeature feat = (SimpleFeature) nt.next();
Geometry point = (Geometry) feat.getAttribute(geomName);
spatialIndex.insert(point.getEnvelopeInternal(), feat);
}
logger.debug("Anzahl Löcher: ");
logger.debug(holes.size());
for (int i = 0; i < holes.size(); i++) {
Polygon p = holes.get(i);
PreparedPolygon ppoly = new PreparedPolygon(p);
double length = 0.0;
SimpleFeature feat = null;
for (final Object o : spatialIndex.query(p
.getEnvelopeInternal())) {
SimpleFeature f = (SimpleFeature) o;
Geometry g = ((Geometry) f.getAttribute(geomName));
if (ppoly.intersects(g)) {
// logger.debug("***************************");