Index: src/uk/me/parabola/mkgmap/reader/osm/SeaGenerator.java =================================================================== --- src/uk/me/parabola/mkgmap/reader/osm/SeaGenerator.java (revision 4600) +++ src/uk/me/parabola/mkgmap/reader/osm/SeaGenerator.java (working copy) @@ -286,6 +286,7 @@ System.err.println(" land-tag=TAG=VAL tag to use for land polygons (default natural=land)"); if (forPrecompSea) return; + System.err.println(" no-sea-sectors disable use of \"sea sectors\""); System.err.println(" extend-sea-sectors extend coastline to reach border"); System.err.println(" close-gaps=NUM close gaps in coastline that are less than this distance (metres)"); @@ -294,7 +295,7 @@ System.err.println(" fbthres=NUM min points contained in a polygon to be flood blocked (default 20)"); System.err.println(" fbratio=NUM min ratio (points/area size) for flood blocking (default 0.5)"); } - + /** * Read the index from stream and populate the index grid. * @param fileStream already opened stream @@ -602,9 +603,9 @@ } // check if the land tags need to be changed - boolean changeLadTag = landTag != null && ("natural".equals(landTag[0]) && !"land".equals(landTag[1])); + boolean changeLandTag = landTag != null && ("natural".equals(landTag[0]) && !"land".equals(landTag[1])); for (Way w : landWays) { - if (changeLadTag) { + if (changeLandTag) { w.deleteTag("natural"); w.addTag(landTag[0], landTag[1]); } @@ -760,7 +761,7 @@ } else if (w.getPoints().size() > 1){ beginMap.put(w.getFirstPoint(), w); } else { - log.info("Discarding coastline way", w.getId(), "because it consists of less than 2 points"); + log.info("Discarding coastline", w.getBasicLogInformation(), "because it consists of less than 2 points"); } } segments.clear(); @@ -785,7 +786,7 @@ // merge the ways and maintain maps and list private static void merge(Map beginMap, List joined, Way w1, Way w2) { - log.info("merging: ", beginMap.size(), w1.getId(), w2.getId()); + log.info("merging:", beginMap.size(), w1.getBasicLogInformation(), "with", w2.getBasicLogInformation()); Way wm; if (FakeIdGenerator.isFakeId(w1.getId())) { wm = w1; @@ -831,7 +832,7 @@ } // clip all shoreline segments - clipShorlineSegments(); + clipShorelineSegments(); if(shoreline.isEmpty()) { // No sea required @@ -877,11 +878,11 @@ createLandPolygons(hitMap); processIslands(seaRelation); - processAntiIslands(true); + checkIslands(true); Way sea = createSeaWay(true); - log.info("sea: ", sea); + log.info("sea:", sea); saver.addWay(sea); if(seaRelation != null) seaRelation.addElement("outer", sea); @@ -888,7 +889,7 @@ } else { // background is land createSeaPolygons(hitMap); - processAntiIslands(false); + checkIslands(false); // generate a land polygon so that the tile's // background colour will match the land colour on the @@ -931,26 +932,59 @@ } /** - * These are bits of sea have been generated as polygons. - * if the tile is also sea based, then check that surrounded by an island - * @param seaRelation if set, add as inner + * Check whether land is enclosed in land or sea within sea. * @param seaBased true if the tile is also sea with land [multi-]polygons */ - private void processAntiIslands(boolean seaBased) { + private void checkIslands(boolean seaBased) { for (Way ai : antiIslands) { - if (seaBased) { - boolean containedByLand = false; - for (Way i : islands) { - if (i.containsPointsOf(ai)) { - containedByLand = true; - break; - } + Way containingLand = null; + Way containingSea = null; + for (Way i : islands) { + if (i.containsPointsOf(ai) && ((containingLand == null) || (containingLand.containsPointsOf(i)))) + containingLand = i; + } + for (Way ai2 : antiIslands) { + if ((ai2 != ai) && ai2.containsPointsOf(ai) && ((containingSea == null) || containingSea.containsPointsOf(ai2))) + containingSea = ai2; + } + if ((containingSea != null) && (containingLand != null)) { + if (containingSea.containsPointsOf(containingLand)) + containingSea = null; + else if (containingLand.containsPointsOf(containingSea)) + containingLand = null; + else { + log.warn("inner sea", ai, "is surrounded by both water", containingSea, "and land", containingLand); + containingSea = null; + containingLand = null; } - if (!containedByLand) { - // found an anti-island that is not contained by land - log.warn("inner sea", ai , "is surrounded by water"); + } + if ((containingLand == null) & (seaBased | (containingSea != null))) + log.error("inner sea", ai, "is surrounded by water", containingSea == null ? "" : containingSea); + } + for (Way i : islands) { + Way containingLand = null; + Way containingSea = null; + for (Way ai : antiIslands) { + if (ai.containsPointsOf(i) && ((containingSea == null) || containingSea.containsPointsOf(ai))) + containingSea = ai; + } + for (Way i2 : islands) { + if ((i2 != i) && i2.containsPointsOf(i) && ((containingLand == null) || (containingLand.containsPointsOf(i2)))) + containingLand = i2; + } + if ((containingSea != null) && (containingLand != null)) { + if (containingSea.containsPointsOf(containingLand)) + containingSea = null; + else if (containingLand.containsPointsOf(containingSea)) + containingLand = null; + else { + log.warn("island", i, "is surrounded by both water", containingSea, "and land", containingLand); + containingSea = null; + containingLand = null; } } + if ((containingSea == null) && (containingLand != null)) + log.error("island", i, "is surrounded by land", containingLand); } } @@ -988,7 +1022,7 @@ /** * Clip the shoreline ways to the bounding box of the map. */ - private void clipShorlineSegments() { + private void clipShorelineSegments() { List toBeRemoved = new ArrayList<>(); List toBeAdded = new ArrayList<>(); for (Way segment : shoreline) { @@ -1110,7 +1144,7 @@ if (Way.clockwise(w.getPoints())) addAsSea(w); else - addAsLand(w); + addAsLand(w); } private void addAsSea(Way w) { @@ -1136,9 +1170,10 @@ private void createLandPolygons(NavigableMap hitMap) { NavigableSet hits = hitMap.navigableKeySet(); while (!hits.isEmpty()) { - Way w = new Way(FakeIdGenerator.makeFakeId()); Double hFirst = hits.first(); Double hStart = hFirst, hEnd; + Way w = new Way(hitMap.get(hFirst).getOriginalId()); + w.markAsGeneratedFrom(hitMap.get(hFirst)); boolean finished = false; do { Way segment = hitMap.get(hStart); @@ -1174,9 +1209,10 @@ private void createSeaPolygons(NavigableMap hitMap) { NavigableSet hits = hitMap.navigableKeySet(); while (!hits.isEmpty()) { - Way w = new Way(FakeIdGenerator.makeFakeId()); Double hFirst = hits.last(); Double hStart = hFirst, hEnd; + Way w = new Way(hitMap.get(hFirst).getOriginalId()); + w.markAsGeneratedFrom(hitMap.get(hFirst)); boolean finished = false; do { Way segment = hitMap.get(hStart); @@ -1245,7 +1281,7 @@ Double hEnd = getEdgeHit(tileBounds, pEnd); if (hStart != null && hEnd != null) { // nice case: both ends touch the boundary - log.debug("hits: ", hStart, hEnd); + log.debug("hits:", hStart, hEnd); hitMap.put(hStart, w); hitMap.put(hEnd, null); // put this for verifyHits which then deletes it } else { @@ -1268,7 +1304,8 @@ */ List points = w.getPoints(); if (allowSeaSectors) { - Way seaOrLand = new Way(FakeIdGenerator.makeFakeId()); + Way seaOrLand = new Way(w.getOriginalId()); + seaOrLand.markAsGeneratedFrom(w); seaOrLand.getPoints().addAll(points); int startLat = pStart.getHighPrecLat(); int startLon = pStart.getHighPrecLon(); @@ -1291,7 +1328,7 @@ } seaOrLand.addPoint(Coord.makeHighPrecCoord(cornerLat, cornerLon)); seaOrLand.addPoint(pStart); - log.info("seaSector: ", generateSeaBackground, startLatIsCorner, Way.clockwise(seaOrLand.getPoints()), seaOrLand); + log.info("seaSector:", generateSeaBackground, startLatIsCorner, Way.clockwise(seaOrLand.getPoints()), seaOrLand.getBasicLogInformation()); } else if (extendSeaSectors) { // join to nearest tile border if (null == hStart) { @@ -1302,7 +1339,7 @@ hEnd = getNextEdgeHit(tileBounds, pEnd); w.getPoints().add(getPoint(tileBounds, hEnd)); } - log.debug("hits (second try): ", hStart, hEnd); + log.debug("hits (second try):", hStart, hEnd); hitMap.put(hStart, w); hitMap.put(hEnd, null); // put this for verifyHits which then deletes it } else { @@ -1328,6 +1365,7 @@ NavigableSet hits = hitMap.navigableKeySet(); Iterator iter = hits.iterator(); int lastStatus = 0, thisStatus; + Double lastHit = 0.0; while (iter.hasNext()) { Double aHit = iter.next(); Way segment = hitMap.get(aHit); @@ -1339,14 +1377,15 @@ thisStatus = +1; } if (thisStatus == lastStatus) - log.error("Adjacent coastlines hit tile edge in same direction", aHit, segment); + log.error("Adjacent coastlines hit tile edge in same direction at", getPoint(tileBounds, lastHit).toDegreeString(), "and", getPoint(tileBounds, aHit).toDegreeString(), segment); lastStatus = thisStatus; + lastHit = aHit; } } // create the point where the shoreline hits the sea bounds private static Coord getPoint(Area a, double edgePos) { - log.info("getPoint: ", a, edgePos); + log.info("getPoint:", a, edgePos); int aMinLongHP = a.getMinLong() << Coord.DELTA_SHIFT; int aMaxLongHP = a.getMaxLong() << Coord.DELTA_SHIFT; int aMinLatHP = a.getMinLat() << Coord.DELTA_SHIFT;