Mengubah objek Shapely Polygon dan MultiPolygon

Apakah ada cara mudah untuk mengubah objek Shapely (yaitu, Polygons dan MultiPolygons) dari satu proyeksi ke yang lain tanpa harus menggali sekitar dan mengekstrak koordinat dengan tangan?

Sebenarnya, saya bahkan tidak peduli apakah mereka objek Shapely pada titik ini, saya hanya ingin melewatkan fitur dan proyeksi, dan mendapatkan serangkaian fitur yang diproyeksikan kembali.

Apakah fungsi semacam ini ada, atau haruskah kode tangan?

Sementara rupanya Shapely tidak memahami sistem koordinat, shapely.ops.transform()dapat melakukannya bersama pyproj. Jika Anda pyproj.Projdapat memahami kedua sistem koordinat Anda, maka itu dapat dibuat menjadi fungsi yang dapat ditransformasi secara rupawan.

Dari dokumen indah :

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

Meskipun bukan solusi Shapely, menggunakan GeoPandas memungkinkan proyeksi yang relatif mudah. Misalnya, jika kita ingin mengonversi shapefile ke 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'})

Jika Anda menggunakan pyproj2, jauh lebih mudah menggunakan Transformer. Ini sebuah contoh:

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

Ini juga jauh lebih cepat, karena pyproj tidak perlu membuat ulang proyeksi untuk setiap titik.