Chuyển đổi các đối tượng Shapely Polygon và MultiPolygon

Có cách nào dễ dàng để chuyển đổi các đối tượng Shapely (cụ thể là Đa giác và MultiPolygons) từ phép chiếu này sang phép chiếu khác mà không phải đào xung quanh và trích xuất tọa độ bằng tay không?

Trên thực tế, tôi thậm chí không quan tâm nếu chúng là các đối tượng Shapely vào thời điểm này, tôi chỉ muốn truyền các tính năng và phép chiếu, và lấy lại một bộ các tính năng bị từ chối.

Là loại chức năng này tồn tại, hoặc nó phải được mã hóa bằng tay?

Mặc dù hoàn toàn không hiểu hệ thống tọa độ, nhưng shapely.ops.transform()có thể làm điều đó cùng với pyproj. Nếu pyproj.Projcó thể hiểu cả hai hệ tọa độ của bạn, thì nó có thể được tạo thành một hàm có thể biến đổi với.

Từ các tài liệu tạo hình :

from functools import partial
import pyproj
from shapely.ops import transform

project = partial(
    pyproj.transform,
    pyproj.Proj(init='epsg:4326'), # source coordinate system
    pyproj.Proj(init='epsg:26913')) # destination coordinate system

g2 = transform(project, g1)  # apply projection

Mặc dù không phải là một giải pháp Shapely, sử dụng GeoPandas cho phép chiếu tương đối đơn giản. Ví dụ: nếu chúng tôi muốn chuyển đổi một shapefile thành ESPG 4326:

import geopandas as gpd

HabModelEnviro = gpd.GeoDataFrame.from_file('data/HabModelEnviro.shp').replace({-999: None})

HabModelEnviroWGS84 = HabModelEnviro.to_crs({'proj':'longlat', 'ellps':'WGS84', 'datum':'WGS84'})

Nếu bạn đang sử dụng pyproj2, việc sử dụng Biến áp sẽ dễ dàng hơn nhiều. Đây là một ví dụ:

import pyproj
from shapely.ops import transform

project = pyproj.Transformer.from_proj(
    pyproj.Proj(init='epsg:4326'), # source coordinate system
    pyproj.Proj(init='epsg:26913')) # destination coordinate system

# g1 is a shapley Polygon

g2 = transform(project.transform, g1)  # apply projection

Điều này cũng nhanh hơn nhiều, vì pyasej không cần phải tạo lại hình chiếu cho mọi điểm.