python - 从 WKT 字符串中提取 XY 范围

Question

使用 python，如何提取 WKT 多边形的 XY 范围？我需要 gdal 兼容格式，例如:minx miny maxx maxy

因此，例如下面的例子，我需要将 wkt 转换为 wkt_extent:

wkt = "Polygon ((366247 156971, 366247 174054, 383331 174054, 383331 156971, 366247 156971))"
# need xy bounding box coordinates of wkt, as shown below
wkt_extent = "366247 156971 383331 174054"

请记住，wkt 并不总是矩形/正方形，并且并不总是像此处所示的示例 wkt 那样按顺时针方向“绘制”。

这里的SRID是27700，单位是米。

Answer 1

假设您的 WKT 作为字符串传递:

def parse_geometry(geometry):
    regex = r'[0-9-\.]+'
    parsed_geom = re.findall(regex, geometry)
    parsed_geom = [float(i) for i in parsed_geom]
    return max(parsed_geom[::2]), min(parsed_geom[::2]), max(parsed_geom[1::2]), min(parsed_geom[1::2])