Δημοσιεύτηκε: 13 Απρ 2012, 01:54
Για τα bindings του proj σε python ένα minimum working example είναι αυτό.
- Κώδικας: Επιλογή όλων
from pyproj import Proj
import numpy as np
to_WGS84 = Proj(init='epsg:4326')
to_EGSA87 = Proj(init='epsg:2100')
x_old = np.array((23.72826, 37.98414))
y_old = np.array((37.98414, 37.98414))
x_new, y_new = to_EGSA87(x_old, y_old)
for x, y in zip(x_new, y_new):
print x, y