Δημοσιεύτηκε: 19 Απρ 2012, 02:00
pmav99 έγραψε:Για τα 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
αυτό δίνει;
- Κώδικας: Επιλογή όλων
476136.794863 4204090.14954
1730940.53067 4297555.98423