def test_convert_002(self): # <node id="583525261" version="6" timestamp="2018-09-20T08:38:47Z" uid="8354582" user="******" changeset="62755031" lat="44.0472726" lon="15.3004193"> x, y = deg2num(44.0472726, 15.3004193, 9) map = num2MapBB(x, y, 9, conversion_type.normal) print("bb_norm z=9 {}".format(map.getinfo())) print("bb {}".format(map.geojson())) print("bb x={},y={},z={}".format(x, y, 9))
def test_convert_ext_z12(self): # sample: 2176-1315-12.osm # <bounds minlat='53.90433816' minlon='11.16210938' maxlat='54.05938789' maxlon='11.42578125'/> tilemap = num2MapBB(2176, 1315, 12, conversion_type.extendet) print("bb_ext z=12 {}".format(tilemap.getinfo())) self.assertAlmostEqual(tilemap.minlat, 53.90433816) self.assertAlmostEqual(tilemap.minlon, 11.16210938) self.assertAlmostEqual(tilemap.maxlat, 54.05938789) self.assertAlmostEqual(tilemap.maxlon, 11.42578125)
def test_convert_ext_z11(self): # sample: 1089-657-11.osm # <bounds minlat='53.85252660' minlon='11.25000000' maxlat='54.16243397' maxlon='11.77734375'/> tilemap = num2MapBB(1089, 657, 11, conversion_type.extendet) print("bb_ext z=12 {}".format(tilemap.getinfo())) self.assertAlmostEqual(tilemap.minlat, 53.85252660) self.assertAlmostEqual(tilemap.minlon, 11.25000000) self.assertAlmostEqual(tilemap.maxlat, 54.16243397) self.assertAlmostEqual(tilemap.maxlon, 11.77734375)
def test_convert_ext_z9(self): # sample 272-164-9.osm # <bounds minlat='53.64463782' minlon='11.07421875' maxlat='54.26522408' maxlon='12.12890625'/> tilemap = num2MapBB(272, 164, 9, conversion_type.extendet) print("bb_ext z=9 {}".format(tilemap.getinfo())) self.assertAlmostEqual(tilemap.minlat, 53.64463782) self.assertAlmostEqual(tilemap.minlon, 11.07421875) self.assertAlmostEqual(tilemap.maxlat, 54.26522408) self.assertAlmostEqual(tilemap.maxlon, 12.12890625)
def test_convert_ext_z10(self): # sample: 544-328-10.osm # <bounds minlat='53.85252660' minlon='11.07421875' maxlat='54.26522408' maxlon='11.77734375'/> tilemap = num2MapBB(544, 328, 10, conversion_type.extendet) print("bb_ext z=12 {}".format(tilemap.getinfo())) self.assertAlmostEqual(tilemap.minlat, 53.85252660) self.assertAlmostEqual(tilemap.minlon, 11.07421875) self.assertAlmostEqual(tilemap.maxlat, 54.26522408) self.assertAlmostEqual(tilemap.maxlon, 11.77734375)
parser.add_argument('-z', dest='z', default="9", help='tile z') args = parser.parse_args() x=int(args.x) y=int(args.y) z=int(args.z) bb = num2MapBB(x, y, z, conversion_type.extendet) link = "http://map.openseamap.org/" link +="?zoom={}".format(z) link += "&lat={}&lon={}".format(bb.lat, bb.lon) link += "&mlat={}&mlon={}".format(bb.lat, bb.lon) link += "&mtext=x{}_y{}_z{}".format(x,y,z) link += "&layers=BFTFFFFFFTF0FFFFFFFFFF" print(link) print("\n\nquery for https://overpass-turbo.eu/") print(bb.GetOverpassQuery()) print("\n\ngeojson info:") print(bb.geojson()())
print("generate config file for command \"osmium extract\": {}".format( in_filename)) print("input file: {}".format(in_filename)) print("output dir: {}".format(args.outdir)) print("x_min={}, x_max={}, y_min={}, y_max={}, z={}".format( x_min_z9, x_max_z9, y_min_z9, y_max_z9, z)) # stage 1 / generate config file with level 9 maplist = list() for y_z9 in range(y_min_z9, y_max_z9 + 1): for x_z9 in range(x_min_z9, x_max_z9 + 1): maplist.append([ x_z9, y_z9, z, num2MapBB(x_z9, y_z9, z, conversion_type.extendet) ]) confdir = args.outdir + "cfg/" ensure_dir(confdir) outfile = confdir + "osmium_extract_z9.cfg" if os.path.exists(outfile): os.remove(outfile) osmdir = args.outdir + "osm-extracts/" ensure_dir(osmdir) write_configfile(outfile, osmdir, maplist) convert_cmd = "osmium extract --overwrite --config {} {}".format( outfile, in_filename)
def test_convert_002(self): tilemap = num2MapBB(272, 164, 9, conversion_type.normal) print("bb_norm z=9 {}".format(tilemap.getinfo()))
def test_convert_ext_z9_max(self): tilemap = num2MapBB(511, 511, 9, conversion_type.extendet) print("bb_ext z=9 {}".format(tilemap.getinfo())) #self.assertAlmostEqual(tilemap.minlat, -90) self.assertAlmostEqual(tilemap.maxlon, 180)