Python 的 Geodaisy 库转换带负号坐标的WKT时的Bug

Geodaisy是一个Python库，用于创建由类型和坐标表示的地理对象，并在各种标准和表示之间进行转换，包括GeoJSON、Well - Known Text和Python的__geo_interface__协议，它被其他地理库使用。

在使用其 converters.wkt_to_geo_interface（）方法时转换带符号坐标的wkt字符串时，发现其转换的结果跟预期的不一样：

例如：

1 2 3

> wkt = "MULTILINESTRING ((3146.2134801800566 2556399.9823166174, -2611.53319329879 2564044.0076796883))" > converters.wkt_to_geo_interface(wkt) # {'type': 'MultiLineString', 'coordinates': (3146.2134801800566, 2556399.9823166174, -2611.53319329879, 2564044.0076796883)}

我们看到 coordinates 的值并不是我们想要的结果：

1	{'type': 'MultiLineString', 'coordinates': [(3146.2134801800566, 2556399.9823166174), (-2611.53319329879, 2564044.0076796883)]}

原因是 wkt_to_geo_interface 方法源码采用正则匹配，并没有考虑带负数坐标的情况（或者根本就不考虑），但是在地方坐标系处理偏移的过程中是肯定有出现负值坐标的情况的。

解决方法：

修改wkt_to_geo_interface方法源码的正则匹配（匹配”-“号和数字）：

def wkt_to_geo_interface(wkt):
    # type: (str) -> dict
    """Converts a WKT string to a geo_interface dictionary."""
    try:
        wkt_type, coords = re.split(r'(?<=[A-Z])\s', wkt)

        geo_type = type_translations[wkt_type]

        # Clean up the strings so they'll covert correctly
        if geo_type in {'Polygon', 'MultiLineString', 'MultiPolygon'}:
            # coords = re.sub(r'(?<=\d)\), \((?=\d)', ')), ((', coords)
            coords = re.sub(r'(?<=[-\d])\), \((?=[-\d])', ')), ((', coords)   # modif by 举个栗子

        # Pairs of coordinates must be enclosed in parentheses
        # coords = re.sub(r'(?<=\d), (?=\d)', '), (', coords)
        coords = re.sub(r'(?<=[-\d]), (?=[-\d])', '), (', coords) # modif by 举个栗子

        # Coordinates within parentheses must be separated by commas
        # coords = re.sub(r'(?<=\d) (?=\d)', ', ', coords)
        coords = re.sub(r'(?<=[-\d]) (?=[-\d])', ', ', coords)    # # modif by 举个栗子

        # Now we can turn the string into a tuple or a tuple of tuples
        coords = literal_eval(coords)

        coords = reformat_coordinates(coords, 'geo_interface')  # type: ignore  # noqa: E501

        # If we only have a simple polygon (no hole), the coordinate array
        # won't be deep enough to satisfy the GeoJSON/geo_interface spec, so
        # we need to enclose it in a list.
        numbers = {float, int}
        if geo_type == 'Polygon' and type(coords[0][0]) in numbers:
            coords = [coords]  # type: ignore
        elif geo_type == 'MultiPolygon' and type(coords[0][0][0]) in numbers:
            coords = [coords]  # type: ignore

    except Exception:
        raise ValueError('{} is not a WKT string'.format(wkt))

    return {'type': geo_type, 'coordinates': coords}