Hi,
it seems that pyproj.Proj does not take into account the "towgs84" optional
parameter.
Here is what I simply tried:
from mpl_toolkits.basemap import *
lonref = -3.
latref = 47.
kwproj = dict(proj="lcc", a=6378249.2, b=6356515., x_0=600000.,
y_0=200000.,
lon_0="2d20'14.025", lat_0="46d48'", lat_1="45d53'56.108",
lat_2="47d41'45.652")
print pyproj.Proj(**kwproj)(lonref, latref)
kwproj.update(towgs84="-168,-60,+320")
print pyproj.Proj(**kwproj)(lonref, latref)
Both transforms give the same result.
Try the same parameters with "cs2cs" changes the results as expected.
Did I miss something?
--
Stephane Raynaud
------------------------------------------------------------------------------
This SF.net email is sponsored by:
SourcForge Community
SourceForge wants to tell your story.
http://p.sf.net/sfu/sf-spreadtheword
_______________________________________________
Matplotlib-users mailing list
Matplotlib-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/matplotlib-users
