wkt convert osmosis .poly ( PolygonFilter input )

Posted 2022-10-25 Hepburn Yang

MULTIPOLYGON WKT 数据地址 21523.txt

Osmosis/Polygon Filter File Format

wkt 转 poly

import sys
import numpy as np

for 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:
                continue
            f.write(f'areaidx\\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 数据