Index: /applications/editors/josm/plugins/smed2/jsearch/build.xml =================================================================== --- /applications/editors/josm/plugins/smed2/jsearch/build.xml (revision 30408) +++ /applications/editors/josm/plugins/smed2/jsearch/build.xml (revision 30408) @@ -0,0 +1,29 @@ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + Index: /applications/editors/josm/plugins/smed2/jsearch/src/jsearch/Extract.java =================================================================== --- /applications/editors/josm/plugins/smed2/jsearch/src/jsearch/Extract.java (revision 30408) +++ /applications/editors/josm/plugins/smed2/jsearch/src/jsearch/Extract.java (revision 30408) @@ -0,0 +1,257 @@ +/* Copyright 2014 Malcolm Herring + * + * This 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, version 3 of the License. + * + * For a copy of the GNU General Public License, see . + */ + +package jsearch; + +import java.io.BufferedReader; +import java.io.FileReader; +import java.io.IOException; +import java.util.ArrayList; +import java.util.HashMap; + +import jsearch.Jsearch.MapBB; + +public class Extract { + + public static ArrayList extractData(String filename, MapBB box) throws IOException { + + HashMap nodes = new HashMap(); + HashMap ways = new HashMap(); + HashMap rels = new HashMap(); + + ArrayList buf = new ArrayList(); + BufferedReader in; + boolean inOsm = false; + boolean inNode = false; + boolean inWay = false; + boolean inRel = false; + + double lat = 0; + double lon = 0; + long id = 0; + ArrayList refs = null; + boolean needed = false; + String ln; + + in = new BufferedReader(new FileReader(filename)); + while ((ln = in.readLine()) != null) { + if (inOsm) { + if (inNode) { + if (ln.contains(" box.minlat) && (lat <= box.maxlat) && (lon > box.minlon) && (lon <= box.maxlon)) { + nodes.put(id, null); + } + if (!ln.contains("/>")) { + inNode = true; + } + } else if (inWay) { + if (ln.contains("(); + needed = false; + if (!ln.contains("/>")) { + inWay = true; + } + } else if (inRel) { + if (ln.contains("(); + needed = false; + if (!ln.contains("/>")) { + inRel = true; + } + } else if (ln.contains("")) { + needed = ways.containsKey(id); + inWay = true; + } + } else if (ln.contains(""); + buf.add(""); + buf.add(String.format("", box.minlat, box.minlon, box.maxlat, box.maxlon)); + while ((ln = in.readLine()) != null) { + if (inOsm) { + if (inNode) { + buf.add(ln); + if (ln.contains("")) { + inNode = true; + } + } + } else if (inWay) { + buf.add(ln); + if (ln.contains("")) { + inWay = true; + } + } + } else if (inRel) { + buf.add(ln); + if (ln.contains("")) { + inRel = true; + } + } + } else if (ln.contains(". + */ + +package jsearch; + +import java.io.BufferedReader; +import java.io.File; +import java.io.FileReader; +import java.io.PrintStream; +import java.util.ArrayList; +import java.util.HashMap; + +public class Jsearch { + public static class MapBB { + public double minlat; + public double minlon; + public double maxlat; + public double maxlon; + } + + public static void main(String[] args) throws Exception { + HashMap cnodes = new HashMap(); + HashMap cways = new HashMap(); + HashMap crels = new HashMap(); + HashMap nnodes = new HashMap(); + HashMap nways = new HashMap(); + HashMap nrels = new HashMap(); + String dir = args[0]; + long id = 0; + HashMap z9s = new HashMap(); + HashMap z10s = new HashMap(); + HashMap z11s = new HashMap(); + HashMap z12s = new HashMap(); + BufferedReader in = new BufferedReader(new FileReader(dir + "diffs")); + String ln; + while ((ln = in.readLine()) != null) { + for (String token : ln.split("[ ]+")) { + if (token.matches("^id=.+")) { + id = Long.parseLong(token.split("[\"\']")[1]); + break; + } + } + if (ln.matches("^<.+")) { + if (ln.contains(".+")) { + if (ln.contains(" buf = new ArrayList(); + while ((ln = in.readLine()) != null) { + if (inOsm) { + if (inNode) { + if (ln.contains("")) { + inNode = true; + } + } else if (inWay) { + if (ln.contains("")) { + inWay = true; + } + } else if (inRel) { + if (ln.contains("")) { + inRel = true; + } + } else if (ln.contains("")) { + inWay = true; + } + } else if (ln.contains("= 0) && (y <= 4095)) + z12s.put((((x < 0) ? 4095 : (x > 4095) ? 0 : x) * 4096) + y, true); + } + } + } + } + } + next = !next; + } while (next); + + for (int t : z9s.keySet()) { + int x = (t / 512) * 8; + int y = (t % 512) * 8; + MapBB bb = new MapBB(); + bb.minlon = tile2lon((x + 4095) % 4096, 12); + bb.maxlon = tile2lon((x + 9) % 4095, 12); + bb.minlat = tile2lat(Math.min((y + 9), 4095), 12); + bb.maxlat = tile2lat(Math.max((y - 1), 0), 12); + ArrayList ext = Extract.extractData(dir + "next.osm", bb); + PrintStream out = new PrintStream(dir + "tmp/" + (t / 512) + "-" + (t % 512) + "-9.osm"); + for (String line : ext) { + out.println(line); + } + out.close(); + } + for (int t : z10s.keySet()) { + int x = (t / 1024) * 4; + int y = (t % 1024) * 4; + MapBB bb = new MapBB(); + bb.minlon = tile2lon((x + 4095) % 4096, 12); + bb.maxlon = tile2lon((x + 5) % 4095, 12); + bb.minlat = tile2lat(Math.min((y + 5), 4095), 12); + bb.maxlat = tile2lat(Math.max((y - 1), 0), 12); + ArrayList ext = Extract.extractData(dir + "tmp/" + ((t / 1024) / 2) + "-" + ((t % 1024) / 2) + "-9.osm", bb); + PrintStream out = new PrintStream(dir + "tmp/" + (t / 1024) + "-" + (t % 1024) + "-10.osm"); + for (String line : ext) { + out.println(line); + } + out.close(); + } + for (int t : z11s.keySet()) { + int x = (t / 2048) * 2; + int y = (t % 2048) * 2; + MapBB bb = new MapBB(); + bb.minlon = tile2lon((x + 4095) % 4096, 12); + bb.maxlon = tile2lon((x + 3) % 4095, 12); + bb.minlat = tile2lat(Math.min((y + 3), 4095), 12); + bb.maxlat = tile2lat(Math.max((y - 1), 0), 12); + ArrayList ext = Extract.extractData(dir + "tmp/" + ((t / 2048) / 2) + "-" + ((t % 2048) / 2) + "-10.osm", bb); + PrintStream out = new PrintStream(dir + "tmp/" + (t / 2048) + "-" + (t % 2048) + "-11.osm"); + for (String line : ext) { + out.println(line); + } + out.close(); + } + + System.exit(0); + } + + static int lon2xtile(double lon, int zoom) { + int xtile = (int) Math.floor((lon + 180) / 360 * (1 << zoom)); + if (xtile < 0) + xtile = 0; + if (xtile >= (1 << zoom)) + xtile = ((1 << zoom) - 1); + return (xtile); + } + + static int lat2ytile(double lat, int zoom) { + int ytile = (int) Math.floor((1 - Math.log(Math.tan(Math.toRadians(lat)) + 1 / Math.cos(Math.toRadians(lat))) / Math.PI) / 2 * (1 << zoom)); + if (ytile < 0) + ytile = 0; + if (ytile >= (1 << zoom)) + ytile = ((1 << zoom) - 1); + return (ytile); + } + + MapBB tile2bb(final int x, final int y, final int zoom) { + MapBB bb = new MapBB(); + bb.maxlat = tile2lat(y, zoom); + bb.minlat = tile2lat(y + 1, zoom); + bb.minlon = tile2lon(x, zoom); + bb.maxlon = tile2lon(x + 1, zoom); + return bb; + } + + static double tile2lon(int x, int z) { + return x / Math.pow(2.0, z) * 360.0 - 180; + } + + static double tile2lat(int y, int z) { + double n = Math.PI - (2.0 * Math.PI * y) / Math.pow(2.0, z); + return Math.toDegrees(Math.atan(Math.sinh(n))); + } +}