MULTIPOLYGON WKT 数据地址 21523.txt

Osmosis/Polygon Filter File Format

wkt 转 poly

import sys
import numpy as npfor ipath in sys.argv[1:]:if not ipath.endswith('.txt'):print(f'skip {ipath}')opath = f'{ipath}.poly'print(opath)with open(ipath) as f:text = f.read()with open(opath, 'w') as f:f.write('geofence\n')cols = text.split('(')for idx, col in enumerate(cols):if len(col) < 100:continuef.write(f'area{idx}\n')coords = np.array([[float(x) for x in lla.split()] for lla in col[:-6].split(',')])for c in coords:f.write(f'    {c[0]} {c[1]}\n')f.write('END\n')# print(coords)f.write('END\n')print(f'wrote to {opath}')

输出 21523.txt.poly

开始过滤数据

 osmosis --read-xml file="msd.osm"--bounding-polygon file="21523.txt.poly" completeWays="yes" clipIncompleteEntities="true" --write-xml file="21523.osm"

从 msd.osm 中以 21523.txt.poly 为 boundary 过滤出了高亮部分的 osm 数据